Source author record

Keith Promislow

Keith Promislow appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Competition and Complexity in Amphiphilic Polymer Morphology

We analyze the competitive evolution of codimension one and two morphologies within the $H^{-1}$ gradient flow of the strong Functionalized Cahn-Hilliard equation. On a slow time scale a sharp hypersurface reduction yields a degenerate Mullins-Sekerka evolution for both codimension one and two hypersurfaces, leading to a geometric flow that depends locally on curvatures couples to the dynamic value of the spatially constant far-field chemical potential. Both codimension one and two morphologies admit two classes of bifurcations, one leads to pearling, a short-wavelength in-plane modulation of interfacial width, the other flips motion by curvature to the locally-ill posed motion against curvature, which leads to fingering instabilities. We present a bifurcation diagram for the morphological competition, and compare our results quantitatively to simulations of the full system and qualitatively to simulations of self-consistent mean field models and laboratory experiments; illuminating the role of the pearling bifurcation in the development of complex network morphologies.

preprint2022arXiv

Undulated Bilayer Interfaces in the Planar Functionalized Cahn-Hilliard Equation

Experiments with diblock co-polymer melts display undulated bilayers that emanate from defects such as triple junctions and endcaps, \cite{batesjain_2004}. Undulated bilayers are characterized by oscillatory perturbations of the bilayer width, which decay on a spatial length scale that is long compared to the bilayer width. We mimic defects within the functionalized Cahn-Hillard free energy by introducing spatially localized inhomogeneities within its parameters. For length parameter $\varepsilon\ll1$, we show that this induces undulated bilayer solutions whose width perturbations decay on an $O\!\left(\varepsilon^{-1/2}\right)$ inner length scale that is long in comparison to the $O(1)$ scale that characterizes the bilayer width.

preprint2020arXiv

Asymptotic Behaviour of Time Stepping Methods for Phase Field Models

Adaptive time stepping methods for metastable dynamics of the Allen Cahn and Cahn Hilliard equations are investigated in the spatially continuous, semi-discrete setting. We analyse the performance of a number of first and second order methods, formally predicting step sizes required to satisfy specified local truncation error $σ$ in the limit of small order parameter $ε\rightarrow 0$ during meta-stable dynamics. The formal predictions are made under stability assumptions that include the preservation of the asymptotic structure of the diffuse interface, a concept we call profile fidelity. In this setting, definite statements about the relative behaviour of time stepping methods can be made. Some methods, including all so-called energy stable methods but also some fully implicit methods, require asymptotically more time steps than others.The formal analysis is confirmed in computational studies. We observe that some provably energy stable methods popular in the literature perform worse than some more standard schemes. We show further that when Backward Euler is applied to meta-stable Allen Cahn dynamics, the energy decay and profile fidelity properties for these discretizations are preserved for much larger time steps than previous analysis would suggest. The results are established asymptotically for general interfaces, with a rigorous proof for radial interfaces. It is shown analytically and computationally that for most reaction terms, Eyre type time stepping performs asymptotically worse due to loss of profile fidelity.

preprint2020arXiv

Codimension One Minimizers of Highly Amphiphilic Mixtures

We present a modified form of the Functionalized Cahn Hilliard (FCH) functional which models highly amphiphilic systems in solvent. A molecule is highly amphiphilic if the energy of a molecule isolated within the bulk solvent molecule is prohibitively high. For such systems once the amphiphilic molecules assemble into a structure it is very rare for a molecule to exchange back into the bulk. The highly amphiphilic FCH functional has a well with limited smoothness and admits compactly supported critical points. In the limit of molecular length epsilon approaches 0, we consider sequences with bounded energy whose support resides within an epsilon-neighborhood of a fixed codimension one interface. We show that the FCH energy is uniformly bounded below, independent of epsilon >0, and identify assumptions on tangential variation of sequences that guarantee the existence of subsequences that converge to a weak solution of a rescaled bilayer profile equation, and show that sequences with limited tangential variation enjoy a lim inf inequality. For fixed codimension one interfaces we construct bounded energy sequences which converge to the bilayer profile and others with larger tangential variation which do not converge to the bilayer profile but whose limiting energy can violate the lim inf inequality, depending upon the energy parameters.

preprint2020arXiv

Gradient Invariance of Slow Energy Descent: Spectral Renormalization and Energy Landscape Techniques

For gradient flows of energies, both spectral renormalization (SRN) and energy landscape (EL) techniques have been used to establish slow motion of orbits near low-energy manifold. We show that both methods are applicable to flows induced by families of gradients and compare the scope and specificity of the results. The SRN techniques capture the flow in a thinner neighborhood of the manifold, affording a leading order representation of the slow flow via as projection of the flow onto the tangent plane of the manifold. The SRN approach requires a spectral gap in the linearization of the full gradient flow about the points on the low-energy manifold. We provide conditions on the choice of gradient under which the spectral gap is preserved, and show that up to reparameterization the slow flow is invariant under these choices of gradients. The EL methods estimate the magnitude of the slow flow, but cannot capture its leading order form. However the EL only requires normal coercivity for the second variation of the energy, and does not require spectral conditions on the linearization of the full flow. It thus applies to a much larger class of gradients of a given energy. We develop conditions under which the assumptions of the SRN method imply the applicability of the EL method, and identify a large family of gradients for which the EL methods apply. In particular we apply both approaches to derive the interaction of multi-pulse solutions within the 1+1D Functionalized Cahn-Hilliard (FCH) gradient flow, deriving gradient invariance for a class of gradients arising from powers of a homogeneous differential operator.

preprint2020arXiv

Manifolds of Amphiphilic Bilayers: Stability up to the Boundary

We consider the mass preserving $L^2$-gradient flow of the strong scaling of the functionalized Cahn Hilliard gradient flow and establish the nonlinear stability of a manifold comprised of quasi-equilibrium bilayer \muckmucks up to the manifold's boundary. In the limit of thin but non-zero interfacial width, $\varepsilon\ll1,$ the bilayer manifold is parameterized by meandering modes that describe the interfacial evolution and "pearling" modes that control the structure of the profile near the interface. The pearling modes are weakly damped and can lead to the dynamic rupture of the interface. Amphiphilic interfaces can lengthen to decrease energy. We introduce an implicitly defined parameterization of the interfacial shape that uncouples this growth from the parameters describing the shape and introduce a nonlinear projection onto the manifold from a surrounding neighborhood. The bilayer manifold has asymptotically large but finite dimension tuned to maximize normal coercivity while preserving the wave-number gap between the meandering and the pearling modes. Modulo a pearling stability assumption, we show that the manifold attracts nearby orbits into a tubular neighborhood about itself so long as the interfacial shape remains sufficiently smooth and far from self-intersection. In a companion paper, arXiv:1907.02196, we identify open sets of initial data whose orbits converge to circular equilibrium after a significant transient, and derive a singularly perturbed interfacial evolution comprised of motion against curvature regularized by an asymptotically weak Willmore term.

preprint2020arXiv

Regularized Curve Lengthening from the Strong FCH Flow

We present a rigorous analysis of the transient evolution of nearly circular bilayer interfaces evolving under the thin interface limit, $\varepsilon\ll1$, of the mass preserving $L^2$-gradient flow of the strong scaling of the functionalized Cahn-Hilliard equation. For a domain $Ω\subset{\mathbb R}^2$ we construct a bilayer manifold with boundary comprised of quasi-equilibrium of the flow and a projection onto the manifold that associates functions $u$ in an $H^2$ tubular neighborhood of the manifold with an interface $Γ$ embedded in $Ω$. These interfaces, and hence the bilayer manifold, are parameterized by a finite but asymptotically large number of degrees of freedom. The manifold contains a unique, up to translation and mass constraint, equilibrium of the gradient flow whose projected interface is circular up to exponentially small corrections. The thin tubular neighborhood is forward invariant under the flow with orbits that ultimately converge to the equilibrium. Projections of these orbits yield an interfacial evolution equivalent at leading order to the regularized curve-lengthening motion characterized by normal motion {\sl against} mean curvature, regularized by a higher order Willmore expression. The curve lengthening is driven by absorption of excess mass from the regions of $Ω$ away from the interface, generically leading to nontrivial dynamics that are ill-posed in the $\varepsilon\to0$ limit.

preprint2016arXiv

Dependence of the dielectric constant of electrolyte solutions on ionic concentration - a microfield approach

We present a novel microfield approach for studying the dependence of the orientational polarization of the water in aqueous electrolyte solutions upon the salt concentration and temperature. The model takes into account the orientation of the solvent dipoles due to the electric field created by ions, and the effect of thermal fluctuations. The model predicts a dielectric functional dependence of the form $\varepsilon(c)=\varepsilon_w-βL(3αc/β),\quadβ=\varepsilon_w-\varepsilon_{\rm ms}$, where $L$ is the Langevin function, $c$ is the salt concentration, $\varepsilon_w$ is the dielectric of pure water, $\varepsilon_{\rm ms}$ is the dielectric of the electrolyte solution at the molten salt limit, and $α$ is the total excess polarization of the ions. The functional form gives a remarkably accurate description of the dielectric constant for a variety of salts and a wide range of concentrations.

preprint2015arXiv

Existence, Bifurcation, and Geometric Evolution of Quasi-Bilayers in the Multicomponent Functionalized Cahn-Hilliard Equation

Multicomponent bilayer structures arise as the ubiquitous plasma membrane in cellular biology and as blends of amphiphilic copolymers used in electrolyte membranes, drug delivery, and emulsion stabilization within the context of synthetic chemistry. We develop the multicomponent functionalized Cahn-Hilliard (mFCH) free energy as a model which allows competition between bilayers with distinct composition and between bilayers and higher codimensional structures, such as co-dimension two filaments and co-dimension three micelles. We investigate the stability and slow geometric evolution of multicomponent bilayer interfaces within the context of gradient flows of the mFCH, addressing the impact of aspect ratio of the lipid/copolymer unit on the intrinsic curvature and the codimensional bifurcation. In particular we derive a Canham-Helfrich sharp interface energy whose intrinsic curvature arises through a Melnikov parameter associated to lipid aspect ratio. We construct asymmetric homoclinic bilayer profiles via a billiard limit potential and show that the dominant co-dimensional bifurcation mechanism is via the layer-by-layer pearling observed experimentally.

preprint2014arXiv

Existence of pearled patterns in the planar Functionalized Cahn-Hilliard equation

The functionalized Cahn-Hilliard (FCH) equation supports planar and circular bilayer interfaces as equilibria which may lose their stability through the pearling bifurcation: a periodic, high-frequency, in-plane modulation of the bilayer thickness. In two spatial dimensions we employ spatial dynamics and a center manifold reduction to reduce the FCH equation to an 8th order ODE system. A normal form analysis and a fixed-point-theorem argument show that the reduced system admits a degenerate 1:1 resonant normal form, from which we deduce that the onset of the pearling bifurcation coincides with the creation of a two-parameter family of pearled equilibria which are periodic in the in-plane direction and exponentially localized in the transverse direction.

preprint2013arXiv

Adiabatic stability under semi-strong interactions: The weakly damped regime

We rigorously derive multi-pulse interaction laws for the semi-strong interactions in a family of singularly-perturbed and weakly-damped reaction-diffusion systems in one space dimension. Most significantly, we show the existence of a manifold of quasi-steady N-pulse solutions and identify a "normal-hyperbolicity" condition which balances the asymptotic weakness of the linear damping against the algebraic evolution rate of the multi-pulses. Our main result is the adiabatic stability of the manifolds subject to this normal hyperbolicity condition. More specifically, the spectrum of the linearization about a fixed N-pulse configuration contains essential spectrum that is asymptotically close to the origin as well as semi-strong eigenvalues which move at leading order as the pulse positions evolve. We characterize the semi-strong eigenvalues in terms of the spectrum of an explicit N by N matrix, and rigorously bound the error between the N-pulse manifold and the evolution of the full system, in a polynomially weighted space, so long as the semi-strong spectrum remains strictly in the left-half complex plane, and the essential spectrum is not too close to the origin.