1
|
Bassani CL, Engel M, Sum AK. Mesomorphology of clathrate hydrates from molecular ordering. J Chem Phys 2024; 160:190901. [PMID: 38767264 DOI: 10.1063/5.0200516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2024] [Accepted: 03/13/2024] [Indexed: 05/22/2024] Open
Abstract
Clathrate hydrates are crystals formed by guest molecules that stabilize cages of hydrogen-bonded water molecules. Whereas thermodynamic equilibrium is well described via the van der Waals and Platteeuw approach, the increasing concerns with global warming and energy transition require extending the knowledge to non-equilibrium conditions in multiphase, sheared systems, in a multiscale framework. Potential macro-applications concern the storage of carbon dioxide in the form of clathrates, and the reduction of hydrate inhibition additives currently required in hydrocarbon production. We evidence porous mesomorphologies as key to bridging the molecular scales to macro-applications of low solubility guests. We discuss the coupling of molecular ordering with the mesoscales, including (i) the emergence of porous patterns as a combined factor from the walk over the free energy landscape and 3D competitive nucleation and growth and (ii) the role of molecular attachment rates in crystallization-diffusion models that allow predicting the timescale of pore sealing. This is a perspective study that discusses the use of discrete models (molecular dynamics) to build continuum models (phase field models, crystallization laws, and transport phenomena) to predict multiscale manifestations at a feasible computational cost. Several advances in correlated fields (ice, polymers, alloys, and nanoparticles) are discussed in the scenario of clathrate hydrates, as well as the challenges and necessary developments to push the field forward.
Collapse
Affiliation(s)
- Carlos L Bassani
- Institute for Multiscale Simulation, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - Michael Engel
- Institute for Multiscale Simulation, Friedrich-Alexander-Universität Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - Amadeu K Sum
- Phases to Flow Laboratory, Chemical and Biological Engineering Department, Colorado School of Mines, Golden, Colorado 80401, USA
| |
Collapse
|
2
|
Frohoff-Hülsmann T, Thiele U. Nonreciprocal Cahn-Hilliard Model Emerges as a Universal Amplitude Equation. PHYSICAL REVIEW LETTERS 2023; 131:107201. [PMID: 37739387 DOI: 10.1103/physrevlett.131.107201] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/29/2022] [Accepted: 08/08/2023] [Indexed: 09/24/2023]
Abstract
Oscillatory behavior is ubiquitous in out-of-equilibrium systems showing spatiotemporal pattern formation. Starting from a linear large-scale oscillatory instability-a conserved-Hopf instability-that naturally occurs in many active systems with two conservation laws, we derive a corresponding amplitude equation. It belongs to a hierarchy of such universal equations for the eight types of instabilities in homogeneous isotropic systems resulting from the combination of three features: large-scale vs small-scale instability, stationary vs oscillatory instability, and instability without and with conservation law(s). The derived universal equation generalizes a phenomenological model of considerable recent interest, namely, the nonreciprocal Cahn-Hilliard model, and may be of a similar relevance for the classification of pattern forming systems as the complex Ginzburg-Landau equation.
Collapse
Affiliation(s)
- Tobias Frohoff-Hülsmann
- Institute of Theoretical Physics, University of Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - Uwe Thiele
- Institute of Theoretical Physics, University of Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
- Center for Nonlinear Science (CeNoS), University of Münster, Corrensstrasse 2, 48149 Münster, Germany
| |
Collapse
|
3
|
Hankins SN, Zhou Y, Lohan DJ, Dede EM. Generative design of large-scale fluid flow structures via steady-state diffusion-based dehomogenization. Sci Rep 2023; 13:14344. [PMID: 37658099 PMCID: PMC10474040 DOI: 10.1038/s41598-023-41316-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Accepted: 08/24/2023] [Indexed: 09/03/2023] Open
Abstract
A computationally efficient dehomogenization technique was developed based on a bioinspired diffusion-based pattern generation algorithm to convert an orientation field into explicit large-scale fluid flow channel structures. Due to the transient nature of diffusion and reaction, most diffusion-based pattern generation models were solved in both time and space. In this work, we remove the temporal dependency and directly solve a steady-state equation. The steady-state Swift-Hohenberg model was selected due to its simplistic form as a single variable equation and intuitive parameter setting for pattern geometry control. Through comparison studies, we demonstrated that the steady-state model can produce statistically equivalent solutions to the transient model with potential computational speedup. This work marks an early foray into the use of steady-state pattern generation models for rapid dehomogenization in multiphysics engineering design applications. To highlight the benefits of this approach, the steady-state model was used to dehomogenize optimized orientation fields for the design of microreactor flow structures involving hundreds of microchannels in combination with a porous gas diffusion layer. A homogenization-based multi-objective optimization routine was used to produce a multi-objective Pareto set that explored the trade-offs between flow resistance and reactant distribution variability. In total, the diffusion-based dehomogenization method enabled the generation of 200 unique and distinctly different microreactor flow channel designs. The proposed dehomogenization approach permits comprehensive exploration of numerous bioinspired solutions capturing the full complexity of the optimization and Swift-Hohenberg design space.
Collapse
Affiliation(s)
- Sarah N Hankins
- Electronics Research Department, Toyota Research Institute of North America, 1555 Woodridge Avenue, Ann Arbor, MI, 48105, USA
| | - Yuqing Zhou
- Electronics Research Department, Toyota Research Institute of North America, 1555 Woodridge Avenue, Ann Arbor, MI, 48105, USA.
| | - Danny J Lohan
- Electronics Research Department, Toyota Research Institute of North America, 1555 Woodridge Avenue, Ann Arbor, MI, 48105, USA
| | - Ercan M Dede
- Electronics Research Department, Toyota Research Institute of North America, 1555 Woodridge Avenue, Ann Arbor, MI, 48105, USA
| |
Collapse
|
4
|
Revell CK, Herrera JA, Lawless C, Lu Y, Kadler KE, Chang J, Jensen OE. Modeling collagen fibril self-assembly from extracellular medium in embryonic tendon. Biophys J 2023; 122:3219-3237. [PMID: 37415335 PMCID: PMC10465709 DOI: 10.1016/j.bpj.2023.07.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 06/24/2023] [Accepted: 07/03/2023] [Indexed: 07/08/2023] Open
Abstract
Collagen is a key structural component of multicellular organisms and is arranged in a highly organized manner. In structural tissues such as tendons, collagen forms bundles of parallel fibers between cells, which appear within a 24-h window between embryonic day 13.5 (E13.5) and E14.5 during mouse embryonic development. Current models assume that the organized structure of collagen requires direct cellular control, whereby cells actively lay down collagen fibrils from cell surfaces. However, such models appear incompatible with the time and length scales of fibril formation. We propose a phase-transition model to account for the rapid development of ordered fibrils in embryonic tendon, reducing reliance on active cellular processes. We develop phase-field crystal simulations of collagen fibrillogenesis in domains derived from electron micrographs of inter-cellular spaces in embryonic tendon and compare results qualitatively and quantitatively to observed patterns of fibril formation. To test the prediction of this phase-transition model that free protomeric collagen should exist in the inter-cellular spaces before the formation of observable fibrils, we use laser-capture microdissection, coupled with mass spectrometry, which demonstrates steadily increasing free collagen in inter-cellular spaces up to E13.5, followed by a rapid reduction of free collagen that coincides with the appearance of less-soluble collagen fibrils. The model and measurements together provide evidence for extracellular self-assembly of collagen fibrils in embryonic mouse tendon, supporting an additional mechanism for rapid collagen fibril formation during embryonic development.
Collapse
Affiliation(s)
- Christopher K Revell
- Department of Mathematics, University of Manchester, Manchester, United Kingdom; Wellcome Centre for Cell-Matrix Research, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Jeremy A Herrera
- Wellcome Centre for Cell-Matrix Research, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Craig Lawless
- Wellcome Centre for Cell-Matrix Research, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Yinhui Lu
- Wellcome Centre for Cell-Matrix Research, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Karl E Kadler
- Wellcome Centre for Cell-Matrix Research, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom.
| | - Joan Chang
- Wellcome Centre for Cell-Matrix Research, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom; Division of Molecular and Cellular Function, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom.
| | - Oliver E Jensen
- Department of Mathematics, University of Manchester, Manchester, United Kingdom; Wellcome Centre for Cell-Matrix Research, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom.
| |
Collapse
|
5
|
Frohoff-Hülsmann T, Holl MP, Knobloch E, Gurevich SV, Thiele U. Stationary broken parity states in active matter models. Phys Rev E 2023; 107:064210. [PMID: 37464596 DOI: 10.1103/physreve.107.064210] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2022] [Accepted: 05/15/2023] [Indexed: 07/20/2023]
Abstract
We demonstrate that several nonvariational continuum models commonly used to describe active matter as well as other active systems exhibit nongeneric behavior: each model supports asymmetric but stationary localized states even in the absence of pinning at heterogeneities. Moreover, such states only begin to drift following a drift-transcritical bifurcation as the activity increases. Asymmetric stationary states should only exist in variational systems, i.e., in models with gradient structure. In other words, such states are expected in passive systems, but not in active systems where the gradient structure of the model is broken by activity. We identify a "spurious" gradient dynamics structure of these models that is responsible for this nongeneric behavior, and determine the types of additional terms that render the models generic, i.e., with asymmetric states that appear via drift-pitchfork bifurcations and are generically moving. We provide detailed illustrations of our results using numerical continuation of resting and steadily drifting states in both generic and nongeneric cases.
Collapse
Affiliation(s)
- Tobias Frohoff-Hülsmann
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - Max Philipp Holl
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - Edgar Knobloch
- Department of Physics, University of California at Berkeley, Berkeley, California 94720, USA
| | - Svetlana V Gurevich
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
- Center for Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstrasse 2, 48149 Münster, Germany
| | - Uwe Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
- Center for Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstrasse 2, 48149 Münster, Germany
- Center for Multiscale Theory and Computation (CMTC), Westfälische Wilhelms-Universität, Corrensstrasse 40, 48149 Münster, Germany
| |
Collapse
|
6
|
Frohoff-Hülsmann T, Thiele U, Pismen LM. Non-reciprocity induces resonances in a two-field Cahn-Hilliard model. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220087. [PMID: 36842986 DOI: 10.1098/rsta.2022.0087] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Accepted: 01/09/2023] [Indexed: 06/18/2023]
Abstract
We consider a non-reciprocally coupled two-field Cahn-Hilliard system that has been shown to allow for oscillatory behaviour and suppression of coarsening. After introducing the model, we first review the linear stability of steady uniform states and show that all instability thresholds are identical to the ones for a corresponding two-species reaction-diffusion system. Next, we consider a specific interaction of linear modes-a 'Hopf-Turing' resonance-and derive the corresponding amplitude equations using a weakly nonlinear approach. We discuss the weakly nonlinear results and finally compare them with fully nonlinear simulations for a specific conserved amended FitzHugh-Nagumo system. We conclude with a discussion of the limitations of the employed weakly nonlinear approach. This article is part of the theme issue 'New trends in pattern formation and nonlinear dynamics of extended systems'.
Collapse
Affiliation(s)
- Tobias Frohoff-Hülsmann
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Str. 9, Münster 48149, Germany
| | - Uwe Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Str. 9, Münster 48149, Germany
- Center for Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstr. 2, Münster 48149, Germany
| | - Len M Pismen
- Department of Chemical Engineering, Technion-Israel Institute of Technology, Haifa 32000, Israel
| |
Collapse
|
7
|
Abstract
Matter self-assembling into layers generates unique properties, including structures of stacked surfaces, directed transport, and compact area maximization that can be highly functionalized in biology and technology. Smectics represent the paradigm of such lamellar materials - they are a state between fluids and solids, characterized by both orientational and partial positional ordering in one layering direction, making them notoriously difficult to model, particularly in confining geometries. We propose a complex tensor order parameter to describe the local degree of lamellar ordering, layer displacement and orientation of the layers for simple, lamellar smectics. The theory accounts for both dislocations and disclinations, by regularizing singularities within defect cores and so remaining continuous everywhere. The ability to describe disclinations and dislocation allows this theory to simulate arrested configurations and inclusion-induced local ordering. This tensorial theory for simple smectics considerably simplifies numerics, facilitating studies on the mesoscopic structure of topologically complex systems.
Collapse
|
8
|
Steinberg AB, Maucher F, Gurevich SV, Thiele U. Exploring bifurcations in Bose-Einstein condensates via phase field crystal models. CHAOS (WOODBURY, N.Y.) 2022; 32:113112. [PMID: 36456347 DOI: 10.1063/5.0101401] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 10/03/2022] [Indexed: 06/17/2023]
Abstract
To facilitate the analysis of pattern formation and the related phase transitions in Bose-Einstein condensates, we present an explicit approximate mapping from the nonlocal Gross-Pitaevskii equation with cubic nonlinearity to a phase field crystal (PFC) model. This approximation is valid close to the superfluid-supersolid phase transition boundary. The simplified PFC model permits the exploration of bifurcations and phase transitions via numerical path continuation employing standard software. While revealing the detailed structure of the bifurcations present in the system, we demonstrate the existence of localized states in the PFC approximation. Finally, we discuss how higher-order nonlinearities change the structure of the bifurcation diagram representing the transitions found in the system.
Collapse
Affiliation(s)
- A B Steinberg
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - F Maucher
- Departament de Física, Universitat de les Illes Balears and IAC-3, Campus UIB, E-07122 Palma de Mallorca, Spain
| | - S V Gurevich
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - U Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| |
Collapse
|
9
|
Clerc MG, Echeverría-Alar S, Tlidi M. Localized states with nontrivial symmetries: Localized labyrinthine patterns. Phys Rev E 2022; 105:L012202. [PMID: 35193233 DOI: 10.1103/physreve.105.l012202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 01/11/2022] [Indexed: 06/14/2023]
Abstract
The formation of self-organized patterns and localized states are ubiquitous in Nature. Localized states containing trivial symmetries such as stripes, hexagons, or squares have been profusely studied. Disordered patterns with nontrivial symmetries such as labyrinthine patterns are observed in different physical contexts. Here we report stable localized disordered patterns in spatially extended dissipative systems. These two- and three-dimensional localized structures consist of an isolated labyrinth embedded in a homogeneous steady state. Their partial bifurcation diagram allows us to explain this phenomenon as a manifestation of a pinning-depinning transition. We illustrate our findings on the Swift-Hohenberg-type of equations and other well-established models for plant ecology, nonlinear optics, and reaction-diffusion systems.
Collapse
Affiliation(s)
- M G Clerc
- Departamento de Física and Millennium Institute for Research in Optics, FCFM, Universidad de Chile, Casilla 487-3, Santiago, Chile
| | - S Echeverría-Alar
- Departamento de Física and Millennium Institute for Research in Optics, FCFM, Universidad de Chile, Casilla 487-3, Santiago, Chile
| | - M Tlidi
- Faculté des Sciences, Université libre de Bruxelles (U.L.B), CP. 231, 1050 Brussels, Belgium
| |
Collapse
|
10
|
Verschueren N, Knobloch E, Uecker H. Localized and extended patterns in the cubic-quintic Swift-Hohenberg equation on a disk. Phys Rev E 2021; 104:014208. [PMID: 34412325 DOI: 10.1103/physreve.104.014208] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Accepted: 06/16/2021] [Indexed: 11/07/2022]
Abstract
Axisymmetric and nonaxisymmetric patterns in the cubic-quintic Swift-Hohenberg equation posed on a disk with Neumann boundary conditions are studied via numerical continuation and bifurcation analysis. Axisymmetric localized solutions in the form of spots and rings known from earlier studies persist and snake in the usual fashion until they begin to interact with the boundary. Depending on parameters, including the disk radius, these states may or may not connect to the branch of domain-filling target states. Secondary instabilities of localized axisymmetric states may create multiarm localized structures that grow and interact with the boundary before broadening into domain-filling states. High azimuthal wave number wall states referred to as daisy states are also found. Secondary bifurcations from these states include localized daisies, i.e., wall states localized in both radius and angle. Depending on parameters, these states may snake much as in the one-dimensional Swift-Hohenberg equation, or invade the interior of the domain, yielding states referred to as worms, or domain-filling stripes.
Collapse
Affiliation(s)
- Nicolás Verschueren
- Physics Department, University of California at Berkeley, Berkeley, California 94720, USA
| | - Edgar Knobloch
- Physics Department, University of California at Berkeley, Berkeley, California 94720, USA
| | - Hannes Uecker
- Institute for Mathematics, Carl von Ossietzky University of Oldenburg, Oldenburg, Germany
| |
Collapse
|
11
|
Vitral E, Leo PH, Viñals J. Phase-field model for a weakly compressible soft layered material: morphological transitions on smectic-isotropic interfaces. SOFT MATTER 2021; 17:6140-6159. [PMID: 34080608 DOI: 10.1039/d1sm00488c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
A coupled phase-field and hydrodynamic model is introduced to describe a two-phase, weakly compressible smectic (layered phase) in contact with an isotropic fluid of different density. A non-conserved smectic order parameter is coupled to a conserved mass density in order to accommodate non-solenoidal flows near the smectic-isotropic boundary arising from density contrast between the two phases. The model aims to describe morphological transitions in smectic thin films under heat treatment, in which arrays of focal conic defects evolve into conical pyramids and concentric rings through curvature dependent evaporation of smectic layers. The model leads to an extended thermodynamic relation at a curved surface that includes its Gaussian curvature, non-classical stresses at the boundary and flows arising from density gradients. The temporal evolution given by the model conserves the overall mass of the liquid crystal while still allowing for the modulated smectic structure to grow or shrink. A numerical solution of the governing equations reveals that pyramidal domains are sculpted at the center of focal conics upon a temperature increase, which display tangential flows at their surface. Other cases investigated include the possible coalescence of two cylindrical stacks of smectic layers, formation of droplets, and the interactions between focal conic domains through flow.
Collapse
Affiliation(s)
- Eduardo Vitral
- Department of Mechanical Engineering, University of Nevada, 1664 N. Virginia St., Reno, NV 89557, USA.
| | - Perry H Leo
- Department of Aerospace Engineering and Mechanics, University of Minnesota, 110 Union St. SE, Minneapolis, MN 55455, USA
| | - Jorge Viñals
- School of Physics and Astronomy, University of Minnesota, 116 Church St. SE, Minneapolis, MN 55455, USA
| |
Collapse
|
12
|
Frohoff-Hülsmann T, Wrembel J, Thiele U. Suppression of coarsening and emergence of oscillatory behavior in a Cahn-Hilliard model with nonvariational coupling. Phys Rev E 2021; 103:042602. [PMID: 34006003 DOI: 10.1103/physreve.103.042602] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2020] [Accepted: 03/04/2021] [Indexed: 12/27/2022]
Abstract
We investigate a generic two-field Cahn-Hilliard model with variational and nonvariational coupling. It describes, for instance, passive and active ternary mixtures, respectively. Already a linear stability analysis of the homogeneous mixed state shows that activity not only allows for the usual large-scale stationary (Cahn-Hilliard) instability of the well-known passive case but also for small-scale stationary (Turing) and large-scale oscillatory (Hopf) instabilities. In consequence of the Turing instability, activity may completely suppress the usual coarsening dynamics. In a fully nonlinear analysis, we first briefly discuss the passive case before focusing on the active case. Bifurcation diagrams and selected direct time simulations are presented that allow us to establish that nonvariational coupling (i) can partially or completely suppress coarsening and (ii) may lead to the emergence of drifting and oscillatory states. Throughout, we emphasize the relevance of conservation laws and related symmetries for the encountered intricate bifurcation behavior.
Collapse
Affiliation(s)
- Tobias Frohoff-Hülsmann
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Str. 9, 48149 Münster, Germany
| | - Jana Wrembel
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Str. 9, 48149 Münster, Germany
| | - Uwe Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Str. 9, 48149 Münster, Germany.,Center for Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstr. 2, 48149 Münster, Germany.,Center for Multiscale Theory and Computation (CMTC), Westfälische Wilhelms-Universität, Corrensstr. 40, 48149 Münster, Germany
| |
Collapse
|
13
|
Ophaus L, Knobloch E, Gurevich SV, Thiele U. Two-dimensional localized states in an active phase-field-crystal model. Phys Rev E 2021; 103:032601. [PMID: 33862772 DOI: 10.1103/physreve.103.032601] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2020] [Accepted: 02/08/2021] [Indexed: 11/07/2022]
Abstract
The active phase-field-crystal (active PFC) model provides a simple microscopic mean field description of crystallization in active systems. It combines the PFC model (or conserved Swift-Hohenberg equation) of colloidal crystallization and aspects of the Toner-Tu theory for self-propelled particles. We employ the active PFC model to study the occurrence of localized and periodic active crystals in two spatial dimensions. Due to the activity, crystalline states can undergo a drift instability and start to travel while keeping their spatial structure. Based on linear stability analyses, time simulations, and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of resting and traveling states. We explore, for instance, how the slanted homoclinic snaking of steady localized states found for the passive PFC model is modified by activity. Morphological phase diagrams showing the regions of existence of various solution types are presented merging the results from all the analysis tools employed. We also study how activity influences the crystal structure with transitions from hexagons to rhombic and stripe patterns. This in-depth analysis of a simple PFC model for active crystals and swarm formation provides a clear general understanding of the observed multistability and associated hysteresis effects, and identifies thresholds for qualitative changes in behavior.
Collapse
Affiliation(s)
- Lukas Ophaus
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany.,Center of Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstrasse 2, 48149 Münster, Germany
| | - Edgar Knobloch
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Svetlana V Gurevich
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany.,Center of Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstrasse 2, 48149 Münster, Germany
| | - Uwe Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany.,Center of Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstrasse 2, 48149 Münster, Germany
| |
Collapse
|
14
|
Chatterjee P, Goldenfeld N. Field-theoretic model for chemotaxis in run and tumble particles. Phys Rev E 2021; 103:032603. [PMID: 33862765 DOI: 10.1103/physreve.103.032603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2020] [Accepted: 02/12/2021] [Indexed: 11/07/2022]
Abstract
In this paper, we develop a field-theoretic description for run and tumble chemotaxis, based on a density-functional description of crystalline materials modified to capture orientational ordering. We show that this framework, with its in-built multiparticle interactions, soft-core repulsion, and elasticity, is ideal for describing continuum collective phases with particle resolution, but on diffusive timescales. We show that our model exhibits particle aggregation in an externally imposed constant attractant field, as is observed for phototactic or thermotactic agents. We also show that this model captures particle aggregation through self-chemotaxis, an important mechanism that aids quorum-dependent cellular interactions.
Collapse
Affiliation(s)
- Purba Chatterjee
- Department of Physics, University of Illinois at Urbana-Champaign, Loomis Laboratory of Physics, 1110 West Green Street, Urbana, Illinois, 61801-3080, USA
| | - Nigel Goldenfeld
- Department of Physics, University of Illinois at Urbana-Champaign, Loomis Laboratory of Physics, 1110 West Green Street, Urbana, Illinois, 61801-3080, USA
| |
Collapse
|
15
|
Holl MP, Archer AJ, Thiele U. Efficient calculation of phase coexistence and phase diagrams: application to a binary phase-field-crystal model. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:115401. [PMID: 33246322 DOI: 10.1088/1361-648x/abce6e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We show that one can employ well-established numerical continuation methods to efficiently calculate the phase diagram for thermodynamic systems described by a suitable free energy functional. In particular, this involves the determination of lines of phase coexistence related to first order phase transitions and the continuation of triple points. To illustrate the method we apply it to a binary phase-field-crystal model for the crystallisation of a mixture of two types of particles. The resulting phase diagram is determined for one- and two-dimensional domains. In the former case it is compared to the diagram obtained from a one-mode approximation. The various observed liquid and crystalline phases and their stable and metastable coexistence are discussed as well as the temperature-dependence of the phase diagrams. This includes the (dis)appearance of critical points and triple points. We also relate bifurcation diagrams for finite-size systems to the thermodynamics of phase transitions in the infinite-size limit.
Collapse
Affiliation(s)
- Max Philipp Holl
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm Klemm Str. 9, 48149 Münster, Germany
| | - Andrew J Archer
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
- Interdisciplinary Centre for Mathematical Modelling, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Uwe Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm Klemm Str. 9, 48149 Münster, Germany
- Center of Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster, Corrensstr. 2, 48149 Münster, Germany
- Center for Multiscale Theory and Computation (CMTC), Westfälische Wilhelms-Universität, Corrensstr. 40, 48149 Münster, Germany
| |
Collapse
|
16
|
Ophaus L, Kirchner J, Gurevich SV, Thiele U. Phase-field-crystal description of active crystallites: Elastic and inelastic collisions. CHAOS (WOODBURY, N.Y.) 2020; 30:123149. [PMID: 33380045 DOI: 10.1063/5.0019426] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/23/2020] [Accepted: 12/07/2020] [Indexed: 06/12/2023]
Abstract
The active Phase-Field-Crystal (aPFC) model combines elements of the Toner-Tu theory for self-propelled particles and the classical Phase-Field-Crystal (PFC) model that describes the transition between liquid and crystalline phases. In the liquid-crystal coexistence region of the PFC model, crystalline clusters exist in the form of localized states that coexist with a homogeneous background. At sufficiently strong activity (related to self-propulsion strength), they start to travel. We employ numerical path continuation and direct time simulations to first investigate the existence regions of different types of localized states in one spatial dimension. The results are summarized in morphological phase diagrams in the parameter plane spanned by activity and mean density. Then we focus on the interaction of traveling localized states, studying their collision behavior. As a result, we distinguish "elastic" and "inelastic" collisions. In the former, localized states recover their properties after a collision, while in the latter, they may completely or partially annihilate, forming resting bound states or various traveling states.
Collapse
Affiliation(s)
- Lukas Ophaus
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - Johannes Kirchner
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - Svetlana V Gurevich
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| | - Uwe Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany
| |
Collapse
|
17
|
Beta C, Gov NS, Yochelis A. Why a Large-Scale Mode Can Be Essential for Understanding Intracellular Actin Waves. Cells 2020; 9:cells9061533. [PMID: 32585983 PMCID: PMC7349605 DOI: 10.3390/cells9061533] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2020] [Revised: 06/16/2020] [Accepted: 06/18/2020] [Indexed: 01/18/2023] Open
Abstract
During the last decade, intracellular actin waves have attracted much attention due to their essential role in various cellular functions, ranging from motility to cytokinesis. Experimental methods have advanced significantly and can capture the dynamics of actin waves over a large range of spatio-temporal scales. However, the corresponding coarse-grained theory mostly avoids the full complexity of this multi-scale phenomenon. In this perspective, we focus on a minimal continuum model of activator-inhibitor type and highlight the qualitative role of mass conservation, which is typically overlooked. Specifically, our interest is to connect between the mathematical mechanisms of pattern formation in the presence of a large-scale mode, due to mass conservation, and distinct behaviors of actin waves.
Collapse
Affiliation(s)
- Carsten Beta
- Institute of Physics and Astronomy, University of Potsdam, 14476 Potsdam, Germany;
| | - Nir S. Gov
- Department of Chemical and Biological Physics, Weizmann Institute of Science, Rehovot 76100, Israel;
| | - Arik Yochelis
- Department of Solar Energy and Environmental Physics, Blaustein Institutes for Desert Research (BIDR), Ben-Gurion University of the Negev, Sede Boqer Campus, Midreshet Ben-Gurion 8499000, Israel
- Department of Physics, Ben-Gurion University of the Negev, Be’er Sheva 8410501, Israel
- Correspondence:
| |
Collapse
|
18
|
Archer AJ, Ratliff DJ, Rucklidge AM, Subramanian P. Deriving phase field crystal theory from dynamical density functional theory: Consequences of the approximations. Phys Rev E 2019; 100:022140. [PMID: 31574721 DOI: 10.1103/physreve.100.022140] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2019] [Indexed: 06/10/2023]
Abstract
Phase field crystal (PFC) theory is extensively used for modeling the phase behavior, structure, thermodynamics, and other related properties of solids. PFC theory can be derived from dynamical density functional theory (DDFT) via a sequence of approximations. Here, we carefully identify all of these approximations and explain the consequences of each. One approximation that is made in standard derivations is to neglect a term of form ∇·[n∇Ln], where n is the scaled density profile and L is a linear operator. We show that this term makes a significant contribution to the stability of the crystal, and that dropping this term from the theory forces another approximation, that of replacing the logarithmic term from the ideal gas contribution to the free energy with its truncated Taylor expansion, to yield a polynomial in n. However, the consequences of doing this are (i) the presence of an additional spinodal in the phase diagram, so the liquid is predicted first to freeze and then to melt again as the density is increased; and (ii) other periodic structures, such as stripes, are erroneously predicted to be thermodynamic equilibrium structures. In general, L consists of a nonlocal convolution involving the pair direct correlation function. A second approximation sometimes made in deriving PFC theory is to replace L with a gradient expansion involving derivatives. We show that this leads to the possibility of the density going to zero, with its logarithm going to -∞ while being balanced by the fourth derivative of the density going to +∞. This subtle singularity leads to solutions failing to exist above a certain value of the average density. We illustrate all of these conclusions with results for a particularly simple model two-dimensional fluid, the generalized exponential model of index 4 (GEM-4), chosen because a DDFT is known to be accurate for this model. The consequences of the subsequent PFC approximations can then be examined. These include the phase diagram being both qualitatively incorrect, in that it has a stripe phase, and quantitatively incorrect (by orders of magnitude) regarding the properties of the crystal phase. Thus, although PFC models are very successful as phenomenological models of crystallization, we find it impossible to derive the PFC as a theory for the (scaled) density distribution when starting from an accurate DDFT, without introducing spurious artifacts. However, we find that making a simple one-mode approximation for the logarithm of the density distribution lnρ(x) rather than for ρ(x) is surprisingly accurate. This approach gives a tantalizing hint that accurate PFC-type theories may instead be derived as theories for the field lnρ(x), rather than for the density profile itself.
Collapse
Affiliation(s)
- Andrew J Archer
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Daniel J Ratliff
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | | | - Priya Subramanian
- School of Mathematics, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
19
|
Gomez H, Bures M, Moure A. A review on computational modelling of phase-transition problems. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2019; 377:20180203. [PMID: 30827214 PMCID: PMC6460062 DOI: 10.1098/rsta.2018.0203] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 12/06/2018] [Indexed: 05/18/2023]
Abstract
Phase-transition problems are ubiquitous in science and engineering. They have been widely studied via theory, experiments and computations. This paper reviews the main challenges associated with computational modelling of phase-transition problems, addressing both model development and numerical discretization of the resulting equations. We focus on classical phase-transition problems, including liquid-solid, gas-liquid and solid-solid transformations. Our review has a strong emphasis on the treatment of interfacial phenomena and the phase-field method. This article is part of the theme issue 'Heterogeneous materials: metastable and non-ergodic internal structures'.
Collapse
|
20
|
Continuation for Thin Film Hydrodynamics and Related Scalar Problems. COMPUTATIONAL METHODS IN APPLIED SCIENCES 2019. [DOI: 10.1007/978-3-319-91494-7_13] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
|
21
|
Ophaus L, Gurevich SV, Thiele U. Resting and traveling localized states in an active phase-field-crystal model. Phys Rev E 2018; 98:022608. [PMID: 30253633 DOI: 10.1103/physreve.98.022608] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2018] [Indexed: 06/08/2023]
Abstract
The conserved Swift-Hohenberg equation (or phase-field-crystal [PFC] model) provides a simple microscopic description of the thermodynamic transition between fluid and crystalline states. Combining it with elements of the Toner-Tu theory for self-propelled particles, Menzel and Löwen [Phys. Rev. Lett. 110, 055702 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.055702] obtained a model for crystallization (swarm formation) in active systems. Here, we study the occurrence of resting and traveling localized states, i.e., crystalline clusters, within the resulting active PFC model. Based on linear stability analyses and numerical continuation of the fully nonlinear states, we present a detailed analysis of the bifurcation structure of periodic and localized, resting and traveling states in a one-dimensional active PFC model. This allows us, for instance, to explore how the slanted homoclinic snaking of steady localized states found for the passive PFC model is amended by activity. A particular focus lies on the onset of motion, where we show that it occurs either through a drift-pitchfork or a drift-transcritical bifurcation. A corresponding general analytical criterion is derived.
Collapse
Affiliation(s)
- Lukas Ophaus
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany, and Center of Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster Corrensstrasse 2, 48149 Münster, Germany
| | - Svetlana V Gurevich
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany, and Center of Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster Corrensstrasse 2, 48149 Münster, Germany
| | - Uwe Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster Wilhelm-Klemm-Strasse 9, 48149 Münster, Germany, and Center of Nonlinear Science (CeNoS), Westfälische Wilhelms-Universität Münster Corrensstrasse 2, 48149 Münster, Germany
| |
Collapse
|
22
|
Yochelis A. Catalytic Membrane Reactor Model as a Laboratory for Pattern Emergence in Reaction-diffusion-advection Media. Isr J Chem 2018. [DOI: 10.1002/ijch.201700145] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Arik Yochelis
- Department of Solar Energy and Environmental Physics, Swiss Institute for Dryland Environmental and Energy Research (SIDEER), Blaustein Institutes for Desert Research (BIDR); Ben-Gurion University of the Negev, Sede Boqer Campus; Midreshet Ben-Gurion 8499000 Israel
- Department of Physics; Ben-Gurion University of the Negev; Be'er Sheva 8410501 Israel
| |
Collapse
|
23
|
Bier S, Gavish N, Uecker H, Yochelis A. From bulk self-assembly to electrical diffuse layer in a continuum approach for ionic liquids: The impact of anion and cation size asymmetry. Phys Rev E 2017; 95:060201. [PMID: 28709248 DOI: 10.1103/physreve.95.060201] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2016] [Indexed: 06/07/2023]
Abstract
Ionic liquids are solvent-free electrolytes, some of which possess an intriguing self-assembly at finite length scale due to Coulombic interactions. Using a continuum framework (based on Onsager's relations), it is shown that bulk nanostructures arise via linear (supercritical) and nonlinear (subcritical) bifurcations (morphological phase transitions), which also directly affect the electrical double layer structure. A Ginzburg-Landau amplitude equation is derived and the bifurcation type is related to model parameters, such as temperature, potential, and interactions. Specifically, the nonlinear bifurcation occurs for geometrically dissimilar ions and, surprisingly, is induced by perturbations on the order of thermal fluctuations. Finally, qualitative insights and comparisons to the experimentally decaying charge layers within the electrical double layer are discussed.
Collapse
Affiliation(s)
- Sariel Bier
- Department of Solar Energy and Environmental Physics, Swiss Institute for Dryland Environmental and Energy Research, Blaustein Institutes for Desert Research (BIDR), Ben-Gurion University of the Negev, Sede Boqer Campus, 8499000 Midreshet Ben-Gurion, Israel
| | - Nir Gavish
- Department of Mathematics, Technion-IIT, 3200003 Haifa, Israel
| | - Hannes Uecker
- Institute for Mathematics, Carl von Ossietzky University of Oldenburg, P.F 2503, 26111 Oldenburg, Germany
| | - Arik Yochelis
- Department of Solar Energy and Environmental Physics, Swiss Institute for Dryland Environmental and Energy Research, Blaustein Institutes for Desert Research (BIDR), Ben-Gurion University of the Negev, Sede Boqer Campus, 8499000 Midreshet Ben-Gurion, Israel
| |
Collapse
|
24
|
Tabbert F, Schelte C, Tlidi M, Gurevich SV. Delay-induced depinning of localized structures in a spatially inhomogeneous Swift-Hohenberg model. Phys Rev E 2017; 95:032213. [PMID: 28415377 DOI: 10.1103/physreve.95.032213] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2016] [Indexed: 06/07/2023]
Abstract
We report on the dynamics of localized structures in an inhomogeneous Swift-Hohenberg model describing pattern formation in the transverse plane of an optical cavity. This real order parameter equation is valid close to the second-order critical point associated with bistability. The optical cavity is illuminated by an inhomogeneous spatial Gaussian pumping beam and subjected to time-delayed feedback. The Gaussian injection beam breaks the translational symmetry of the system by exerting an attracting force on the localized structure. We show that the localized structure can be pinned to the center of the inhomogeneity, suppressing the delay-induced drift bifurcation that has been reported in the particular case where the injection is homogeneous, assuming a continuous wave operation. Under an inhomogeneous spatial pumping beam, we perform the stability analysis of localized solutions to identify different instability regimes induced by time-delayed feedback. In particular, we predict the formation of two-arm spirals, as well as oscillating and depinning dynamics caused by the interplay of an attracting inhomogeneity and destabilizing time-delayed feedback. The transition from oscillating to depinning solutions is investigated by means of numerical continuation techniques. Analytically, we use an order parameter approach to derive a normal form of the delay-induced Hopf bifurcation leading to an oscillating solution. Additionally we model the interplay of an attracting inhomogeneity and destabilizing time delay by describing the localized solution as an overdamped particle in a potential well generated by the inhomogeneity. In this case, the time-delayed feedback acts as a driving force. Comparing results from the later approach with the full Swift-Hohenberg model, we show that the approach not only provides an instructive description of the depinning dynamics, but also is numerically accurate throughout most of the parameter regime.
Collapse
Affiliation(s)
- Felix Tabbert
- Institute for Theoretical Physics, University of Münster, Wilhelm-Klemm-Strasse 9, D-48149 Münster, Germany
| | - Christian Schelte
- Institute for Theoretical Physics, University of Münster, Wilhelm-Klemm-Strasse 9, D-48149 Münster, Germany
| | - Mustapha Tlidi
- Faculté des Sciences, Université Libre de Bruxelles, Campus Plaine, C.P. 231, Brussels B-1050, Belgium
| | - Svetlana V Gurevich
- Institute for Theoretical Physics, University of Münster, Wilhelm-Klemm-Strasse 9, D-48149 Münster, Germany
- Center for Nonlinear Science (CeNoS), University of Münster, Corrensstrasse 2, D-48149 Münster, Germany
| |
Collapse
|
25
|
Gustave F, Radwell N, McIntyre C, Toomey JP, Kane DM, Barland S, Firth WJ, Oppo GL, Ackemann T. Observation of Mode-Locked Spatial Laser Solitons. PHYSICAL REVIEW LETTERS 2017; 118:044102. [PMID: 28186807 DOI: 10.1103/physrevlett.118.044102] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/27/2016] [Indexed: 06/06/2023]
Abstract
A stable nonlinear wave packet, self-localized in all three dimensions, is an intriguing and much sought after object in nonlinear science in general and in nonlinear photonics in particular. We report on the experimental observation of mode-locked spatial laser solitons in a vertical-cavity surface-emitting laser with frequency-selective feedback from an external cavity. These spontaneously emerging and long-term stable spatiotemporal structures have a pulse length shorter than the cavity round-trip time and may pave the way to completely independent cavity light bullets.
Collapse
Affiliation(s)
- F Gustave
- Université Côte d'Azur, CNRS, 06560 Valbonne, France
| | - N Radwell
- SUPA and Department of Physics and Astronomy, University of Glasgow, Glasgow G12 8QQ, United Kingdom
| | - C McIntyre
- SUPA and Department of Physics, University of Strathclyde, Glasgow G4 ONG, Scotland, United Kingdom
| | - J P Toomey
- MQ Photonics Research Centre, Department of Physics & Astronomy, Macquarie University, Sydney 2109, Australia
| | - D M Kane
- MQ Photonics Research Centre, Department of Physics & Astronomy, Macquarie University, Sydney 2109, Australia
| | - S Barland
- Université Côte d'Azur, CNRS, 06560 Valbonne, France
| | - W J Firth
- SUPA and Department of Physics, University of Strathclyde, Glasgow G4 ONG, Scotland, United Kingdom
| | - G-L Oppo
- SUPA and Department of Physics, University of Strathclyde, Glasgow G4 ONG, Scotland, United Kingdom
| | - T Ackemann
- SUPA and Department of Physics, University of Strathclyde, Glasgow G4 ONG, Scotland, United Kingdom
| |
Collapse
|
26
|
Beltrame P, Makhoul M, Joelson M. Deterministic particle transport in a ratchet flow. Phys Rev E 2016; 93:012208. [PMID: 26871074 DOI: 10.1103/physreve.93.012208] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2013] [Indexed: 11/07/2022]
Abstract
This study is motivated by the issue of the pumping of particle through a periodic modulated channel. We focus on a simplified deterministic model of small inertia particles within the Stokes flow framework that we call "ratchet flow." A path-following method is employed in the parameter space in order to retrace the scenario which from bounded periodic solutions leads to particle transport. Depending on whether the magnitude of the particle drag is moderate or large, two main transport mechanisms are identified in which the role of the parity symmetry of the flow differs. For large drag, transport is induced by flow asymmetry, while for moderate drag, since the full transport solution bifurcation structure already exists for symmetric settings, flow asymmetry only makes the transport effective. We analyzed the scenarios of current reversals for each mechanism as well as the role of synchronization. In particular we show that, for large drag, the particle drift is similar to phase slip in a synchronization problem.
Collapse
Affiliation(s)
- Philippe Beltrame
- Université d'Avignon et des Pays de Vaucluse-INRA, UMR1114 EMMAH, F-84914 Avignon, France
| | - Mounia Makhoul
- Université d'Avignon et des Pays de Vaucluse-INRA, UMR1114 EMMAH, F-84914 Avignon, France
| | - Maminirina Joelson
- Université d'Avignon et des Pays de Vaucluse-INRA, UMR1114 EMMAH, F-84914 Avignon, France
| |
Collapse
|
27
|
Praetorius S, Voigt A. A Navier-Stokes phase-field crystal model for colloidal suspensions. J Chem Phys 2015; 142:154904. [DOI: 10.1063/1.4918559] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Simon Praetorius
- Institute of Scientific Computing, Technische Universität Dresden, D-01062 Dresden, Germany
| | - Axel Voigt
- Institute of Scientific Computing, Technische Universität Dresden, D-01062 Dresden, Germany
| |
Collapse
|
28
|
Schüler D, Alonso S, Torcini A, Bär M. Spatio-temporal dynamics induced by competing instabilities in two asymmetrically coupled nonlinear evolution equations. CHAOS (WOODBURY, N.Y.) 2014; 24:043142. [PMID: 25554062 DOI: 10.1063/1.4905017] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Pattern formation often occurs in spatially extended physical, biological, and chemical systems due to an instability of the homogeneous steady state. The type of the instability usually prescribes the resulting spatio-temporal patterns and their characteristic length scales. However, patterns resulting from the simultaneous occurrence of instabilities cannot be expected to be simple superposition of the patterns associated with the considered instabilities. To address this issue, we design two simple models composed by two asymmetrically coupled equations of non-conserved (Swift-Hohenberg equations) or conserved (Cahn-Hilliard equations) order parameters with different characteristic wave lengths. The patterns arising in these systems range from coexisting static patterns of different wavelengths to traveling waves. A linear stability analysis allows to derive a two parameter phase diagram for the studied models, in particular, revealing for the Swift-Hohenberg equations, a co-dimension two bifurcation point of Turing and wave instability and a region of coexistence of stationary and traveling patterns. The nonlinear dynamics of the coupled evolution equations is investigated by performing accurate numerical simulations. These reveal more complex patterns, ranging from traveling waves with embedded Turing patterns domains to spatio-temporal chaos, and a wide hysteretic region, where waves or Turing patterns coexist. For the coupled Cahn-Hilliard equations the presence of a weak coupling is sufficient to arrest the coarsening process and to lead to the emergence of purely periodic patterns. The final states are characterized by domains with a characteristic length, which diverges logarithmically with the coupling amplitude.
Collapse
Affiliation(s)
- D Schüler
- Physikalisch-Technische Bundesanstalt, Abbestrasse 2-12, 10587 Berlin, Germany
| | - S Alonso
- Physikalisch-Technische Bundesanstalt, Abbestrasse 2-12, 10587 Berlin, Germany
| | - A Torcini
- CNR-Consiglio Nazionale delle Ricerche, Istituto dei Sistemi Complessi - Via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - M Bär
- Physikalisch-Technische Bundesanstalt, Abbestrasse 2-12, 10587 Berlin, Germany
| |
Collapse
|
29
|
Archer AJ, Walters MC, Thiele U, Knobloch E. Solidification in soft-core fluids: Disordered solids from fast solidification fronts. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:042404. [PMID: 25375507 DOI: 10.1103/physreve.90.042404] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2014] [Indexed: 06/04/2023]
Abstract
Using dynamical density functional theory we calculate the speed of solidification fronts advancing into a quenched two-dimensional model fluid of soft-core particles. We find that solidification fronts can advance via two different mechanisms, depending on the depth of the quench. For shallow quenches, the front propagation is via a nonlinear mechanism. For deep quenches, front propagation is governed by a linear mechanism and in this regime we are able to determine the front speed via a marginal stability analysis. We find that the density modulations generated behind the advancing front have a characteristic scale that differs from the wavelength of the density modulation in thermodynamic equilibrium, i.e., the spacing between the crystal planes in an equilibrium crystal. This leads to the subsequent development of disorder in the solids that are formed. In a one-component fluid, the particles are able to rearrange to form a well-ordered crystal, with few defects. However, solidification fronts in a binary mixture exhibiting crystalline phases with square and hexagonal ordering generate solids that are unable to rearrange after the passage of the solidification front and a significant amount of disorder remains in the system.
Collapse
Affiliation(s)
- A J Archer
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - M C Walters
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - U Thiele
- Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm Klemm Str. 9, D-48149 Münster, Germany and Center of Nonlinear Science (CeNoS), Westfälische Wilhelms Universität Münster, Corrensstr. 2, 48149 Münster, Germany
| | - E Knobloch
- Department of Physics, University of California at Berkeley, Berkeley, California 94720, USA
| |
Collapse
|
30
|
Thiele U. Patterned deposition at moving contact lines. Adv Colloid Interface Sci 2014; 206:399-413. [PMID: 24331374 DOI: 10.1016/j.cis.2013.11.002] [Citation(s) in RCA: 63] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2013] [Revised: 11/07/2013] [Accepted: 11/07/2013] [Indexed: 01/09/2023]
Abstract
When a simple or complex liquid recedes from a smooth solid substrate it often leaves a homogeneous or structured deposit behind. In the case of a receding non-volatile pure liquid the deposit might be a liquid film or an arrangement of droplets depending on the receding speed of the meniscus and the wetting properties of the system. For complex liquids with volatile components as, e.g., polymer solutions and particle or surfactant suspensions, the deposit might be a homogeneous or structured layer of solute--with structures ranging from line patterns that can be orthogonal or parallel to the receding contact line via hexagonal or square arrangements of drops to complicated hierarchical structures. We review a number of recent experiments and modelling approaches with a particular focus on mesoscopic hydrodynamic long-wave models. The conclusion highlights open question and speculates about future developments.
Collapse
Affiliation(s)
- Uwe Thiele
- Department of Mathematical Sciences, Loughborough University, Loughborough, Leicestershire LE11 3TU, UK; Institut für Theoretische Physik, Westfälische Wilhelms-Universität Münster, Wilhelm Klemm Str. 9, D-48149 Münster, Germany.
| |
Collapse
|
31
|
Affiliation(s)
- Shai Cohen
- Raymond
and Beverly Sackler
School of Physics and Astronomy, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
| | - David Andelman
- Raymond
and Beverly Sackler
School of Physics and Astronomy, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
| |
Collapse
|