Source-linked AI summary
A Thermodynamically Consistent Model for Multicomponent Vesicles
Wangbo Luo, Zhonghua Qiao, Yanxiang Zhao
TL;DR
Prior multicomponent-vesicle models included mechanochemical coupling, while this work formulates membrane deformation and protein segregation through a coupled free energy. It derives an L2 gradient flow and stabilized alternating ETD schemes, proving discrete energy dissipation under regularity assumptions. Two- and three-dimensional experiments validate dissipation and illustrate segregation and deformation.
Problem
Earlier work localized Ohta-Kawasaki energy and coupled it to a force-balance equation, whereas this paper targets a thermodynamically consistent variational coupling of membrane deformation and protein segregation.
Source: §1 · p. 2Method
The authors derive a coupled L2 gradient flow from a free energy combining protein-dependent bending, surface tension, volume penalty, and membrane-associated Ohta-Kawasaki energy, then construct stabilized alternating ETD1 and Strang-ETDRK2 schemes.
Source: §1 · p. 3, §2 · p. 4, §6 · p. 26Results
The alternating ETD1 and Strang-ETDRK2 schemes satisfy discrete energy dissipation under regularity assumptions, and two- and three-dimensional experiments validate dissipation while illustrating protein segregation and membrane deformation.
Source: §1 · p. 3, §4.4 · pp. 20–21Takeaways & Limitations
The coupled free-energy and gradient-flow structure provides a framework for studying energy minimizers and well-posedness of the coupled dynamics.
Source: §1 · p. 2The discrete energy-dissipation analysis depends on regularity assumptions, and the uniform bounds those assumptions require are not proved in this work.
Source: §6 · p. 26
Abstract
from arXiv · showhide
We develop a thermodynamically consistent diffuse-interface model for multicomponent membranes. The proposed model is derived from a coupled free energy functional that incorporates protein-dependent bending elasticity, diffuse surface tension, a volume penalty, and a membrane-associated Ohta--Kawasaki energy. Applying the Onsager variational principle, we derive a coupled $L^2$ gradient flow system and its energy dissipation law. We then construct a stabilized alternating ETD1 scheme and a stabilized alternating ETDRK2 scheme with Strang-type composition (alternating Strang-ETDRK2). To the best of our knowledge, the proposed alternating Strang-ETDRK2 scheme has not previously been developed and analyzed for coupled phase field systems. We further prove the discrete energy dissipation for both schemes under some regularity assumptions on the numerical solutions. Numerical experiments in two and three dimensions validate the discrete energy dissipation law, and present the protein segregation and membrane deformation produced by the proposed model.
1. Introduction
The paper formulates membrane deformation and protein segregation as a coupled variational model, motivated by observed interactions between membrane composition and geometry. It also develops stabilized alternating ETD schemes and proves their discrete energy dissipation under regularity assumptions.
1. Introduction
- Membrane domains have different thicknesses and mechanical properties, while curvature and line tension can produce buds, necks, tubes, and periodic shapes.Protein composition can modify bending rigidity, and curvature gradients can redistribute or localize softer and stiffer membrane areas.Source: §1 · p. 1
- Diffuse-interface methods model membrane deformation and topological changes without explicitly tracking the moving surface.Prior phase-field formulations include multicomponent membranes and two-component membranes on fixed computational domains.Source: §1 · p. 2
- The Ohta-Kawasaki energy models protein segregation through competition between local surface energy, which favors larger domains, and long-range repulsion, which favors smaller domains.The resulting competition produces microphase separation.Source: §1 · p. 2
- The model couples membrane deformation and protein segregation through a free energy, with protein composition affecting bending rigidity and membrane geometry affecting protein segregation.The authors contrast this variational coupling with their earlier mechanochemical model and note that it offers a framework for future analysis of minimizers and well-posedness.Source: §1 · pp. 2–3
- Prior ETD work addressed stiff phase-field and coupled systems, but nonlinear coupling in multicomponent vesicles makes scheme construction and stability analysis challenging.The paper identifies its alternating Strang-ETDRK2 scheme as not previously studied for the coupled multicomponent-membrane phase-field system.Source: §1 · p. 3
- The authors construct stabilized alternating ETD1 and Strang-ETDRK2 schemes and prove discrete energy dissipation for both under regularity assumptions on numerical solutions.Two- and three-dimensional experiments validate discrete dissipation and illustrate protein segregation and membrane deformation.Source: §1 · p. 3
2. Model Setting and the Coupled Free Energy Functional
The model combines membrane mechanics and protein segregation in a coupled free energy defined using membrane and protein phase fields. Its terms represent protein-dependent bending, surface tension and volume penalty, and membrane-associated Ohta-Kawasaki interactions.
2. Model Setting and the Coupled Free Energy Functional
- The membrane phase field ϕ marks the vesicle interior and exterior, while the protein phase field u marks protein-rich and protein-poor regions.Both variables take values near 1 in their respective regions and near 0 outside them.Source: §2 · p. 3
- The total energy combines protein-dependent bending, surface tension and enclosed-volume penalty, and membrane-associated Ohta-Kawasaki energy.The Ohta-Kawasaki term consists of local, long-range, and soft protein-mass components.Source: §2 · p. 4
- Bending rigidity depends on protein density, taking values κ+ in protein-rich regions and κ− in protein-poor regions.The phase-dependent rigidity function is chosen so its first and second derivatives vanish at the pure-region values.Source: §2 · p. 4
- The surface-tension term and volume penalty use coefficient λsurf, prescribed volume V0, and penalty coefficient Kvol.The volume penalty favors the prescribed volume but does not impose exact conservation.Source: §2 · pp. 4–5
- The membrane-associated Ohta-Kawasaki energy localizes protein interactions near the diffuse membrane and balances short-range and inverse-Laplacian long-range interactions.Its soft protein-mass penalty favors a prescribed protein fraction without enforcing exact conservation.Source: §2 · pp. 4–5
- In the sharp-interface limit, the local Ohta-Kawasaki term approximates line tension up to normalization, while quadratic potential extensions provide bounded derivatives for dissipation analysis.The local and long-range interactions compete to produce microphase separation.Source: §2 · pp. 4–5
3. L2 Gradient Flow from the Onsager Variational Principle
The authors derive the coupled evolution equations by applying the Onsager variational principle to the total free energy and a dissipation potential. The resulting L2 gradient flow has a continuous energy dissipation law.
3. L2 Gradient Flow from the Onsager Variational Principle
- The model uses the rate space V = L2(Ω) × L2(Ω) and obtains the coupled evolution equations by taking variations of the total energy along smooth perturbations.The first variation is represented using the L2 Riesz representation.Source: §3 · p. 6
- The Onsager dissipation potential uses a symmetric positive-definite operator K, chosen as a diagonal operator from the inverse mobilities.The mobility constants Mϕ and Mu are positive.Source: §3 · p. 6
- At each fixed state, minimizing the Rayleighian selects the evolution rate, and positive mobilities make its minimizer unique.The Rayleighian is strictly convex and coercive.Source: §3 · p. 6
- The resulting coupled L2 gradient flow satisfies a continuous energy dissipation law.The dynamics are nonconservative L2 gradient flow dynamics.Source: §3 · pp. 6–7
4. Exponential Time Differencing for the Coupled System
This section develops stabilized alternating ETD1 and Strang-ETDRK2 schemes for the coupled system and establishes their discrete energy dissipation. For Strang-ETDRK2, dissipation holds at every positive time step under the stated regularity and stabilization assumptions.
4.1. Stabilized Linear-nonlinear Splitting for the Coupled System
- The stabilized splitting uses a constant reference rigidity for protein-dependent bending rigidity to retain a constant-coefficient, self-adjoint, positive-definite linear operator.The φ- and u-equations are split using stabilization constants α1, α2 > 0 and β1, β2 > 0, respectively.Source: §4.1 · p. 8
4.2. Spatial Discretization by Second-order Finite Difference
- The spatial scheme uses second-order finite differences on a periodic grid, with discrete linear operators that are self-adjoint and positive-definite.The stated positivity follows from −Δh ⪰ 0 and positive κref, α1, α2, β1, and β2.Source: §4.2 · pp. 9, 11
4.3. Temporal Discretization by Stabilized Alternating ETD Schemes
- Because traditional alternating ETDRK2 may not attain second-order accuracy for the coupled system, the authors propose a Strang-type alternating ETDRK2 composition.Its local truncation error is O(τ^3), and for a fixed spatial grid, standard stability assumptions give second-order temporal convergence.Source: §4.3 · pp. 12, 14
4.4. Energy Dissipation of the Stabilized Alternating ETD Schemes
- Under Assumption 4.1 and the stabilization conditions, the alternating ETD1 method satisfies discrete energy dissipation for every time step τ > 0.The assumptions concern regularity of the numerical solutions; the theorem states the dissipation result for all positive τ.Source: §4.4 · p. 20
- The proof applies Lemmas 4.1 and 4.2 to the scheme’s three steps and combines estimates (4.44), (4.45), and (4.46) to reach its conclusion.Source: §4.4 · p. 21
5. Numerical Experiments
Numerical experiments validate discrete energy dissipation and show how long-range interaction strength, bending rigidities, and the membrane-associated energy scale affect protein domains and membrane shape. Three-dimensional simulations produce several distinct membrane morphologies through coupled protein segregation and deformation.
5. Numerical Experiments
- Simulations use periodic Cartesian grids, second-order finite differences, and fast Cosine transforms; the reported 2D and 3D time step is 10^-3.The grids are 256 × 256 in 2D and 128 × 128 × 128 in 3D.Source: §5 · p. 22
5.1. Discrete Energy Dissipation of the ETD Schemes
- The discrete free energy decreases during protein relaxation on a fixed membrane and decreases further after coupled membrane evolution begins.The membrane is held fixed until t = 200; the experiments use alternating ETD1 and Strang-ETDRK2 schemes.Source: §5.1 · p. 22, Figure 5.1 · p. 23
5.2. Effect of the Long-range Interaction Strength γ
- Increasing γ from 14000 to 25000 yields eight to eleven protein-rich subdomains, with stronger long-range interaction favoring smaller, nearly uniformly spaced domains.Figure 5.2 compares random initial distributions with the corresponding equilibrium patterns.Source: §5.2 · p. 23, Figure 5.2 · p. 24
5.3. Effects of λu and the Bending Rigidities κ+ and κ−
- Changing bending rigidities changes which protein phase occupies high-curvature regions and can alter membrane morphology; varying λu produces small changes in vesicle size and curvature.In the tested two-dimensional cases, λu's line-tension effect is less apparent because the two phases interact at only two boundary points.Source: §5.3 · p. 24, Figure 5.3 · p. 25
5.4. Representative Three-dimensional Patterns
- Three-dimensional simulations produce completely phase-separated, dumbbell-shape, pear-shape, and multi-bump membranes through coupled protein segregation and membrane deformation.The experiments vary protein fractions, bending rigidities, and protein-domain number and location.Source: §5.4 · pp. 25–26, Figure 5.4 · p. 26
6. Conclusion
The conclusion summarizes a thermodynamically consistent model and energy-dissipating schemes for coupled membrane deformation and protein segregation, alongside numerical demonstrations of parameter effects and three-dimensional morphologies. It identifies regularity assumptions and several analytical and modeling extensions as future work.
6. Conclusion
- The model couples deformation and protein segregation through a free energy, and the alternating ETD1 and Strang-ETDRK2 schemes satisfy proved discrete energy dissipation laws.The energy combines protein-dependent diffuse bending, surface tension, an enclosed-volume penalty, and membrane-associated Ohta-Kawasaki energy.Source: §6 · p. 26
- Numerical experiments validate discrete energy dissipation and show parameter-dependent membrane morphologies, including phase-separated, dumbbell-shape, pear-shape, and multi-bump forms.Source: §6 · p. 26
- The dissipation analysis depends on regularity assumptions for numerical solutions; convergence analysis, well-posedness, sharp-interface limits, hydrodynamics, and adhesion remain future directions.The authors propose investigating whether smooth initial data can yield the required regularity.Source: §6 · pp. 26–27