Since Minkowskian space-time does not follow Riemannian geometry, we need a small adjustment of definitions. For Minkowski space-time, which fits to pseudo-Riemannian geometry, we need definitions of arbitrary base pairs of ambient space, even though the definitions look exactly the same as for Riemannian geometry embedded in Euclidean ambient space [26–28]. The summarized relationships about Riemannian geometry embedded in Euclidean space are given in tensor calculus text books [26, 27].

Combination of three ordinary dimensions with the single time dimension forms a four-dimensional manifold and represents Minkowski space-time. In this framework Minkowski four-dimensional space-time is the mathematical model of physical space in which Einsteins general theory is formulated. Minkowski space is independent of the inertial frame of reference and is a consequence of the postulates of special relativity [27, 29].

Euclidean space is the flat analog of Riemannian geometry while Minkowski space is considered as the flat analog of curved space-time, which is known in mathematics as pseudo-Riemannian geometry. Considerations of four-dimensional space-time make embedded moving manifolds three dimensional, where parametric time t, describing the motion of manifolds, may not have anything to do with proper time τ used in general relativity.

To briefly describe Minkowskian space-time, let us refer to arbitrary coordinates X^α, α = 0, …, 3, where the position vector R is expressed in coordinates as R = R(X^α). Bold letters throughout the manuscript designate vectors. Latin letters in indexes indicate surface related tensors. Greek letters in indexes show tensors related to the ambient space. All equations are fully tensorial and follow the Einstein summation convention.

Suppose that Sⁱ (i = 0, 1, 2) are the surface coordinates of the moving manifold S (Figure 1). Coordinates Sⁱ, X^α are arbitrarily chosen so that sufficient differentiation is achieved in both space and parametric time. The surface equation in ambient coordinates can be written as X^α = X^α(t, Sⁱ) and the position vector can be expressed as

Covariant bases for the ambient space are introduced as X_α = ∂_αR, where ∂α=∂/∂Xα. The covariant metric tensor is the dot product of covariant bases

The contravariant metric tensor is defined as the matrix inverse of the covariant metric tensor, so that ηανηνβ=δβα, where δβα is the Kronecker delta. From definition (3) it follows that η₀₀ = X₀·X₀ and consequently if for Minkowskian space-time, the space like signature is set (−1, +1, +1, +1), then X₀ = (i, 0, 0, 0)⁴. Therefore, vector components are complex numbers in general. As far as the ambient space is set to be Minkowskian, the covariant bases are linearly independent, so that the square root of the negative metric tensor determinant is unit -|η..|=1. Furthermore, the Christoffel symbols given by

vanish and the equality between partial and curvilinear derivatives follows ∂_α = ∇_α. In Minkowski space-time (later space) the ∂_α partial derivative and ∇_α curvilinear derivative are the same. Everywhere in calculations we use ∂ letter for the ambient space derivative and keep in mind that when referring to Minkowski space the derivative has index in Greek letters and, in that case, it is the same as partial derivative. When indexes are mixed Greek and Latin letters the last statement, as is shown below, does not hold in general.

Now let's discuss tensors on the embedded surface with arbitrary coordinates Sⁱ, where i = 0, 1, 2. Latin indexes throughout the text are used exclusively for curved surfaces and curvilinear derivative ∇_i is no longer the same as the partial derivative ∂i=∂/∂Si. Similar to the bases of ambient space, covariant bases of an embedded manifold are defined as S_i = ∂_iR and the covariant surface metric tensor is the dot product of the covariant surface bases:

The definition (4) dictates that the surface is three dimensional pseudo Riemannian manifold, because ambient space is four dimensional Minkowskian space and the surface in four manifold is three manifold.

Analogically to space metric tensor g^ij the contravariant surface metric tensor is the matrix inverse of the covariant one g_ij. The matrix inverse nature of covariant-contravariant metrics gives possibility to raise and lower indexes of tensors defined on the manifold. The surface Christoffel symbols are given by

and along with Christoffel symbols of the ambient space provide all the necessary tools for covariant derivatives to be defined as tensors with mixed space/surface indexes:

where ηiγ is the shift tensor which reciprocally shifts space bases to surface bases, as well as space metric to surface metric. For instance, Si=ηiαXα and

The metrilinic property ∇_ig_mn = 0 of the surface metric tensor is a direct consequence of (4, 5) definitions, therefore S_m · ∇_iS_n = 0. The S_m and ∇_iS_n vectors are orthogonal, so that ∇_iS_n must be parallel to the N surface normal

where B_ij is the tensorial coefficient of the (6) relationship and is generally referred to as the symmetric curvature tensor. The trace of the curvature tensor with upper and lower indexes is the mean curvature and its determinant is the Gaussian curvature. It is well-known that a surface with constant Gaussian curvature is a sphere, consequently a sphere can be expressed as:

When the constant becomes null the surface becomes either a plane or a cylinder. Equation (7) is the expression of constant mean curvature (CMC) surfaces in general. Finding the curvature tensor defines the way of finding covariant derivatives of surface base vectors and so (6, 7) provide the way of finding surface base vectors which indirectly leads to the identification of the surface.

After defining the metric tensor for ambient space η_μν (3) and the metric tensor for a moving surface g_ij (4), we now proceed with a brief review of surface velocity, t explicit (parametric) time derivative of surface tensors and time differentiation theorems for the surface/space integrals. The original definitions of time derivatives for moving surfaces were given in Hadamard [30] and recently extended in Grinfeld [26].

For the definition of surface velocity we need to define ambient coordinate velocity V^α first and to show that the coordinate velocity is the α component of the surface velocity. Indeed, by the velocity definition

taking into account (2), R the position vector is tracking the material point coordinate Sⁱ. Therefore, by the partial time differentiation of (2) and definition of ambient base vectors, we find that V the surface velocity is

Consequently V^α is the ambient component of the surface velocity. According to (9), the normal component of the surface velocity is the dot product with the surface normal

The normal component C of the surface velocity is generally referred to as an interface velocity and is invariant in contrast with the coordinate velocity V^α. Its sign depends on a choice of the normal. The projection of the surface velocity on the tangent space (Figure 2) [28] is tangential velocity and can be expressed as

Graphical illustrations of coordinate velocity V^α, interface velocity C and tangential velocity Vⁱ are given in Figure 2. The surface velocity V can be expressed as

There is a clear geometric interpretation of the interface velocity [26, 28]. Let the surface at two nearby moments of time t and t + Δt be S_t, S_{t + Δt}. Suppose that A is a point on the S_t surface and the corresponding point B, belonging to S_{t + Δt}, has the same surface coordinate as A (Figure 3), then AB ≈ VΔt. Let P be the point where the unit normal N ∈ S_t intersects the surface S_t+Δt. Then for small enough Δt, the angle ∠APB ≈ π/2 and AP ≈ V · NΔt. Therefore, C can be defined as

and can be interpreted as the instantaneous velocity of the surface in the normal direction. It is worth mentioning that the sign of the interface velocity depends on the choice of the normal. Although, C is a scalar, it is called interface velocity because the normal direction is implied.

In this section we briefly explain the concept behind the invariant time derivative for scalar and tensor fields defined on moving manifolds, even though these concepts are already given [26]. Suppose that invariant tensor field T is defined on the manifold at all time. To define the time invariant derivative of the tensor field it is necessary to capture the rate of change of T in the normal direction. The physical explanation of why the deformations along the normal direction are so important is given below when discussing integrals. This is similar to how C measures the rate of deformation in the normal direction. For a given point A ∈ S_t, find the point B ∈ S_t+Δt and P the intersection of S_t+Δt and the straight line orthogonal to S_t (Figure 3). Then the geometrically intuitive definition dictates that

Because (13) is entirely geometric, it must be free from choice of a reference frame. Therefore, it is invariant. On the other hand, from the geometric construction it follows that

T(B) is related to T(P) because B,P are nearby points and are situated on the S_t+Δt surface B, P ∈ S_{t + Δt}. Then

since ∇_iT shows rate of change in the tensor field along the surface and ΔtVⁱ indicates the directed distance BP. After a few lines of algebra, taking into account Equations (14, 15) in (13), we find

Generalization of (16) to any arbitrary tensors with mixed space and surface indexes is given by the formula

where Christoffel symbol Γ·mn for moving surfaces is Γ·mn=∇mVn-CBmn. The derivative commutes with contraction, satisfies sum, product and chain rules, is metrinilic with respect to the ambient metrics and does not commute with the surface derivative [26]. Also from (13) it is clear that the invariant time derivative applied to the time independent scalar vanishes.

Time differentiation of surface and space integrals have a central role in evaluation of the principle of least action. Dependence of time variation of the potential energy on the geometry becomes rigorously clarified from these theorems. For any scalar field T = T(t, Sⁱ) defined on a Minkwoskian domain Ω with boundary S manifold evolving with the interface velocity C, the evolution of the space integral and surface integral for closed compact manifolds are given by the formulas

The first term in the integral represents the rate of change of the tensor field, while the second term shows changes in the geometry. It therefore properly takes into account the convective and advective terms due to volume motion. We are not going to reproduce proof of these theorems here⁵, but instead we give intuitive explanation of why only interface velocity has a role and tangential velocities do not appear in the integration. If the surface velocity has no interface velocity and has only tangential components, then the tangent velocity translates each point to it's neighboring ones and does not add new area and volume to the surface and space. Therefore, it provokes rotational movement of the material object and can be excluded from the integration. This statement becomes obvious for one dimensional motion. If the material point is moving along some trajectory, then the velocity is tangential to the curve. However, motion of the material point along the curve can be understood as the motion of the curve embedded in the plane. If the curve has no interface velocity, then it only slides in the ambient plane without changing the local length⁶.

In this section we provide several theorems, that will be directly used to deduce equations of motions. The first such theorem is the general Gauss theorem about integration, which gives the rule for the reciprocal transfer of space integral to surface integral. For a domain Ω in Minkowski space with the boundary S, for any sufficiently smooth tensor field T^α, the generalized Gauss theorem reads

The proof is simple if one uses the Voss-Weyl formula to deduce the theorem. For any sufficiently smooth tensor field in Minkowski space, the Voss-Weyl formula [26] reads

Using (21) in the right part of (20) and the designation η = −|η_..|, we have

where dX^α = dX⁰dX¹dX²dX³. This term is subject to the Gausss theorem in the arithmetic space. Since arithmetic space and Minkowski space, which is pseudo-Euclidean, can be related to Cartesian coordinates, Minkowski space can be identified as an arithmetic one and the Gauss theorem for the arithmetic space can be used. Thus, using unity of the Minkowski space metric tensor determinant one may prove that⁷

where g = |g_..|. This proves that the generalized Gauss's theorem holds for pseudo-Riemannian manifolds embedded in Minkowski space.

The next step is to provide short proofs for Weingarten's and Thomas' formulas by using the relation between the surface derivative and the interface velocity.

Weingarten's formula expresses the surface covariant derivative of the surface normal in the product of the shift and mixed curvature tensors. Proof follows from the definition NαNα=1, from which we find Nα∇iNα=0. On the other hand

If we apply the covariant derivative to (22) and take into account that from (6) Bji=Nα∇jηiα then by the product rule we find

Let's contract both sides of (23) with η^iβ and take into account the commonly used relationship in tensor calculus NαNβ+ηiαηβi=δβα, then we find

Since the second term of the last equality vanishes, we get (24), also known as Weingarten's formula.

Now we turn to the Thomas formula which allows calculation of the invariant time derivative of the surface normal. Indeed, using the invariant time derivative formula for the surface base vector [26]

and dotting both sides of (25) with N, N·∇·Si=N·N∇iC and using the product rule, taking into account that N · S_i = 0, we find ∇iC=-∇·NSi, therefore

Equation (26) is generally referred to as the Thomas formula.

Since we have all mathematical preliminaries in hand we can proceed with derivation of master equations of motion. To derive the equations we apply the calculus of moving surfaces to the motion of compact and closed manifolds in an electromagnetic field. In this step we only discuss free motion of the single closed surface, where “single” surface means boundary of the single material body and “free” means contact with environment is ignored⁸. The interaction with the environment can be incorporated into the equations later on⁹. The surface is treated as a continuum medium of material particles (points), where charge and mass distribution is heterogeneous. The boundary of the body is the surface with a surface mass density ρ and a surface charge density q. The surface can be semipermeable at some material points, meaning the charge can flow through the surface. Interaction between the material points is exclusively electromagnetic, as the mass of each material particle is set to be infinitely small compared to unit charges. The ambient space is set to be Minkowskian, the body is four dimensional and has the surface boundary of three dimensional manifold. Electromagnetic interaction between the material particles and the heterogeneous distribution of charges throughout the object induces motion of the surface and the potential energy of the interaction can be modeled as

where the electromagnetic tensor F_αβ is the combination of the electric and magnetic fields in a covariant antisymmetric tensor [25, 29]. The electromagnetic covariant four-potential is a covariant four vector A_· = (−φ/c, a) composed of the φ electric potential and the a magnetic potential. Contravariant four current J^· = (cQ, j) is the contravariant four vector combining j electric current density and Q the charge density, c is the speed of light and μ₀ is the magnetic permeability of the vacuum. The Minkowski space metric tensor signature is set to be space-like (− + + +) throughout the paper. This formulation is fully relativistic though it can be easily simplified for non-relativistic cases. Raising and lowering the indexes is performed by the Minkowski metric η_αβ. The relation between the four potentials and the electromagnetic tensor is given by

As far as the boundary of the material body is a moving three manifold, the surface kinetic energy with variable surface mass density ρ and surface velocity V is

Subtraction of the potential energy (27) from the kinetic energy (29) leads to the system Lagrangian

where S is the boundary of Ω. Hamilton's least action principle [31] for the given Lagrangian (30) reads

For proper evaluation of the (31) Lagrangian we start from the simplest term first, the potential energy. Since (27) is the space integral by theorem (18) we have

According to (32) determination of variation of potential energy is calculated from the time differential of the space integrand. Following standard algebraic manipulations for classical electrodynamics, we find

where u=(1/4μ0)FμνFμν+AμJμ and the fact that u is a function of A_α and ∂_βA_α and at the boundary condition ∂A_α/∂t = 0 the last term vanishes. It is easy to show that,

To calculate the last integrand (33), we take into account the definition (28) and note that the covariant electromagnetic tensor can be obtained by lowering indexes in contravariant tensor Fαβ=ηγαηκβFγκ. The electromagnetic tensor is antisymmetric F^αβ = −F^βα, so that

Taking into account (33–35) in (32) we find the variation of the potential energy

Now we turn to the calculation of the kinetic energy variation. To deduce the variation for the kinetic energy let's define the generalization of conservation of mass law first. The variation of the surface mass density must be so that dm/dt = 0, where

is the surface mass with ρ surface mass density. Since we discuss compact closed manifolds the boundary conditions v=niVi=0 dictate that a pass integral along any curve across the surface must vanish. This statement formally taking into consideration (37), can be rewritten as

where n is a normal of the curve that lies in the tangent space, v is the velocity of the γ curve. Since the last integral from (38) mast be identical to zero for any integrand, one immediately finds the generalization of conservation of mass law

Incidentally, an equation for the surface charge conservation can analogously be deduced and it has exactly the same form. The equation (39) was also reported in Grinfeld [10]. To calculate the variation of the kinetic energy we use (19, 29, 39) and after a few lines of algebra, we find

Here we used the fact that at the end of variations the surface reaches the stationary point and therefore, by the Gauss theorem integral for ∇i(ρViV2/2) converted to line integral, vanishes [as we used it already in (38)]. To deduce the final form of equations of motion we decompose the dot product in the integral (40) into normal and tangential components. After a few lines of algebraic manipulations, we find

Using Weingartens formula (24), the metrilinic property of the Minkowski space base vectors ∇_iX_α = 0 and the definition of the surface normal N=NαXα, the last equation of (41) transforms

Taking into account (12) and its covariant and invariant time derivatives in (42), we find

Continuing algebraic manipulations using the formula for the surface derivative of the interface velocity (25), Thomas formula (26), and the definition of the curvature tensor (6) in (43), yield

Dotting (44) on V and combining it with (40) the last derivation reveals the variation of the kinetic energy

where the first part is the normal component and the second part is the tangent component of the dot product. Combination of (36, 45) with (31) reveals

To find the final form of the equations of motion we separate the dot product of the space integrand from (46) into normal and tangential components. Let the vector F with contravariant α component be

where F and Fi are normal and tangential components of F, by analogy we have for ∂A/∂t four vector partial time derivative

where A,Ai are the normal and tangential components of the partial time derivative of the four vector potential. Using the definitions (47, 48) the dot product of the two vectors is

Since equation (46) must hold for every V, F, ∂A/∂t vector, the normal and tangential components of the dot product must be equal so that taking into account (47–49) in (46), we find

After applying the Gauss theorem to the surface integrals in (50), the surface integrals are converted to a space integral so that one gets

To summarize (39, 50–52), equations of moving manifolds in an electromagnetic field read

Equations (53) are the master equations of motions.

A case that deserves some attention is the homogeneous symmetrical surface. In that case the only nonzero allowed “force” is FA≠0 and FiAi=0. This leads to significant simplification of the third equation from (53) and the second equation can be analytically solved for homogeneous, equilibrium surfaces as we have done for micelles [28]. When FA≠0 and FiAi≠0 then motion of the surface induces swimming of the body. The case FA=0 and FiAi≠0, as shown below, simplifies to the Euler equation for dynamic fluid free motion and to the Navier-Stokes equation or to magneto-hydrodynamic (MHD) equations if one takes into account interactions with the environment.

Equations (53) are correct for freely moving manifolds of the body in a vacuum. Generalization can be trivially achieved if instead of the electromagnetic tensor F_αβ one proposes the electromagnetic stress energy tensor T^αβ, which is related to the electromagnetic tensor by the relationship

For objects in matter the electromagnetic tensor F^μν in (47, 53) is replaced by the electric displacement tensor D^μν and by the magnetization-polarization tensor M^μν so that

The charge density Q and four current J become the sum of bound and free charges and of bound and free four currents, respectively. The electric displacement tensor, magnetization tensor, free and bound charges/currents can be modeled differently for different problems, therefore the general equations (53) can be modified as needed.

To link the above formulated problem with real physical surfaces it is necessary to do some modeling. To begin let's illustrate macromolecules¹⁰ as two dimensional fluid manifolds with the thickness of variable mass density Figure 4. Even though molecular surfaces are three manifold in Minkowski space, in some cases¹¹ it can be modeled as moving two manifolds in Euclidean space. The surface is considered to be semipermeable against partial charges and water molecules. The permeability defines the surface mass density as a variable and the volume charge also becomes variable. The variability of charge and mass densities is properly taken into account in the equations of motion (53).

Let's model a bio-macromolecular surface as a Gaussian map contoured at 2 Å to 8 Å resolution. Figure 4 shows Gaussian maps for the protein (Figures 4A,B) and for the DNA (Figures 4C,D). Ω is the space inside the macromolecules and the boundary of the space is the surface S. S_i base vectors are defined in the tangent space of the Gaussian map. S_ij is the metric tensor of the map. These are illustrations of surfaces as two-manifolds embedded in Euclidian space and are only true for non-relativistic representations, therefore they do not show the shape of three-manifolds in Minkowski space-time. Figures 4A,C show a Gaussian map of the polypeptide main chain of a protein and of a polynucleotide double helical DNA dodecamer, respectively¹². Figures 4B,D show thickness variations, captured by surface mass density, of the modeled surfaces for the protein and DNA. Light gray is the Gaussian map at 5 Å resolution while the blue surface indicates a more detailed surface contoured at 2 Å resolution. Thickness variation can be induced by diffusion of solvent molecules at solvent accessible sites; e.g., sites marked by water molecules obtained from crystal structures as illustrated in the Figures 4A,C (red and white sticks), or by thermal fluctuation of amino acid sidechains. In all these cases, the surface thickness variation, captured by ρ surface mass density, is in the range of angstrom to nanometer. This range is higher for micelles, cell membranes, fluid films etc. If the system is in aqueous solution then the surface motion is determined by so called hydrophobic-hydrophilic interactions.

As we already stated in the introduction, hydrophobic and hydrophilic interactions incorporate dispersive interactions throughout the molecules, mainly related to electrostatics and electrodynamics (Van der Waals forces), induced by permanent (water molecules) or induced dipoles (dipole-dipole interactions) and possibly quadrupole-quadrupole interactions (for instance stacking or London forces) plus ionic interactions (Coulomb forces) [20]. The hydrophobic effect can be considered as synonymous with dispersive interactivity with water molecules and the hydrophilic one as synonymous with polar interactivity with water molecules [14, 19, 20]. All these interactions have one common feature and can be unified as electromagnetic interaction's dependence on interacting bodies' geometries, where by geometries we mean shape of the objects' surfaces. To model potential energy we note that on the scale of hydrophobic-hydrophilic interactions, which usually occurs at nanometer distances [19, 20], no interactions other than electromagnetic forces are available. An electromagnetic field is set up by dipole moments of water molecules and partial charges of molecules. In other words, we have a closed, smooth manifold in aqueous solution where charge and water molecules could migrate through the surface Figure 4. The surface can be of mixed nature (hydrophobic, hydrophilic, or both) with randomly distributed polar or non-polar groups and can be compressible, continuously deformable and permeable against water and ionic charges. At the nanometer scale, for small masses, potential energy can be electromagnetic only. Therefore, we have potential energy density constructed from the electromagnetic tensor plus the term related to variation of charges as it is defined in (27). Even though modeling of potential energy as electromagnetic interaction energy is fairly clear, the dependence of these interactions on the object's geometry is not. The geometry dependence becomes visible only after the complete formulation of the equations of motion (53).

To demonstrate effectiveness of (53) let's discuss free motion of two manifolds embedded in three dimensional Euclidean space for the stationary surface in an electrostatic field. We have the following conditions: V = 0 stationary surface in electrostatic field where a = 0, j = 0, A_· = (−φ/c, 0), J^· = (cQ, 0) and ∂ = (0, ∂_x, ∂_y, ∂_z). Then from second equation of (53) with the condition (46), we find

Taking into account the definition of electromagnetic tensor and that we consider the electrostatic field, the partial derivative of the electromagnetic tensor in (56) is (1/c)∂βEβ and therefore

By the definition of the electric field E^β = −∂^βφ and c2μ0=1/ϵ0 so that (57) transforms to

The equation (58) is generally known as the Poisson-Boltzmann equation in vacuum and was proposed to describe the distribution of the electric potential in the direction of the normal to a charged surface [34–36].

Here we demonstrate that the Poisson-Boltzmann equation is a special case and can be obtained from the equations of motion (53) for stationary surfaces in an electrostatic field. To support this statement we have generated electrostatic field lines using the Adaptive Poisson-Boltzmann Solver (APBS) [37] software for the protein [32] and the DNA [33] (Figure 5). As can be seen from Figure 5 field lines follow the surface normal as expected from the second equation of (53) for the V = 0 stationary case.

In this subsection we demonstrate that, the equations of motion simplify to Maxwell equations for stationary interfaces C = 0 and massless ρ = 0 three manifolds embedded in Minkowski space. Indeed, from the second and the third equations of (53), taking into account that in the stationary case the second term in (32) vanishes, we find

Adding (59) to (60) and taking into account (47) and (49) one obtains

(61) must hold for any partial time derivative of the four vector potential, therefore

and the Maxwell equations with the source in the vacuum follow.

This is a somewhat unexpected result: any three manifold with stationary interface C = 0 and with massless surface mass density ρ = 0¹³, satisfies Maxwell equations. However, the photon is arguably the only massless particle that satisfies the Maxwell equation, therefore the photon can be interpreted as a stationary interface three manifold embedded in Minkowski space with vanishing surface mass density.

In this section we simplify the equations of motion using physical arguments and demonstrate that the equation system (53) yields the Euler equation for a dynamic fluid for some simplified cases. Let's propose that a moving fluid has a planar surface B_ij = 0 with stationary interface C = 0 and is embedded in Euclidean three space. Then simplifications of (53) lead to a system of equations of motion

The first equation of (63) is the continuity equation for the surface mass density and is conservation of mass at the flat space; the second one yields that normal component of the dot product F · (∂A/∂t) vanishes. To simplify the last equation of (63) we note that the total 'force' acting on the volume is equal to the integral -∫SpdS of the total pressure p, taken over the boundary (surface) of the volume. Applying Gauss' theorem to the surface integral by taking into account that pressure across the surface acts in a normal direction so that it can be written as p=pαNα, then

On the other hand FiAi is a cause of the gradient of the tangential velocity and the tangential gradient of the pressure, therefore

Taking into consideration (64, 65) in (63) and applying Gauss' theorem to the space integral, we find

According to Weingarten's formula (24) N^α is invariant vs the surface derivative for flat manifolds and therefore, can be taken into the surface covariant derivative, so that ViNα∇ipα=Vi∇ip¹⁴. Then (66) after subtracting V_i yields

Taking into account that for flat surfaces Christoffel symbols vanish and ∇_j = ∂_j one immediately recognizes that the last equation (67) becomes the classical Euler equation of fluid dynamics.

As stated above the equations of motion (53) are formulated for freely moving manifolds; i.e., interaction with the environment is ignored and matter is set to be a vacuum. Though it can be trivially generalized for the matter and then simplified, instead of giving Euler equation, will lead to the more complete Navier-Stokes equation and or magnetohydrodynamic equations. For instance, in matter, according to (55), the electromagnetic tensor becomes the sum of the electric displacement and magnetization tensors. Therefore, in (67), instead of a pure pressure gradient we will have an additive term coming from the magnetic field so that (67) will transform to the ideal magneto hydrodynamic equation.

Analogously, if interaction with an environment is taken into account, then instead of a single surface we have two surfaces at the surface/environment interface and the Lagrangian (30) is split into two kinetic energy terms, one for surface and another one for the environmental interface. All these will occur as additive terms in the third equation of (53) so that equation (67) will transform to the Navier-Stokes equation.

Let's answer the question: what is the shape of micelles formed from lipid molecules when they are in thermodynamic equilibrium with solvent. Lipids have hydrophilic heads and hydrophobic tails, so that in solutions they tend to form a surface with heads on one side and tails on the other. Since the tails disperse the water molecules, the surface made is closed and has some given volume. Such structures are called micelles [38]. Since lipids form a homogeneous surface, in equilibrium conditions we must have

and FiAi=0. Usually the speed of micelle interface motion is in the range of nm/ns and, therefore there is no necessity of discussion of a relativistic formalism so that the surface is two dimensional and the space is Euclidean. The surface dynamic is slow, magnetic field is much smaller then electric field B² < < E² and the potential energy becomes

Using the first law of thermodynamics, (69) can be modeled as a volume integral from the surface pressure [28] and

On the other hand, taking into account the conditions (68, 69), the total potential energy of the surface can be modeled as

Taking into account (71), the system Lagrangian becomes the same as it is in (1) and its variation leads to the equation

(72) was first reported in Grinfeld [10]. Using (69, 70, 72) in the equations of motion (53), after simple algebra we find

When the homogeneous surface, such as a micelle, is in equilibrium with the environment then the solution of the (73)¹⁵ is

From equation (74) the generalized Young-Laplace relation which connects the surface pressure to the curvature and the surface tension immediately follows. (74) dictates that the homogeneous surfaces in equilibrium with environment adopt a shape with constant mean curvatures (CMC), which explains the well anticipated lamellar, cylindrical and spherical shapes of micelles. This is another unexpected and surprisingly simple solution to the equations of motion (53).

Equations of motion (53) further simplify for two dimensional surfaces. In non-relativistic framework the space is three dimensional Euclidean (α = 1, 2, 3 and V⁰ = 0 limit), the surface is two-dimensional Riemannian (i = 1, 2) and the potential energy becomes

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 same formalism as in (65, 66, 69, 70) into account we find

Taking (76, 77) along with that p=ϵ02E2-12μ0B2+φQ-aj into account we end up with the following equations of motion for two dimensional surfaces:

Alternative way of deducing (78–80) without using (53) is given in Svintradze [28]. In equations of motions for two dimensional surfaces (78–80) we mention that only first equation (78) is the same as dynamic fluid film equations [10]. (78) captures surface mass density variation during the manifold motion, second equations (79) shows how the surface moves in normal direction and third equation (80) shows how the surface moves in tangent directions. Soup films, water droplets and for all two dimensional surfaces which satisfies preconditions: (1) the surface is homogeneous, (2) surface potential energy density σ is time invariable, and (3) the surface is in thermodynamic equilibrium with environment can be solved exactly the same way as it is done in sub-section 4.4 with the solution (74).

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.