Showing 1–50 of 52 results for author: de Pablo, J J

Thermodynamic and Transport Properties of Binary Mixtures of Polyethylene and Higher n-Alkanes from Physics-Informed and Machine-Learned Models

Authors: Maria Ley-Flores, Riccardo Alessandri, Sean Najmi, Michele Valsecchi, George Jackson, Amparo Galindo, LaShanda Korley, Dionisios G. Vlachos, Juan J. de Pablo

Abstract: The thermodynamics and transport properties of polymeric materials are essential for the design of reactors and for the development of polymer deconstruction processes. Existing property prediction tools such as correlations based on entropy scaling, kinetic gas theory, and free-volume model are inadequate for polymers. In this paper, we introduce a data-driven model for polyolefins based on data… ▽ More The thermodynamics and transport properties of polymeric materials are essential for the design of reactors and for the development of polymer deconstruction processes. Existing property prediction tools such as correlations based on entropy scaling, kinetic gas theory, and free-volume model are inadequate for polymers. In this paper, we introduce a data-driven model for polyolefins based on data from molecular dynamics simulations that can accurately predict the transport properties of polyethylenes and their binary mixtures with higher n-alkanes across a range of temperatures, pressures, concentrations, and oligomer molecular weights. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Numerical Study of Cleavable Bond-Modified Polyethylene for Circular Polymer Design

Authors: Maria Ley-Flores, Archit Chabbi, Riccardo Alessandri, Sam Marsden, Isabella Vettese, Stuart J. Rowan, Juan J. de Pablo

Abstract: There is considerable interest in designing new polymeric materials with built-in mechanisms for recycling. In this work, we present a systematic exploration of several polyethylene-based polymers that contain strategically incorporated cleavable bonds that are susceptible to chemical degradation. We consider ten distinct telechelic functionalities across varying chain lengths: ester, aromatic est… ▽ More There is considerable interest in designing new polymeric materials with built-in mechanisms for recycling. In this work, we present a systematic exploration of several polyethylene-based polymers that contain strategically incorporated cleavable bonds that are susceptible to chemical degradation. We consider ten distinct telechelic functionalities across varying chain lengths: ester, aromatic ester, anhydride, carbonate, urethane, siloxane, oxalate, acetal, urea, and amide linkages. Our study is focused on the effects of cleavable bonds on density, diffusion, enthalpy, and melting temperatures. We also examine the crystallization process. These are key properties for polymer processing and, by comparing them to those of linear polyethylene, we establish a platform that will help guide future synthetic efforts towards the design of sustainable plastics with performance comparable to that of traditional, high-performance polyolefin-based materials. △ Less

Submitted 16 April, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Synthetic active liquid crystals powered by acoustic waves

Authors: Andrey Sokolov, Jaideep Katuri, Juan J. de Pablo, Alexey Snezhko

Abstract: Active nematics are materials composed of mobile, elongated particles that can transform energy from the environment into a mechanical motion. Current experimental realizations of the active nematics are of biological origin and include cell layers, suspensions of elongated bacteria in liquid crystal, and combinations of bio-filaments with molecular motors. Here, we report the realization of a ful… ▽ More Active nematics are materials composed of mobile, elongated particles that can transform energy from the environment into a mechanical motion. Current experimental realizations of the active nematics are of biological origin and include cell layers, suspensions of elongated bacteria in liquid crystal, and combinations of bio-filaments with molecular motors. Here, we report the realization of a fully synthetic active nematic system comprised of a lyotropic chromonic liquid crystal energized by ultrasonic waves. This artificial active liquid crystal is free from biological degradation and variability, exhibits stable material properties, and enables precise and rapid activity control over a significantly extended range. We demonstrate that the energy of the acoustic field is converted into microscopic extensile stresses disrupting long-range nematic order and giving rise to an undulation instability, development of active turbulence, and proliferation of topological defects. We reveal the emergence of unconventional free-standing persistent vortices in the nematic director field at high activity levels. The results provide a foundation for the design of externally energized active nematic fluids with stable material properties and tunable topological defects dynamics crucial for the realization of reconfigurable microfluidic systems. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Light-responsive active particles in a thermotropic liquid crystal

Authors: Antonio Tavera-Vázquez, Andrés Córdoba, Sam Rubin, Vincenzo Vitelli, Juan J. de Pablo

Abstract: The development of synthetic microswimmers has advanced our understanding of the fundamental self-propelled mechanisms of living systems. However, there are scarce studies at the microscale within highly structured anisotropic media, such as bacteria or cellular receptors that swim in concentrated solutions of filamentous proteins or lipids with viscoelastic properties. Synthetic liquid crystals (… ▽ More The development of synthetic microswimmers has advanced our understanding of the fundamental self-propelled mechanisms of living systems. However, there are scarce studies at the microscale within highly structured anisotropic media, such as bacteria or cellular receptors that swim in concentrated solutions of filamentous proteins or lipids with viscoelastic properties. Synthetic liquid crystals (LCs) have the potential to serve as biomimetic surrogates to study the structure and dynamics of living systems. Nevertheless, studies on thermotropic LCs have mainly focused on electro and magneto-phoretic effects, with a few others on diffusiophoresis or light-driven distortions of the LC nematic director. To the best of our knowledge, here we report self-thermophoretic experiments on thermotropic LCs for the first time. Our system consists of 2D confined Janus particles in 5CB with homeotropic anchoring on the particle and LC cell surfaces. The Janus particles include a conductive titanium coating that, upon exposure to an LED source, is heated and induces a local steady nematic-isotropic phase transition, leading to the self-propulsion of the particles orthogonally to the LC director. The trajectories of the Janus particles were tracked at different intensities of the applied light. A model is developed to describe the mean-squared displacement of a Janus particle suspended in a nematic LC. The model assumes that the Janus particle feels the LC as a continuum with the anisotropic viscosity of the bulk nematic phase. Moreover, the viscoelasticity of the LC is also considered. The model describes the experimental data well, and the fitting parameter related to the magnitude of the swimming force increases with the intensity of the applied light. Moreover, our approach represents an important step for developing a platform for highly structured anisotropic active materials. △ Less

Submitted 6 May, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 15 pages, 6 figures

Rheology of bidisperse non-Brownian suspensions

Authors: Abhinendra Singh, Christopher Ness, Abhishek K. Sharma, Juan J. de Pablo, Heinrich M. Jaeger

Abstract: We study the rheology of bidisperse non-Brownian suspensions using particle-based simulation, mapping the viscosity as a function of the size ratio of the species, their relative abundance, and the overall solid content. The variation of the viscosity with applied stress exhibits shear thickening phenomenology irrespective of composition, though the stress-dependent limiting solids fraction govern… ▽ More We study the rheology of bidisperse non-Brownian suspensions using particle-based simulation, mapping the viscosity as a function of the size ratio of the species, their relative abundance, and the overall solid content. The variation of the viscosity with applied stress exhibits shear thickening phenomenology irrespective of composition, though the stress-dependent limiting solids fraction governing the viscosity and its divergence point are non-monotonic in the mixing ratio. Contact force data demonstrate an asymmetric exchange in dominant stress contribution from large-large to small-small particle contacts as the mixing ratio of the species evolves. Combining a prior model for shear thickening with one for composition-dependent jamming, we obtain a full description of the rheology of bidisperse non-Brownian suspensions capable of predicting effects such as the viscosity reduction observed upon adding small particle fines to a suspension of large particles. △ Less

Submitted 31 May, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 5 pages, 4 figures

LCPOM: Precise Reconstruction of Polarized Optical Microscopy Images of Liquid Crystals

Authors: Chuqiao Chen, Viviana Palacio-Betancur, Sepideh Norouzi, Pablo F. Zubieta Rico, Monirosadat Sadati, Stuart J. Rowan, Juan J. de Pablo

Abstract: When viewed with a cross-polarized optical microscope (POM), liquid crystals display interference colors and complex patterns that depend on the material's microscopic orientation. That orientation can be manipulated by application of external fields, which provides the basis for applications in optical display and sensing technologies. The color patterns themselves have a high information content… ▽ More When viewed with a cross-polarized optical microscope (POM), liquid crystals display interference colors and complex patterns that depend on the material's microscopic orientation. That orientation can be manipulated by application of external fields, which provides the basis for applications in optical display and sensing technologies. The color patterns themselves have a high information content. Traditionally, however, calculations of the optical appearance of liquid crystals have been performed by assuming that a single-wavelength light source is employed, and reported in a monochromatic scale. In this work, the original Jones matrix method is extended to calculate the colored images that arise when a liquid crystal is exposed to a multi-wavelength source. By accounting for the material properties, the visible light spectrum and the CIE color matching functions, we demonstrate that the proposed approach produces colored POM images that are in quantitative agreement with experimental data. Results are presented for a variety of systems, including radial, bipolar, and cholesteric droplets, where results of simulations are compared to experimental microscopy images. The effects of droplet size, topological defect structure, and droplet orientation are examined systematically. The technique introduced here generates images that can be directly compared to experiments, thereby facilitating machine learning efforts aimed at interpreting LC microscopy images, and paving the way for the inverse design of materials capable of producing specific internal microstructures in response to external stimuli. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 12 pages, 5 figures (main text). 6 pages, 6 figures (appendices)

Structure and Dynamics of Hybrid Colloid-Polyelectrolyte Coacervates: Insights from Molecular Simulations

Authors: Boyuan Yu, Heyi Liang, Paul F. Nealey, Matthew Tirrell, Artem M. Rumyantsev, Juan J. de Pablo

Abstract: Electrostatic interactions in polymeric systems are responsible for a wide range of liquid-liquid phase transitions that are of importance for biology and materials science. Such transitions are referred to as complex coacervation, and recent studies have sought to understand the underlying physics and chemistry. Most theoretical and simulation efforts to date have focused on oppositely charged li… ▽ More Electrostatic interactions in polymeric systems are responsible for a wide range of liquid-liquid phase transitions that are of importance for biology and materials science. Such transitions are referred to as complex coacervation, and recent studies have sought to understand the underlying physics and chemistry. Most theoretical and simulation efforts to date have focused on oppositely charged linear polyelectrolytes, which adopt nearly ideal-coil conformations in the condensed phase. However, when one of the coacervate components is a globular protein, a better model of complexation should replace one of the species with a spherical charged particle or colloid. In this work, we perform coarse-grained simulations of colloid-polyelectrolyte coacervation using a spherical model for the colloid. Simulation results indicate that the electroneutral cell of the resulting (hybrid) coacervates consists of a polyelectrolyte layer adsorbed on the colloid. Power laws for the structure and the density of the condensed phase, which are extracted from simulations, are found to be consistent with the adsorption-based scaling theory of coacervation. The coacervates remain amorphous (disordered) at a moderate colloid charge, $Q$, while an intra-coacervate colloidal crystal is formed above a certain threshold, at $Q > Q^{*}$. In the disordered coacervate, if $Q$ is sufficiently low, colloids diffuse as neutral non-sticky nanoparticles in the semidilute polymer solution. For higher $Q$, adsorption is strong and colloids become effectively sticky. Our findings are relevant for the coacervation of polyelectrolytes with proteins, spherical micelles of ionic surfactants, and solid organic or inorganic nanoparticles. △ Less

Submitted 31 May, 2023; originally announced June 2023.

Journal ref: Macromolecules 56 (2023) 7256-7270

Scattering Evidence of Positional Charge Correlations in Polyelectrolyte Complexes

Authors: Yan Fang, Artem M. Rumyantsev, Angelika E. Neitzel, Heyi Liang, William T. Heller, Paul F. Nealey, Matthew V. Tirrel, Juan J. de Pablo

Abstract: Polyelectrolyte complexation plays an important role in materials science and biology. The internal structure of the resultant polyelectrolyte complex (PEC) phase dictates properties such as physical state, response to external stimuli, and dynamics. Small-angle scattering experiments with X-rays and neutrons have revealed structural similarities between PECs and semidilute solutions of neutral po… ▽ More Polyelectrolyte complexation plays an important role in materials science and biology. The internal structure of the resultant polyelectrolyte complex (PEC) phase dictates properties such as physical state, response to external stimuli, and dynamics. Small-angle scattering experiments with X-rays and neutrons have revealed structural similarities between PECs and semidilute solutions of neutral polymers, where the total scattering function exhibits an Ornstein-Zernike form. In spite of consensus among different theoretical predictions, the existence of positional correlations between polyanion and polycation charges has not been confirmed experimentally. Here, we present small-angle neutron scattering profiles where the polycation scattering length density is matched to that of the solvent to extract positional correlations among anionic monomers. The polyanion scattering functions exhibit a peak at the inverse polymer screening radius of Coulomb interactions, $q^{*} \approx 0.2 Å^{-1}$. This peak, attributed to Coulomb repulsions between the fragments of polyanions and their attractions to polycations, is even more pronounced in the calculated charge scattering function that quantifies positional correlations of all polymer charges within the PEC. Screening of electrostatic interactions by adding salt leads to the gradual disappearance of this correlation peak, and the scattering functions regain an Ornstein-Zernike form. Experimental scattering results are consistent with those calculated from the random phase approximation, a scaling analysis, and molecular simulations. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Journal ref: Proc. Natl. Acad. Sci. U.S.A. 120(2023) e2302151120

Acousto-microfluidic Control of Liquid Crystals

Authors: Gustavo A. Vasquez-Montoya, Tadej Emersic, Noe Atzin, Antonio Tavera-Vazquez, Ali Mozaffari, Rui Zhang, Orlando Guzman, Alexey Snezhko, Paul F. Nealey, Juan J. de Pablo

Abstract: The optical properties of liquid crystals serve as the basis for display, diagnostic, and sensing technologies. Such properties are generally controlled by relying on electric fields. In this work, we investigate the effects of microfluidic flows and acoustic fields on the molecular orientation and the corresponding optical response of nematic liquid crystals. Several previously unknown structures… ▽ More The optical properties of liquid crystals serve as the basis for display, diagnostic, and sensing technologies. Such properties are generally controlled by relying on electric fields. In this work, we investigate the effects of microfluidic flows and acoustic fields on the molecular orientation and the corresponding optical response of nematic liquid crystals. Several previously unknown structures are identified, which are rationalized in terms of a state diagram as a function of the strengths of the flow and the acoustic field. The new structures are interpreted by relying on calculations with a free energy functional expressed in terms of the tensorial order parameter, using continuum theory simulations in the Landau-de Gennes framework. Taken together, the findings presented here offer promise for the development of new systems based on combinations of sound, flow, and confinement. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Measuring response functions of active materials from data

Authors: Mehdi Molaei, Steven A. Redford, Wen-hung Chou, Danielle Scheff, Juan J. de Pablo, Patrick W. Oakes, Margaret L. Gardel

Abstract: From flocks of birds to biomolecular assemblies, systems in which many individual components independently consume energy to perform mechanical work exhibit a wide array of striking behaviors. Methods to quantify the dynamics of these so called active systems generally aim to extract important length or time scales from experimental fields. Because such methods focus on extracting scalar values, t… ▽ More From flocks of birds to biomolecular assemblies, systems in which many individual components independently consume energy to perform mechanical work exhibit a wide array of striking behaviors. Methods to quantify the dynamics of these so called active systems generally aim to extract important length or time scales from experimental fields. Because such methods focus on extracting scalar values, they do not wring maximal information from experimental data. We introduce a method to overcome these limitations. We extend the framework of correlation functions by taking into account the internal headings of displacement fields. The functions we construct represent the material response to specific types of active perturbation within the system. Utilizing these response functions we query the material response of disparate active systems composed of actin filaments and myosin motors, from model fluids to living cells. We show we can extract critical length scales from the turbulent flows of an active nematic, anticipate contractility in an active gel, distinguish viscous from viscoelastic dissipation, and even differentiate modes of contractility in living cells. These examples underscore the vast utility of this method which measures response functions from experimental observations of complex active systems. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Dynamics of active defects on the anisotropic surface of an ellipsoidal droplet

Authors: Martina Clairand, Ali Mozaffari, Jérôme Hardoüin, Rui Zhang, Claire Doré, Jordi Ignés-Mullol, Francesc Sagués, Juan J. de Pablo, Teresa Lopez-Leon

Abstract: Cells are fundamental building blocks of living organisms displaying an array of shapes, morphologies, and textures that encode specific functions and physical behaviors. Elucidating the rules of this code remains a challenge. In this work, we create biomimetic structural building blocks by coating ellipsoidal droplets of a smectic liquid crystal with a protein-based active cytoskeletal gel, thus… ▽ More Cells are fundamental building blocks of living organisms displaying an array of shapes, morphologies, and textures that encode specific functions and physical behaviors. Elucidating the rules of this code remains a challenge. In this work, we create biomimetic structural building blocks by coating ellipsoidal droplets of a smectic liquid crystal with a protein-based active cytoskeletal gel, thus obtaining core-shell structures. By exploiting the patterned texture and anisotropic shape of the smectic core, we were able to mold the complex nematodynamics of the interfacial active material and identify new time-dependent states where topological defects periodically oscillate between rotational and translational regimes. Our nemato-hydrodynamic simulations of active nematics demonstrate that, beyond topology and activity, the dynamics of the active material are profoundly influenced by the local curvature and smectic texture of the droplet, as well as by external hydrodynamic forces. These results illustrate how the incorporation of these constraints into active nematic shells orchestrates remarkable spatio-temporal motifs, offering critical new insights into biological processes and providing compelling prospects for designing bio-inspired micro-machines. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 16 pages, 8 figures

Generation and propagation of solitary waves in nematic liquid crystals

Authors: Xingzhou Tang, Ali Mozaffari, Noe Atzin, Soumik Das, Nicholas L. Abbott, Juan J. de Pablo

Abstract: Solitons in nematic liquid crystals offer intriguing opportunities for transport and sensing in microfluidic systems. Little is known about the elementary conditions that are needed to create solitons in nematic materials. In this work, theory, simulations and experiments are used to study the generation and propagation of solitary waves (or "solitons") in nematic liquid crystals upon the applicat… ▽ More Solitons in nematic liquid crystals offer intriguing opportunities for transport and sensing in microfluidic systems. Little is known about the elementary conditions that are needed to create solitons in nematic materials. In this work, theory, simulations and experiments are used to study the generation and propagation of solitary waves (or "solitons") in nematic liquid crystals upon the application of an alternating current (AC) electric field. We find that these solitary waves exhibit "butterfly"-like or "bullet"-like structures that travel in the direction perpendicular to the applied electric field. Such structures propagate over long distances without losing their initial shape. The theoretical model adopted here serves to identify some of the key requirements that are needed to generate solitons in the absence of electrostatic interactions. These include surface imperfections that introduce a twist in the director, unequal elastic constants, and negative anisotropic dielectric permittivity. The results of simulations are shown to be in good agreement with our own experimental observations, serving to establish the validity of the theoretical concepts advanced in this work. △ Less

Submitted 2 November, 2022; originally announced November 2022.

A minimal model of solitons in nematic liquid crystals

Authors: Noe Atzin, Ali Mozaffari, Xingzhou Tang, Soumik Das, Nicholas L. Abbott, Juan J. de Pablo

Abstract: Solitons in liquid crystals have generated considerable interest. Several hypotheses of varying complexity have been advanced to explain how they emerge, and a consensus has not emerged yet about the underlying forces responsible for their formation or their structure. In this work, we present a minimal model for soliton structures in achiral nematic liquid crystals, which reveals the key requirem… ▽ More Solitons in liquid crystals have generated considerable interest. Several hypotheses of varying complexity have been advanced to explain how they emerge, and a consensus has not emerged yet about the underlying forces responsible for their formation or their structure. In this work, we present a minimal model for soliton structures in achiral nematic liquid crystals, which reveals the key requirements needed to generate traveling solitons in the absence of added charges. These include a surface imperfection or inhomogeneity capable of producing a twist, flexoelectricity, dielectric contrast, and an applied AC electric field that can couple to the director's orientation. Our proposed model is based on a tensorial representation of a confined liquid crystal, and it predicts the formation of "butterfly" structures, quadrupolar in character, in regions of a slit channel where the director is twisted by the surface imperfection. As the applied electric field is increased, solitons (or "bullets") become detached from the wings of the butterfly, which then rapidly propagate throughout the system. The main observations that emerge from the model, including the formation and structure of butterflies, bullets, and stripes, as well as the role of surface imperfections and the strength of the applied field, are consistent with our own experimental findings presented here for nematic LCs confined between two chemically treated parallel plates. △ Less

Submitted 16 October, 2022; originally announced October 2022.

Prediction of Electronic Properties of Radical-Containing Polymers at Coarse-Grained Resolutions

Authors: Riccardo Alessandri, Juan J. de Pablo

Abstract: The properties of soft electronic materials depend on the coupling of electronic and conformational degrees of freedom over a wide range of spatiotemporal scales. Description of such properties requires multiscale approaches capable of, at the same time, accessing electronic properties and sampling the conformational space of soft materials. This could in principle be realized by connecting the co… ▽ More The properties of soft electronic materials depend on the coupling of electronic and conformational degrees of freedom over a wide range of spatiotemporal scales. Description of such properties requires multiscale approaches capable of, at the same time, accessing electronic properties and sampling the conformational space of soft materials. This could in principle be realized by connecting the coarse-grained (CG) methodologies required for adequate conformational sampling to conformationally-averaged electronic property distributions via backmapping to atomistic-resolution level models and repeated quantum-chemical calculations. Computational demands of such approaches, however, have hindered their application in high-throughput computer-aided soft materials discovery. Here, we present a method that, combining machine learning and CG techniques, can replace traditional backmapping-based approaches without sacrificing accuracy. We illustrate the method for an emerging class of soft electronic materials, namely non-conjugated, radical-containing polymers, promising materials for all-organic energy storage. Supervised machine learning models are trained to learn the dependence of electronic properties on polymer conformation at CG resolutions. We then parametrize CG models that retain electronic structure information, simulate CG condensed phases, and predict the electronic properties of such phases solely from the CG degrees of freedom. We validate our method by comparing it against a full backmapping-based approach, and find good agreement between both methods. This work demonstrates the potential of the proposed method to accelerate multiscale workflows, and provides a framework for the development of CG models that retain electronic structure information. △ Less

Submitted 24 April, 2023; v1 submitted 5 September, 2022; originally announced September 2022.

Current rectification by nanoparticles in bipolar nanopores

Authors: Andrés Córdoba, Joan M. Montes de Oca, Johnson Dhanasekaran, Seth B. Darling, Juan J. de Pablo

Abstract: Bipolar nanochannels comprising two domains of positively and negatively charged walls along the pore axis are known to rectify current when exposed to an electric potential bias. We find that addition of charged nanoparticles can increase rectification considerably, by approximately one order of magnitude. Two bipolar channel geometries are considered here; their behavior is examined at rest and… ▽ More Bipolar nanochannels comprising two domains of positively and negatively charged walls along the pore axis are known to rectify current when exposed to an electric potential bias. We find that addition of charged nanoparticles can increase rectification considerably, by approximately one order of magnitude. Two bipolar channel geometries are considered here; their behavior is examined at rest and under the influence of a negative bias and a positive bias, respectively. We do so by relying on a molecular-level model of the electrolyte solution in the channels. The large increase in current rectification can be explained by the inherent electric field that charged nanoparticles generate within the channel. This effect is found to be largely dependent on the pore's geometry, its charge distribution, and the sign of the nanoparticles' charge, thereby offering new opportunities for design of engineered nanopore membrane-nanoparticle systems for energy storage. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 30 pages, 4 figures

Systematic Modification of Functionality in Disordered Elastic Networks Through Free Energy Surface Tailoring

Authors: Dan Mendels, Fabian Byléhn, Timothy W. Sirk, Juan J. de Pablo

Abstract: Advances in manufacturing and characterization of complex molecular systems have created a need for new methods for design at molecular length scales. Emerging approaches are increasingly relying on the use of Artificial Intelligence (AI), and the training of AI models on large data libraries. This paradigm shift has led to successful applications, but shortcomings related to interpretability and… ▽ More Advances in manufacturing and characterization of complex molecular systems have created a need for new methods for design at molecular length scales. Emerging approaches are increasingly relying on the use of Artificial Intelligence (AI), and the training of AI models on large data libraries. This paradigm shift has led to successful applications, but shortcomings related to interpretability and generalizability continue to pose challenges. Here, we explore an alternative paradigm in which AI is combined with physics-based considerations for molecular and materials engineering. Specifically, collective variables, akin to those used in enhanced sampled simulations, are constructed using a machine learning (ML) model trained on data gathered from a single system. Through the ML-constructed collective variables, it becomes possible to identify critical molecular interactions in the system of interest, the modulation of which enables a systematic tailoring of the system's free energy landscape. To explore the efficacy of the proposed approach, we use it to engineer allosteric regulation, and uniaxial strain fluctuations in a complex disordered elastic network. Its successful application in these two cases provides insights regarding how functionality is governed in systems characterized by extensive connectivity, and points to its potential for design of complex molecular systems. △ Less

Submitted 15 August, 2022; v1 submitted 8 July, 2022; originally announced July 2022.

Catapulting of topological defects through elasticity bands in active nematics

Authors: Nitin Kumar, Rui Zhang, Steven A. Redford, Juan J. de Pablo, Margaret L. Gardel

Abstract: Active materials are those in which individual, uncoordinated local stresses drive the material out of equilibrium on a global scale. Examples of such assemblies can be seen across scales from schools of fish to the cellular cytoskeleton and underpin many important biological processes. Synthetic experiments that recapitulate the essential features of such active systems have been the object of st… ▽ More Active materials are those in which individual, uncoordinated local stresses drive the material out of equilibrium on a global scale. Examples of such assemblies can be seen across scales from schools of fish to the cellular cytoskeleton and underpin many important biological processes. Synthetic experiments that recapitulate the essential features of such active systems have been the object of study for decades as their simple rules allow us to elucidate the physical underpinnings of collective motion. One system of particular interest has been active nematic liquid crystals (LCs). Because of their well understood passive physics, LCs provide a rich platform to interrogate the effects of active stress. The flows and steady state structures that emerge in an active LCs have been understood to result from a competition between nematic elasticity and the local activity. However most investigations of such phenomena consider only the magnitude of the elastic resistance and not its peculiarities. Here we investigate a nematic liquid crystal and selectively change the ratio of the material's splay and bend elasticities. We show that increases in the nematic's bend elasticity specifically drives the material into an exotic steady state where elongated regions of acute bend distortion or "elasticity bands" dominate the structure and dynamics. We show that these bands strongly influence defect dynamics, including the rapid motion or "catapulting" along the disintegration of one of these bands thus converting bend distortion into defect transport. Thus, we report a novel dynamical state resultant from the competition between nematic elasticity and active stress. △ Less

Submitted 31 March, 2022; originally announced April 2022.

Journal ref: Soft Matter 18, 5271-5281 (2022)

Transformation between elastic dipoles, quadrupoles, octupoles and hexadecapoles driven by surfactant self-assembly in nematic emulsion

Authors: Bohdan Senyuk, Ali Mozaffari, Kevin Crust, Rui Zhang, Juan J. de Pablo, Ivan I. Smalyukh

Abstract: Emulsions comprising isotropic fluid drops within a nematic host are of interest for applications ranging from biodetection to smart windows, which rely on changes of molecular alignment structures around the drops in response to chemical, thermal, electric and other stimuli. We show that absorption or desorption of trace amounts of common surfactants can drive continuous transformations of elasti… ▽ More Emulsions comprising isotropic fluid drops within a nematic host are of interest for applications ranging from biodetection to smart windows, which rely on changes of molecular alignment structures around the drops in response to chemical, thermal, electric and other stimuli. We show that absorption or desorption of trace amounts of common surfactants can drive continuous transformations of elastic multipoles induced by the droplets within the uniformly aligned nematic host. Out-of-equilibrium dynamics of director structures emerge from a controlled self-assembly or desorption of different surfactants at the drop-nematic interfaces, with ensuing forward and reverse transformations between elastic dipoles, quadrupoles, octupoles and hexadecapoles. We characterize inter-transformations of droplet-induced surface and bulk defects, probe elastic pair interactions and discuss emergent prospects for fundamental science and applications of the reconfigurable nematic emulsions. △ Less

Submitted 6 January, 2022; originally announced January 2022.

Stability and molecular pathways to the formation of spin defects in silicon carbide

Authors: Elizabeth M. Y. Lee, Alvin Yu, Juan J. de Pablo, Giulia Galli

Abstract: Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate--the divacancy complex (VV) in silic… ▽ More Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate--the divacancy complex (VV) in silicon carbide (SiC). Using atomistic models, enhanced sampling simulations, and density functional theory calculations, we find that VV formation is a thermally activated process that competes with the conversion of silicon ($V_{Si}$) to carbon monovacancies ($V_{C}$), and that VV reorientation can occur without dissociation. We also find that increasing the concentration of $V_{Si}$ relative to $V_{C}$ favors the formation of divacancies. Moreover, we identify pathways to create spin defects consisting of antisite-double vacancy complexes and determine their electronic properties. The detailed view of the mechanisms that underpin the formation and dynamics of spin defects presented here may facilitate the realization of qubits in an industrially relevant material. △ Less

Submitted 14 September, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: 23 pages (main), 16 pages (supporting info), 6 main figures

Collective Variables for Free Energy Surface Tailoring -- Understanding and Modifying Functionality in Systems Dominated by Rare Events

Authors: Dan Mendels, Juan J. de Pablo

Abstract: We introduce a method for elucidating and modifying the functionality of systems dominated by rare events that relies on the automated tuning of their underlying free energy surface. The proposed approach seeks to construct collective variables (CVs) that encode the essential information regarding the rare events of the system of interest. The appropriate CVs are identified using Harmonic Linear D… ▽ More We introduce a method for elucidating and modifying the functionality of systems dominated by rare events that relies on the automated tuning of their underlying free energy surface. The proposed approach seeks to construct collective variables (CVs) that encode the essential information regarding the rare events of the system of interest. The appropriate CVs are identified using Harmonic Linear Discriminant Analysis (HLDA), a machine-learning based method that is trained solely on data collected from short ordinary simulations in the relevant metastable states of the system. Utilizing the interpretable form of the resulting CVs, the critical interaction potentials that determine the system's rare transitions are identified and purposely modified to tailor the free energy surface in manner that alters functionality as desired. The applicability of the method is illustrated in the context of three different systems, thereby demonstrating that thermodynamic and kinetic properties can be tractably modified with little to no prior knowledge or intuition. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Stress-activated Constraints in Dense Suspension Rheology

Authors: Abhinendra Singh, Grayson L. Jackson, Michael van der Naald, Juan J. de Pablo, Heinrich M. Jaeger

Abstract: Dispersing small particles in a liquid can produce surprising behaviors when the solids fraction becomes large: rapid shearing drives these systems out of equilibrium and can lead to dramatic increases in viscosity (shear-thickening) or even solidification (shear jamming). These phenomena occur above a characteristic onset stress when particles are forced into frictional contact. Here we show via… ▽ More Dispersing small particles in a liquid can produce surprising behaviors when the solids fraction becomes large: rapid shearing drives these systems out of equilibrium and can lead to dramatic increases in viscosity (shear-thickening) or even solidification (shear jamming). These phenomena occur above a characteristic onset stress when particles are forced into frictional contact. Here we show via simulations how this can be understood within a framework that abstracts details of the forces acting at particle-particle contacts into general stress-activated constraints on relative particle movement. We find that focusing on just two constraints, affecting sliding and rolling at contact, can reproduce the experimentally observed shear thickening behavior quantitatively, despite widely different particle properties, surface chemistries, and suspending fluids. Within this framework parameters such as coefficients of sliding and rolling friction can each be viewed as a proxy for one or more forces of different physical or chemical origin, while the parameter magnitudes indicate the relative importance of the associated constraint. In this way, a new link is established that connects features observable in macroscale rheological measurements to classes of constraints arising from micro- or nano-scale properties. △ Less

Submitted 22 August, 2021; originally announced August 2021.

Comments: 7 Pages, 4 Figures

Defect Spirograph: Dynamical Behavior of Defects in Spatially Patterned Active Nematics

Authors: Ali Mozaffari, Rui Zhang, Noe Atzin, Juan J. de Pablo

Abstract: Topological defects in active liquid crystals can be confined by introducing gradients of activity. Here, we examine the dynamical behavior of two defects confined by a sharp gradient of activity that separates an active circular region and a surrounding passive nematic material. Continuum simulations are used to explain how the interplay among energy injection into the system, hydrodynamic intera… ▽ More Topological defects in active liquid crystals can be confined by introducing gradients of activity. Here, we examine the dynamical behavior of two defects confined by a sharp gradient of activity that separates an active circular region and a surrounding passive nematic material. Continuum simulations are used to explain how the interplay among energy injection into the system, hydrodynamic interactions, and frictional forces governs the dynamics of topologically required self-propelling $+1/2$ defects. Our findings are rationalized in terms of a phase diagram for the dynamical response of defects in terms of activity and frictional damping strength. Different regions of the underlying phase diagram correspond to distinct dynamical modes, namely immobile defects (ID), steady rotation of defects (SR), bouncing defects (TB), bouncing-cruising defects (BC), dancing defects (DA), and multiple defects with irregular dynamics (MD). These dynamic states raise the prospect of generating synchronized defect arrays for microfluidic applications. △ Less

Submitted 22 May, 2021; originally announced May 2021.

Poisson's ratio and angle bending in spring networks

Authors: Nidhi Pashine, Daniel R. Reid, Meng Shen, Juan J. de Pablo, Sidney R. Nagel

Abstract: The Poisson's ratio of a spring network system has been shown to depend not only on the geometry but also on the relative strength of angle-bending forces in comparison to the bond-compression forces in the system. Here we derive the very simple analytic result that in systems where the spring interaction strength is equal to the bond-reorientation interaction, the Poisson's ratio identically goes… ▽ More The Poisson's ratio of a spring network system has been shown to depend not only on the geometry but also on the relative strength of angle-bending forces in comparison to the bond-compression forces in the system. Here we derive the very simple analytic result that in systems where the spring interaction strength is equal to the bond-reorientation interaction, the Poisson's ratio identically goes to zero and is independent of the network geometry. △ Less

Submitted 7 April, 2021; originally announced April 2021.

Comments: 3 pages, 2 figures

Machine learning active-nematic hydrodynamics

Authors: Jonathan Colen, Ming Han, Rui Zhang, Steven A. Redford, Linnea M. Lemma, Link Morgan, Paul V. Ruijgrok, Raymond Adkins, Zev Bryant, Zvonimir Dogic, Margaret L. Gardel, Juan J. De Pablo, Vincenzo Vitelli

Abstract: Hydrodynamic theories effectively describe many-body systems out of equilibrium in terms of a few macroscopic parameters. However, such hydrodynamic parameters are difficult to derive from microscopics. Seldom is this challenge more apparent than in active matter where the energy cascade mechanisms responsible for autonomous large-scale dynamics are poorly understood. Here, we use active nematics… ▽ More Hydrodynamic theories effectively describe many-body systems out of equilibrium in terms of a few macroscopic parameters. However, such hydrodynamic parameters are difficult to derive from microscopics. Seldom is this challenge more apparent than in active matter where the energy cascade mechanisms responsible for autonomous large-scale dynamics are poorly understood. Here, we use active nematics to demonstrate that neural networks can extract the spatio-temporal variation of hydrodynamic parameters directly from experiments. Our algorithms analyze microtubule-kinesin and actin-myosin experiments as computer vision problems. Unlike existing methods, neural networks can determine how multiple parameters such as activity and elastic constants vary with ATP and motor concentration. In addition, we can forecast the evolution of these chaotic many-body systems solely from image-sequences of their past by combining autoencoder and recurrent networks with residual architecture. Our study paves the way for artificial-intelligence characterization and control of coupled chaotic fields in diverse physical and biological systems even when no knowledge of the underlying dynamics exists. △ Less

Submitted 23 June, 2020; originally announced June 2020.

Comments: SI Movie 1: https://www.youtube.com/watch?v=9WzIT7OG9pY SI Movie 2: https://youtu.be/Trc4RyU7-dw SI Movie 3: https://youtu.be/Epm_P_EakH8

Shear thickening and jamming of dense suspensions: the "roll" of friction

Authors: Abhinendra Singh, Christopher Ness, Ryohei Seto, Juan J de Pablo, Heinrich M Jaeger

Abstract: Particle-based simulations of discontinuous shear thickening (DST) and shear jamming (SJ) suspensions are used to study the role of stress-activated constraints, with an emphasis on resistance to gear-like rolling. Rolling friction decreases the volume fraction required for DST and SJ, in quantitative agreement with real-life suspensions with adhesive surface chemistries and "rough" particle shape… ▽ More Particle-based simulations of discontinuous shear thickening (DST) and shear jamming (SJ) suspensions are used to study the role of stress-activated constraints, with an emphasis on resistance to gear-like rolling. Rolling friction decreases the volume fraction required for DST and SJ, in quantitative agreement with real-life suspensions with adhesive surface chemistries and "rough" particle shapes. It sets a distinct structure of the frictional force network compared to only sliding friction, and from a dynamical perspective leads to an increase in the velocity correlation length, in part responsible for the increased viscosity. The physics of rolling friction is thus a key element in achieving a comprehensive understanding of strongly shear-thickening materials. △ Less

Submitted 7 May, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 5 pages, 4 figures. +Supplemental Material--- --- Supplementary Material also consists of the following movies: ---https://youtu.be/wgrfN7jiu1g (Force network in 3D at volume fraction 0.45, at high stress limit, with and without rolling friction) --- https://youtu.be/SMzTB7CsRsY (Force network in 2D at packing fraction 0.7, at high stress limit, with and without rolling friction)

Journal ref: Phys. Rev. Lett. 124, 248005 (2020)

Structure and dynamics of hydrodynamically interacting finite-size Brownian particles in a spherical cavity: spheres and cylinders

Authors: Jiyuan Li, Xikai Jiang, Abhinendra Singh, Olle G. Heinonen, Juan P. Hernández-Ortiz, Juan J. de Pablo

Abstract: The structure and dynamics of confined suspensions of particles of arbitrary shape is of interest in multiple disciplines, from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cyl… ▽ More The structure and dynamics of confined suspensions of particles of arbitrary shape is of interest in multiple disciplines, from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an Immersed-Boundary General geometry Ewald-like method to capture lubrication and long-range hydrodynamics, and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than $10\%$, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion towards the walls, and a sub-diffusive regime --caused by crowding-- in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells. △ Less

Submitted 12 February, 2020; originally announced February 2020.

Comments: 10 pages, 7 Figures+SI

Journal ref: J. Chem. Phys. 152, 204109 (2020)

Shape induced segregation and anomalous particle transport under spherical confinement

Authors: Abhinendra Singh, Jiyuan Li, Xikai Jiang, Juan P. Hernández-Ortiz, Heinrich M. Jaeger, Juan J. de Pablo

Abstract: Colloid or nanoparticle mobility under confinement is of central importance to a wide range of physical and biological processes. Here, we introduce a minimal model of particles in a hydrodynamic continuum to examine how particle shape and concentration affect the transport of particles in spherical confinement. Specifically, an immersed boundary-General geometry Ewald-like approach is adopted to… ▽ More Colloid or nanoparticle mobility under confinement is of central importance to a wide range of physical and biological processes. Here, we introduce a minimal model of particles in a hydrodynamic continuum to examine how particle shape and concentration affect the transport of particles in spherical confinement. Specifically, an immersed boundary-General geometry Ewald-like approach is adopted to simulate the dynamics of spheres and cylinders under the influence of short-and long-range fluctuating hydrodynamic interactions with appropriate non-slip conditions at the confining walls. An efficient $\it{O(N)}$ parallel finite element algorithm is used, thereby allowing simulations at high concentrations, while a Chebyshev polynomial approximation is implemented in order to satisfy the fluctuation-dissipation theorem. A concentration-dependent anomalous diffusion is observed for suspended particles. It is found that introducing cylinders in a background of spheres, i.e. particles with a simple degree of anisotropy, has a pronounced influence on the structure and dynamics of the particles. First, increasing the fraction of cylinders induces a particle segregation effect, where spheres are pushed towards the wall and cylinders remain near the center of the cavity. This segregation leads to lower mobility for the spheres relative to that encountered in a system of pure spheres at the same volume fraction. Second, the diffusive-to-anomalous transition and the degree of anomaly--quantified by the power-law exponent in the mean square displacement vs. time relation-both increase as the fraction of cylinders becomes larger. These findings are of relevance for studies of diffusion in the cytoplasm, where proteins exhibit a distribution of size and shapes that could lead to some of the effects identified in the simulations reported here. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: 9 pages, 5 figures

Journal ref: Physics of Fluids 32, 053307 (2020)

Sculpting stable structures in pure liquids

Authors: Tadej Emeršič, Rui Zhang, Žiga Kos, Simon Čopar, Natan Osterman, Juan J. de Pablo, Uroš Tkalec

Abstract: Pure liquids in thermodynamic equilibrium are structurally homogeneous. In liquid crystals, flow and light pulses are used to create reconfigurable domains with polar order. Moreover, through careful engineering of concerted microfluidic flows and localized opto-thermal fields, it is possible to achieve complete control over the nucleation, growth, and shape of such domains. Experiments, theory, a… ▽ More Pure liquids in thermodynamic equilibrium are structurally homogeneous. In liquid crystals, flow and light pulses are used to create reconfigurable domains with polar order. Moreover, through careful engineering of concerted microfluidic flows and localized opto-thermal fields, it is possible to achieve complete control over the nucleation, growth, and shape of such domains. Experiments, theory, and simulations indicate that the resulting structures can be stabilized indefinitely, provided the liquids are maintained in a controlled non-equilibrium state. The resulting sculpted liquids could find applications in microfluidic devices for selective encapsulation of solutes and particles into optically active compartments that interact with external stimuli. △ Less

Submitted 2 January, 2020; originally announced January 2020.

Journal ref: Science Advances, vol. 5, eaav4283 (2019)

Structuring Stress for Active Materials Control

Authors: Rui Zhang, Steven A. Redford, Paul V. Ruijgrok, Nitin Kumar, Ali Mozaffari, Sasha Zemsky, Aaron R. Dinner, Vincenzo Vitelli, Zev Bryant, Margaret L. Gardel, Juan J. de Pablo

Abstract: Active materials are capable of converting free energy into mechanical work to produce autonomous motion, and exhibit striking collective dynamics that biology relies on for essential functions. Controlling those dynamics and transport in synthetic systems has been particularly challenging. Here, we introduce the concept of spatially structured activity as a means to control and manipulate transpo… ▽ More Active materials are capable of converting free energy into mechanical work to produce autonomous motion, and exhibit striking collective dynamics that biology relies on for essential functions. Controlling those dynamics and transport in synthetic systems has been particularly challenging. Here, we introduce the concept of spatially structured activity as a means to control and manipulate transport in active nematic liquid crystals consisting of actin filaments and light-sensitive myosin motors. Simulations and experiments are used to demonstrate that topological defects can be generated at will, and then constrained to move along specified trajectories, by inducing local stresses in an otherwise passive material. These results provide a foundation for design of autonomous and reconfigurable microfluidic systems where transport is controlled by modulating activity with light. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Ideal isotropic auxetic networks from random networks

Authors: Daniel R. Reid, Nidhi Pashine, Alec S. Bowen, Sidney R. Nagel, Juan J. de Pablo

Abstract: Auxetic materials are characterized by a negative Poisson's ratio, $\mathrmν$. As the Poisson's ratio becomes negative and approaches the lower isotropic mechanical limit of $\mathrm{ν= -1}$, materials show enhanced resistance to impact and shear, making them suitable for applications ranging from robotics to impact mitigation. Past experimental efforts aimed at reaching the $\mathrm{ν= -1}$ limit… ▽ More Auxetic materials are characterized by a negative Poisson's ratio, $\mathrmν$. As the Poisson's ratio becomes negative and approaches the lower isotropic mechanical limit of $\mathrm{ν= -1}$, materials show enhanced resistance to impact and shear, making them suitable for applications ranging from robotics to impact mitigation. Past experimental efforts aimed at reaching the $\mathrm{ν= -1}$ limit have resulted in highly anisotropic materials, which show a negative Poisson's ratio only when subjected to deformations along specific directions. Isotropic designs have only attained moderately auxetic behavior, or have led to structures that cannot be manufactured in 3D. Here, we present a design strategy to create isotropic structures from disordered networks that leads to Poisson's ratios as low as $\mathrm{ν= -0.98}$. The materials conceived through this approach are successfully fabricated in the laboratory and behave as predicted. The Poisson's ratio $\mathrmν$ is found to depend on network structure and bond strengths; this sheds light on the structural motifs that lead to auxetic behavior. The ideas introduced here can be generalized to 3D, a wide range of materials, and a spectrum of length scales, thereby providing a general platform that could impact technology. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: 16 pages, 6 figures

Tunable structure and dynamics of active liquid crystals

Authors: Nitin Kumar, Rui Zhang, Juan J. de Pablo, Margaret L. Gardel

Abstract: Active materials are capable of converting free energy into directional motion, giving rise to striking dynamical phenomena. Developing a general understanding of their structure in relation to the underlying non-equilibrium physics would provide a route towards control of their dynamic behavior and pave the way for potential applications. The active system considered here consists of a quasi-two-… ▽ More Active materials are capable of converting free energy into directional motion, giving rise to striking dynamical phenomena. Developing a general understanding of their structure in relation to the underlying non-equilibrium physics would provide a route towards control of their dynamic behavior and pave the way for potential applications. The active system considered here consists of a quasi-two-dimensional sheet of short ($\approx$ 1 $μ$m) actin filaments driven by myosin-II motors. By adopting a concerted theoretical and experimental strategy, new insights are gained into the non-equilibrium properties of active nematics over a wide range of internal activity levels. In particular, it is shown that topological defect interactions can be led to transition from attractive to repulsive as a function of initial defect separation and relative orientation. Furthermore, by examining the +1/2 defect morphology as a function of activity, it is found that the apparent elastic properties of the system (the ratio of bend-to-splay elastic moduli) are altered considerably by increased activity, leading to an effectively lower bend elasticity. At high levels of activity, the topological defects that decorate the material exhibit a liquid-like structure, and adopt preferred orientations depending on their topological charge. Taken together, these results suggest that it should be possible to tune internal stresses in active nematic systems with the goal of designing out-of-equilibrium structures with engineered dynamic responses. △ Less

Submitted 16 October, 2018; originally announced October 2018.

Journal ref: Science Advances 4(10), eaat7779 (2018)

Evolutionary strategy for inverse charge measurements of dielectric particles

Authors: Xikai Jiang, Jiyuan Li, Victor Lee, Heinrich M. Jaeger, Olle G. Heinonen, Juan J. de Pablo

Abstract: We report a computational strategy to obtain the charges of individual dielectric particles from experimental observation of their interactions as a function of time. This strategy uses evolutionary optimization to minimize the difference between trajectories extracted from experiment and simulated trajectories based on many-particle force fields. The force fields include both Coulombic interactio… ▽ More We report a computational strategy to obtain the charges of individual dielectric particles from experimental observation of their interactions as a function of time. This strategy uses evolutionary optimization to minimize the difference between trajectories extracted from experiment and simulated trajectories based on many-particle force fields. The force fields include both Coulombic interactions and dielectric polarization effects that arise due to particle-particle charge mismatch and particle-environment dielectric contrast. The strategy was applied to systems of free falling charged granular particles in vacuum, where electrostatic interactions are the only driving forces that influence the particles' motion. We show that when the particles' initial positions and velocities are known, the optimizer requires only an initial and final particle configuration of a short trajectory in order to accurately infer the particles' charges; when the initial velocities are unknown and only the initial positions are given, the optimizer can learn from multiple frames along the trajectory to determine the particles' initial velocities and charges. While the results presented here offer a proof-of-concept demonstration of the proposed ideas, the proposed strategy could be extended to more complex systems of electrostatically charged granular matter. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Auxetic metamaterials from disordered networks

Authors: Daniel R Reid, Nidhi Pashine, Justin M Wozniak, Heinrich M Jaeger, Andrea J Liu, Sidney R Nagel, Juan J de Pablo

Abstract: Recent theoretical work suggests that systematic pruning of disordered networks consisting of nodes connected by springs can lead to materials that exhibit a host of unusual mechanical properties. In particular, global properties such as the Poisson's ratio or local responses related to deformation can be precisely altered. Tunable mechanical responses would be useful in areas ranging from impact… ▽ More Recent theoretical work suggests that systematic pruning of disordered networks consisting of nodes connected by springs can lead to materials that exhibit a host of unusual mechanical properties. In particular, global properties such as the Poisson's ratio or local responses related to deformation can be precisely altered. Tunable mechanical responses would be useful in areas ranging from impact mitigation to robotics and, more generally, for creation of metamaterials with engineered properties. However, experimental attempts to create auxetic materials based on pruning-based theoretical ideas have not been successful. Here we introduce a new and more realistic model of the networks, which incorporates angle-bending forces and the appropriate experimental boundary conditions. A sequential pruning strategy of select bonds in this model is then devised and implemented that enables engineering of specific mechanical behaviors upon deformation, both in the linear and non-linear regimes. In particular, it is shown that the Poisson's ratio can be tuned to arbitrary values. The model and concepts discussed here are validated by preparing physical realizations of the networks designed in this manner, which are produced by laser cutting two-dimensional sheets and are found to behave as predicted. Furthermore, by relying on optimization algorithms, we exploit the networks' susceptibility to tuning to design networks that posses a distribution of stiffer and more compliant bonds, and whose auxetic behavior is even greater than that of homogeneous networks. Taken together, the findings reported here serve to establish that pruned networks represent a promising platform for the creation of novel mechanical metamaterials. △ Less

Submitted 2 October, 2017; originally announced October 2017.

Interplay of Structure, Elasticity and Dynamics in Actin-Based Nematic Materials

Authors: Rui Zhang, Nitin Kumar, Jennifer Ross, Margaret L. Gardel, Juan J. de Pablo

Abstract: Achieving control and tunability of lyotropic materials has been a long-standing goal of liquid crystal research. Here we show that the elasticity of a liquid crystal system consisting of a dense suspension of semiflexible biopolymers can be manipulated over a relatively wide range of elastic moduli. Specifically, thin films of actin filaments are assembled at an oil-water interface. At sufficient… ▽ More Achieving control and tunability of lyotropic materials has been a long-standing goal of liquid crystal research. Here we show that the elasticity of a liquid crystal system consisting of a dense suspension of semiflexible biopolymers can be manipulated over a relatively wide range of elastic moduli. Specifically, thin films of actin filaments are assembled at an oil-water interface. At sufficiently high concentrations, one observes the formation of a nematic phase riddled with $\pm1/2$ topological defects, characteristic of a two-dimensional nematic system. As the average filament length increases, the defect morphology transitions from a U-shape into a V-shape, indicating the relative increase of the material's bend over splay modulus. Furthermore, through the sparse addition of rigid microtubule filaments, one can further control the liquid crystal elasticity. We show how the material's bend constant can be raised linearly as a function of microtubule filament density, and present a simple means to extract absolute values of the elastic moduli from purely optical observations. Finally, we demonstrate that it is possible to predict not only the static structure of the material, including its topological defects, but also the evolution of the system into dynamically arrested states. Despite the non-equilibrium nature of the material, our continuum model, which couples structure and hydrodynamics, is able to capture the annihilation and movement of defects over long time scales. Thus, we have experimentally realized a lyotropic liquid crystal system that can be truly engineered, with tunable mechanical properties, and a theoretical framework to capture its structure, mechanics, and dynamics. △ Less

Submitted 14 September, 2017; originally announced September 2017.

Journal ref: PNAS 115(2), E-124-E133 (2018)

Low-temperature anomalies of a vapor deposited glass

Authors: Beatriz Seoane, Daniel R. Reid, Juan J. de Pablo, Francesco Zamponi

Abstract: We investigate the low temperature properties of two-dimensional Lennard-Jones glass films, prepared in silico both by liquid cooling and by physical vapor deposition. We identify deep in the solid phase a crossover temperature $T^*$, at which slow dynamics and enhanced heterogeneity emerge. Around $T^*$, localized defects become visible, leading to vibrational anomalies as compared to standard so… ▽ More We investigate the low temperature properties of two-dimensional Lennard-Jones glass films, prepared in silico both by liquid cooling and by physical vapor deposition. We identify deep in the solid phase a crossover temperature $T^*$, at which slow dynamics and enhanced heterogeneity emerge. Around $T^*$, localized defects become visible, leading to vibrational anomalies as compared to standard solids. We find that on average, $T^*$ decreases in samples with lower inherent structure energy, suggesting that such anomalies will be suppressed in ultra-stable glass films, prepared both by very slow liquid cooling and vapor deposition. △ Less

Submitted 12 January, 2018; v1 submitted 14 September, 2017; originally announced September 2017.

Comments: 10 pages including appendices, 8 figures. Version accepted for Physical Review Materials

Journal ref: Phys. Rev. Materials 2, 015602 (2018)

Controlling mechanical properties of a polymeric glass

Authors: George J. Papakonstantopoulos, Juan J. de Pablo

Abstract: In this work we use molecular simulations to examine methods of controlling mechanical properties of polymeric glass materials such as elastic moduli, mechanical heterogeneity as well as their glass transition temperature. We study filled and unfilled polymers and examine the effect of particle size, volume fraction and polymer-particle interactions. We identify a relationship between mobility and… ▽ More In this work we use molecular simulations to examine methods of controlling mechanical properties of polymeric glass materials such as elastic moduli, mechanical heterogeneity as well as their glass transition temperature. We study filled and unfilled polymers and examine the effect of particle size, volume fraction and polymer-particle interactions. We identify a relationship between mobility and dynamic heterogeneity with elastic moduli and glass transition temperature. △ Less

Submitted 12 October, 2016; originally announced October 2016.

Mesoscale Structure of Chiral Nematic Shells

Authors: Ye Zhou, Ashley Guo, Rui Zhang, Julio C. Armas-Perez, Jose A. Martinez-Gonzalez, Mohammad Rahimi, Monirosadat Sadati, Juan J. de Pablo

Abstract: There is considerable interest in understanding and controlling topological defects in nematic liquid crystals (LCs). Confinement, in the form of droplets, has been particularly effective in that regard. Here, we employ the Landau-de Gennes method to explore the geometrical frustration of nematic order in shell geometries, and focus on chiral materials. By varying the chirality and thickness in un… ▽ More There is considerable interest in understanding and controlling topological defects in nematic liquid crystals (LCs). Confinement, in the form of droplets, has been particularly effective in that regard. Here, we employ the Landau-de Gennes method to explore the geometrical frustration of nematic order in shell geometries, and focus on chiral materials. By varying the chirality and thickness in uniform shells, we construct a phase diagram that includes tetravalent structures, bipolar structures (BS), bent structures and radial spherical structures (RSS). It is found that, in uniform shells, the BS-to-RSS structural transition, in response to both chirality and shell geometry, is accompanied by an abrupt change of defect positions, implying a potential use for chiral nematic shells as sensors. Moreover, we investigate thickness heterogeneity in shells and demonstrate that non-chiral and chiral nematic shells exhibit distinct equilibrium positions of their inner core that are governed by shell chirality c. △ Less

Submitted 27 July, 2016; originally announced July 2016.

Structural Transitions in Fibers of Bent-Core Liquid Crystals from Field-Theory Monte Carlo Simulations

Authors: N. Atzin, O. Guzmán, J. J. de Pablo

Abstract: Fibers of bent-core liquid crystals present an internal structure of a rolled smectic layer and can be used as optical waveguides. We used a field-theoretical Monte Carlo simulation to analyze the internal configuration of such fibers as a function of the radial coordinate and to study their equilibrium sates. In contrast to previous studies, we analyzed the fully nonlinear model proposed in [Bail… ▽ More Fibers of bent-core liquid crystals present an internal structure of a rolled smectic layer and can be used as optical waveguides. We used a field-theoretical Monte Carlo simulation to analyze the internal configuration of such fibers as a function of the radial coordinate and to study their equilibrium sates. In contrast to previous studies, we analyzed the fully nonlinear model proposed in [Bailey et al., Phys. Rev. E, 2007, 75, 031701] and revised in [Pérez-Ortiz et al., Phys. Rev. E, 2011, 84, 011701]. We found that, due to the non-differentiable character of such model, the Euler-Lagrange equations are not able to find all equilibrium states. Our Monte Carlo procedure identified both differentiable and non-differentiable equilibria and first order transitions between them. In all cases, the equilibrium states show inhomogenous configurations that display a boundary layer. This methodology can by applied to other models of liquid crystals that have more degrees of freedom, including those with non-differentiable minima. The equilibrium structures found can be used as inputs to model the transmission of light along the liquid crystal fibers. △ Less

Submitted 3 August, 2016; v1 submitted 21 July, 2016; originally announced July 2016.

Comments: Second version. 12 pages, 7 figures in PDF format. Added references to typical values of spontaneous polarization of LC phases of fully organic and metallorganic bent-core mesogens

Nanoscopic Characterization of DNA within Hydrophobic Pores:Thermodynamics and Kinetics

Authors: Fernando J. A. L. Cruz, Juan J. de Pablo, José P. B. Mota

Abstract: The energetic and transport properties of a double-stranded DNA dodecamer encapsulated in hydrophobic carbon nanotubes are probed employing two limiting nanotube diameters, D=4 nm and D=3 nm, corresponding to (51,0) and (40,0) zig-zag topologies, respectively. It is observed that the thermodynamically spontaneous encapsulation in the 4 nm nanopore (40 kJ/mol) is annihilated when the solid diameter… ▽ More The energetic and transport properties of a double-stranded DNA dodecamer encapsulated in hydrophobic carbon nanotubes are probed employing two limiting nanotube diameters, D=4 nm and D=3 nm, corresponding to (51,0) and (40,0) zig-zag topologies, respectively. It is observed that the thermodynamically spontaneous encapsulation in the 4 nm nanopore (40 kJ/mol) is annihilated when the solid diameter narrows down to 3 nm, and that the confined DNA termini directly contact the hydrophobic walls with no solvent slab in-between. During the initial moments after confinement (2-3 ns),the biomolecule translocates along the nano pore's inner volume according to Fick's law (t) with a self-diffusion coefficient D=1.713 x 10-9m2/s, after which molecular diffusion assumes a single-file type mechanism (t1/2). As expected, diffusion is anisotropic, with the pore main axis as the preferred direction, but an in-depth analysis shows that the instantaneous velocity probabilities are essentially identical along the x, y and z directions. The 3D velocity histogram shows a maximum probability located at v=30.8 m/s, twice the observed velocity for a single-stranded three nucleotide DNA encapsulated in comparable armchair geometries (v=16.7 m/s, D=1.36-1.89 nm). Because precise physiological conditions (310 K, [NaCl] = 134 mM) are employed throughout, the present study establishes a landmark for the development of next generation in vivo drug delivery technologies based on carbon nanotubes as encapsulation agents. △ Less

Submitted 4 May, 2016; originally announced May 2016.

Comments: arXiv admin note: text overlap with arXiv:1605.01317

Journal ref: Biochem. Eng. J. 104, 2015, 41-47

Endohedral confinement of a DNA dodecamer onto pristine carbon nanotubes and the stability of the canonical B form

Authors: Fernando J. A. L. Cruz, Juan J. de Pablo, José P. B. Mota

Abstract: Although carbon nanotubes are potential candidates for DNA encapsulation and subsequent delivery of biological payloads to living cells, the thermodynamical spontaneity of DNA encapsulation under physiological conditions is still a matter of debate. Using enhanced sampling techniques, we show for the first time that, given a sufficiently large carbon nanotube, the confinement of a double-stranded… ▽ More Although carbon nanotubes are potential candidates for DNA encapsulation and subsequent delivery of biological payloads to living cells, the thermodynamical spontaneity of DNA encapsulation under physiological conditions is still a matter of debate. Using enhanced sampling techniques, we show for the first time that, given a sufficiently large carbon nanotube, the confinement of a double-stranded DNA segment, 5'-D(*CP*GP*CP*GP*AP*AP*TP*TP*CP*GP*CP*G)-3', is thermodynamically favourable under physiological environments (134 mM, 310 K, 1 bar), leading to DNA-nanotube hybrids with lower free energy than the unconfined biomolecule. A diameter threshold of 3 nm is established below which encapsulation is inhibited. The confined DNA segment maintains its translational mobility and exhibits the main geometrical features of the canonical B form. To accommodate itself within the nanopore, the DNA end-to-end length increases from 3.85 nm up to approximately 4.1 nm, due to a 0.3 nm elastic expansion of the strand termini. The canonical Watson-Crick H-bond network is essentially conserved throughout encapsulation, showing that the contact between the DNA segment and the hydrophobic carbon walls results in minor rearrangements of the nucleotides H-bonding. The results obtained here are paramount to the usage of carbon nanotubes as encapsulation media for next generation drug delivery technologies. △ Less

Submitted 4 May, 2016; originally announced May 2016.

Journal ref: J. Chem. Phys. 140, 225103 (2014)

Monte Carlo simulation study of diblock copolymer self assembly

Authors: George J. Papakonstantopoulos, Kostas Ch. Daoulas, Marcus Muller, Juan J. de Pablo

Abstract: A technique is presented which maps the parameters of a bead spring model, using the Flory Huggins theory, to a specific experimental system. By keeping only necessary details, for the description of these systems, the mapping procedure turns into an estimation of a few characteristic parameters. An asset of this technique is that it is simple to apply and captures the behavior of block copolymer… ▽ More A technique is presented which maps the parameters of a bead spring model, using the Flory Huggins theory, to a specific experimental system. By keeping only necessary details, for the description of these systems, the mapping procedure turns into an estimation of a few characteristic parameters. An asset of this technique is that it is simple to apply and captures the behavior of block copolymer phase separation. In our study this mapping technique is utilized in conjunction with a Monte Carlo (MC) algorithm to perform simulations on block copolymer systems. The microphase separation is investigated in the bulk and under confinement, on unpatterned and patterned surfaces. △ Less

Submitted 18 April, 2016; originally announced April 2016.

Inherent Structure Energy is a Good Indicator of Molecular Mobility in Glasses

Authors: Julian Helfferich, Ivan Lyubimov, Daniel Reid, Juan J. de Pablo

Abstract: Glasses produced via physical vapor deposition can display greater kinetic stability and lower enthalpy than glasses prepared by liquid cooling. While the reduced enthalpy has often been used as a measure of the stability, it is not obvious whether dynamic measures of stability provide the same view. Here, we study dynamics in vapor-deposited and liquid-cooled glass films using molecular simulatio… ▽ More Glasses produced via physical vapor deposition can display greater kinetic stability and lower enthalpy than glasses prepared by liquid cooling. While the reduced enthalpy has often been used as a measure of the stability, it is not obvious whether dynamic measures of stability provide the same view. Here, we study dynamics in vapor-deposited and liquid-cooled glass films using molecular simulations of a bead-spring polymer model as well as a Lennard-Jones binary mixture in two and three dimensions. We confirm that the dynamics in vapor-deposited glasses is indeed slower than in ordinary glasses. We further show that the inherent structure energy is a good reporter of local dynamics, and that aged systems and glasses prepared by cooling at progressively slower rates exhibit the same behavior as vapor-deposited materials when they both have the same inherent structure energy. These findings suggest that the stability inferred from measurements of the energy is also manifested in dynamic observables, and they strengthen the view that vapor deposition processes provide an effective strategy for creation of stable glasses. △ Less

Submitted 4 April, 2016; originally announced April 2016.

Comments: 10 pages, 8 figures, to be submitted to Soft Matter

Turning Statistical Physics Models Into Materials Design Engines

Authors: Marc Z. Miskin, Gurdaman S. Khaira, Juan J. de Pablo, Heinrich M. Jaeger

Abstract: Despite the success statistical physics has enjoyed at predicting the properties of materials for given parameters, the inverse problem, identifying which material parameters produce given, desired properties, is only beginning to be addressed. Recently, several methods have emerged across disciplines that draw upon optimization and simulation to create computer programs that tailor material respo… ▽ More Despite the success statistical physics has enjoyed at predicting the properties of materials for given parameters, the inverse problem, identifying which material parameters produce given, desired properties, is only beginning to be addressed. Recently, several methods have emerged across disciplines that draw upon optimization and simulation to create computer programs that tailor material responses to specified behaviors. However, so far the methods developed either involve black-box techniques, in which the optimizer operates without explicit knowledge of the material's configuration space, or they require carefully tuned algorithms with applicability limited to a narrow subclass of materials. Here we introduce a formalism that can generate optimizers automatically by extending statistical mechanics into the realm of design. The strength of this new approach lies in its capability to transform statistical models that describe materials into optimizers to tailor them. By comparing against standard black-box optimization methods, we demonstrate how optimizers generated by this formalism can be faster and more effective, while remaining straightforward to implement. The scope of our approach includes new possibilities for solving a variety of complex optimization and design problems concerning materials both in and out of equilibrium. △ Less

Submitted 19 October, 2015; originally announced October 2015.

Self-consistent description of electrokinetic phenomena in particle-based simulations

Authors: Juan P. Hernandez-Ortiz, Juan J. de Pablo

Abstract: A new computational method is presented for study suspensions of charged soft particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are inclu… ▽ More A new computational method is presented for study suspensions of charged soft particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are included through a Green's function formalism. An Ewald-like split is adopted in order to satisfy arbitrary boundary conditions for the Stokeslet and Poisson Green functions, thereby providing a formalism that is applicable to any geometry and that can be extended to deformable objects. The convection-diffusion equation for the continuum ions is solved simultaneously considering Nernst-Planck diffusion. The method can be applied to systems at equilibrium and far from equilibrium. Its applicability is demonstrated in the context of electrokinetic motion, where it is shown that the ionic clouds associated with individual particles can be severely altered by the flow and concentration, leading to intriguing cooperative effects. △ Less

Submitted 12 April, 2015; originally announced April 2015.

Comments: 5 pages, 4 figures

Optimizing directed self-assembled morphology

Authors: Jian Qin, Gurdaman S. Khaira, Yongrui Su, Grant P. Garner, Marc Miskin, Heinrich M. Jaeger, Juan J. de Pablo

Abstract: Directed assembly of block polymers is rapidly becoming a viable strategy for lithographic patterning of nanoscopic features. One of the key attributes of directed assembly is that an underlying chemical or topographic substrate pattern used to direct assembly need not exhibit a direct correspondence with the sought after block polymer morphology, and past work has largely relied on trial-and-erro… ▽ More Directed assembly of block polymers is rapidly becoming a viable strategy for lithographic patterning of nanoscopic features. One of the key attributes of directed assembly is that an underlying chemical or topographic substrate pattern used to direct assembly need not exhibit a direct correspondence with the sought after block polymer morphology, and past work has largely relied on trial-and-error approaches to design appropriate patterns. In this work, a computational evolutionary strategy is proposed to solve this optimization problem. By combining the Cahn-Hilliard equation, which is used to find the equilibrium morphology, and the covariance-matrix evolutionary strategy, which is used to optimize the combined outcome of particular substrate-copolymer combinations, we arrive at an efficient method for design of substrates leading to non-trivial, desirable outcomes. △ Less

Submitted 2 August, 2013; originally announced August 2013.

Comments: 9 pages, 8 figures

Model Vapor-Deposited Glasses: Growth Front and Composition Effects

Authors: Ivan Lyubimov, Mark D. Ediger, Juan J. de Pablo

Abstract: A growing body of experimental work indicates that physical vapor deposition provides an effective route for preparation of stable glasses, whose properties correspond in some cases to those expected for glasses that have been aged for thousands of years. In this work, model binary glasses are prepared in a process inspired by physical vapor deposition, in which particles are sequentially added to… ▽ More A growing body of experimental work indicates that physical vapor deposition provides an effective route for preparation of stable glasses, whose properties correspond in some cases to those expected for glasses that have been aged for thousands of years. In this work, model binary glasses are prepared in a process inspired by physical vapor deposition, in which particles are sequentially added to the free surface of a growing film. The resulting glasses are shown to be more stable than those prepared by gradual cooling from the liquid phase. However, it is also shown that the composition of the resulting glass, which is difficult to control in physical vapor deposition simulations of thin films, plays a significant role on the physical characteristics of the material. That composition dependence leads to a re- evaluation of previous results from simulations of thinner films than those considered here, where the equivalent age of the corresponding glasses was overestimated. The simulations presented in this work, which correspond to films that are approximately 38 molecular diameters thick, also enable analysis of the devitrification mechanism of vapor-deposited glasses. Consistent with experiments, it is found that this mechanism consists of a mobility front that propagates from the free interface into the interior of the films. Eliminating surface mobility eliminates this route of transformation into the supercooled liquid. △ Less

Submitted 19 June, 2013; originally announced June 2013.

Comments: 11 pages, 6 figures. Submitted to J.Chem.Phys

Density of states of a binary Lennard-Jones Glass

Authors: Roland Faller, Juan J. de Pablo

Abstract: We calculate the density of states of a binary Lennard-Jones glass using a recently proposed Monte Carlo algorithm. Unlike traditional molecular simulation approaches, the algorithm samples distinct configurations according to self-consistent estimates of the density of states, thereby giving rise to uniform internal-energy histograms. The method is applied to simulate the equilibrium, low-tempe… ▽ More We calculate the density of states of a binary Lennard-Jones glass using a recently proposed Monte Carlo algorithm. Unlike traditional molecular simulation approaches, the algorithm samples distinct configurations according to self-consistent estimates of the density of states, thereby giving rise to uniform internal-energy histograms. The method is applied to simulate the equilibrium, low-temperature thermodynamic properties of a widely studied glass former consisting of a binary mixture of Lennard-Jones particles. We show how a density-of-states algorithm can be combined with particle identity swaps and configurational bias techniques to study that system. Results are presented for the energy and entropy below the mode coupling temperature. △ Less

Submitted 28 May, 2003; originally announced May 2003.

Comments: 6 pages, 3 figures, accepted by J Chem Phys

Journal ref: J Chem Phys 119(13), 6883-6894 (2003)

Local Elastic Constants in Thin Films of an FCC Crystal

Authors: K. Van Workum, J. J. de Pablo

Abstract: In this work we present a formalism for the calculation of the local elastic constants in inhomogeneous systems based on a method of planes. Unlike previous work, this formalism does not require the partitioning of the system into a set of finite volumes over which average elastic constants are calculated. Results for the calculation of the local elastic constants of a nearest neighbor Lennard-J… ▽ More In this work we present a formalism for the calculation of the local elastic constants in inhomogeneous systems based on a method of planes. Unlike previous work, this formalism does not require the partitioning of the system into a set of finite volumes over which average elastic constants are calculated. Results for the calculation of the local elastic constants of a nearest neighbor Lennard-Jones fcc crystal in the bulk and in a thin film are presented. The local constants are calculated at exact planes of the (001) face of the crystal. The average elastic constants of the bulk system are also computed and are consistent with the local constants. Additionally we present the local stress profiles in the thin film when a small uniaxial strain is applied. The resulting stress profile compares favorably with the stress profile predicted via the local elastic constants. The surface melting of a model for argon for which experimental and simulation data are available is also studied within the framework of this formalism. △ Less

Submitted 11 October, 2002; originally announced October 2002.

Comments: 7 pages, 10 figures

Report number: UW-CHE-02-10-11

Potential of Mean Force between a Spherical Particle Suspended in a Nematic Liquid Crystal and a Substrate

Authors: Evelina B. Kim, Roland Faller, Qiliang Yan, Nicholas L. Abbott, Juan J. de Pablo

Abstract: We consider a system where a spherical particle is suspended in a nematic liquid crystal confined between two walls. We calculate the liquid-crystal mediated potential of mean force between the sphere and a substrate by means of Monte Carlo simulations. Three methods are used: a traditional Monte Carlo approach, umbrella sampling, and a novel technique that combines canonical expanded ensemble s… ▽ More We consider a system where a spherical particle is suspended in a nematic liquid crystal confined between two walls. We calculate the liquid-crystal mediated potential of mean force between the sphere and a substrate by means of Monte Carlo simulations. Three methods are used: a traditional Monte Carlo approach, umbrella sampling, and a novel technique that combines canonical expanded ensemble simulations with a recently proposed density of states formalism. The latter method offers clear advantages in that it ensures good sampling of phase space without prior knowledge of the energy landscape of the system. The resulting potential of mean force, computed as a function of the normal distance between the sphere and a surface, suggests that the sphere is attracted to the surface, even in the absence of attractive molecular interactions. △ Less

Submitted 1 August, 2002; originally announced August 2002.

Comments: 11 pages, 9 figures, accepted by J Chem Phys

Journal ref: J Chem Phys 117(16) 7781-7787 (2002)

Density of States Monte Carlo Method for Simulation of Fluids

Authors: Qiliang Yan, Roland Faller, Juan J. de Pablo

Abstract: A Monte Carlo method based on a density-of-states sampling is proposed for study of arbitrary statistical mechanical ensembles in a continuum. A random walk in the two-dimensional space of particle number and energy is used to estimate the density of states of the system; this density of states is continuously updated as the random walk visits individual states. The validity and usefulness of th… ▽ More A Monte Carlo method based on a density-of-states sampling is proposed for study of arbitrary statistical mechanical ensembles in a continuum. A random walk in the two-dimensional space of particle number and energy is used to estimate the density of states of the system; this density of states is continuously updated as the random walk visits individual states. The validity and usefulness of the method are demonstrated by applying it to the simulation of a Lennard-Jones fluid. Results for its thermodynamic properties, including the vapor-liquid phase coexistence curve, are shown to be in good agreement with high-accuracy literature data. △ Less

Submitted 25 January, 2002; originally announced January 2002.

Comments: 5 pages, 3 figures, accepted by J Chem Phys

Journal ref: J. Chem. Phys. 116(20) 8745-8749 (2002)