Front. Phys. Frontiers in Physics Front. Phys. 2296-424X Frontiers Media S.A. 10.3389/fphy.2018.00136 Physics Original Research Closed, Two Dimensional Surface Dynamics Svintradze David V. * School of Health Sciences, University of Georgia, Tbilisi, Georgia

Edited by: Carla M. A. Pinto, Instituto Superior de Engenharia do Porto (ISEP), Portugal

Reviewed by: Q. H. Liu, Hunan University, China; Daniel Luiz Nedel, Universidade Federal da Integraçăo Latino-Americana, Brazil

 *Correspondence: David V. Svintradze dsvintra@yahoo.com; d.svintradze@ug.edu.ge

This article was submitted to Mathematical Physics, a section of the journal Frontiers in Physics

 06 12 2018 2018 6 136 06 07 2018 21 11 2018 Copyright © 2018 Svintradze. 2018 Svintradze

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

We present dynamic equations for two dimensional closed surfaces and analytically solve it for some simplified cases. We derive final equations for surface normal motions by two different ways. The solution of the equations of motions in normal direction indicates that any closed, two dimensional, homogeneous surface with time invariable surface energy density adopts constant mean curvature shape when it comes in equilibrium with environment. In addition, we show that the shape equation is an approximate solution to our equation of motion in the normal direction and is valid for stationary or near to stationary shapes. As an example, we apply the formalism to analyze equilibrium shapes of micelles and explain why they adopt spherical, lamellar, and cylindrical shapes. Theoretical calculation for micellar optimal radius is in good agreement with all atom simulations and experiments.

 shape dynamics differential geometry membrane dynamics micelles fluid film dynamics 1. Introduction

Biological systems exhibit a variety of morphologies and experience large shape deformations during a motion. Such “choreography” of shape motility is characteristic not only for all living organisms and cells [1] but also for proteins, nucleic acids, and to all biomacromolecules in general. Shape motility, which is a motion of two-dimensional surfaces, may be a result of active (by consuming energy) or passive (without consuming energy) processes. The time scale for shape dynamics may vary from slow (nanometer per nanoseconds) to very fast (nanometer per femtosecond) [2, 3]. Slowly moving surfaces are considered as overdamped systems. An example is cell motility. In that case one may use well developed the Helfrich formalism to describe the motion, where free energy expansion in powers of the curvature tensor gives coarse-grained modeling of membranes [4]. However, while the formalism [4] are applicable to slowly moving surfaces they are not applicable to fast moving surfaces, where biomolecules could be fitted. Surface dynamics for proteins or DNA [2, 3] may reach nm/fs range. So that surfaces may be represented as virtual three dimensional pseudo Riemannian manifolds. We derived fully generic equations of motions for three manifolds [5], but purposefully omitted lengthy discussion about motion of two-dimensional surfaces, which is a topic for this paper.

Currently, significant progress on fluidic models of membrane dynamics has already been made. The role of geometric constraints in self-assembly have been elucidated by linking together thermodynamics, interaction free energies, and geometry [6, 7]. The Helfrich formalism provides the foundation for a purely differential geometric approach whereby the membrane surface potential energy density is considered as a functional of the static curvature [4], see also review papers [8–10]. The model has been improved by adding force and torque balance equations [11, 12]. Specific dynamical equations accounting for bending as well as electrodynamic effects have also been reported [13–15]. Furthermore, active membrane theories, by including external forces, have extended the understanding of passive membranes [16–19].

Among the remarkable aspects of the fluid lipid membranes deduced from the large body of theoretical works [8–10], is that the physical behavior of a membrane on the length scale of its own thickness, in first approximation, can be expressed by a purely geometric Hamiltonian [4, 20, 21]. Associated Euler-Lagrange equations [22, 23], so called “shape equations,” are fourth order partial nonlinear differential equations, and finding a general analytical solution is typically difficult, even though it has been analytically [24] and numerically solved for some specific [25–32] and general cases [33, 34].

In fluid dynamics, material particles can be treated as a vertex of geometric figure and virtual layers as surfaces and equations of motion for such surfaces can be searched. We refer to the formalism as differentially variational surfaces (DVS) (or DVS formalism) [5].

In this paper, we propose different approach to the “shape choreography” problem. We use DVS formalism, tensor calculus of moving surfaces and the first law of thermodynamics to derive the final equation for the closed 2D surface dynamics (later on referred as a surface) and to solve it analytically for the equilibrium case. In other words, we derive generic equations of motions for closed two-dimensional surfaces and without any a priori symmetric assumptions, we show that constant mean curvature shapes are equilibrium solutions. In contrast to the Young-Laplace law these solutions, are universally correct descriptions of capillary surfaces as well as molecular surfaces. In addition, our equations of motions (20–25) are generic and exact. It advances our understanding of fluid dynamics because generalizes ideal magneto-hydrodynamic and Naiver-Stokes equations [5] and in contrast to Navier-Stokes, as we demonstrate in this paper, are trivially solvable for equilibrium shapes. To demonstrate the validity of these equations and their analytical solutions we apply them to micelles. Within our formalism it becomes simple task to show micelles lamellar, cylindrical, spherical shapes, and assert their optimal spherical radius.

For clarity, we shall give brief description of micelles and their structures. A micelle consists of monolayer of lipid molecules containing hydrophilic head and hydrophobic tail. These amphiphilic molecules, in aqueous environment, aggregate spontaneously into a monomolecular layer held together due to a hydrophobic effect [35, 36] (see also 5, 37–40) by weak non-covalent forces [41]. They form flexible, diverse surfaces with different topology, but remarkably in thermodynamic equilibrium conditions they are spherical, lamellar (plane), or cylindrical in shape.

 2. Methods

In the section we provide basics of tensor calculus for moving surfaces and summarize the theorems we used directly or indirectly to derive equations for two-dimensional surface dynamics. Differential geometry preliminaries we used here are available in tensor calculus textbook [42] and in our work [5].

2.1. Basics of Differential Geometry

Let Si (i = 1, 2) be the coordinates of the moving surface S and Xα be the coordinates in the ambient Euclidean space (Figure 1). The surface is smooth enough for sufficient differentiability in space-time and it is expressed as Xα = Xα(t, Si) in the ambient coordinates. The surface equation for the position vector R→ is1:

R→=R→(Xα)=R→(t,Si)

Graphical illustration of the arbitrary surface and its' local tangent plane. S→1,S→2,N→ are local tangent plane base vectors and local surface normal, respectively. X→1,X→2,X→3 are arbitrary base vectors of the ambient Euclidean space and R→=R→(X)=R→(t,S) is radius vector of the point. V→ is arbitrary surface velocity and C, V1, V2 display projection of the velocity to the N→,S→1,S→2 directions, respectively.

A dot product of the covariant bases X→α=∂αR→ is a metric tensor:

Xαβ=X→αX→β

where ∂α=∂/∂Xα. The matrix inverse of the covariant metric forms the contravariant one: XαβXβγ=δγα (the Kronecker delta δγα). In the ambient Euclidean space, due to mutual linear independence of the base vectors, the Christoffel symbols Γβγα=X→α·∂βX→γ are zeros and, therefore, covariant/contravariant derivatives become partial ones ∇α = ∂α.

The surface base vectors Si→=∂iR→ and metric tensors Sij,Sij are similarly defined

Sij=S→i·S→j

where ∂i=∂/∂Si and SabSbc=δac. The surface Christoffel symbols given by Γjki=S→i·∂jS→k are not zeros and covariant/contravariant derivatives differ from partial ones ∇i ≠ ∂i. The space/surface Christoffel symbols form the basic concept for defining curvilinear derivatives on mixed space/surface tensors:

∇iTβkαj=∂iTβkαj+XiγΓγναTβkνj-XiγΓγβμTμkαj+ΓimjTβkαm-ΓikmTβmαj

where Xiγ is the shift tensor, so that

S→i=XiαX→α and Sij=S→i·S→j=XiαX→αXjβX→β=XiαXjβXαβ. Note that in (4) the Christoffel symbols with Greek indexes are zeros.

Using (2, 4), one may directly prove metrilinic property of the surface metric tensor ∇iSmn = 0, from where follows S→m·∇iS→n=0, meaning that S→m⊥∇iS→n are orthogonal vectors, therefore:

∇iS→j=N→Bij

where N→ is an unit surface normal and Bij is the symmetric curvature tensor. The trace of the mixed curvature tensor is the mean curvature Bii and the determinant is the Gaussian curvature. A particular case:

Bii=λ

where λ is some non-zero constant is a sphere if the surface is closed. According to (5, 6), finding the curvature tensor non-directly implies identification of the surface.

 2.2. Basics of Tensor Calculus for Moving Surfaces

All Equations written above are generally true for moving surfaces. We now turn to a brief review of definitions of coordinate velocity Vα, interface velocity C (which is the same as normal velocity), tangent velocity Vi (Figure 1), time derivatives of the surface tensors and the space/surface integrals. The original definitions can be found in textbooks [42, 43].

Suppose that Vα is the coordinate velocity defined as:

Vα=∂Xα∂t

Then, given that the position vector R→ (1) is tracking the coordinate particle Si the surface velocity can be written as:

V→=∂R→(t,Si)∂t=∂R→∂Xα∂Xα(t,Si)∂t=VαX→α

Therefore, Vα is ambient component of the surface velocity V→. Taking into account (8) the surface normal velocity C can be represented as:

C=V→·N→=VαX→αNβX→β=VαNβδβα=VαNα

C is called the interface velocity, is invariant and its sign depends on the choice of the surface normal. The surface velocity projection on the tangent plane (Figure 1) is Vi the tangential velocity:

Vi=VαXαi

Taking (9, 10) into account one may write the surface velocity as V→=CN→+ViS→i.

Geometric interpretation of the interface velocity C and of the curvilinear time derivative ∇˙ applied to invariant field F. A is arbitrary chosen point so that it lays on F(St) ∈ St curve and B is its' corresponding point on the St+Δt surface. P is the point where St surface normal, applied on the point A, intersects the surface St+Δt. By the geometric construction, for small enough Δt → 0, ∠APB → π/2, AB→≈V→Δt, and AP≈V→N→Δt. On other hand, by the same geometric construction the field F in the point B can be estimated as F(B) ≈ F(A)+Δt∂F/∂t, while from viewpoint of the St+Δt surface the F(B) value can be estimated as F(P)+ΔtVi∇iF, where ∇iF shows rate of change in F ∈ St+Δt and along the BP ≈ ΔtVi.

The surface velocity in general and the interface velocity in particular has clear geometric interpretation. Suppose that the surfaces at the two nearby moments of time t, t+Δt are St, St+Δt. Then, for any arbitrarily chosen A ∈ St and the corresponding point B ∈ St+Δt the vector AB→ is approximately the surface velocity times increment of time: AB→≈V→Δt (Figure 2). If the P is the point where the surface normal N→∈St intersects the surface St+Δt, then for infinitely small Δt, ∠APB → π/2 and AP→V→·N→Δt, therefore:

C=limΔt→0APΔt

The interface velocity C is the instantaneous velocity of the interface in the normal direction.

 2.3. Invariant Time Differentiation

Among the key definitions in calculus for moving surfaces, perhaps one of the most important is the invariant time derivative ∇˙. As we have already stated, invariant time derivative is already well defined in the literature [42, 43]. In this paragraph, we just give geometrically intuitive definition.

Let take smooth functional F defined on the surface F ∈ St. The invariant time derivative has to be defined so that: (1) it is free from choice of a reference frame (invariance) and (2) properly takes into account deformations along the normal. We explain importance of normal deformations in integration subsection. This can be achieved by following geometric construction: for an arbitrary chosen point A from the St surface A ∈ St let us find points B ∈ St+Δt and P so that the B is the intersection of the surface velocity to the St+Δt and the P is intersection of the surface normal to the St+Δt (Figure 2). Then, the invariant time derivative ∇˙ can be defined as:

∇˙F=limΔt→0F(P)-F(A)Δt

(12) is geometric, therefore must be an invariant. From the geometric construction (12) one can estimate value of F in point B, so that

F(B)≈F(A)+Δt∂F∂t

B, P ∈ St+Δt are nearby points, therefore according to the definition of covariant derivative F(B) can be estimated as:

F(B)≈F(P)+ΔtVi∇iF

since ∇iF shows rate of change in F ∈ St+Δt and Δt·Vi indicates BP. Determining F(A), F(P) values from (13, 14) and putting it in (12), gives

∇˙F=∂F∂t-Vi∇iF

The definition (15) extends for space/surface mixed tensors by the following formula:

∇˙Tβjαi=∂Tβjαi∂t-Vk∇kTβjαi+VγΓγμαTβjμi-VγΓγβμTμjαi+Γ˙kiTβjαk-Γ˙jkTβkαi

where analog of the Christoffel symbol Γ˙ji for moving surfaces is defined by the formula Γ˙ji=∇jVi-CBji. The derivative (16) satisfies: commutativity with the contraction, product and chain rules, metrilinic properties against the ambient metric and does not commute with the surface derivative [42]. According to (12) the invariant time derivative vanishes when applied to time independent scalars.

 2.4. Time Differentiation of Integrals

The calculus of moving surfaces are effective due to two fundamental theorems about taking time derivatives of space and surface integrals [42]. In evaluation of the least action principle of the Lagrangian there is a central role for time differentiation of the surface and the space integrals, from where the geometry dependence is rigorously clarified.

Let the scalar field F = F(t, S) be defined on an Euclidean domain Ω which has the surface boundary S and the surface evolves with the interface velocity C while volume moves (expands or shrinks). Analogically suppose that the closed surface evolves and F = F(t, S) is again some scalar functional defined on the evolving surface S. Then, there are theorems for taking derivative of space and surface integrals from F integrand:

ddt∫ΩFdΩ=∫Ω∂F∂tdΩ+∫SCFdS ddt∫SFdS=∫S∇˙FdS-∫SCFBiidS

On the right hand side of (17, 18) the first terms indicate how the scalar field changes and the second terms show how geometry evolves. There are rigorous mathematical proofs of these formulas in the tensor calculus textbooks; we do not reproduce them here. Instead, we intuitively show why only interface velocity has to be taken into account. Rigorous mathematical proof follows from fundamental theorem of calculus

ddt∫ab(t)F(t,x)dx=∫ab(t)∂F(t,x)∂tdx+b′(t)F(t,b(t))

In the case of volume integral or surface integral it can be shown that b′(t) is replaced by interface velocity C.

Intuitive explanation is pretty simple. Propose there is no interface velocity then closed surface velocity only has tangent component. For each given time tangent velocity (if there is no interface velocity) translates each point to its neighboring point and therefore, does not add new area to the closed surface (or new volume to the closed space, or new length to the closed curve). As so, tangential velocity just induces rotational movement (or uniform translational motion) of the object and can be excluded from additive terms in the integration. Perhaps, it is easier to understand this statement for one dimensional motion. Let's assume that material point is moving along some trajectory (some closed curve or loop), then, in each point, the velocity of the material point is tangential to the curve. Now one can translate this motion into the motion of the closed curve where the loop has only tangential velocity. In this aspect, the embedded loop only rotates (uniformly translates in the plane) without changing the length locally, therefore tangential velocity of the curve does not add new length to the curve (same is true for open curve with fixed ends).

 3. General Equations of Surface Motions

Fully non-restrained and exact equations for moving three-dimensional surfaces in electromagnetic field, when the interaction with an ambient environment is ignored, reads

∇˙ρ+∇i(ρVi)=ρCBii ∂α(ρVα(∇˙C+2Vi∇iC+ViVjBij)             −Vα(14μ0FμνFμν+AμJμ))=fa ∫SρVi(∇˙Vi+Vj∇jVi-C∇iC-CVjBji)dS=∫ΩfiaidΩ

where ρ is the surface mass density, Vα, Vi are coordinate and tangential components of the surface velocity, C is interface velocity, α = 0, 1, 2, 3 for Minkowski four-dimensional space-time ambient space, i = 0, 1, 2 for pseudo-Riemannian manifold (surface), Bij is the surface curvature tensor, Fμν is electromagnetic tensor, Fα=Jα-∂βFβα, Jα is α component of J→=(Jα) four current, f, fi are normal and tangential components of F→=(Fα), a, ai are the normal and the tangential components of the ∂A→/∂t (here A→=(Aα) is four vector), S, Ω stand for the surface and the space integrals, respectively. The exact derivation of (20–22) is given in our work [5], we do not reproduce derivation of this set in this paper, rather just mention that first one is the consequence of mass conservation, second and third equations come from minimum action principle of a Lagrangian and imply motion in normal direction (21) and in tangent direction (22). For two dimensional surface dynamics, Minkowskian space becomes Euclidean, so that α = 1, 2, 3 and the surface is two-dimensional Riemannian manifold i = 1, 2. So that, after modeling the potential energy as a negative volume integral of the internal pressure and inclusion interaction with an environment, (20–22) further simplifies as

∇˙ρ+∇i(ρVi)=ρCBii ∂α(ρVα(∇˙C+2Vi∇iC+ViVjBij)            +Vα(P++Π))=−Vα∂α(P++Π) ρVi(∇˙Vi+Vj∇jVi-C∇iC-CVjBji)=0

where P+, Π are internal hydrodynamic and osmotic pressures, respectively. Derivation of (20–22) can be found in Svintradze [5]. We derive (23-25) in Appendix section. It is noteworthy that from the last equations set only the second equation (24) differs from the dynamic fluid film equations [42, 44]

ρ(∇˙C+2Vi∇iC+ViVjBij)=σBii

where σ is the surface tension. Equation (26) is only valid when the surface can be described with time invariable surface tension [42, 44], meaning that the surface is homogeneous and the surface tension is constant, while (24) does not have that restriction. Using (26) in (24) and taking into account that in equilibrium processes internal pressure is the same as external pressure, one gets exactly the same equation of motion in normal direction (39) as we get from using the first law of thermodynamics (see below).

∂α(σVαBii+(P++Π)Vα)=-(∂αP++∂αΠ)Vα

It is worth of mentioning that (23–25) also follows from (20–22) if one applies same formalism as it is given in (30-32). Indeed, the space is 3D Euclidean so that α = 1, 2, 3 and the surface is 2D Riemannian (i = 1, 2) for relatively slowly moving surfaces and the potential energy could be written as:

U=∫Ω(14μ0FμνFμν+AμJμ)=∫Ω(-ϵ02E2+1μ0B2-qφ+A→J→)dΩ

where E→,B→ are electric and magnetic fields and q,φ,A→,J→ are charge density, electric potential, magnetic vector potential and current density vector respectively. Using (30–32) formalism into account, we find

dU=-(P++Π)dΩ=(-ϵ02E2+1μ0B2-qφ+A→J→)dΩ

Taking into account (29) and that the pressure comes from the normal force applied to the surface, we find fa=-Vα∂α(P++Π) and in tangent direction fiai=0, then (20–22) becomes (23–25). Electromagnetic potential energy can be generalized if one takes into account environment, which enters in energy terms as bound and free charges and electric/magnetic fields are replaced by polarization and magnetization vectors [5].

 4. Results and Discussion 4.1. General Assumptions

In this section we apply basics of thermodynamics and fundamental theorems of calculus of moving surfaces to demonstrate shortest derivation of the equation, describing motion of homogeneous, closed two dimensional surface with time invariable surface tension at normal direction (27). We consider the system consisted of aqueous media with the formed closed surface in it (Figure 3). The system is isolated with constant temperature and there is no absorbed or dissipated heat on the surface; in other words, a process is adiabatic. According to the first law of thermodynamic, as far as there is no dissipated or absorbed heat, the change of the internal energy of the surface must be

dE=δW

where δW is infinitesimal work done on the subsystem and dE is infinitesimal change of the internal energy. Because the temperature of the system is constant, the differential of the subsystems' internal energy can be remodeled as

dE=dU

where U is the total potential energy of the surface. By the definition the elementary work done on the subsystem is

δW=(P-+Π)dΩ

where, P−, Π are external hydrodynamic and osmotic pressures applied on the surface by the surroundings correspondingly and Ω is the volume that surface encloses with boundary of S surface area. Let's propose that the surface is homogeneous (i.e., material particles are homogeneously distributed on the surface) so that the total potential energy is integration of the potential energy per unit area over the surface, then

dU=σdS

where σ is the potential energy per unit area and is called surface tension in the paper. As far as we discuss simplest case of the system consisted of aqueous medium and single closed surface, we can suggest that the surface tension is not time variable. Using (30–33) after few lines of algebra, we fined

∫SσdS=∫Ω(P-+Π)dΩ

Graphical illustration of the isolated system containing aqueous solution. Water molecules are represented as red and white sticks. The system boundary is shown as white faces with black edges. The subsystem-micelle is closed surface, blue blob in the center of the system.

Assuming the surface is moving so that (34) stays valid for any time variations, then time differentiation of the left side must be equal to the time differentiation of the right integral. As far as on the right hand side we have space integral, time differentiation can be taken into the integral, using general theorems for differentiation of space and surface integrals (17–18), so that the integration theorem for the space integral which takes into account the volume motion holds, therefore:

ddt∫Ω(P-+Π)dΩ=∫Ω(∂αP-+∂αΠ)∂Xα∂tdΩ+∫SC(P-+Π)dS

To calculate a time derivative of the surface integral we have to take into account the theorem about time differentiation of the surface integral (18), from which follows that for the time invariable surface tension

ddt∫SσdS=∫S-σCBiidS

Where C=VαNα is interface velocity, Nα is α component of the surface normal and V = ∂Xα/∂t is coordinate velocity, Xα is general coordinate and Bii is the trace of the mixed curvature tensor generally known as mean curvature. After few lines of algebra putting (34–36) together, we find

∫S(σCBii+C(P-+Π))dS=-∫Ω(∂αP-+∂αΠ)VαdΩ

Generalized Gauss theorem converts the surface integral of the left hand side of (37) into space integral, so that

∫SNαVα(σBii+P−+Π)dS=∫Ω∂α(σVαBii                                      +(P−+Π)Vα)dΩ

Combination of (37) and (38) immediately gives equation of motion for surface in normal direction

∂α(σVαBii+(P-+Π)Vα)=-(∂αP-+∂αΠ)Vα

For equilibrium processes internal and external pressures are identical P− = P+, so that (39) becomes identical to the equation of motion in normal direction observed from master equations (23–27). Also, we should note that (39) is only valid for motion of the homogeneous surfaces with time invariable surface tension at the normal direction, therefore, it does not display any deformation in tangent directions. Equation (39) further simplifies when the surface comes in equilibrium with the solvent where divergence of the surface velocity ∂αVα (stationary interface) along with ∂P/∂t (where P = P−+Π) vanishes, then the solution to (39), taking into account the condition (35), becomes

Bii=-Pσ

The result (40) shows that the solution is constant mean curvatures (CMC) surfaces. Such CMC are rare and can be many if one relaxes the condition we restricted to the system. We consider isolated system where the surface is closed subsystem, these two preconditions mathematically mean that the surface we discuss is compact embedded surface in ℝ3. According to A. D. Alexandrov uniqueness theorem, if a compact surface embedded in ℝ3 has constant non-zero mean curvature then it is a sphere [45]. Correspondingly the solution (40) is a sphere (as far as we have compact two-manifold in the Euclidean space). When

Pσ≠0

the surface is spheroid (or a cylinder if one relaxes compactness restriction making the cylinder infinitely long) and becomes plane (again when compactness argument is relaxed) or other zero mean curvature shape when compactness argument is not relaxed but contour of the surface remains fixed. This surprisingly simple and elegant derivation explains all the shapes surfaces can adopt in aqueous solution at equilibrium conditions2. If the compactness condition is relaxed then (40) predicts that in addition to cylinder and plane all other CMC surfaces are also equilibrium shapes for moving surfaces. Taking into account that the surface tension in general can be a function of many variables, such as Gaussian curvature, bending rigidity, spontaneous curvature, molecules concentration, geometry of surfactant molecules and etc., then (39) may predict possible deformations of differently shaped surfaces and their wide range of static shapes. In fact, if considered that the surface tension, which is defined as potential energy per unit area, can be a function of mean curvature σ=σ(Bii), then Taylor expansion of σ(Bii) naturally rises all additional terms. These generalizations and temperature fluctuations can be included in the equations, but it is not scope of this paper and should be addressed separately. One may even propose σ as time independent the Helfrich Hamiltonian and then (39, 40) will become equation of static shapes for homogeneous surfaces with time invariable surface tension. In fact we clarify this statement in the following subsection.

 4.2. Relevance to the Shape Equation

Here we give brief comparison of our equations of motions to the Helfrich formalism and show that the shape equation is a stationary approximate solution to the equation of the normal motion (39).

Let us propose that the system is overdamped, therefore variations in kinetic energy is much smaller than variations in the potential energy. Throughout of our equations we model potential energy as

U=∫ΩP-dΩ

Assume that the Hamiltonian of the system H=∫SρV2/2dS+∫ΩP-dΩ is modeled as the free energy on the surface ∫S(σ0+k/2(Bii-K0)2+kKG)dS, where σ0 is classically defined constant surface tension, k is bending rigidity constant, K0 is spontaneous curvature and KG is the Gaussian curvature. This interpretation of the free energy is called spontaneous curvature model [8, 31]. As far as the system is overdamped, the variation of the kinetic energy is infinitely small, therefore in the first approximation δ∫SρV2/2dS≈0 and we have

H=∫SρV22dS+∫ΩP-dΩ=∫S(σ0+k/2(Bii-K0)2+kKG)dSδ∫ΩP-dΩ=δ∫S(σ0+k/2(Bii-K0)2+kKG)dS

For evaluation of the Hamiltonian variation, we note that according to the time differentiation of the space integrals (17) δ∫ΩP-dΩ=∫Ω∂P-/∂tdΩ+∫SCP-dS. On the other hand, the variation of the Helfrich free energy, which can be restricted to the normal variations [8], is

δ∫S(σ0+k/2(Bii−K0)2+kKG)dS=∫Sδ(σ0+k/2(Bii−K0)2+kKG)δNδNδtdS=∫S(σ0Bii−k(∇a∇aBii−12(Bii−K0)[(Bii−K0)Bii−2(Bii)2+4KG]))CdS

where N is the normal coordinate and thus the interface velocity becomes C = ∂N/∂t. For stationary or near to stationary shapes ∂P−/∂t ≈ 0, taking into account (42, 43) and the variation of the ∫ΩP-dΩ, we find

∫S(σ0Bii−k(∇a∇aBii−12(Bii−K0)[(Bii−K0)Bii−2(Bii)2+4KG]))CdS=∫SCP−dS

The last Equation (44) has to hold for any interface velocity C, therefore

σ0Bii−k(∇a∇aBii−12(Bii−K0)[(Bii−K0)Bii−2(Bii)2+4KG])=P−

(45) is exactly the shape Equation [4, 8, 22, 23, 31]. Note that the shape Equation (45) can be obtained from the equation of motion in the normal direction (39) for near to stationary case. Indeed, according to (39) taking into account that the surface pressure can be modeled from (44, 45), we find

∂α(σBiiVα+(σ0Bii−k(∇a∇aBii−12(Bii−K0)                   [(Bii−K0)Bii−2(Bii)2+4KG]))Vα)=−∂P−∂t

Taking into account that for stationary shapes ∂P-/∂t=0,∂αVα=0, the last equation has trivial solution

σ0,Bii−k(∇a∇aBii−12(Bii−K0)[(Bii−K0)Bii−2(Bii)2+4KG])=P0

where σ0,,P0 are some constant surface tension and pressure respectively. These calculations formally establish that the shape equation is an approximate solution of our equations of motions constrained by the conditions: time invariable surface tension and stationary interfaces; and is valid only for overdamped systems. Note that the analytic solution (40) is exact and only valid for the surfaces in equilibrium with environment, while the shape equation is the approximate solution for stationary surfaces. It is no surprise that the stationary shapes are geometrically richer than equilibrium shapes. Equilibrium shapes can only adopt constant mean curvature.

 4.3. Physical Application, Micelle

We can put Equation (39) and its solution (40) under the test for homogeneous micellar surface equilibrated with the aqueous solution. Based on (40) we can calculate minimal value of a micelle radius. The value of the trace of the mixed curvature tensor for a sphere is

Bii=-2R

where R is a radius.

Let's calculate value of the surface pressure when the micelle still can exist. Lipids in a micelle are confined in the surface by hydrophobic interactions with average energy in the range of hydrogen bonding. As far as values of hydrogen bonding energy are somewhat uncertain in the literature, by the first approximation we take average energy for the hydrogen bonding energy interval and assign it to the lipid molecule. Low boundary of the interval (minimum energy) for XH···Y hydrogen bond is about 1 kJ per mol (CH···C unit) and high boundary is about 161 kJ per mol (FH···F unit), the low and high values are taken according to references [46, 47]. Therefore, average energy is about (1+161)/2 = 81kJ/mol ≈ 13·10−20J. To estimate hydrogen bonding energy per molecule with the undefined shape (lipid molecule) in the first approximation is to assign average energy to it and consider the spherical shape with the gyration radius. Of course it is low level approximation, but even such rough calculations produce reasonable results. After all these rough estimations the pressure to move one lipid from the surface, in order to induce critical deformations of the surface, is about average energy per the average volume of the lipid molecule

P≈3·13·10-204πrG3≈3.1·107N/m2

where 4πrG3/3 is the estimated volume of a lipid molecule considered as sphere with the gyration radius rG ≈ 1nm. On the other hand, the surface tension of a fluid monolayer at optimal packing of the lipids is about σ ≈ 3·10−2N/m [7, 48, 49], using these and (48, 49) in (40) the estimated micelle radius is

R≈2·3·10-23.1·107=19.3±0.1Å

These calculations put the minimum radius in nanometer scale and is in very good agreement with experimental as well as computational frameworks [50, 51]. To further validate the (50) result, we ran a CHARMM based Micelle Builder simulation [52, 53] for 100 phospholipid molecules (DHPC lipids). The simulation result (Figure 4) generated a spherical micelle with diameter 38.5±0.1Å. These calculations indeed indicate that even such rough estimations produce reasonable accuracy.

Simulated three dimensional coordinates of the micelle in aqueous solution display sphere with diameter 38.5Å. (Left) dihexanoylphosphatidylcholine (DHPC phospholipids) are modeled as orange balls. (Right) Gaussian mapping at contour resolution 8Å of the micelle shows spherical structure.

To get more convincing estimations it is necessary to take into account that neither lipids are spherical nor hydrophobic interactions per lipid are average energy of single hydrogen bond. In the second approximation lipids are no longer undefined spheres, but have well defined surfactant geometry. The Hydrophobic energy is no longer average energy of single hydrogen bond, but is 1 kJ per mol per −CH2− unit. In all atom simulations we used dihexanoylphosphatidylcholine (DHPC) lipid molecule having 12−CH2− units (Figure 5) per hydrophobic tail, so hydrophobic energy is about 12kJ/mol ≈ 1.99·10−20J. Accurate calculation of the lipid molecule volume using cavity, channel and cleft volume calculator [54], gives the volume estimation of about 894Å3. Using this value, one gets

P≈1.99·10-200.894·10-27≈2.22·107N/m2

On other hand, using the same surface tension of a fluid monolayer at the optimal packing of the lipids, one gets R = 27 ± 0.1Å. All atom simulation also generates spherical structure with diameter 54 ± 0.1Å (Figure 5). There is still some uncertainty in this estimation because we assigned 1 kJ/mol energy per −CH2− unit and we based on references data [46, 47], while in other literature it is mentioned that the hydrophobic interactions are about 4 kJ/mol per −CH2− unit [55]. In our opinion, this discrepancy can be resolved if one calculates hydrophobic energy based on the potential energy

U=-∫Ωϵ02ECH22dΩ

where E→CH2 is electric field per −CH2−, ϵ0 is dielectric constant in the vacuum and Ω stands for the volume of the lipid molecules. Equation (52) directly emerges from FμνFμν term written in the equations of motion (20–22 and 28–29). For electrostatics

U=∫Ω14μ0FμνFμνdΩ=-∫Ωϵ02ECH22dΩ

so one should go to the scrutiny of calculating electric field for each −CH2− units, then take a sum of the electric field and square it (we are not going to do it in this paper). Also, one may ask why the hydrophilic interaction energy is not taken into account in these calculations. Hydrophilic head of the lipid molecule is in contact with water molecules so there is no work needed to drag it in aqueous solution from the lipids layer. Therefore, hydrophilic interaction energy can be neglected. The most work goes on overcoming hydrophobic interactions between lipid tails.

All atom simulation of DHPC micelle. (A) The figure shows a geometry of the DHPC surfactant molecule used in simulation and gives parametric description of volume, surface area, sphericity, and effective radius. (B) Indicates atomistic simulation result contoured by Gaussian map and the diameter of the micelle, measured by PyMol. The diameter of the simulated micelle appears to be 54.0Å with the uncertainty of the measurement 0.1Å.

 5. Conclusion

We have presented a framework for the analysis of two dimensional surface dynamics (identified as micelle in the text) using first law of thermodynamics and calculus of moving surfaces. In final equations of normal motion (39,27) we assume that a surface is homogeneous and has time invariable surface tension. However, the general equations (23–25) do not have these constrains and indicate arbitrary motion along normal deformation, as well as into tangent directions, but are analytically more complex. The solution to the normal equations of motion in equilibrium conditions are surprisingly simple and display all possible equilibrium shapes. We applied the formalism to estimate micelle optimal radius and compared estimations to all atom simulations. Even for low-level approximations, we found remarkable agreement between theoretically calculated radius and one obtained from atomistic simulations and from experiments. One can readily apply the theory to any closed surfaces; such are vesicles, membranes, water droplets or soap films.

As a final remark, even though the analytic solution (40) looks like generalized Young-Laplace law, the difference is obvious. Bii is a trace of mixed curvature tensor, known as mean curvature, and when the mean curvature is constant, it defines whole class of constant mean curvature (CMC) surfaces. Generalized the Young-Laplace law is a priori formulated for spherical morphologies and therefore in some particular cases can be obtained from (40) constant mean curvature shapes. The condition for holding the particular case is a compactness. However, the compactness argument can be relaxed in our derivation if the considered system is set to be much larger than the subsystem. Therefore, the solution (40) effectively predicts formation of all CMC surfaces while Young-Laplace law is correct for spherical structures alone. Also in derivation of Young-Laplace relation one of cornerstone idea is suggestion of spherical symmetries, while our derivation is free of symmetries and explains why CMC surfaces are such abundant shapes in nature, observable even on molecular levels. In fact, according to the results, any homogeneous closed surface with time invariable surface tension adopts CMC shape when it comes in equilibrium with environment.

 Author Contributions

The author confirms being the sole contributor of this work and has approved it for publication.

Conflict of Interest Statement

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

We thank Dr. Frank Julicher (MPIPKS) and Dr. Erwin Frey (LMU), for stimulating discussions. The paper in it's current form was initiated at Aspen Center for Physics, which is supported by National Science Foundation grant PHY-1607611 and was partially supported by a grant from the Simons Foundation.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2018.00136/full#supplementary-material

 References Cahalan MD Parker I. Choreography of cell motility and interaction dynamics imaged by two-photon microscopy in lymphoid organs. Annu Rev Immunol. (2008) 26:585–626. 10.1146/annurev.immunol.24.021605.09062018173372 Pal SK Peon J Bagchi B Zewail AH. Biological water:â femtosecond dynamics of macromolecular hydration. J Phys Chem B (2002) 106:12376–95. 10.1021/jp0213506 Wan C Fiebig T Schiemann O Barton JK Zewail AH. Femtosecond direct observation of charge transfer between bases in DNA. Proc Natl Acad Sci USA. (2000) 97:14052–5. 10.1073/pnas.25048329711106376 Helfrich W. Elastic properties of lipid bilayers: theory and possible experiments. Z Naturforsch C (1973) 28:693–703. 4273690 Svintradze DV. Moving manifolds in electromagnetic fields. Front Phys. (2017) 5:37. 10.3389/fphy.2017.00037 Israelachvili JN Mitchell DJ Ninham BW. Theory of self-assembly of hydrocarbon amphiphiles into micelles and bilayers. J Chem Soc Faraday Trans 2 (1976) 72:1525–68. Israelachvili JN Mitchell DJ Ninham BW. Theory of self-assembly of lipid bilayers and vesicles. Biochim Biophys Acta (1977) 470:185–201. 911827 Deserno M. Fluid lipid membranes: from differential geometry to curvature stresses. Chem Phys Lipids (2015) 185:11–45. 10.1016/j.chemphyslip.2014.05.00124835737 Lipowsky R Sackmann E. Structure and Dynamics of Membranes. Amsterdam: Elsevier (1995). Seifert U. Configurations of fluid membranes and vesicles. Adv Phys. (1997) 46:13–137. Capovilla R Guven J. Stresses in lipid membranes. J Phys A (2002) 35:6233. 10.1088/0305-4470/35/30/302 Fournier JB. On the stress and torque tensors in fluid membranes. Soft Matter. (2007) 3:883–8. 10.1039/B701952A Scriven LE. Dynamics of a fluid interface equation of motion for newtonian surface fluids. Chem Eng Sci. (1960) 12:98–108. Arroyo M DeSimone A. Relaxation dynamics of fluid membranes. Phys Rev E (2009) 79:031915. 10.1103/PhysRevE.79.03191519391979 Gao LT Feng XQ Yin YJ Gao H. An electromechanical liquid crystal model of vesicles. J Mech Phys Solids (2008) 56:2844–62. 10.1016/j.jmps.2008.04.006 Salbreux G Jülicher F. Mechanics of active surfaces. Phys Rev E (2017) 96:032404. 10.1103/PhysRevE.96.03240429346890 Ramaswamy S Toner J Prost J. Nonequilibrium fluctuations, traveling waves, and instabilities in active membranes. Phys Rev Lett. (2000) 84:3494–7. 10.1103/PhysRevLett.84.349411019123 Chen HY. Internal states of active inclusions and the dynamics of an active membrane. Phys Rev Lett. (2004) 92:168101. 10.1103/PhysRevLett.92.16810115169263 Gov N. Membrane undulations driven by force fluctuations of active proteins. Phys Rev Lett. (2004) 93:268104. 10.1103/PhysRevLett.93.26810415698026 Canham PB. The minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell. J Theor Biol. (1970) 26:61–81. 5411112 Evans EA. Bending resistance and chemically induced moments in membrane bilayers. Biophys J. (1974) 14:923–31. 10.1016/S0006-3495(74)85959-X4429770 Zhong-Can OY Helfrich W. Instability and deformation of a spherical vesicle by pressure. Phys Rev Lett. (1987) 59:2486–8. 10.1103/PhysRevLett.59.248610035563 Zhong-Can OY Helfrich W. Bending energy of vesicle membranes: general expressions for the first, second, and third variation of the shape energy and applications to spheres and cylinders. Phys Rev A (1989) 39:5280–8. 9901091 Mladenov IAM Djondjorov PA Hadzhilazova MT Vassilev VM. Equilibrium configurations of lipid bilayer membranes and carbon nanostructures. Commun Theor Phys. (2013) 59:213–28. 10.1088/0253-6102/59/2/14 Svetina S Žekš B. Membrane bending energy and shape determination of phospholipid vesicles and red blood cells. Eur Biophys J. (1989) 17:101–11. 2766997 Seifert U Lipowsky R. Adhesion of vesicles. Phys Rev A (1990) 42:4768–71. 9904586 Lipowsky R. The conformation of membranes. Nature (1991) 349:475–81. 1992351 Seifert U Berndl K Lipowsky R. Shape transformations of vesicles: phase diagram for spontaneous-curvature and bilayer-coupling models. Phys Rev A (1991) 44:1182–202. 9906067 Jülicher F Lipowsky R. Domain-induced budding of vesicles. Phys Rev Lett. (1993) 70:2964–7. 10053698 Jülicher F Lipowsky R. Shape transformations of vesicles with intramembrane domains. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics (1996) 53:2670–83. 9964554 Jülicher F Seifert U. Shape equations for axisymmetric vesicles: a clarification. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics (1994) 49:4728–31. 9961774 Miao L Seifert U Wortis M Döbereiner HG. Budding transitions of fluid-bilayer vesicles: the effect of area-difference elasticity. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics (1994) 49:5389-407. 9961866 Heinrich V Svetina S Žekš B. Nonaxisymmetric vesicle shapes in a generalized bilayer-couple model and the transition between oblate and prolate axisymmetric shapes. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics (1993) 48:3112–23. 9960950 Kralj-Iglič V Svetina S Žekš B. The existence of non-axisymmetric bilayer vesicle shapes predicted by the bilayer couple model. Eur Biophys J. (1993) 22:97–103. 8359147 Chandler D. Interfaces and the driving force of hydrophobic assembly. Nature (2005) 437:640–7. 10.1038/nature0416216193038 Leikin S Parsegian VA Rau DC Rand RP. Hydration forces. Annu Rev Phys Chem. (1993) 44:369–95. 8257560 Svintradze DV. Hydrophobic and hydrophilic interactions. Biophys J. (2010) 98:43a–4a. 10.1016/j.bpj.2009.12.250 Svintradze DV. Moving macromolecular surfaces under hydrophobic/hydrophilic stress. Biophys J. (2015) 108:512a. 10.1016/j.bpj.2014.11.2807 Svintradze DV. Cell motility and growth factors according to differentially variational surfaces. Biophys J. (2016) 110:623a. 10.1016/j.bpj.2015.11.3342 Svintradze DV. Geometric diversity of living organisms and viruses. Biophys J. (2017) 112:309a. 10.1016/j.bpj.2016.11.1676 Tanford C. The Hydrophobic Effect Formation of Micelles and Biological Membranes. New York, NY: Wiley-Interscience (1973). Grinfeld P. Introduction to Tensor Analysis and the Calculus of Moving Surfaces. New York, NY: Springer (2010). Hadamard J. Mmoire sur le problme danalyse relatif lquilibre des plaques elastiques encastres. Oeuvres, Hermann, Tome 2; 1968. Grinfeld P. Exact nonlinear equations for fluid films and proper adaptations of conservation theorems from classical hydrodynamics. J Geom Symm Phys. (2009) 16:1–21. Alexandrov AD. Uniqueness theorem for surfaces in the large. Leningrad Univ 13, 19 (1958), 58. Am Math Soc Trans (Series 2). (1958) 21:412–6. Larson J McMahon T. Gas-phase bihalide and pseudobihalide ions. An ion cyclotron resonance determination of hydrogen bond energies in XHY-species (X, Y= F, Cl, Br, CN). Inorg Chem. (1984) 23:2029–33. Emsley J. Very strong hydrogen bonding. Chem Soc Rev. (1980) 9:91–124. Jähnig F. What is the surface tension of a lipid bilayer membrane? Biophys J. (1996) 71:1348–9. 8874009 Israelachvili JN. Intermolecular and Surface Forces: Revised Third Edition. Amsterdam: Academic Press (2011). Feller SE Zhang Y Pastor RW. Computer simulation of liquid/liquid interfaces. II. Surface tension-area dependence of a bilayer and monolayer. J Chem Phys. (1995) 103:10267–76. Egberts E Marrink SJ Berendsen HJ. Molecular dynamics simulation of a phospholipid membrane. Eur Biophys J. (1994) 22:423–36. 8149924 Jo S Kim T Iyer VG Im W. CHARMM-GUI: a web-based graphical user interface for CHARMM. J Comput Chem. (2008) 29:1859–65. 10.1002/jcc.2094518351591 Cheng X Jo S Lee HS Klauda JB Im W. CHARMM-GUI micelle builder for pure/mixed micelle and protein/micelle complex systems. J Chem Inf Model. (2013) 53:2171–80. 10.1021/ci400268423865552 Voss NR Gerstein M. 3V: cavity, channel and cleft volume calculator and extractor. Nucleic Acids Res. (2010) 38:W555–62. 10.1093/nar/gkq39520478824 Liu AL. Advances in Planar Lipids Bilayers and Liposomes, Vol. 4. Amsterdam: Academic Press; Elsevier (2011).

1Latin letters in indexes are related to the surface tensors and Greek letters are related to space tensors. Equations throughout the text are tensorial and repeated upper and lower indexes indicate Einstein summation convention, for instance XaYa=ΣXaYa.

2Even though we set environment as aqueous, it enters into equations as osmotic pressure term, which due to a generality of arguments can be anything. Therefore, as a medium one may pick any liquid or gas.