Search

A method of simulating the drying process of a soft matter solution with an implicit solvent model by moving the liquid-vapor interface is applied to various solution films and droplets. For a solution of a polymer and nanoparticles, we observe "polymer-on-top" stratification, similar to that found previously with an explicit solvent model. Furthermore, "polymer-on-top" is found even when the nanoparticle size is smaller than the radius of gyration of the polymer chains. For a suspension droplet of a bidisperse mixture of nanoparticles, we show that core-shell clusters of nanoparticles can be obtained via the "small-on-outside" stratification mechanism at fast evaporation rates. "Large-on-outside" stratification and uniform particle distribution are also observed when the evaporation rate is reduced. Polymeric particles with various morphologies, including Janus spheres, core-shell particles, and patchy particles, are produced from drying droplets of polymer solutions by combining fast evaporation with a controlled interaction between the polymers and the liquid-vapor interface. Our results validate the applicability of the moving interface method to a wide range of drying systems. The limitations of the method are pointed out and cautions are provided to potential practitioners on cases where the method might fail., Comment: 14 pages, 7 figures

Large scale molecular dynamics simulations are used to study drying suspensions of a binary mixture of large and small particles in explicit and implicit solvents. The solvent is first modeled explicitly and then mapped to a uniform viscous medium by matching the diffusion coefficients and the pair correlation functions of the particles. `Small-on-top' stratification of the particles, with an enrichment of the smaller ones at the receding liquid-vapor interface during drying, is observed in both models under the same drying conditions. With the implicit solvent model, we are able to model much thicker films and study the effect of the initial film thickness on the final distribution of particles in the dry film. Our results show that the degree of stratification is controlled by the P\'{e}clet number defined using the initial film thickness as the characteristic length scale. When the P\'{e}clet numbers of large and small particles are much larger than 1, the degree of `small-on-top' stratification is first enhanced and then weakens as the P\'{e}clet numbers are increased., Comment: 10 pages, 7 figures

A potential strategy for controlling stratification in a drying suspension of bidisperse particles is studied using molecular dynamics simulations. When the suspension is maintained at a constant temperature during fast drying, it can exhibit "small-on-top" stratification with an accumulation (depletion) of smaller (larger) particles in the top region of the drying film, consistent with the prediction of current theories based on diffusiophoresis. However, when only the region near the substrate is thermalized at a constant temperature, a negative temperature gradient develops in the suspension because of evaporative cooling at the liquid-vapor interface. Since the associated thermophoresis is stronger for larger nanoparticles, a higher fraction of larger nanoparticles migrate to the top of the drying film at fast evaporation rates. As a result, stratification is converted to "large-on-top". Very strong "small-on-top" stratification can be produced with a positive thermal gradient in the drying suspension. Here we explore a way to produce a positive thermal gradient by thermalizing the vapor at a temperature higher than that of the solvent. Possible experimental approaches to realize various thermal gradients in a suspension undergoing solvent evaporation, and thus to produce different stratification states in the drying film, are suggested., Comment: 10 pages, 4 figures, supporting information (17 pages, 15 figures)

We study the meniscus on the outside of a small spherical particle with radius $R$ at a liquid-vapor interface. The liquid is confined in a cylindrical container with a finite radius $L$ and has a contact angle $\pi/2$ at the container surface. The center of the particle is placed at various heights along the central axis of the container. By varying $L$, we are able to systematically study the crossover of the meniscus from nanometer to macroscopic scales. The meniscus rise or depression on the particle is found to grow as $\ln (2L/R)$ when $R\ll L\ll \kappa^{-1}$ with $\kappa^{-1}$ being the capillary length and saturate to a value predicted by the Derjaguin-James formula when $R \ll \kappa^{-1} \ll L$. The capillary force on the particle exhibits a linear dependence on the particle's displacement from its equilibrium position at the interface when the displacement is small. The associated spring constant is found to be $2\pi\gamma\ln^{-1} (2L/R)$ for $L\ll \kappa^{-1}$ and saturates to $2\pi\gamma\ln^{-1} (3.7\kappa^{-1}/R)$ for $L\gg \kappa^{-1}$. At nanometer scales, we perform molecular dynamics simulations of the described geometry and the results agree well with the predictions of the macroscopic theory of capillarity. At micrometer to macroscopic scales, comparison to experiments by Anachkov \textit{et al.} [Soft Matter {\bf 12}, 7632 (2016)] shows that the finite span of a liquid-vapor or liquid-liquid interface needs to be considered to interpret experimental data collected with $L \sim \kappa^{-1}$., Comment: 16 pages, 11 figures

Large scale molecular dynamics simulations for bidisperse nanoparticle suspensions with an explicit solvent are used to investigate the effects of evaporation rates and volume fractions on the nanoparticle distribution during drying. Our results show that ``small-on-top'' stratification can occur when ${\rm Pe}_s\phi_s \gtrsim c$ with $c\sim 1$, where ${\rm Pe}_s$ is the P\'{e}clet number and $\phi_s$ is the volume fraction of the smaller particles. This threshold of ${\rm Pe}_s\phi_s$ for ``small-on-top'' is larger by a factor of $\sim\alpha^2$ than the prediction of the model treating solvent as an implicit viscous background, where $\alpha$ is the size ratio between the large and small particles. Our simulations further show that when the evaporation rate of the solvent is reduced, the ``small-on-top'' stratification can be enhanced, which is not predicted by existing theories. This unexpected behavior is explained with thermophoresis associated with a positive gradient of solvent density caused by evaporative cooling at the liquid-vapor interface. For ultrafast evaporation the gradient is large and drives the nanoparticles toward the liquid/vapor interface. This phoretic effect is stronger for larger nanoparticles and consequently the ``small-on-top'' stratification becomes more distinct when the evaporation rate is slower (but not too slow such that a uniform distribution of nanoparticles in the drying film is produced), as thermophoresis that favors larger particles on the top is mitigated. A similar effect can lead to ``large-on-top'' stratification for ${\rm Pe}_s\phi_s$ above the threshold when ${\rm Pe}_s$ is large but $\phi_s$ is small. Our results reveal the importance of including the solvent explicitly when modeling evaporation-induced particle separation and organization and point to the important role of density gradients brought about by ultrafast evaporation., Comment: 11 pages, 4 figures, supporting information (19 pages, 19 figures)

We systematically study the meniscus on the outside of a small circular cylinder vertically immersed in a liquid bath in a cylindrical container that is coaxial with the cylinder. The cylinder has a radius $R$ much smaller than the capillary length, $\kappa^{-1}$, and the container radius, $L$, is varied from a small value comparable to $R$ to $\infty$. In the limit of $L \ll \kappa^{-1}$, we analytically solve the general Young-Laplace equation governing the meniscus profile and show that the meniscus height, $\Delta h$, scales approximately with $R\ln (L/R)$. In the opposite limit where $L \gg \kappa^{-1}$, $\Delta h$ becomes independent of $L$ and scales with $R\ln (\kappa^{-1}/R)$. We implement a numerical scheme to solve the general Young-Laplace equation for an arbitrary $L$ and demonstrate the crossover of the meniscus profile between these two limits. The crossover region has been determined to be roughly $0.4\kappa^{-1} \lesssim L \lesssim 4\kappa^{-1}$. An approximate analytical expression has been found for $\Delta h$, enabling its accurate prediction at any values of $L$ that ranges from microscopic to macroscopic scales., Comment: 9 pages, 6 figures

We study the interplay between the electron-phonon (e-ph) and on-site electron-electron (e-e) interactions in a three-orbital Hubbard-Holstein model on an extended one-dimensional lattice using determinant quantum Monte Carlo. For weak e-e and e-ph interactions, we observe a competition between an orbital-selective Mott phase (OSMP) and a (multicomponent) charge-density-wave (CDW) insulating phase, with an intermediate metallic phase located between them. For large e-e and e-ph couplings, the OSMP and CDW phases persist, while the metallic phase develops short-range orbital correlations and becomes insulating when both the e-e and e-ph interactions are large but comparable. Many of our conclusions are in line with those drawn from a prior dynamical mean field theory study of the two-orbital Hubbard-Holstein model [Phys. Rev. B 95, 12112(R) (2017)] in infinite dimension, suggesting that the competition between the e-ph and e-e interactions in multiorbital Hubbard-Holstein models leads to rich physics, regardless of the dimension of the system., Comment: 12 pages and 11 figures