Logarithmic finite-size effects on interfacial free energies: phenomenological theory and Monte Carlo studies

Interfacial Forces in Free-Standing Layers of Melted Polyethylene, from Critical to Nanoscopic Thicknesses

Molecular dynamics simulations of ultrathin free-standing layers made of melted (373.15–673.15 K) polyethylene chains, which exhibit a lower melting temperature (compared to the bulk value), were carried out to investigate the dominant pressure forces that shape the conformation of chains at the interfacial and bulk liquid regions. We investigated layer thicknesses, tL, from the critical limit of mechanical stability up to lengths of tens of nm and found a normal distribution of bonds dominated by slightly stretched chains across the entire layer, even at large temperatures. In the bulk region, the contribution of bond vibrations to pressure was one order of magnitude larger than the contributions from interchain interactions, which changed from cohesive to noncohesive at larger temperatures just at a transition temperature that was found to be close to the experimentally derived onset temperature for thermal stability. The interchain interactions produced noncohesive interfacial regions at all temperatures in both directions (normal and lateral to the surface layer). Predictions for the value of the surface tension, γ, were consistent with experimental results and were independent of tL. However, the real interfacial thickness—measured from the outermost part of the interface up to the point where γ reached its maximum value—was found to be dependent on tL, located at a distance of 62 Å from the Gibbs dividing surface in the largest layer studied (1568 chains or 313,600 bins); this was ~4 times the length of the interfacial thickness measured in the density profiles.

Interfacial stiffness of nematic-smectic B interface in Gay-Berne liquid crystals using capillary wave theory

The interfacial stiffness for nematic-smectic B (nm-smB) interface in a liquid crystalline (LC) material is calculated using Capillary Wave Theory (CWT) and molecular dynamics simulations. The Gay-Berne (GB) pair potential with parameters κ, κ', μ, and ν equal to 3, 5, 2, and 1 is used to model the LC material. Using a smart three-step recipe, we have obtained an nm-smB phase coexistence in our simulations where the nm and smB directors are nearly parallel to each other and perpendicular to the interface normal. The density profiles are used to compute the nm-smB coexisting density range, the interfacial width, and its position. The smectic phase is differentiated from the nematic phase by using the local bond order parameter (q₆q₆), which has helped us to demonstrate that the interface is indeed rough. Finally, the interfacial stiffness of the nm-smB interface is computed by following the CWT analysis and is found to be γ̃_nm-smB=0.39861k_BT/σ_ee ²=0.04429/σ_ss ², where σ_ee and σ_ss are the length and diameter of the GB LC particles.

Computation of the solid-liquid interfacial free energy in hard spheres by means of thermodynamic integration

We used a thermodynamic integration scheme, which is specifically designed for disordered systems, to compute the interfacial free energy of the solid-liquid interface in the hard-sphere model. We separated the bulk contribution to the total free energy from the interface contribution, performed a finite-size scaling analysis, and obtained for the (100)-interface γ=0.591(11)k_{B}Tσ^{-2}.

Phase Behavior of Poly(ethylene oxide) in Room Temperature Ionic Liquids: A Molecular Simulation and Deep Neural Network Study

The phase behavior of polymers in room temperature ionic liquids is a topic of considerable interest. In this work we study the phase diagram of poly(ethylene oxide) in four imidazolium ionic liquids (ILs) using molecular simulation. We develop united atom models for 1-butyl-2,3-dimethylimidazolium ([BMMIM]), 1-ethyl-2,3-dimethylimidazolium ([EMMIM]), and 1-ethyl-3-methylimidazolium ([EMIM]) in an analogous fashion to previously developed models for 1-butyl-3-methylimidazolium ([BMIM]) and tetrafluoroborate ([BF₄]) using symmetry-adapted perturbation theory. At high temperatures we obtain the coexistence concentrations using an interface method where the polymer and IL are simulated in a large elongated box, and an interface between coexisting phases is formed. At lower temperatures we use a deep neural network (DNN) method. The input descriptors for the DNN are the cohesive energy of mixing, the volume change of mixing, and the coordination numbers between cation and polymer, all of which are obtained from simulations of mixed systems at a series of temperatures. The DNN is trained by using the phase-separated systems at high temperatures and a mixed phase at low temperatures. The method predicts a lower critical solution temperature which decreases as the alkyl chain length on the cation is decreased, consistent with experiment. The simulations show that methylation of the cation has little effect on the phase diagram. This is in contrast to what is seen in experiments but could be because the polymer chains in the simulations are too short. At low temperatures the chains display two conformational motifs, namely a crown ether conformation and a ring conformation, each of which can wrap the chain around a single cation. This provides the entropic penalty for mixing and a reason for demixing as the temperature is raised. Such conformations might not be possible for longer chains. The combination of data-driven techniques and molecular simulation shows promise in the study of the phase behavior and physical properties of complex fluids.

Molecular Structure and Modeling of Water-Air and Ice-Air Interfaces Monitored by Sum-Frequency Generation

From a glass of water to glaciers in Antarctica, water-air and ice-air interfaces are abundant on Earth. Molecular-level structure and dynamics at these interfaces are key for understanding many chemical/physical/atmospheric processes including the slipperiness of ice surfaces, the surface tension of water, and evaporation/sublimation of water. Sum-frequency generation (SFG) spectroscopy is a powerful tool to probe the molecular-level structure of these interfaces because SFG can specifically probe the topmost interfacial water molecules separately from the bulk and is sensitive to molecular conformation. Nevertheless, experimental SFG has several limitations. For example, SFG cannot provide information on the depth of the interface and how the orientation of the molecules varies with distance from the surface. By combining the SFG spectroscopy with simulation techniques, one can directly compare the experimental data with the simulated SFG spectra, allowing us to unveil the molecular-level structure of water-air and ice-air interfaces. Here, we present an overview of the different simulation protocols available for SFG spectra calculations. We systematically compare the SFG spectra computed with different approaches, revealing the advantages and disadvantages of the different methods. Furthermore, we account for the findings through combined SFG experiments and simulations and provide future challenges for SFG experiments and simulations at different aqueous interfaces.

Reducing uncertainty in simulation estimates of the surface tension through a two-scale finite-size analysis: thicker is better

Recent simulation studies of the surface tension γ, and other properties of thin free-standing films, have revealed unexpected finite size effects in which the variance of the properties vary monotonically with the in-plane width of the films, complicating the extrapolation of estimates of film properties to the thermodynamic limit. We carried out molecular dynamics simulations to determine the origin of this phenomenon, and to address the practical problem of developing a more reliable methodology for estimating γ in the thermodynamic limit. We find that there are two distinct finite size effects that must be addressed in a finite size analysis of γ in thin films. The first finite size scale is the in-plane width of the films and the second scale is the simulation cell size in the transverse direction. Increasing the first scale enhances fluctuations in γ, measured by the standard deviation of their distribution, while increasing the second reduces γ fluctuations due to a corresponding increased 'freedom' of the film to fluctuate out of plane. We find that using progressively large simulation cells in the transverse direction, while keeping the film width fixed to an extent in which the full bulk liquid zone is developed, allows us to obtain a smooth extrapolation to the thermodynamic limit, enabling a reduction of the γ uncertainty to a magnitude on the order of 1% for systems having a reasonably large size, i.e., O (1 μm).

A simulation method for the phase diagram of complex fluid mixtures

The phase behavior of complex fluid mixtures is of continuing interest, but obtaining the phase diagram from computer simulations can be challenging. In the Gibbs ensemble method, for example, each of the coexisting phases is simulated in a different cell, and ensuring the equality of chemical potentials of all components requires the transfer of molecules from one cell to the other. For complex fluids such as polymers, successful insertions are rare. An alternative method is to simulate both coexisting phases in a single simulation cell, with an interface between them. The challenge here is that the interface position moves during the simulation, making it difficult to determine the concentration profile and coexisting concentrations. In this work, we propose a new method for single cell simulations that uses a spatial concentration autocorrelation function to (spatially) align instantaneous concentration profiles from different snapshots. This allows one to obtain average concentration profiles and hence the coexisting concentrations. We test the method by calculating the phase diagrams of two systems: the Widom-Rowlinson model and the symmetric blends of freely jointed polymer molecules for which phase diagrams from conventional methods are available. Excellent agreement is found, except in the neighborhood of the critical point where the interface is broad and finite size effects are important. The method is easy to implement and readily applied to any mixture of complex fluids.

Unexpected finite size effects in interfacial systems: Why bigger is not always better-Increase in uncertainty of surface tension with bulk phase width

We present an unexpected finite size effect affecting interfacial molecular simulations that is proportional to the width-to-surface-area ratio of the bulk phase L_l/A. This finite size effect has a significant impact on the variance of surface tension values calculated using the virial summation method. A theoretical derivation of the origin of the effect is proposed, giving a new insight into the importance of optimising system dimensions in interfacial simulations. We demonstrate the consequences of this finite size effect via a new way to estimate the surface energetic and entropic properties of simulated air-liquid interfaces. Our method is based on macroscopic thermodynamic theory and involves comparing the internal energies of systems with varying dimensions. We present the testing of these methods using simulations of the TIP4P/2005 water forcefield and a Lennard-Jones fluid model of argon. Finally, we provide suggestions of additional situations, in which this finite size effect is expected to be significant, as well as possible ways to avoid its impact.

Heterogeneous nucleation of a droplet pinned at a chemically inhomogeneous substrate: A simulation study of the two-dimensional Ising case

Heterogeneous nucleation is studied by Monte Carlo simulations and phenomenological theory, using the two-dimensional lattice gas model with suitable boundary fields. A chemical inhomogeneity of length b at one boundary favors the liquid phase, while elsewhere the vapor is favored. Switching on the bulk field H_b favoring the liquid, nucleation and growth of the liquid phase starting from the region of the chemical inhomogeneity are analyzed. Three regimes occur: for small fields, H_b<H_b^crit, the critical droplet radius is so large that a critical droplet having the contact angle θ_c required by Young's equation in the region of the chemical inhomogeneity does not yet "fit" there since the baseline length of the circle-cut sphere droplet would exceed b. For H_b^crit<H_b<H_b^*, such droplets fit inside the inhomogeneity and are indeed found in simulations with large enough observation times, but these droplets remain pinned to the chemical inhomogeneity when their baseline has grown to the length b. Assuming that these pinned droplets have a circle cut shape and effective contact angles θ_eff in the regime θ_c < θ_eff < π/2, the density excess due to these droplets can be predicted and is found to be in reasonable agreement with the simulation results. On general grounds, one can predict that the effective contact angle θ_eff and the excess density of the droplets, scaled by b, are functions of the product bH_b but do not depend on both variables separately. Since the free energy barrier for the "depinning" of the droplet (i.e., growth of θ_eff to π - θ_c) vanishes when θ_eff approaches π/2, in practice only angles θ_eff up to about θ_eff^max≃70° were observed. For larger fields (H_b>H_b^*), the droplets nucleated at the chemical inhomogeneity grow to the full system size. While the relaxation time for the growth scales as τ_G∝H_b^-1, the nucleation time τ_N scales as lnτ_N∝H_b^-1. However, the prefactor in the latter relation, as evaluated for our simulations results, is not in accord with an extension of the Volmer-Turnbull theory to two-dimensions, when the theoretical contact angle θ_c is used.

Fluctuations near the liquid–liquid transition in a model of silica

Molecular dynamics simulations reveal anomalous small-angle scattering and liquid–liquid phase separation in an ionic model of silica.

Free-energy barriers for crystal nucleation from fluid phases

Monte Carlo simulations of crystal nuclei coexisting with the fluid phase in thermal equilibrium in finite volumes are presented and analyzed, for fluid densities from dense melts to the vapor. Generalizing the lever rule for two-phase coexistence in the canonical ensemble to finite volume, "measurements" of the nucleus volume together with the pressure and chemical potential of the surrounding fluid allows us to extract the surface free energy of the nucleus. Neither the knowledge of the (in general nonspherical) nucleus shape nor of the angle-dependent interface tension is required for this task. The feasibility of the approach is demonstrated for a variant of the Asakura-Oosawa model for colloid-polymer mixtures, which form face-centered cubic colloidal crystals. For a polymer to colloid size ratio of 0.15, the colloid packing fraction in the fluid phase can be varied from melt values to zero by the variation of an effective attractive potential between the colloids. It is found that the approximation of spherical crystal nuclei often underestimates actual nucleation barriers significantly. Nucleation barriers are found to scale as ΔF^{*}=(4π/3)^{1/3}γ[over ¯](V^{*})^{2/3}+const with the nucleus volume V^{*}, and the effective surface tension γ[over ¯] that accounts implicitly for the nonspherical shape can be precisely estimated.

Rapid ordering of block copolymer thin films

Block-copolymers self-assemble into diverse morphologies, where nanoscale order can be finely tuned via block architecture and processing conditions. However, the ultimate usage of these materials in real-world applications may be hampered by the extremely long thermal annealing times-hours or days-required to achieve good order. Here, we provide an overview of the fundamentals of block-copolymer self-assembly kinetics, and review the techniques that have been demonstrated to influence, and enhance, these ordering kinetics. We discuss the inherent tradeoffs between oven annealing, solvent annealing, microwave annealing, zone annealing, and other directed self-assembly methods; including an assessment of spatial and temporal characteristics. We also review both real-space and reciprocal-space analysis techniques for quantifying order in these systems.

Variance-reduced estimator of the connected two-point function in the presence of a broken Z(2)-symmetry

The exchange or geometric cluster algorithm allows us to define a variance-reduced estimator of the connected two-point function in the presence of a broken Z(2)-symmetry. We present numerical tests for the improved Blume-Capel model on the simple-cubic lattice. We perform simulations for the critical isotherm, the low-temperature phase at vanishing external field, and, for comparison, also the high-temperature phase. For the connected two-point function, a substantial reduction of the variance can be obtained, allowing us to compute the correlation length ξ with high precision. Based on these results, estimates for various universal amplitude ratios that characterize the universality class of the three-dimensional Ising model are computed.

Negative Interfacial Tension in Phase-Separated Active Brownian Particles

We study numerically a model for active suspensions of self-propelled repulsive particles, for which a stable phase separation into a dilute and a dense phase is observed. We exploit the fact that for nonsquare boxes a stable "slab" configuration is reached, in which interfaces align with the shorter box edge. Evaluating a recent proposal for an intensive active swimming pressure, we demonstrate that the excess stress within the interface separating both phases is negative. The occurrence of a negative tension together with stable phase separation is a genuine nonequilibrium effect that is rationalized in terms of a positive stiffness, the estimate of which agrees excellently with the numerical data. Our results challenge effective thermodynamic descriptions and mappings of active Brownian particles onto passive pair potentials with attractions.

Latent Alignment in Pathway-Dependent Ordering of Block Copolymer Thin Films

Block copolymers spontaneously form well-defined nanoscale morphologies during thermal annealing. Yet, the structures one obtains can be influenced by nonequilibrium effects, including processing history or pathway-dependent assembly. Here, we explore various pathways for ordering of block copolymer thin films, using oven-annealing, as well as newly disclosed methods for rapid photothermal annealing and photothermal shearing. We report the discovery of an efficient pathway for ordering self-assembled films: ultrarapid shearing of as-cast films induces "latent alignment" in the disordered morphology. Subsequent thermal processing can then develop this directly into a uniaxially aligned morphology with low defect density. This deeper understanding of pathway-dependence may have broad implications in self-assembly.

Crystal-liquid interfacial free energy of hard spheres via a thermodynamic integration scheme

The hard-sphere crystal-liquid interfacial free energy γcl is determined from molecular dynamics simulations using a thermodynamic integration (TI) scheme. The advantage of this TI scheme compared to previous methods is to successfully circumvent hysteresis effects due to the movement of the crystal-liquid interface. This is accomplished by the use of extremely-short-range and impenetrable Gaussian flat walls that prevent the drift of the interface while imposing a negligible free-energy penalty. We find that it is crucial to analyze finite-size effects in order to obtain reliable estimates of γcl in the thermodynamic limit.

Calculation of the absolute free energy of a smectic-A phase

In this paper, we provide a scheme to compute the absolute free energy of a smectic-A phase via the "indirect method." The state of interest is connected through a three-step reversible path to a reference state. This state consists of a low-density layer of rods coupled to two external fields maintaining these rods close to the layer's plane and oriented preferably normal to the layer. The low-density free energy of the reference state can be computed on the basis of the relevant second virial coefficients between two rods coupled to the two external fields. We apply this technique to the Gay-Berne potential for calamitics with a parameter set leading to stable isotropic (I), nematic (N), smectic-A (SmA), and crystal (Cr) phases. We locate the I-SmA phase transition at low pressure and the sequence of phase transitions I-N-SmA along higher-pressure isobars and we establish the location of the I-N-SmA triple point. Close to this triple point, we show that the N-SmA transition is clearly first order. Our results are compared to the coexistence lines of the approximate phase diagram elucidated by de Miguel et al. [J. Chem. Phys. 121, 11183 (2004)] established through the direct observation of the sequence of phase transitions occurring along isobars under heating or cooling sequences of runs. Finally, we discuss the potential of our technique in studying similar transitions observed on layered phases under confinement.