H.H. Wensink,G.J. Vroege
Abstract:The isotropic-nematic phase behavior of length polydisperse hard rods with arbitrary length distributions is calculated. Within a numerical treatment of the polydisperse Onsager model using the Gaussian trial function Ansatz we determine the onset of isotropic-nematic phase separation, coming from a dilute isotropic phase and a dense nematic phase. We focus on parent systems whose lengths can be described by either a Schulz or a `fat-tailed' log-normal distribution with appropriate lower and upper cutoff lengths. In both cases, very strong fractionation effects are observed for parent polydispersities larger than roughly 50 %. In these regimes, the isotropic and nematic phases are completely dominated by respectively the shortest and the longest rods in the system. Moreover, for the log-normal case, we predict triphasic isotropic-nematic-nematic equilibria to occur above a certain threshold polydispersity. By investigating the properties of the coexisting phases across the coexistence region for a particular set of cutoff lengths we explicitly show that the region of stable triphasic equilibria does not extend up to very large parent polydispersities but closes off at a consolute point located not far above the threshold polydispersity. The experimental relevance of the phenomenon is discussed.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to study the behavior of the isotropic - nematic phase transition in a system of length - polydisperse hard rods. Specifically, the author calculates the starting point of phase separation in a hard - rod system with an arbitrary length distribution during the transition from a sparse isotropic phase to a dense nematic phase by numerically processing the polydisperse Onsager model.
#### Main problems include:
1. **The influence of polydispersity on phase behavior**: The author focuses on systems whose length distributions can be described by the Schulz distribution or the "fat - tailed" log - normal distribution. They specifically explore the strong fractionation effects, i.e., the phenomenon that the shortest and the longest rods dominate the isotropic and nematic phases respectively, when the polydispersity of the parent system is greater than about 50%.
2. **Three - phase equilibrium phenomenon**: For the case of the log - normal distribution, the author predicts that when the polydispersity exceeds a certain threshold, an isotropic - nematic - nematic triphasic equilibria may occur. By studying the properties of the coexisting phases in the coexistence region, the author clearly shows that this stable three - phase equilibrium region does not extend to a very large parent polydispersity, but closes near the threshold polydispersity.
3. **Experimental relevance**: The author discusses the experimental relevance of these phenomena, especially how to interpret the influence of polydispersity in experimental results.
### Research methods
To study these problems, the author uses the Gaussian trial function Ansatz to calculate the isotropic - nematic phase behavior of the hard - rod system. The advantages of this method are:
- All integrals required for the analytic treatment of the monodisperse Onsager model can be processed, so only numerical integration in the polydisperse case needs to be considered.
- The Gaussian ODF (orientation distribution function) can better describe the highly ordered nematic phase and is more accurate compared to the P2 approximation, especially in the presence of extremely long rods.
### Results and conclusions
Through a detailed analysis of the numerical results of different parent distributions (Schulz distribution and log - normal distribution), the author finds that:
- When the polydispersity of the parent system is large, the nematic phase is mainly composed of the longest rods, while the isotropic phase is mainly composed of the shortest rods.
- For the log - normal distribution, there indeed exists a specific polydispersity threshold, and when this threshold is exceeded, the three - phase equilibrium phenomenon occurs, but the stable region of this three - phase equilibrium is limited and only exists in the range close to the threshold.
In conclusion, through numerical simulation and theoretical analysis, this paper reveals the complex phase behavior in polydisperse hard - rod systems and provides a theoretical basis for understanding experimental results.