Finsler geometry modeling and Monte Carlo study of liquid crystal elastomers under electric fields

Numerical study of anisotropic diffusion in Turing patterns based on Finsler geometry modeling

We numerically study the anisotropic Turing patterns (TPs) of an activator-inhibitor system described by the reaction-diffusion (RD) equation of Turing, focusing on anisotropic diffusion using the Finsler geometry (FG) modeling technique. In FG modeling, the diffusion coefficients are dynamically generated to be direction dependent owing to an internal degree of freedom (IDOF) and its interaction with the activator and inhibitor. Because of this dynamical diffusion coefficient, FG modeling of the RD equation sharply contrasts with the standard numerical technique in which direction-dependent coefficients are manually assumed. To find the solution of the RD equations in FG modeling, we use a hybrid numerical technique combining the Metropolis Monte Carlo method for IDOF updates and discrete RD equations for steady-state configurations of the activator-inhibitor variables. We find that the newly introduced IDOF and its interaction are a possible origin of spontaneously emergent anisotropic patterns of living organisms, such as zebra and fishes. Moreover, the IDOF makes TPs controllable by external conditions if the IDOF is identified with the direction of cell diffusion accompanied by thermal fluctuations.

Numerical Methods in Studies of Liquid Crystal Elastomers

Liquid crystal elastomers (LCEs) are a type of material with specific features of polymers and of liquid crystals. They exhibit interesting behaviors, i.e., they are able to change their physical properties when met with external stimuli, including heat, light, electric, and magnetic fields. This behavior makes LCEs a suitable candidate for a variety of applications, including, but not limited to, artificial muscles, optical devices, microscopy and imaging systems, biosensor devices, and optimization of solar energy collectors. Due to the wide range of applicability, numerical models are needed not only to further our understanding of the underlining mechanics governing LCE behavior, but also to enable the predictive modeling of their behavior under different circumstances for different applications. Given that several mainstream methods are used for LCE modeling, viz. finite element method, Monte Carlo and molecular dynamics, and the growing interest and reliance on computer modeling for predicting the opto-mechanical behavior of complex structures in real world applications, there is a need to gain a better understanding regarding their strengths and weaknesses so that the best method can be utilized for the specific application at hand. Therefore, this investigation aims to not only to present a multitude of examples on numerical studies conducted on LCEs, but also attempts at offering a concise categorization of different methods based on the desired application to act as a guide for current and future research in this field.

Parallel Tempering Monte Carlo Studies of Phase Transition of Free Boundary Planar Surfaces

We numerically study surface models defined on hexagonal disks with a free boundary. 2D surface models for planar surfaces have recently attracted interest due to the engineering applications of functional materials such as graphene and its composite with polymers. These 2D composite meta-materials are strongly influenced by external stimuli such as thermal fluctuations if they are sufficiently thin. For this reason, it is very interesting to study the shape stability/instability of thin 2D materials against thermal fluctuations. In this paper, we study three types of surface models including Landau-Ginzburg (LG) and Helfirch-Polyakov models defined on triangulated hexagonal disks using the parallel tempering Monte Carlo simulation technique. We find that the planar surfaces undergo a first-order transition between the smooth and crumpled phases in the LG model and continuous transitions in the other two models. The first-order transition is relatively weak compared to the transition on spherical surfaces already reported. The continuous nature of the transition is consistent with the reported results, although the transitions are stronger than that of the reported ones.