Structure of the hard ellipsoid fluid

Computational approach for structure generation of anisotropic particles (CASGAP) with targeted distributions of particle design and orientational order

The macroscopic properties of materials are governed by their microscopic structure which depends on the materials' composition (i.e., building blocks) and processing conditions. In many classes of synthetic, bioinspired, or natural soft and/or nanomaterials, one can find structural anisotropy in the microscopic structure due to anisotropic building blocks and/or anisotropic domains formed through the processing conditions. Experimental characterization and complementary physics-based or data-driven modeling of materials' structural anisotropy are critical for understanding structure-property relationships and enabling targeted design of materials with desired macroscopic properties. In this pursuit, to interpret experimentally obtained characterization results (e.g., scattering profiles) of soft materials with structural anisotropy using data-driven computational approaches, there is a need for creating real space three-dimensional structures of the designer soft materials with realistic physical features (e.g., dispersity in building block sizes) and anisotropy (i.e., aspect ratios of the building blocks, their orientational and positional order). These real space structures can then be used to compute and complement experimentally obtained characterization results or be used as initial configurations for physics-based simulations/calculations that can then provide training data for machine learning models. To address this need, we present a new computational approach called CASGAP - Computational Approach for Structure Generation of Anisotropic Particles - for generating any desired three dimensional real-space structure of anisotropic building blocks (modeled as particles) adhering to target distributions of particle shape, size, and positional and orientational order.

Electrochemical Impedance Characterization of Blood Cell Suspensions. Part 1: Basic Theory and Application to Two-Phase Systems

Electrochemical impedance spectra of composite materials contain information on the topological arrangement, volume fraction, and shape of particles, as well as the dielectric properties of the matrix and particles. The objective of this study is to investigate how these parameters affect the dielectric spectrum and what reliable information can be extracted from experimental data. The main attention was focused on systems with dielectric behavior similar to that of human blood. Mostly plasma and erythrocytes determine the dielectric properties of whole blood. Erythrocytes suspended in plasma can be considered as three-phase systems with single-shelled particles. A theoretical approach based on the effective medium theory is developed for calculating the effective permittivity and conductivity of three-phase composites at a wide frequency range (from 0 to 1 GHz). A finite-difference method is applied to model three-dimensional periodic structures. A special case of two-phase materials is used to demonstrate the influence of the shape and arrangement of particles on dielectric properties. Theoretical and numerical approaches are applied to two-phase composites with spherical, spheroidal and biconcave particles and are compared with each other and with published data. It is shown that two-phase composites exhibit only β-dispersion. In contrast to the quasi-static limit, the wide-bandwidth impedance spectroscopy makes it possible to distinguish between disordered and regular arrangements of spheroidal and biconcave particles. The results can be used to analyze the dielectric properties of blood, which is very promising for various medical applications. This study of two-phase composites can be further extended to three-phase composites.

Strong effect of weak charging in suspensions of anisotropic colloids

Suspensions of hard colloidal particles frequently serve as model systems in studies on fundamental aspects of phase transitions. But often colloidal particles that are considered as "hard" are in fact weakly charged. If the colloids are spherical, weak charging has only a weak effect on the structural properties of the suspension, which can be easily corrected for. However, this does not hold for anisotropic particles. We introduce a model for the interaction potential between charged ellipsoids of revolution (spheroids) based on the Derjaguin approximation of Debye-Hückel theory and present a computer simulation study on aspects of the system's structural properties and phase behaviour. In line with previous experimental observations, we find that even a weak surface charge has a strong impact on the correlation functions. A likewise strong impact is seen on the phase behaviour, in particular, we find stable cubatic order in suspensions of oblate ellipsoids.

Computer simulations of biaxial nematics

Biaxial nematic (N(b)) liquid crystals are a fascinating condensed matter phase that has baffled, for more than thirty years, scientists engaged in the challenge of demonstrating its actual existence, and which has only recently been experimentally found. During this period computer simulations of model N(b) have played an important role, both in providing the basic physical properties to be expected from these systems, and in giving clues about the molecular features essential for the thermodynamic stability of N(b) phases. However, simulation studies are expected to be even more crucial in the future for unravelling the structural features of biaxial mesogens at the molecular level, and for helping in the design and optimization of devices towards the technological deployment of N(b) materials. This review article gives an overview of the simulation work performed so far, and relying on the recent experimental findings, focuses on the still unanswered questions which will determine the future challenges in the field.

A new anisotropic soft-core model for the simulation of liquid crystal mesophases

A new anisotropic soft-core model is presented, which is suitable for the rapid simulation of liquid crystal mesophases. The potential is based on a soft spherocylinder, which can be easily tuned to favor different liquid crystal mesophases. The soft-core nature of the potential makes it suitable for long-time step molecular dynamics or dissipative particle dynamics simulations, particularly as a reference model for mesogens or as an anisotropic solvent for use in combination with atomistic models. Results are presented for two variants of the new potential, which show different mesophase behaviors. Variants of the potential can also be linked together to produce more complicated molecular structures. Here, as an example, results are provided for a model multipedal liquid crystal, which has eight liquid crystalline groups linked to a central core via semiflexible chains. Here, despite the complexity of molecular structure, the model succeeds in showing the spontaneous formation of a liquid crystal phase. The results also demonstrate that there is a very strong coupling between the internal structure of the multipedal mesogen and the molecular order of the phase, with the mesogen spontaneously undergoing major structural rearrangement at the transition to the liquid crystal phase.

Direct correlation functions of binary mixtures of hard Gaussian overlap molecules

We study the direct correlation function (DCF) of a classical fluid mixture of nonspherical molecules. The components of the mixture are two types of hard ellipsoidal molecules with different elongations, interacting through the hard Gaussian overlap (HGO) model. Two different approaches are used to calculate the DCFs of this fluid, and the results are compared. Here, the Pynn approximation [J. Chem. Phys. 60, 4579 (1974)] is extended to calculate the DCF of the binary mixtures of HGO molecules, then we use a formalism based on the weighted density functional theory introduced by Chamoux and Perera [J. Chem. Phys. 104, 1493 (1996)]. These results are fairly in agreement with each other. The pressure of this system is also calculated using the Fourier zero components of the DCF. The results are in agreement with the Monte Carlo molecular simulation.

Stokes–Einstein relations for a square-well fluid

A Stokes-Einstein relation, relating the shear viscosity eta to the self-diffusion coefficient D, is constructed for a classical fluid subject to an effective two-body intermolecular force, derived from a square-well potential, undergoing dynamics as described by a Smoluchowski equation for pair diffusion. The time correlation functions for eta and 1D are separated into contributions from delta function, hard-sphere forces, and from delta function, square-well soft forces. Furthermore, D is separated into its two- and three-body time correlation functions, and eta into its two- to four-body terms. D shows activated diffusion, as in Arrhenius behavior, and on the level of two-body dynamics, the Deta product adheres to the Stokes-Einstein relation, subject to a small correction for potential softness. Three-body time correlation functions increase D, whereas three- and four-body correlation functions in eta are partially offsetting. The deviation of Deta product from the Stokes-Einstein law arises from the three-body time correlations functions in D.