Capillary waves

From SklogWiki
Jump to navigation Jump to search

Thermal capillary waves[edit]

Thermal capillary waves are oscillations of an interface which are thermal in origin. These take place at the molecular level, where only the contribution due to surface tension is relevant. Capillary wave theory is a classic account of how thermal fluctuations distort an interface (Ref. 1). It starts from some intrinsic surface that is distorted. In the Monge representation, the surface is given as Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z=h(x,y)} . An increase in area of the surface causes a proportional increase of energy:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_\mathrm{st}= \sigma \iint dx\, dy\ \sqrt{1+\left( \frac{dh}{dx} \right)^2+\left( \frac{dh}{dy} \right)^2} -1 }

for small values of the derivatives (surfaces not too rough):

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_\mathrm{st} \approx \frac{\sigma}{2} \iint dx\, dy\ \left[ \left( \frac{dh}{dx} \right)^2+\left( \frac{dh}{dy} \right)^2 \right]. }

A Fourier analysis treatment begins by writing the intrinsic surface as an infinite sum of normal modes:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h(x,y)= \sum_\vec{q} a_\vec{q} e^{i\vec{q}\vec{r}}.}

Since normal modes are orthogonal, the energy is easily expressible as a sum of terms . Each term of the sum is quadratic in the amplitude; hence equipartition holds, according to standard classical statistical mechanics, and the mean energy of each mode will be Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle k_B T/2} . Surprisingly, this result leads to a divergent surface (the width of the interface is bound to diverge with its area) (Ref 2). This divergence is nevertheless very mild; even for displacements on the order of meters, the deviation of the surface is comparable to the size of the molecules. Moreover, the introduction of an external field removes this divergence: the action of gravity is sufficient to keep the width fluctuation on the order of one molecular diameter for areas larger than about 1 mm2 (Ref. 2). The action of gravity is taken into account by integrating the potential energy density due to gravity, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho g z} from a reference height to the position of the surface, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z=h(x,y)} :

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_\mathrm{g}= \iint dx\, dy\, \int_0^h dz \rho g z = \frac{\rho g}{2} \int dx\, dy\, h^2.}

(For simplicity, one neglects the density of the gas above, which is often acceptable; otherwise, instead of the density the difference in densities appears).

Recently, a procedure has been proposed to obtain a molecular intrinsic surface from simulation data (Ref. 3), the intrinsic sampling method. The density profiles obtained from this surface are, in general, quite different from the usual mean density profiles.

Gravity-capillary waves[edit]

These are ordinary waves excited in an interface, such as ripples on a water surface. Their dispersion relation reads, for waves on the interface between two fluids of infinite depth:

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega} is the angular frequency, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle g} the acceleration due to gravity, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma} the surface tension, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho} and the mass density of the two fluids (Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho > \rho^'} ) and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle k} the is wavenumber.

Derivation[edit]

This is a sketch of the derivation of the general dispersion relation, see Ref. 4 for a more detailed description. The problem is unfortunately a bit complex. As Richard Feynman put it (Ref. 6):

"...[water waves], which are easily seen by everyone and which are used as an example of waves in elementary courses... are the worst possible example... they have all the complications that waves can have"

Defining the problem[edit]

Three contributions to the energy are involved: the surface tension, gravity, and hydrodynamics. The parts due to surface tension (again the derivatives are taken to be small) and gravity are exactly as above. The new contribution involves the kinetic energy of the fluid:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T= \frac{\rho}{2} \iint dx\, dy\, \int_{-\infty}^h dz v^2,}

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v} is the module of the velocity field . (Again, we are neglecting the flow of the gas above for simplicity.)

Wave solutions[edit]

Let us suppose the surface of the liquid is described by a traveling plane wave:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h(x,y,t)=\eta(t)e^{i\vec{q}\vec{r}},}

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta(t)=\exp[i\omega t]} and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{q}=(q_x,q_y)} is a two dimensional wave number vector, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{r}=(x,y)} being the horizontal position. We may take without loss of generality:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h(x,y,t)=\eta(t)e^{i q x}.}

In this case it is easy to perform the integrations involved in the expressions for the energies. The integration over Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x} can taken over a period of oscillation Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \lambda=2\pi/q} , then multiplied by the number of oscillations in our very large (in principle, infinite) system: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_x / \lambda} . The integration over trivially yields Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle L_y} . Performing the integrations, keeping in mind that only the real part of complex numbers is to be taken as physical, one finds:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_\mathrm{g} = \frac{A}{2} \frac{\rho g}{2} \eta^2,}
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_\mathrm{st} = \frac{A}{2} \frac{\sigma}{2} q^2 \eta^2,}

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle A=L_x\times L_y} is the area of the system.

To tackle the kinetic energy, suppose the fluid is incompressible and its flow is irrotational (often, sensible approximations) - the flow will then be potential: , and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \phi} is a potential (scalar field) which must satisfy Laplace's equation Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \nabla^2\phi=0} . If we try try separation of variables with the potential:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \phi(x,y,z,t)=\xi(t) f(z) e^{i q x},}

with some function of time Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \xi(t)} , and some function of vertical component (height) . Laplace's equation then requires on the later

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f''(z)= q^2 f(z) .}

This equation can be solved with the proper boundary conditions: first, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{v}} must vanish well below the surface (in the "deep water" case, which is the one we consider, otherwise a more general relation holds, which is also well known in oceanography). Therefore

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(z) = a \, \exp(|q| z) } ,

with some constant Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a} . The less trivial condition is the proper matching between and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h} : the potential field must correspond to a velocity field that is adjusted to the movement of the surface:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v_z (z=h) =\partial h/\partial t} .

(Actually, this is the linearized version of a more general expression, see below.)

This implies that:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \xi(t)=\eta(t)'} , and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f'(z=h) = 1, }

so that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(z) = \exp( -|q|(h-z))/|q| } .


We may now find the velocity field, , which shows the well-known circles: the elements of fluid undergo circular motion in the Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x,z} plane, with the circles getting smaller at deeper levels. The displacement of a fluid element is given by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \partial\vec{\psi}/\partial t= \vec{v}} , and is plotted in Figure 1.

Figure 1. Displacement of particle (snapshot)

For the kinetic energy, we need Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v^2=|\nabla\phi|^2} , which is:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v^2= (\eta')^2 e^{ -2 |q|(h-z)}, }

with no dependence on or Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle y} ; the other integration provides:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T= \frac{\rho A }{2|q|} ( \eta' )^2 \int_{-\infty}^h e^{ -2 |q|(h-z)} = \frac{A}{2} \frac{\rho }{2|q|} ( \eta' )^2 . }

The problem is thus specified by just a potential energy involving the square of Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \eta(t)} and a kinetic energy involving the square of its time derivative: a regular harmonic oscillator. In particular:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E= \frac{A}{2} \left[ \left( \rho g + {\sigma} q^2 \right) \frac{\eta^2}{2}+ \frac{\rho}{q} \frac{(\eta')^2}{2} \right]. }

Identifying the oscillator's "spring constant" Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \kappa = \rho g + {\sigma} q^2 } , and its "mass" Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle m= \rho / q} , the oscillation frequency must be given by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega^2=\kappa/m} , which results in:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \omega^2=g q+\frac{\sigma}{\rho} q^3,}

the same dispersion as above if Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho'} is neglected.

Alternative derivation[edit]

In Reference 7 the dispersion relation is derived in a somewhat different manner. The same assumptions are made regarding the fluid (it is inviscid, irrotational, and incompressible), so Laplace's Equation is to be satisfied: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \nabla^2\phi=0} , with Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \vec{v}=\nabla\phi} . The boundary conditions, on the other hand, are sufficient to solve the problem.

One boundary condition is the requirement that the surface of the liquid, defined by Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z=h(x,y;t)} follows the velocity field:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial h}{\partial t}+v_x \frac{\partial h}{\partial x}+v_y \frac{\partial h}{\partial y}= v_z . }

A simpler condition follows from linearization: Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \partial h /\partial t =v_z } , as in the previous derivation. There is an additional boundary condition at the bottom of the fluid, which we take here as Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v_z=0} when Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle z\rightarrow -\infty} .

In the same fashion as above, we seek surface wave solutions, of the form Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h(x,y,t)=a e^{i (qx-\omega t)}} . We may guess a solution of the form

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \phi=-i\omega h(x,y;t) f(z).}

This first condition implies Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f'(z=h)=1} . Together with Laplace's equation, this leads to a function

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f=(1/q) \exp(q(z-h)). }

(see Ref 8 for a discussion on when Laplace's equation admits wave solutions.)


The other surface boundary condition is a Bernoulli equation, stating that the pressure just below the surface, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p_-} , must equal the saturation pressure of coexistence, minus a contribution due to the surface:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p_-=p - \rho\left[ \frac{\partial \phi}{\partial t}+ \frac{1}{2} v^2+ gh \right] . }

The linearized condition is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial \phi}{\partial t}+gh = \frac{p-p_-}{\rho} }

The connection with the curvature of the surface can be introduced by Young's equation for the pressure drop across a curved interface, whose linearized form is:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p-p_-=\sigma\left(\frac{\partial^2 h}{\partial x^2}+\frac{\partial^2 h}{\partial y^2}\right), }

where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \sigma} is the surface tension. The linearized condition is finally

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \frac{\partial \phi}{\partial t}+gh = \frac{\sigma}{\rho} \left(\frac{\partial^2 h}{\partial x^2}+\frac{\partial^2 h}{\partial y^2}\right). }

This second condition, when applied to the surface wave above, establishes that Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle f(z=h)=(g+\sigma/\rho q^2)/\omega^2} .

For the two conditions to hold, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle 1/q} must equal Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle (g+\sigma/\rho q^2)/\omega^2} , which is precisely the same dispersion relation as the one above.

This derivation makes clear the assumptions introduced. In particular, the linearization will only hold for smooth waves, the ones for which the wave amplitude, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle a} is smaller than the wavelength. Mathematically, the limit is Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle q a \ll 1} . For ocean waves, this happens when waves approach the shore and the amplitude grows (in this limit, a bottom boundary condition Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle v_z (z=-H) =0} must be employed, and waves are not dispersive, see Ref 7.)

References[edit]

  1. F. P. Buff, R. A. Lovett, and F. H. Stillinger, Jr. "Interfacial density profile for fluids in the critical region" Physical Review Letters 15 pp. 621-623 (1965)
  2. J. S. Rowlinson and B. Widom "Molecular Theory of Capillarity". Dover 2002 (originally: Oxford University Press 1982) ISBN 0486425444
  3. E. Chacón and P. Tarazona "Intrinsic profiles beyond the capillary wave theory: A Monte Carlo study", Physical Review Letters 91 166103 (2003)
  4. Samuel A. Safran "Statistical thermodynamics of surfaces, interfaces, and membranes" Addison-Wesley 1994 ISBN 9780813340791
  5. P. Tarazona, R. Checa, and E. Chacón "Critical Analysis of the Density Functional Theory Prediction of Enhanced Capillary Waves", Physical Review Letters 99 196101 (2007)
  6. R.P. Feynman, R.B. Leighton, and M. Sands "The Feynman lectures on physics" Addison-Wesley 1963. Section 51-4. ISBN 0201021153
  7. W.K. Melville "Surface, gravity and capillary waves", in "Encyclopedia of Ocean Sciences", Eds: Steve A. Thorpe and Karl K. Turekian. Elsevier 2001, page 2916. ISBN 978-0-12-227430-5
  8. F Behroozi "Fluid viscosity and the attenuation of surface waves: a derivation based on conservation of energy", European Journal of Physics 25 115 (2004)

External links[edit]