Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations

Many objects of nanoscopic dimensions involve fluid-tethered chains interfaces. These systems are of interest for basic science and several applications, in particular for design of nanodevices for specific purposes. We review recent developments of theoretical methods in this area of research and i...

Повний опис

Збережено в:

Бібліографічні деталі
Дата:	2014
Автори:	Sokołowski, S., Ilnytskyi, Ja., Pizio, O.
Формат:	Стаття
Мова:	English
Опубліковано:	Інститут фізики конденсованих систем НАН України 2014
Назва видання:	Condensed Matter Physics
Онлайн доступ:	http://dspace.nbuv.gov.ua/handle/123456789/153450
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations / S. Sokołowski, Ja. Ilnytskyi, O. Pizio // Condensed Matter Physics. — 2014. — Т. 17, № 1. — С. 12601:1-34. — Бібліогр.: 215 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-153450
record_format	dspace
spelling	irk-123456789-1534502019-06-15T01:25:47Z Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations Sokołowski, S. Ilnytskyi, Ja. Pizio, O. Many objects of nanoscopic dimensions involve fluid-tethered chains interfaces. These systems are of interest for basic science and several applications, in particular for design of nanodevices for specific purposes. We review recent developments of theoretical methods in this area of research and in particular of density functional (DF) approaches, which provide important insights into microscopic properties of such interfaces. The theories allow to describe dependence of adsorption, wettability, solvation forces and electric interfacial phenomena on thermodynamic states and on characteristics of tethered chains. Computer simulations for the problems in question are overviewed as well. Theoretical results are discussed in relation to simulation results and some experimental observations. Iнтерфейси мiж флюїдом та приєднаними полiмерними ланцюжками зустрiчаються у численних наноскопiчних об’єктах. Дослiдження їх властивостей є цiкавим як з точки зору фундаментальних наукових знань, так i для низки застосувань, зокрема для проектування нанопристроїв iз спецiальним призначенням. В оглядi подано сучасний стан досягнень теоретичних пiдходiв в цiй областi (насамперед методiв, що базуються на функцiоналi густини), якi просунули наше розумiння мiкроскопiчних властивостей таких iнтерфейсiв. Цi теоретичнi дослiдження дозволяють описати залежностi адсорбцiйних, змочувальних та сольватацiйних сил та особливостi електро-iнтерфейсних явищ вiд термодинамiчного стану та характеристик приєднаних ланцюжкiв. Зроблено також огляд комп’ютерних симуляцiй таких систем, якi, як i теоретичнi дослiдження, спiвставлено з наявними експериментальними даними. 2014 Article Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations / S. Sokołowski, Ja. Ilnytskyi, O. Pizio // Condensed Matter Physics. — 2014. — Т. 17, № 1. — С. 12601:1-34. — Бібліогр.: 215 назв. — англ. 1607-324X arXiv:1403.1491 DOI:10.5488/CMP.17.12601 PACS: 61.20.-p, 61.25.-f, 61.25.he, 61.25.Em, 68.15.+e, 68.43.-h, 68.47.Pe, 78.30.cd http://dspace.nbuv.gov.ua/handle/123456789/153450 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	English
description	Many objects of nanoscopic dimensions involve fluid-tethered chains interfaces. These systems are of interest for basic science and several applications, in particular for design of nanodevices for specific purposes. We review recent developments of theoretical methods in this area of research and in particular of density functional (DF) approaches, which provide important insights into microscopic properties of such interfaces. The theories allow to describe dependence of adsorption, wettability, solvation forces and electric interfacial phenomena on thermodynamic states and on characteristics of tethered chains. Computer simulations for the problems in question are overviewed as well. Theoretical results are discussed in relation to simulation results and some experimental observations.
format	Article
author	Sokołowski, S. Ilnytskyi, Ja. Pizio, O.
spellingShingle	Sokołowski, S. Ilnytskyi, Ja. Pizio, O. Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations Condensed Matter Physics
author_facet	Sokołowski, S. Ilnytskyi, Ja. Pizio, O.
author_sort	Sokołowski, S.
title	Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations
title_short	Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations
title_full	Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations
title_fullStr	Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations
title_full_unstemmed	Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations
title_sort	description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations
publisher	Інститут фізики конденсованих систем НАН України
publishDate	2014
url	http://dspace.nbuv.gov.ua/handle/123456789/153450
citation_txt	Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations / S. Sokołowski, Ja. Ilnytskyi, O. Pizio // Condensed Matter Physics. — 2014. — Т. 17, № 1. — С. 12601:1-34. — Бібліогр.: 215 назв. — англ.
series	Condensed Matter Physics
work_keys_str_mv	AT sokołowskis descriptionoffluidtetheredchainsinterfacesadvancesindensityfunctionaltheoriesandofflatticecomputersimulations AT ilnytskyija descriptionoffluidtetheredchainsinterfacesadvancesindensityfunctionaltheoriesandofflatticecomputersimulations AT pizioo descriptionoffluidtetheredchainsinterfacesadvancesindensityfunctionaltheoriesandofflatticecomputersimulations
first_indexed	2025-07-14T04:36:47Z
last_indexed	2025-07-14T04:36:47Z
_version_	1837595684339974144
fulltext	Condensed Matter Physics, 2014, Vol. 17, No 1, 12601: 1–34 DOI: 10.5488/CMP.17.12601 http://www.icmp.lviv.ua/journal Review Article Description of interfaces of fluid-tethered chains: advances in density functional theories and off-lattice computer simulations S. Sokołowski1 , J. Ilnytskyi2, O. Pizio3¤ 1 Department for the Modelling of Physico-Chemical Processes, Maria Curie-Sklodowska University, 33 Gliniana St., Lublin, Poland 2 Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, 1 Svientsitskii St., 79011 Lviv, Ukraine 3 Instituto de Química, Universidad Nacional Autonoma de México, Circuito Exterior, Ciudad Universitaria, 04510 México, D.F., México Received May 28, 2013, in final form October 24, 2013 Many objects of nanoscopic dimensions involve fluid-tethered chain interfaces. These systems are of interest for basic science and for several applications, in particular for design of nanodevices for specific purposes. We review recent developments of theoretical methods in this area of research and in particular of density functional (DF) approaches, which provide important insights into microscopic properties of such interfaces. The theories permit to describe the dependence of adsorption, wettability, solvation forces and electric interfacial phenomena on thermodynamic states and on characteristics of tethered chains. Computer simulations for the problems in question are overviewed as well. Theoretical results are discussed in relation to simulation results and to some experimental observations. Key words: tethered chains, polymer brushes, fluid, adsorption, wetting, electrolyte solutions, differential capacitance PACS: 61.20.-p, 61.25.-f, 61.25.he, 61.25.Em, 68.15.+e, 68.43.-h, 68.47.Pe, 78.30.cd 1. Introduction Systems involving interfaces formed by polymers tethered onto a solid surface and a fluid contacting them are widespread in nature and are the subject of many studies. One of the reasons is that the grafted layers of chains find several practical applications, for example in efficient colloid stabilization, protec- tion and prevention of surfaces from undesirable adsorptions, manufacturing of functional surfaces and nanoscale fluidic channels, and many others. Interesting properties of such systems have also stimulated the development of several methods of synthesis for preparation of assemblies of end-grafted polymers. However, up to now, the development of new preparation methods is mainly based on empirical knowl- edge or intuition. From the very beginning, experimental studies of tethered chains are accompanied by the develop- ment of theoretical methods of their description. They involve thermodynamic and structural properties of systems containing uncharged and charged species under the effect of various external factors includ- ing confinement at nanoscale as well. The purpose of this article is to present a brief review of the statistical mechanical DFmethods and off- lattice computer simulations of classical fluids in contact with the walls modified by tethered polymers, ¤E-mail: oapizio@gmail.com © S. Sokołowski, J. Ilnytskyi, O. Pizio, 2014 12601-1 http://dx.doi.org/10.5488/CMP.17.12601 http://www.icmp.lviv.ua/journal S. Sokołowski, J. Ilnytskyi, O. Pizio to discuss some recent achievements in this area of research, as well as to provide some prospects for fur- ther investigations. After introductory remarks, in the next section (section 2), we briefly outline different types of theoretical approaches used for the description of tethered chains-fluid interfaces. We beginwith the scaling theories, but our emphasis is on the theories that are based on the microscopic description and the use of the intermolecular potentials as the input, i.e., on the integral equations, self-consistent field and density functional theories. Particular attention is paid to the DF approaches. In section 3 we discuss one version of the density functional theory (DFT) that has been successful in describing several phenomena occurring within tethered chains-fluid interfaces. Some illustrative examples are given in section 4. Section 5 extends the approach of section 4 to the case of systems having electrostatic inter- actions. We discuss the theory and then, in section 6, present some applications concerning structural, thermodynamic and electric properties of the systems with tethered chains. Mesoscopic simulations by the dissipative particle dynamics method and the (semi-) atomistic simulations using the Monte Carlo (MC) and Molecular Dynamics (MD) techniques are reviewed in sections 7.1 and 7.2, respectively. The manuscript ends with a short summary in section 8. 2. Overview of theoretical methods Essential features of the systems involving end-tethered polymers onto a solid and in contact with a fluid can be captured by modelling at atomistic or semi-atomistic (at groups of atoms) scale by using contemporary theories for interfaces between a solid and a complex fluid. A link between fundamentals and applications can be reached by combining theoretical (classical and quantum) tools, intensive com- puter simulations and laboratory experiments. Current studies at a simple level of modelling put focus on the development of a theory for various phenomena occurring in the systems of our interest. Computer simulations provide data to confront theoretical predictions and simulation results. The most important, they serve as predictive tools to guide elaboration of interfaces with desired functional properties. Theoretical tools to investigate systems with tethered chains are similar to the methods designed to study nonuniform polymeric systems. They can be classified into several groups by using different crite- ria. One can divide them according to the way in which translational degrees of freedom are considered. Namely, either lattice or off-lattice models can be distinguished. On the other hand, the classification can be made according to theoretical approaches used to describe a model. One can distinguish the theo- ries based on the application of scaling arguments [1, 2], integral equations approaches, self-consistent field (SCF) as well as the DF methods. The scaling approaches provide qualitatively correct trends of the system behavior over a variety of conditions. The integral equations and the SCF approaches are more ample in the sense that they provide information (e.g., details on the structure of interfaces) which scal- ing approaches do not. Perhaps, the most accurate are the DFTs that besides an insight into the structure also provide an elegant way to study phase transitions in the systems. Finally, one can classify the models into the ones with or without electrostatic interactions. The systems involving tethered molecules have been also extensively studied using computer simulation methods. 2.1. Scaling approaches. A few comments A system to which the below described theoretical methods apply, is a set of long flexible polymer chains which are anchored on a grafting substrate through a special single end group. In the majority of cases in question, the substrate is taken as flat. The degree of polymerization determines the “chain length” M . On the other hand, the important parameter of the system is the grafting density ½C, i.e., the fraction of grafted sites per unit area. Usually, this system is put in contact with a solvent. These three parameters (M , ½C and solvent quality) determine several regimes of behavior of polymer layer hav- ing distinct static (structural and thermodynamic) and dynamic properties. The systems of this type are intrinsically characterized by a set of mesoscopic length scales which makes the polymer layers under different conditions amenable to the scaling analysis. Of course, scaling relation imply that all character- istic lenghts are much larger than any atomistic scales, e.g., the range over which the grafting substrate induces density oscillations or the persistence length. 12601-2 Fluid-grafted chains interfaces It is commonly assumed that all other segments of chain molecules, besides grafted end group, are not attracted to the surface. Within the scaling theories, the structure of a grafted polymer layer is character- ized by the average height and the values of gyration radius along selected axis of grafted molecules. If the grafting density is very low, different coils do not overlap and practically do not interact. Each chain occupies roughly a half-sphere with a radius comparable to the Flory radius for a polymer coil in a good solvent, such that values of gyration radius along selected axis in the z-direction normal to the surface Rgz and in the directions perpendicular to the z-axis Rg? differ only by their prefactor in the power law relating size to the chain length [1, 2], Rgz / Rg?/M º (1) with º ¼ 3/5, M is the number of segments forming a chain. This behavior is the so-called “mushroom regime”. For denser grafted layers, the chains overlap close to the grafting surface, and are bound to stretch away from the surface, like the bristles in a brush. This situation is quite different from a typical behavior of polymer chains in a solution. If the grafting density ½C is high (the polymer brush regime) several char- acteristic length scales appear. According to the scaling approach by Alexander [3] each grafted polymer chain is considered as a linear array of M 0 “blobs”, M 0 ÆM/nblob , where nblob is the number of chain segments per blob (only segements of a single chain are inside a blob), nblob / » 1/º g (»g Æ ½ ¡1/2 C is the average distance between the grafting points or the blob diameter). The theory by Alexander, see reference [3] for details, yields the following scaling law for the tethered polymer layer height, h, h/ ½ (1/º¡1)/2 C M ¼ ½ 1/3 C M . (2) This is a very important result showing that the brush height at a chosen fixed grafted density, is pro- portional to M , while the unstretched chain dimension, Rg, grows slower. Consequently, the properties of chains under such conditions may be expected to differ much from unstretched chains in a solution. This scaling argument implies a step-like density profile of the segments. It is difficult to expect a such situation for more realistic models. In fact, entropy favors a nonuniform blob picture. A nonuniform stretching of chains can occur, some chains can turn backwards and some segments can be located anywhere within the brush (of course, there is a maximum distance from the surface within which the segments can be found). Due to nonuni- formity, with an increasing distance from the surface, »g grows from the value ½¡1/2 C assumed by Alexan- der, to larger values. In particular, for semi-dilute grafted layers, grafted at z Æ 0, as documented by de Gennes [2], the density profile »g(z) Æ »g(0)[½s(z)/½s(0)℄ [º/(1¡3º)℄, is related to the segment local density ½s(z), that is assumed to be parabolic, ½s(z) Æ ½s(0)[1¡ (z/h) 2 ℄ in the self-consistent field theory commented below. The scaling approaches do not take into account the interactions in the system explicitly. Consequently, specific correlations between particular species of importance in various phenomena cannot be captured in such framework. 2.2. Integral equations Development of theoretical approaches in which intermolecular potentials are the input to the the- ory is particularly important for discussing phase transitions. One of the first theoretical approaches to describe uniform, as well as nonuniform polymeric systems that explicitly takes into account interaction potentials between all the species was based on the integral equation theory of fluids. Common integral equation theories of atomic fluids were extended to describe polymers by introducing the so-called Poly- mer Reference Interaction Site Model (PRISM) [5–7] that generalized the Reference Interaction Site Model (RISM) by Chandler [8]. For the case of a fluid of homopolymers built ofM indentical sites (segments), the PRISM integral equation in the Fourier space reads, � h(k)Æ � !(k) � (k)[ � !(k)Å½ M � h(k)℄, (3) where ½ M is the site density (½ M ÆM½, with ½ being the density of molecules), � h(k) is the Fourier trans- form of the site-site total correlation function, � (k) is the Fourier transform of the site-site direct correla- tion function. Since all the sites along a chain molecule are equivalent, it follows that � h(k)Æ � h ®,° (k) and 12601-3 S. Sokołowski, J. Ilnytskyi, O. Pizio � (k)Æ � ®,° (k). Moreover, � !(k)Æ 1 M M X ®,°Æ1 � ! ®,° (k), (4) is the single chain structure factor consisting of Fourier transforms of the site-site intramolecular corre- lation functions. ! ®,° (k). A generalization of equations (3)–(4) to the cases of mixtures of molecules with different types of sites is available in the literature [5]. After introducing closure approximations, the system of the PRISM integral equations can be solved numerically, yielding the correlation functions that contain a complete information about the microscopic structure of the molecular system. By computing the integrals of the correlation functions, one can obtain thermodynamic properties. The PRISM-based approaches have been used to study a type of model systems involving small solid particles of different geometry (e.g., spherical), coated with polymer chains, see e.g., closely related works [7, 9–11]. One recent interesting example is due to Jayaraman and Schweizer [12–14] who used the PRISM theory to study model polymer nanocomposites consisting of nanoparticles (called “fillers”) with a few grafted chains embedded in a polymer “matrix” (solvent or melt). Coated nanoparticles with a relatively small and controlled number of grafted polymers in this type of mixtures are of interest for various ap- plications. Both the microscopic structure and the phase behavior of such complex systems are sensitive to the effective interactions between species. Those can be tuned by chemical means: synthetic routes of laboratory preparation determine the attraction strength between bare nanoparticles, tether grafting density and matrix chain length. Other important factors are the total packing density and mixture com- position, along with the diameter ratio between nanoparticles and polymerizing monomers belonging to grafted chains and to the polymer solvent. We would just like to mention three scenarios leading to the development of different structures. In general terms, there is a competition between the assumed nanoparticle-nanoparticle attraction, steric repulsion, i.e., excluded volume effects between grafted tethers, and filler-homopolymer correla- tions. The matrix subsystem can induce a depletion attraction between fillers. Under certain circum- stances, the tether-tether repulsion can sterically “screen” the attraction between nanoparticles and if the depletion type entropic attraction is weak, then the fillers form a well-dispersed fluid of a filler in the homopolymer matrix. If the intertether repulsion is weak and the attraction between fillers is strong enough, then the macroscopic phase separation can be induced by depletion matrix effects in addition to bare attraction between nanoparticles. In this case, the matrix-rich and the tethered filler-rich coexisting phases can be formed. Still another, “intermediate” physical situation occurs if nanoparticle-nanoparticle interactions are rather strong and a sufficient steric stabilization due to tether repulsion exists to preclude macrophase demixing, then aggregate formation of microphase-type ordering can take place. So far the PRISM theories have not been applied to describe situations where the grafting occurs on planar surfaces, or, in general, on the surfaces with linear dimension much larger than the size of a polymer segment. 2.3. Self-consistent field approaches One of the popular approaches for the systems of our interest is the self-consistent field theory (SCFT), developed by Edwards [15], DiMarzio [16], Vrij [17] and Helfand and Tagami [18, 19]. The key idea behind SCFT is to approximate an ensemble of interacting polymers by the system of mutually non-interacting polymer molecules in an effective, position dependent field. The SCFT methods have been successful in predicting adsorption of polymers, tethered polymers, polymer blends and copolymers, but they are intrinsically incapable of capturing fine details of the system structure at short-range scale [20–22]. Due to space limitations, we would like to restrict the discussion of the method to the statement of the problem rather than to the enumeration of several principal results. For illustrative purposes, we begin with a general statistical mechanical arguments and consider the Gaussian chain model. The system is a blend of N ¸ kinds of linear homopolymers. The length of the homopolymers of a given kind ¸ isM ¸ . For Gaussian chains, the potential energy of the system is [23] ¯Vpot Æ X ¸® M ¸ ¡1 X iÆ1 3 2l 2 ¸ ¯ ¯ ¯ r ¸ ®,iÅ1 ¡r ¸ ®,i ¯ ¯ ¯ 2 Å ¯ 2 X ¸®i,¸ 0 ® 0 i 0 u ¸¸ 0 ³ r ¸ ®,i ¡r ¸ 0 ® 0 ,i 0 ´ Å¯ X ®¸i v ¸ ³ r ¸ ®,i ´ , (5) 12601-4 Fluid-grafted chains interfaces the sum over ®¸i is over all the species, over all the particles of the given species and over all the sites of a molecule of a given kind and l ¸ is the corresponding Kuhn length; ¯ Æ 1/kT , u ¸¸ 0 and v ¸ are the segment-segment interaction potential and the external field, respectively. The grand partition function for the model in question is, ¥Æ Y ¸ 1 X n ¸ Æ0 exp � ¯ X ¸ ¹ ¸ n ¸ ! ¡ 3/2¼l 2 ¸ ¢ 3(M ¸ ¡1)/2 Q ¸ h 2 n ¸ (n ¸ !)¤ 3n ¸ M ¸ ¸ i ( Y ¸ n ¸ Y ®Æ1 M ¸ Y iÆ1 Z dr ¸ ®,i ) exp ¡ ¡¯Vpot ¢ , (6) where ¹ ¸ is the chemical potential per chain of species ¸, and ¤ ¸ Æ h/(2¼m ¸ kT ) 1/2 is the thermal de Broglie length for monomers of species ¸. It is convenient to define the effective field functionW ¸ (r) via the relation Z dr X ¸ � ½ ¸ (r)W ¸ (r)Æ¯ X ¸ ¹ ¸ n ¸ ¡¯ X ¸®i v ³ r ¸ ®i ´ , (7) where the density operator for the species ¸ has been introduced � ½ ¸ (r)Æ X ®i ± ³ r¡r ¸ ®,i ´ . (8) Then, the partition function (6) can be rewritten as ¥ Æ Y ¸ 1 X n ¸ Æ0 " Y ¸ 2 n ¸ (n ¸ !)¤ 3n ¸ M ¸ ¸ # ¡1 � 3 2¼l 2 ¸ ! 3(M ¸ ¡1)/2 £ ( Y ¸ n ¸ Y ®Æ1 M ¸ Y iÆ1 Z dr ¸ ®,i ) exp " ¡¯VpotÅ Z dr X ¸ � ½ ¸ (r)W ¸ (r) # . (9) The segment density profile is obtained from ¥ by using the definition ½ ¸ (r)Æ h � ½ ¸ (r)i Æ ± ln¥ ±W ¸ (r) . (10) Then, the free energy follows from the Legendre transformation [23]. In general terms, the density-field relation (10) can be inverted analytically for an ideal system to yield ±W ¸ (r) as a function of segment density profiles ½ ¸ (r). Such a task can be realized for the Gaussian chain ideal model that permits to reach the principal equations of the self-consistent field method. Tang and Freed [24] have chosen the model with a vanishing pair interaction, u ¸¸ 0 Æ 0, as the ideal system such that the potential energy given by equation (5) contains solely the connectivity of chains contribution as well as the term coming from an arbitrary external field. The partition function of the model factorizes as follows: ¥ideal Æ Y ¸ ¥ ¸ ideal Æ Y ¸ 1 X n ¸ Æ0 Q n ¸ ¸ n ¸ ! Æ Y ¸ exp(Q ¸ ), (11) whereQ ¸ denotes the partition function for a single Gaussian chain belonging to the species ¸ under the effect of the arbitrary external fieldW ¸ . The density profile of segments can be easily obtained, ½ ¸ (r)Æ ± ln¥ ±W ¸ (r) Æ ¡ 3/2¼l 2 ¸ ¢ 3(M ¸ ¡1)/2 2¤ 3M ¸ ¸ � M ¸ Y iÆ1 Z dr ¸ i !" X i ± ³ r¡r ¸ i ´ # £exp " ¡ M ¸ ¡1 X i 3 2l 2 ¸ ¯ ¯ ¯ r ¸ iÅ1 ¡r ¸ i ¯ ¯ ¯ 2 Å M ¸ X i W ¸ ³ r ¸ i ´ # (12) such that Z dr½ ¸ (r)ÆM ¸ Q ¸ . (13) 12601-5 S. Sokołowski, J. Ilnytskyi, O. Pizio It is more convenient to develop the following derivation presented below using a continuous chain notation, or, in other words, substituting sums over monomer indices by integrals. The single chain par- tition function,Q (the subscript ¸ is dropped to simplify the notation) reads Q Æ ¡ 1/2¤ 3M ¢ Z drdr 0 G ¡ r,r 0 ;M ,0j[W ℄ ¢ , (14) where G(r,r 0 ;M ,0j[W ℄) is the unnormalized end-vector distribution function, given as G(r,r 0 ;M ,0j[W ℄) Æ r(M )Ær Z r(0)Ær 0 D[r(¿)℄exp 2 4 ¡ 3 2l 2 M Z 0 d¿jr(¿)j 2 Å M Z 0 d¿W [r(¿)℄ 3 5 . (15) The chain contour variable ¿ changes between 0 and M . In general, the end-vector distribution function G(r,r 0 ;¿,¿ 0 j[W ℄) gives the number of chain configurations that start at r0,¿0 and end at r,¿. The path in- tegral representation of the end-vector distribution function given by equation (15) may be alternatively transformed into the partial differential equation · � �¿ ¡ l 2 6 r 2 r ¡W (r) ¸ G(r,r 0 ;¿,¿ 0 j[W ℄)Æ ±(r¡r 0 )±(¿¡¿ 0 ) (16) that describes a diffusion in an external field W (r). This type of equations serves as the starting point in several research works performed by using the SCFT, see e.g., [9, 25]. Additional boundary conditions can be used for a set of specific problems that include confinement effects. It is worth mentioning that the inversion of a field-density relation can be explored in a wider context for systems of interacting chain molecules. The present simple illustration intrinsically leads us to the explanation of the density functional approaches. 2.4. Density functional approaches. Generalities The powerful group of methods to describe nonuniform complex fluids is built based on the ideas of DFT for simple fluids [26]. Perhaps, the first formulation of the DFT for polyatomic fluids was given by Chandler et al. [27, 28]. They extended the DFT of nonuniform simple fluids to the systems composed of polyatomic species, considered in the framework of the RISM approach [29, 30]. By using Legendre trans- forms, these authors proved the existence of a free energy functional with the local densities referring to the interaction sites rather than to molecular coordinates. The methodology proposed by Chandler et al. [27, 28] retains mathematical simplicity of the traditional theory for atomic fluids, but it requires the “site-site” direct correlation function as an input. Another class of DFT approaches for polymeric fluids was proposed by Woodward [31], who used the concept of a weighted density approximation, common in the DFTs for simple fluids [32]. Woodward applied the theory to a fluid of particles made of tangentially bonded hard spheres and borrowed the construction of weighted densities from DFT for hard spheres. He obtained the expression for a nonlo- cal free energy functional that depends on the segment density only. This approximation was also used to investigate a mixture of polymers and hard spheres in slit-like pores. The approach by Woodward was extended to polydisperse polymeric systems [33] with the Schulz-Flory mass distribution, and to in- finitely long polymers [34]. At present it is clear that that the weighing procedure of this method is crude and requires improvements. However, a weakness of the theory is due to the presence of an adjustable parameter to account for complex correlation effects. An important group of DF approaches has been constructed by using the thermodynamic pertur- bation theory (TPT) by Wertheim [35–37]. The earliest theory of this type was proposed by Kierlik and Rosinberg [38–41]. This approach applies the intramolecular distribution function for pairs of bonded atoms to evaluate the functional of the free energy of the system and a more sophisticated procedure for obtaining the weighted densities. The theory by Kierlik and Rosinberg was refined by Yu and Wu [42], who introduced a new concept for the treatment of intramolecular degrees of freedom. The novelty is in the use of the fundamental measure theory to accurately evaluate the free energy of a hard sphere reference system and pair cavity 12601-6 Fluid-grafted chains interfaces distribution function at contact necessary in the TPT formulation. The approach developed by Yu and Wuwas used in several theoretical studies of inhomogeneous polymeric systems including surface phase transitions [43], phase equilibria in confined geometry [44–46] and effective interactions [47]. The Yu and Wu functional permitted to derive an analytical expression for the polymer-wall interfacial tension [48], to study star polymers [49–51], rod-polymer mixtures [52], and hyperbranched polymers [53]. The ap- proach by Yu and Wu has been widely used to study the systems involving tethered chains [54–66]. Making a short overview of the DF approaches used in the description of nonuniform fluids involving chain molecules, one should mention a quite accurate TPT-based DFT approach by Chapman et al. [67– 69]. It concerns the associating fluids that can form different complexes and chains, in particular. The DF formulation for associating monomers is supplemented by appropriate mass action laws. The limit of complete association is imposed at the final stage of the procedure, in contrast to the method by Yu and Wu. On the other hand, associative approach is useful to describe the bonding of chain molecules with a solid surface. It should be noted that there exists a link between DFmethods and the SCFT for systems of polyatomic molecules. Several DFTs can be recast to the SCFT form [23, 24] with a different degree of difficulty. Generally, the so-called DFT-SCFT approaches rely on a combination of computer simulations of a single chain in an effective position dependent field that results from nonlocal DFT [70–75]. The advantage of the DFT-SCFT approaches is that practically all global and local properties of polymers at interfaces, such as the center-of-mass profile or the end-to-end vector profile, can be calculated and analyzed [76–79]. In particular, those methods can be used to take into account the shrinking or the swelling polymers at the surface. So far, such calculations have been carried out only for bulk systems using PRISM integral equation approaches [80–82]. Quite recently, Bryk and MacDowell proposed a new hybrid self-consistent field DF approach to de- scribe the solvation effects for polymers at interfaces [83]. Using DFT of Kierlik and Rosinberg [38–41], they showed that when the second-order terms of Wertheim thermodynamic perturbation approach (TPT2) [35–37] are taken into account, the DFT-TPT2 is capable of capturing the effects of solvation of tethered chains. Theoretical predictions were tested with respect to Monte Carlo simulations for moder- ate chain lengths [83]. The results for the end-to-end distances of chain molecules in the solvent in the bulk phase and at the solid surface were in a reasonable agreement with simulations. 3. Density functional theory Aswe have alreadymentioned in the introductory section, different formulations of DFT are available in the theory of inhomogeneous complex and polymeric fluids. In this text, however, we restrict ourselves to solely one, frequently used version, for illustrative purposes. It contains elements of other approaches and in the course of presentation we will address the issues related to other DFTs. Let us consider the fluid (indexed as 1) in contact with the surface covered with the end-grafted chains, P , [54–63] (an extension to fluidmixtures is straightforward). The chain consists ofM tangentially bonded segments. Although other models involving different segments can be considered, we restrict our attention to the case where all the segments are identical. The connectivity of segments within a chain is provided by the bonding potential, VB, defined as follows: exp[¡¯VB(R)℄Æ M¡1 Y iÆ1 ± ¡ jr iÅ1 ¡r i j¡¾ (P ) ¢ . 4¼ ¡ ¾ (P ) ¢ 2 , (17) where R ´ (r 1 ,r 2 , ¢ ¢ ¢ ,r M ) is the vector specifying the positions of segments and ¾ (P ) is the diameter of the segments. It is assumed that all the segments and fluid molecules interact with each other and with solvent molecules by means of (12–6) Lennard-Jones potential, u (kl ) (r )Æ 4" (kl ) · ³ ¾ (kl ) /r ´ 12 ¡ ³ ¾ (kl ) /r ´ 6 ¸ , (18) where "(kl ) is the strength of interactions between the species k and l and ¾ (kl ) Æ ¡ ¾ (k) Å¾ (l ) ¢ /2, (k, l Æ P,1). 12601-7 S. Sokołowski, J. Ilnytskyi, O. Pizio The chains are pinned to the surface. The interaction of their first segments with the surface is given by exp h ¡¯v (P ) 1 (z) i ÆC± ¡ z¡¾ (P ) /2 ¢ , (19) where z is the distance from the surface and C is a constant. The remaining segments of the chains, j Æ 2,3, . . . ,M , as well as the fluid molecules interact with the surface via a hard-wall potential, v (k) j (z)Æ ½ 1, z Ç¾ (k) /2, 0, otherwise, (20) or, depending on the model used, via the Lennard-Jones (9–3) potential, v (k) j (z)Æ " (k) s · ³ z (k) 0 /z ´ 9 ¡ ³ z (k) 0 /z ´ 3 ¸ . (21) In the above equation, " (k) s is the strength of interaction and z 0 Æ ¾ (k) /2. The geometry of the system can be changed. One can study a slit-like pore, then the potentials (19), (20) and (21) must be modified, cf. [57, 58]. In order to present the theory, let us introduce the notation ½(P ) (R) and ½(1) (r) for the density profiles of chains and fluid, respectively, and define the segment densities ½ (P ) s j and the total segment density of chains, ½(P ) s , via the relation, ½ (P ) s (r)Æ M X jÆ1 ½ (P ) s j (r)Æ M X jÆ1 Z dR±(r¡r j )½ (P ) (R) . (22) To simplify the equations, we use the notation ½(1) s (r)´ ½ (1) s1 (r)´ ½ (1) (r). The system is considered in the grand canonical ensemble. The thermodynamic potential, , is as follows: Æ F £ ½ (P ) (R),½ (1) (r) ¤ Å Z dR½ (1) ( r ) h v (1) 1 (r)¡¹ i Å M X jÆ1 Z ½ (P ) s j (r)v (P ) j (z), (23) where F £ ½ (P ) (R),½ (1) (r) ¤ is the Helmholtz free energy functional and ¹ is the chemical potential of species 1. The functional is evaluated under the constraint 1 A Z dr½ (P ) s j (r)Æ ½C , (24) where A is the surface area and ½C is the surface density of grafted chains. The free energy functional is divided into ideal, Fid, and the excess, Fex, terms. The excess terms are due to the hard-sphere interactions (hs), due to the connectivity of segments of chain molecules (c) and due to attractive interactions between all the species (att), F Æ Fid ÅFex, Fex Æ Fhs ÅFc ÅFatt. The configurational ideal part of the free energy functional is as follows: !!! ¯Fid Æ ¯ Z dR½ (P ) (R)VB(R)Å Z dR½ (P ) (R) h ln ¡ ½ (P ) (R) ¢ ¡1 i Å Z dr½ (1) (r) h ln ¡ ½ (1) (r) ¢ ¡1 i . (25) The excess free energy, due to hard-sphere interactions, Fhs Æ R ©hs(r)dr, results from the “white bear” version of the fundamental measure theory. Then, we have [42] ©hs Æ¡n0 ln(1¡n 3 )Å n 1 n 2 ¡n 1 ¢n 2 1¡n 3 Ån 3 2 (1¡» 2 ) 3 n 3 Å (1¡n 3 ) 2 ln(1¡n 3 ) 36¼n 2 3 (1¡n 3 ) 2 . (26) In the above, »(r)Æ jn 2 (r)j/n 2 (r). The definitions of scalar, n ® , for ®Æ 0,1,2,3 and vector, n ® , for ®Æ 1,2 weighted densities are given in reference [42]. 12601-8 Fluid-grafted chains interfaces The contribution Fc Æ R ©c(r)dr is evaluated from Wertheim’s first-order TPT [35–37] © (c) Æ 1¡M M n (P ) 0 ³ (P ) ln h y (P ) (hs) ¡ ¾ (P ) ¢ i , (27) where ³ (P ) Æ 1¡n (P ) 2 ¢n (P ) 2 ±¡ n (P ) 2 ¢ 2 and the contact value of the cavity function of hard spheres, y (P ) (hs) , results from Boublik-Mansoori-Carnahan-Starling-Leland equation of state, y (P ) hs (¾ (P ) )Æ 1 1¡n 3 Å n 2 ¾ (P ) ³ 4(1¡n 3 ) 2 Å ¡ n 2 ¾ (P ) ¢ 2 ³ 72(1¡n 3 ) 3 , (28) with ³Æ 1¡n 2 ¢n 2 /(n 2 ) 2. Finally, the attractive contribution is given by the mean-field approximation Fatt Æ 1 2 X k ,lÆ1,P Z drdr 0 ½ (k) s (r)½ (l ) s (r 0 )u (kl ) att (jr¡r 0 j), (29) where u(kl ) att (jr¡r 0 j) is defined according to Weeks-Chandler-Andersen scheme, u (kl ) att (r )Æ ½ ¡" (kl ) , r É 2 1/6 ¾ (kl ) , u (kl ) (r ), r È 2 1/6 ¾ (kl ) , (30) and the diameter is assumed d (i ) Æ¾ (i ) , (i ÆP,1). The equilibrium density profiles are obtained by minimizing under the constraint (24), ± ±½ (P ) R) Æ 0, ± ±½ (1) (r) Æ 0. (31) The final equations for the segment density profile and for the fluid density profile can be found in refer- ences [42, 57, 58]. The bulk fluid in equilibrium with the system of interest is the fluid of Lennard-Jones particles of type 1. The configurational chemical potential ¹ in terms of the bulk density ½b of species 1 is ¯¹Æ ln½bÅ¯¹(hs)Å¯½b Z u (11) att (r )dr, (32) where ¹ (hs) is the excess chemical potential of hard spheres. The knowledge of the density profiles is crucial for the evaluation of thermodynamic properties. The density profiles provide the excess adsorption isotherms of fluid species, ¡ex, ¡ex Æ Z dz £ ½ (1) (z)¡½b ¤ , (33) where the limits of integration are determined by the geometry of the system. The knowledge of the density profiles yields the value of and thus permits for the evaluation of the solvation force ¯ f /(2A)Æ¡ 1 A �¯ �H Æ Z dz½ (1) (z) �v (1) 1 (z) �z Å M X jÆ1 Z dz½ (P ) s j (z) �v (P ) j (z) �z . (34) The solvation force between the surfaces is a measurable property. It comes out from the experiments performed by using a surface force apparatus, see e.g., [84]. Our final comment concerns another important property. Upon fluid adsorption, the average height of tethered chains changes in a nontrivial way. In order to describe this property, we use the normalized first moment, h Æ a Z dzz½ (P ) s (z) . Z dz½ (P ) s (z). (35) Depending on the physical argument, the coefficient a is taken equal to 8/3 (cf. [85, 86]) or a Æ 2 [87]. The theory presented above permits for several extensions. In particular, it is possible to consider either tethered copolymers [88] or chain layers of different species. Flexibility of bonding between seg- ments can be modified as well, assuming, for example, preferential angles for bonding between certain groups of segments [89]. It is also possible to consider the branched polymers and to change the identity of segments (and their sequence). 12601-9 S. Sokołowski, J. Ilnytskyi, O. Pizio 4. Applications of the density functional approaches In this section, we would like to focus on a few results of DFTs. In general, one needs to be convinced that the obtained density profiles and the values of are accurate for different values of parameters describing the system. The profiles yield the excess adsorption and the values of allow for the investi- gation of the phase behavior. We proceed to a discussion of density profiles. The configuration of end-grafted polymers in a good solvent has been well documented. Specifically, the SCF theory by Milner et al. [90] based on the ideas of Alexander [3] and de Gennes [1], lead to conclusions that at a moderate surface grafting density, the chains exhibit parabolic density profiles. This result has been supported by simulation data. We focus on exploring the possibilities of DFTs and the DFT by Chapman et al. [69], denominated as iSAFT (interfacial statistical associating fluid theory). As we have alreadymentioned, this approach slightly differs from the theory described above by the treatment of bonding of chain molecules to the solid surface. The first illustration [69] shows the structure of athermal chains with M Æ 100 hard-sphere segments tethered to a hard wall. The density profile is parabolic (left hand panel of figure 1). At a high grafting density, the deviations from a parabolic shape are not big. In general, this DFT overestimates the density profile values close to the wall and underestimates the profile further from the wall. However, the dif- ference between simulation data [91] and theory is not large. In addition, Chapman et al. [69] showed that the DFT yields predictions similar to the SCF approach [92]. The right hand panel of figure 1 shows Figure 1. Left hand panel: segment density profiles of chains ofM Æ 100 segments tethered to a hardwall at a grafting density ½C¾ 2=0.03 (squares), 0.07 (triangles), and 0.1 (circles). Right hand panel: segment density profiles for a grafting density ½C¾ 2 Æ 0.1 in the presence of implicit solvent (circles) or explicit polymer solvent with M f Æ 2 (triangles up), M f Æ 5 (squares), or M f Æ 10 (triangles down). Symbols are the simulation results [91], and the curves result from iSAFT. The total segment density of fluid is ½b¾ 3 Æ 0.682. Reprinted with permission from reference [69]. Copyright AIP. the density profiles of chains tethered to a planar surface in the presence of free polymer solution of short chains resulting from iSAFT theory and compares them with simulation results and with the pro- files for an implicit solvent. The density profiles for the models with explicit solvent significantly differ from the implicit solvent case. The segments of solvent chains interact via a purely repulsive potential. Consequently, the observed effect is entropic. In an explicit solvent, the solvent molecules compress the chains and the layer height is reduced compared to an implicit solvent Nc kind of linear homopolymers case. The total solvent density is high, see the figure caption. The compression increases with an increase of the number of segments of the solvent, showing that penetration of longer (M Æ 10) free chains into the tethered layer decreases, compared to shorter chains. Thus, the chain solvent effect is to reduce the tethered layer height and to make it denser in the vicinity of the solid surface. The observations from the DFT are in agreement with computer simulations [91]. 12601-10 Fluid-grafted chains interfaces 4.1. Phase transitions After evaluating the density profiles, one can try to describe the thermodynamics of adsorption and to analyze the phase behavior. Such problems were studied by Borówko et al. [62], who considered a monoatomic fluid in contact with modified and non-modified solid surfaces. The fluid-wall interaction was given by the Lennard-Jones (9,3) potential. The surface was modified by short (5 segments) chain molecules. The model assumed the same interactions " Æ " (kl ) between all species, the diameters of seg- ments and of fluid particles were the same and "s/" Æ 8; "s Æ " (P ) s Æ " (1) s . The DF approach described in section 3 was used. 0.01 0.015 ρ b * 0 2 4 Γ∗ 0.9 0.92 0.95 0 2 4 6 8 10 z * 0 1 2 3 ρ(1 )∗ (z )0.008655 0.08655 0.00915 1 2 3 4 z * 0 0.05 0.1 ρ(P )* s (z ) 0.02 0.04 ρ∗ b 0.8 1 T * ρ C * = 0.02 T * = a ρ∗ b = b T * =0.87 T * =0.87 ρ∗ C =0.02 ρ∗ C =0.01 ρ∗ C =0c d Figure 2. (Color online) (a): Adsorption isotherms on a modified surface with the grafting density ½¤ C Æ 0.02 at different temperatures. Dotted vertical lines denote the bulk liquid-vapor transition. Dashed ver- tical lines denote the equilibrium thin-thick film transition. (b) and (c): Fluid (b) and segment (c) local densities at T¤ Æ 0.87 and for bulk fluid densities given in the (b). (d): The surface phase diagram. Solid line is the bulk phase diagram, dashed lines are the wetting diagrams for ½¤ C Æ 0.01, 0.2 and for non- modified surface, ½¤ C Æ 0. A set of adsorption isotherms, at different reduced temperatures, T ¤ Æ kT /", is shown in figure 2 (a). At high temperatures, the adsorption isotherms, ¡¤ Æ ¡¾2, monotonously increase with bulk density, ½ ¤ b Æ ½b¾ 3. However, at lower temperatures, they exhibit hysteresis loops. The first-order phase transition point at each temperature was found from the crossing of two branches of the thermodynamic potential, , as functions of chemical potential of the fluid. Dashed lines in figure 2 (a) show the bulk densities at the transition points. The changes in themicroscopic structure upon transition can be seen on the density profiles, see figure 2 (b) for T ¤ Æ 0.87. One observes a transition from a thin to a thick fluid film that is accompanied by changes of the tethered layer structure [figure 2 (c)]. It is evident that the jump of the fluid film thickness occurs within the tethered polymer layer. The changes of the structure of tethered chains are shown in [figure 2 (c)]. The growth of the fluid film causes the swelling of tethered chains. In order to characterize these phenomena it is necessary to construct the phase diagram. Numeri- cal calculations revealed that for the system without tethered chains, the wetting temperature is T ¤w Æ kTw/"Æ 1.03, whereas the surface critical temperature is T ¤sc Æ kTsc/"Æ 1.125. Upon introducing even a small amount of chains, one observes a shift of both characteristic temperatures to lower values. At the grafting density ½¤C Æ ½c¾ 2 Æ 0.01, the values of T ¤w and T ¤sc are T ¤ w Æ 0.83, and T ¤sc Æ 0.92, respectively, see figure 2 (d). Similarly to the changes of the wetting properties, one can investigate layering transitions. The lay- ering transitions are usually discussed for fluids in contact with strongly adsorbing surfaces. The aim is 12601-11 S. Sokołowski, J. Ilnytskyi, O. Pizio to explore how tethered chains change the scenarios for layering transitions. In particular, whether it is possible to suppress the layering of a fluid or to eliminate some transitions. However, for some systems, one can also observe new layering transitions as a consequence of chemical modification of a weakly attractive surface by strongly attractive chains. The important factors effecting the changes of the phase behavior are the differences in the interaction energies between the bare surface and the segments of chains and the entropic effects due to the changes in the structure of chains. The parameters such as the number of segments in the chains and grafting density are of importance as well. These factors permit to tune the phase behavior in terms of the formation or destruction of consecutive layers. In the case of adsorption in pores, the modification of the pore walls effects the capillary condensation. In this case, the ratio between the chain length and the pore width becomes a new important parameter. The chains can also enhance or inhibit layering transitions that may preceed the capillary condensation. All these problems were discussed in references [58, 59]. If the binding of the terminal segments of chains to the walls of a slit-like pore is not strong, a new kind of phase transition can take place. This transition is connected with the change of the symmetry of the distribution of chains and fluid molecules inside the pore. Upon the adsorption of a fluid, a part of weakly tethered chains can be “wiped out” from one wall and bonded to the opposite one. Consequently, the distribution of the chains inside the pore becomes non-symmetric. As a result, the symmetry of the distribution of fluid molecules with respect to the pore center is also broken. A further increase of the fluid density may induce a reentrant behavior manifested in the recovery of the symmetry. Such changes of symmetry may occur as first- or second-order transitions [64]. The symmetry breaking transitions are also observed for slit-like pores modified with chains whose two ends are permanently attached to the opposing walls, provided that the length of chains is suffi- ciently larger than the pore width [65, 66]. Such chains are often called pillars. In this case, the inner segments of the chain molecules can move towards one of the pore walls. This problem was discussed in references [59, 62]. The phase transitions, associated with the changes of the symmetry of local densities can be observed if the chains are much longer than the pore width. For chains of the length comparable to the pore width, such transitions do not occur [65]. Moreover, the symmetry changes are possible if the surface density of the chains is moderate. In the limiting case of a very small grafting density, the phase transformations in modified and non-modified pores are similar. In another limiting case of a very high grafting density, all phase transitions within the confined fluid disappear. The values of the fluid-wall energy parame- ters are also important. In the case of strongly attractive adsorbing potential, the symmetry breaking is linked with the layering (that takes place at one wall). For weaker adsorbing potentials, the growth of a non-symmetric film at one wall resembles the growth of a film after the prewetting jump. Finally, for a very weak adsorbing potential, the symmetry breaking does not occur. In the case of a strongly ad- sorbing surface, the first-order symmetry change transitions end up at the tricritical points from which the second-order transition lines depart. There exists a triple point between the “return-to-symmetry” transition and the capillary condensation envelope. In the case of weaker adsorbing potentials, however, the second-order symmetry change line joins the capillary condensation part of the phase diagram at the critical end point. In such cases, a symmetry recovery at lower temperatures takes place during capillary condensation. The phase behavior becomes even more complicated if the adsorption of a binary mixture inside the modified pore is considered. If the mixture exhibits a demixing in the bulk phase, then demixing tran- sitions inside the pore, or even within particular layers formed upon the layering transitions, compete with the symmetry breaking. This competition may lead to complex phase diagrams that include critical points, lower and upper critical end points, tricritical points, ¸-lines connected with demixing as well as ¸-lines connected with symmetry breaking and symmetry recovery transitions. Some of these problems were discussed in reference [66]. 4.2. Solvation force There are several works focused on the study of the effective interaction between surfaces coated with chains. In particular, Wu et al. [54] studied the solvation force between two surfaces covered with telechelic chains. Telechelic chains are chain molecules with certain functional groups at both ends. 12601-12 Fluid-grafted chains interfaces Thus, the end segments of a polymer can bind to a single surface simultaneously to form a loop, or to the opposite surfaces to form a bridge-like structure. The backbone segments in between the two termi- nal groups are chemically inert. Experimental and theoretical studies show that the interaction between highly stretched telechelic chains is mostly repulsive, in close similarity to the behavior observed for the system of two layers of singly-tethered chains. If the two layers of telechelic chains are put almost at contact, i.e., if their separation distance is close to the double height of an isolated chain, a weak attrac- tive force is observed. The origin and the magnitude of this attraction has been discussed by different authors [93–95]. Similarly to the DFT by Chapman et al. [67–69], the approach by Wu et al. [54] involves Figure 3. Solvation force between two surfaces mediated by telechelic polymers: (a) chain lengthM Æ 11, and (b) M Æ 32, (M is the number of segments). The inset plots are magnified views near the classical contact, i.e., when the separation between the two surfaces is about twice the brush thickness. Reprinted with permission from reference [54]. Copyright ACS. a parameter describing the strength of an attractive interaction between terminal segments and the sur- face, surface-adhesive energy, "Aw. In figure 3, taken from [54], the solvation force between two identical surfaces covered by telechelic chains is shown. The dependence of the solvation force on the value of "Aw and on the number of segments,M , was studied. The segment packing fraction in the bulk phase of chains was rather low, ½¤ b Æ 0.05. If the surface-adhesive energy increases from "Aw/kT Æ 1, to 8, 12 and 15, the solvation force exhibits changes from a short-range attraction to a long-range repulsion. At the lowest en- ergy "Aw/kT , the solvation force is similar to that for nonadsorbing polymers, "Aw/kT Æ 0; the attractive short-range part of the force is due to the entropic depletion of chains from the surface. At the highest value of "Aw/kT , the repulsive force arises from the interaction between tethered telechelic chains. The strength of repulsion depends on the energy "Aw/kT . The magnified view of the solvation force shown in the inset indicates that in a nonadsorbing regime (small "Aw/kT ), the theory leads to entropic depletion, as well as to a weak, secondary repulsion between surfaces. If the polymers are strongly adsorbed, the theory predicts the existence of a weak attraction between telechelic chains at contact. Moreover, the location of the attraction shows that the chains are highly stretched. 12601-13 S. Sokołowski, J. Ilnytskyi, O. Pizio DF observations are in agreement with experimental findings and other theories. If the chain length increases, similar conclusions concerning the behavior of the solvation force preserve. A more detailed discussion can be found in reference [54]. To summarize our illustrations in this part, it is worth men- tioning that similar successful developments have been used for the description of inhomogenous fluid mixtures, tethered chains of copolymer type and other related models with nonelectrostatic interparticle interactions. We proceed now to another class of systems, namely to the models with long-range Coulomb interaction. 5. Electrostatic density functionals An ample set of problems of interest for laboratory and for applications involves fluids of charged particles, as well as tethered polyelectrolytes and polyampholytes. In the case of tethered polyelectrolyte chains, the presence of long-range electrostatic interactions makes the properties of such systems much more complicated than those corresponding to neutral grafted chains. Since the pioneering works by Pin- cus [96], Balastre et al. [97] and Borisov et al. [98], numerous theoretical and experimental studies have been reported concerning the structure and surface properties of charged grafted chains. Similarly to the case of uncharged chains, the theoretical investigations were primarily based on the scaling analyses and the self-consistent mean-field theories that all rely on a number of simplifications [96, 98–102]. If the number and positions of charges at each chain are fixed, theoretical and experimental investigations elucidate the existence of different brush regimes that depend on the grafting density, degree of disso- ciation, and ionic strength. However, more recent experiments by Wu et al. [103] and Currie et al. [104] have shown that the predicted scaling relations do not necessarily hold throughout the entire range of the grafting densities and ionic strengths. Recently, Witte et al. [105] developed a new version of a SCF approach from the grand canonical partition function. The theory properly accounts for the local nature of the charge equilibrium and captures the basic behaviors of grafted polyelectrolytes. In agreement with recent experiments, the theory predicts that the scaling of grafted chains height with grafting density can be qualitatively different at intermediate chain lengths from that predicted by scaling arguments. This difference was attributed to the relative strength of electrostatic interactions compared to finite segment size packing constraints. On the other hand, the trend of decreasing chains height with an increasing grafting density predicted by the scaling analysis is recovered for large molecular weight polymers im- mersed in a solution of a very weak ionic strength. One should also mention here the SCF approach [106] aimed at predicting the structural properties of polyelectrolytes grafted via one of their ends to solid surfaces which explicitly incorporates the acid-base equilibrium responsible for the charge regulation of the polyelectrolyte groups, as well as the conforma- tions, size, shape, and charge distributions of all the molecular species present. The rest of that theory was compared with experimental observations for the height of the layer as a function of ionic strength for different chain molecular weight and grafting densities, and a good agreement was found. Similarly to the case of uncharged chains, DFT provides an alternative method to traditional scaling andmean-field theories by taking into account both short- and long-range correlations in a self-consistent manner [107, 108]. In order to illustrate the current developments of the DFT for this type of systems, we would like to briefly consider the adsorption of an electrolyte solution on the solid surface modified by grafted chain molecules containing charged segments. Commonly, an electrolyte solution is modelled in the framework of the primitivemodel, i.e., the fluid is considered as amixture of hard-spheres interacting via Coulomb potential. The same form of the potential is assumed for the interaction between the fluid species and charged segments of the tethered chains, u (´®) i j (r )Æ 8 > > < > > : 1, r Ç¾, e 2 Z (´) i Z (®) j 4¼²² 0 1 r , r È¾, (36) where the superscripts ´,®ÆP,1,2, and Z® i is the valency of ions or of the charged segments. The solvent is treated as a continuum with a given relative dielectric permittivity, ² (² 0 is the permittivity of the vacuum). 12601-14 Fluid-grafted chains interfaces The interaction between the ions and the wall, as well as between the charged segments of the chains and the wall, contains the electrostatic term, ¯v (®) j ,el (z)Æ¡2¼lBQZ (®) j z , (37) in addition to non-electrostatic interaction. Here, Qe is the surface charge density on the wall and lB Æ e 2 /(4¼kT ²² 0 ) is the Bjerrum length. The bulk densities of ionic species of the solution must satisfy the electroneutrality condition Z (1) 1 ½ (1) b Å Z (2) 1 ½ (2) b Æ 0. The electroneutrality holds for the entire system and the charge of tethered chains is neutralized by charges of counterions. The distribution of all charges in the system satisfies the Poisson equation, r 2 ª(z)Æ¡ 1 ²² 0 q(z), (38) whereª(z) is the electrostatic potential. The solution of this equation requires boundary conditions. For a fluid in contact with a single (modified) wall at z Æ 0, we can specify the value of the electrostatic poten- tial at the wall, setting V 0 Æª(z Æ 0). The wall charge follows then from the electroneutrality condition, QeÅ Z dz q(z)Æ 0, (39) where q(z)/e Æ P M jÆ1 Z (P ) j ½ (P ) j (z)Å P ´Æ1,2 Z (´) 1 ½ (´) (z). One can work either under the condition of a con- stantQ or V 0 . Similarly to the previous section, the principal issue is in the determination of the density profiles of fluid species and of tethered chains segments by minimization of the grand thermodynamic potential. In the present problem, the excess grand potential is given by Æ F Å Z dr M X jÆ1 v (C) (z j )½ (C) sj (z)Å X ´Æ1,2 Z dr h v (´) 1 (z)½ (´) (z)¡¹ ´ i Å Z drq(z)ª(z). (40) In contrast to section 3 the free energy F also contains the term due to the coupling between electrostatic and hard sphere interactions [109, 110]. It is given by Fe/kT Æ¡ 1 2 X ®,´,i , j Z drdr 0 ¢ ¯ (®´) i j (r,r 0 )¢½ ®,i (r)¢½ ´, j (r 0 ). (41) The summation in equation (41) is carried out over cations, anions and the segments (i , j Æ 1,2, . . . ,M) of the chains, ¢½ P,i (r)Æ ½ (P ) i (r), ¢½ 1,i (r)Æ ½ (1) (r)¡½ (1) b and ¢½ 2,i (r)Æ ½ (2) (r)¡½ (2) b . The functions ¢ ¯ (®´) i j are short-range parts of electrostatic two-particle direct correlation functions, that can be taken from the analytic solution of the mean spherical approximation for a restricted primitive model of electrolyte solutions and weighted correlation approach by Wang et al. [109, 110]. Another difference between the systems involving charged and uncharged segments is that the free energy term describing the connectivity of chains is modified. According toWertheim’s theory [35–37] for chains built of tangentially jointed hard spheres, the free energy contribution Fc (cf. section. 3) involves the contact value of the cavity function of the reference system of hard spheres, y (P ) hs . For uncharged systems, a key assumption of the DFT is that the contact value of the segment-segment cavity function is represented by that corresponding to a fluid composed of hard spheres, y (P ) hs . In the case of chains built of charged segments, an approximation for the cavity function includes the effects of screened segment- segment electrostatic interactions and takes the form [111, 112], y j , jÅ1 (¾)Æ y (P ) hs (¾)exp 2 4 ¡ 1 T ¤ Z (P ) j Z (P ) jÅ1 ¡ 2¡¾Å¡ 2 ¾ 2 ¢ (1Å¡¾) 2 3 5 . (42) In the above, ¡Æ ¡p 1Å2·¾¡1 ¢ /2¾, T ¤ Æ¾/lB is the “electrostatic” reduced temperature and · 2 Æ (4¼lB) " M X jÆ1 ³ Z (P ) j ´ 2 n (P ) 0 j (z)Å X ´Æ1,2 ³ Z (´) 1 ´ 2 n (´) 0 (z) # (43) 12601-15 S. Sokołowski, J. Ilnytskyi, O. Pizio is the square of the “local” screening parameter [109, 110]. The application of the above approximation for chains with different charge sequences under different thermodynamic conditions leads to qualitatively correct results [112]. 6. Selected applications of the electrostatic density functional approach In general, tethered polyelectrolytes respond to external stimuli such as the properties of electrolyte solutions, density, ionic strength, temperature, see, references [113–116], as well external electric field. There exist a vast literature dedicated to experimental, simulation and theoretical studies of such re- sponsive chains. Among various external stimuli, electric field is the one of much interest because it can induce changes of the tethered layer height [114, 117, 118] and consequently yields changes of, e.g., permeability, transport properties and electric conductance of nanochannels [119]. Upon putting an ionic fluid in contact with the modified surface, the height of the tethered poly- electrolyte changes depending on different factors. A set of results concerning these trends from the DF approach by Jiang et al. [55] is given in figure 4. In this particular case, the polyelectrolyte grafted at the Figure 4. (Color online) Left hand panel. Dependence of the tethered layer thickness on salt concentra- tion for chain length M Æ 20, 50, 100. The reduced grafting density is ½C¾ 2 Æ 0.01. Dashed lines are the predictions of the scaling laws. Right hand panel. Dependence of the tethered layer thickness on the graft- ing density in the osmotic regime (½b¾ 3 Æ 10 ¡4) and in the salted-brush regime (½b¾ 3 Æ 0.1). The chain length is M Æ 50, dashed lines are the predictions of the scaling laws. Reprinted with permission from reference [55]. Copyright ACS. uncharged solid surface consists of univalent negatively charged spherical segments, whereas the fluid (salt) is a mixture of cations and anions of the same valency. For simplicity, the diameters of ions and segments are assumed to be equal. At a given chain length, the height of polyelectrolyte layer practically does not depend on the salt density in the interval of low salt densities. This is the so-called osmotic-brush regime and this type of behavior was first established theoretically, see e.g., [96] and confirmed experi- mentally [97]. By contrast, at higher salt concentrations, the results shown in the figure are in accordance with the salted-brush regime, in which the brush thickness obeys the power law, established previously in references [96, 98]. In the unifying picture for a wide range of salt concentrations, one can say that poly- electrolyte chains undergo a smooth transition from “osmotic sticks” to salted coils. The DFT reproduces the previously established scaling laws for long chains. Moreover, for shorter chains, e.g., forM Æ 20, the DFT provides very good estimates for the tethered layer height on salt density in close similarity to other approaches. Another result concerns the dependence of the height of chains on the grafting density. Although the grafting density in the previous example was low, this value still yields the brush rather than mushroom- type structure, as it can be deduced from the values in the figure. The grafting density parameter effects the height of the chains in a nontrivial manner. Namely, in the osmotic regime, the height of the chains 12601-16 Fluid-grafted chains interfaces increases slightly with an increase of grafting density showing deviation from the scaling law behavior. Such a growth is in accordance with computer simulations and recent experimental findings [120, 121]. This “nonlinear osmotic” regime (i.e., when the concentration of counterions inside the tethered chains is lower than in the bulk fluid) is well described by DFT. Moreover, the slope of the curve changes with the grafting density, or, in other words, there is no simple scaling relation between two variables of the figure at low salt concentration. In the “salted-brush” regime, the height of the chains grows slowly until the separation between the grafting points of polyions becomes small, such that at a higher grafting density, the height of the chains growth obeys a quasi-power law. In fact, the emerging picture is richer and gives the details that are not intrinsically provided by simple scaling laws. Previous examples have concerned the polyelectrolytes tethered at uncharged surfaces. However, there is an increasing interest in developing nanofluidic systems that can control various properties, in- cluding transport of ions andmolecules in aqueous solutions in charged nanochannels. The next example shows the effects of non-zero electrostatic potential at the surface on the distribution of chain segments and of ions in the interfacial region. One would expect that the electric double layer will be substantially affected by the presence of tethered chains. Our presentation concerns the primitivemodel of electrolyte. Moreover, we assume that tethered chains consist of uncharged segments of the same diameter as fluid ions. Two walls covered with tethered chains constitute a slit-like pore. We apply the theory briefly de- scribed in the first part of this section. The walls are charged, and the charge is determined by the elec- trostatic potential on the wall V ¤ Æ eV 0 /kT . One of the properties of interest is the differential capacitance. It is defined as the derivative of the surface charge with respect to the applied electrostatic potential, V ¤, C ¤ d Æ dQ ¤ dV ¤ . (44) Recently there has been much interest in the dependence of a differential capacitance on the pore width [122–126]. This interest has been inspired by a possibility of developing new electric double layer ca- pacitors. In particular, it is quite challenging for the theory to provide interpretation of experimental observations concerning an anomalous growth of capacitance with a decreasing pore width for very nar- row pores. A few examples of the dependence of the reduced differential capacitance on pore width are given in figure 5. 1 2 3 4 5 6 H * 0.02 0.04 0.06 0.08 C * d a 1 2 3 4 5 6 H * 0.02 0.04 0.06 C * d ρ* C =0.025, V * =V PZC ρ∗ C =0.025, V * =3 ρ* C =0.025, V * =V PZC ρ* C =0.025, V * =3 ρ* C =0.05, V * =3 P3M3 P5M5 b Figure 5. (Color online) (a): The dependence of C¤ d on the slit-like pore width, the pore walls are mod- ified by uncharged chains (N6), V ¤ Æ 0.5. The parameters are: ½¤ b Æ 0.3, T¤ Æ 0.15. The curves are for grafting density ½¤ C Æ 0, =– black curve and circles, 0.05— red curve and squares, 0.1— blue curve and triangles. (b): The dependence of C¤ d on the pore width. The line decorated with circles is for ½¤ C Æ 0.05, all remaining lines — for ½¤ C Æ 0.025. The calculations are for the brushes P3M3 and P5M5 and the line nomenclature is explained in the figure. Reprinted with permission from reference [128]. Copyright AIP. 12601-17 S. Sokołowski, J. Ilnytskyi, O. Pizio In the absence of tethered chains, the dependence of C¤ d on the reduced pore width, H¤ Æ H/¾, ex- hibits oscillations, with two well pronounced maxima at H¤ ¼ 1.3 and 3.3, (V ¤ Æ 0.5) figure 5 (a). The maxima are separated by a trough. In the presence of tethered chains made of six uncharged segments (this system is abbreviated as N6), the magnitude of oscillations decreases, but at a low grafting density, ½ ¤ C Æ 0.05, the capacitance curve still preserves the shape similar to the case without chains. At higher grafting density, ½¤C Æ 0.1, the first maximum is washed out and only the reminiscence of the second maximum preserves. Even a relatively small amount of chains substantially alters the electric capacitance in narrow pores compared to a pore with bare walls. Since the segments are uncharged, the observed effect of chains is due to excluded volume effects. Figure 5 (b) shows the results for the pores modified by grafted polyam- pholytes, i.e., copolymers consisting of positively and negatively charged segments (in some models un- charged groups are considered as well, see e.g., [127]). The segments with charges of different sign can be distributed randomly or can form certain patterns. In particular, the charges of the same sign can be arranged in blocks. The results displayed here have been obtained for the chains whose net charge is zero, i.e., the number of positively and negatively charged segments within a chain is the same. Two polyampholyte chains were considered [128], namely built of 6 and of 10 segments, whose 3 or 5 initial segments bear positive unit charges, while the following segments bear negative unit charges. These sys- tems are abbreviated as P3M3 and P5M5, respectively. The majority of the results in figure 5 (b) are for a low grafting density, ½¤C Æ 0.025. Higher grafting density suppresses oscillations, similarly to the case of uncharged segments. The electrostatic potential at the walls is low because pronounced oscillations appear only at low voltages, as in the case of non-modified pore walls [58]. Dashed lines in figure 5 (b) were obtained at V ¤ Æ 3. Solid lines, however, were evaluated at V ¤ ÆVPZC, the potential of zero charge. Note that for the system involving polyampholytes, the value of the potential of zero charge differs from zero, even for chain particles whose net charge is zero [130]. Of course, the value of VPZC changes with the pore width. In the case of non-modified pores, the first (and the most significant) maximum appears at the distance close to H ¤ ¼ 1.3. The presence of a brush significantly diminishes this maximum. The maximum still preserves for the brush P3M3, but vanishes for the brush of longer chains, P5M5. In the case of a pore with non-modified walls [129], the the distribution of charges inside the pores, and, consequently, the dependence of Q¤ on V ¤, results from the possibility of the development of a number of layers of ions at both walls and from the interference between the layers formed at the op- posing walls. For pores having walls covered with brushes whose segments are charged, there exists a constant amount of charges inside the pore, due to charged segments. Their distribution is not “free”, but is restricted by the existing bonds between the segments. This “background” charge distribution reduces the freedom of changes of the distribution of ions and, consequently, variations of Q¤ with V ¤ become smaller. This is reflected by a smaller (in comparison to non-modified pores) amplitude of the function C ¤ d vs. H¤. Another interesting question is the effect of the brush on the dependence of the double layer capac- itance on the voltage (or the charge) at the wall. It is known that for non-modified surfaces, the shape of this curve changes from a single-hump (“bell-like”) to a double hump (“camel like”) if the bulk den- sity changes. For non-modified surfaces, the appearance of single, as well as double humps was found in Monte Carlo simulations by Fedorov et al. [131, 132], and in DF calculations by Wu et al. [133] and [134]. This behavior was experimentally observed in the systems that involve surfaces free from grafted species [135, 136]. The calculations carried out in reference [128] have revealed that for narrow pores modified with grafted polyampholytes, even three humps can appear on the dependence of the capacitance on the elec- trostatic potential on the wall and that the shape of the capacitance curve depends on the grafting den- sity. In general, in the presence of tethered species in slit-like pores, the shape of C¤ d (Q ¤ ) changes from a camel-like to a bell-like, if the bulk RPM fluid density changes from low to moderate and high. The value of crossover bulk density depends on the grafting density of the chains, their length and on the charge sequence along the chains. It is possible to alter the shape of C¤ d (Q ¤ ) curve from double to single hump by changing the grafting density, if the bulk fluid density is adequately chosen. The symmetry of C¤ d (Q ¤ ) curves with respect to Q¤ Æ 0, observed for several model brushes of uncharged segments, is destroyed if one explores the brushes made of the chains having charged segments. To summarize, these few examples show the capability of the DF approaches to describe the models 12601-18 Fluid-grafted chains interfaces of complex systems comprising a solid surface with end-tethered polyelectrolytes or polyampholytes and electrolyte solutions. There is much room for improving the modelling and theory. One of the principal problems is an adequate description of the aqueous solvent and its dielectric properties, as well as to reasonably capture the screening of ion-solvent interactions. The inclusion of orientation degrees of free- dom of species into the DFTs is, therefore, necessary. In general, we expect that the systems described above will find their place in the research aimed at optimizing the energy density of supercapacitors, see e.g., [137]. 7. Computer simulations Computer simulation is considered as an intermediate method bridging up the experimental and the- oretical studies. Theoretical DF approaches have already been discussed above. Concerning the experi- mental studies, one should note that recent technological progress in the fabrication of modified surfaces involving tethered chains is not adequately followed by the advances in themethods for their experimen- tal analysis [138]. The existing techniques such as optical, X-ray and neutron scattering, operate on the meso- andmicro-length scales and cannot be informative enough for the systems rapidly changing on the nano-scale (as in the case of nano-patterned surfaces). On the other hand, the Atomic Force Microscopy (AFM) can be invasive [139], and can lead to non-negligible changes of the material structure [138]. Other difficulties are connected with a precise control of polymer synthesis (problems with polydispersity, de- fects, purity, etc.). The benefit of computer simulations is that they provide access to the positions (and in some methods to the velocities) of individual particles enabling a high resolution insight into the structure and inter- nal dynamics of the particles constituting the material. As a result, the relevant microscopic mechanisms of particular phenomena can be clarified or even predicted in the way inaccessible for theory and the experiment. This level of resolution is achieved at the expense of simplification of the interaction poten- tials (often being empirical) and of a finite size of the simulated system. Most simulations of the systems involving tethered chains can be split into two categories [140–142]: the (semi-)atomistic (the particle represents an atom or a group of atoms involving hydrogens, e.g., CH 2 ) and mesoscopic (the particle represents a fragment of a polymer chain or a collection of solvent molecules). Both approaches can be thought of as mutually complementary. While the (semi-) atomistic modelling allows one to address the real chemistry of the system (and the effects like hydrogen bonding, entanglements etc.), it lacks a pos- sibility to cover large length- and time-scales of the phenomena under study. The mesoscopic modelling is capable of doing that, but at the expense of introducing a certain level of coarse-graining focused on chemical (in)compatibility of fragments and the effects of their microphase separation. For further de- tails of various particle-based simulation techniques for the systems involving tethered chains see the reviews [140–144]. 7.1. Mesoscopic modelling The mesoscopic level simulations have been mainly based on the Brownian [145–149] and dissipative particle dynamics (DPD) [150–153] methods and we will restrict our review to the latter approach only. DPD is an off-lattice mesoscopic simulation technique which involves a set of particles moving in continuous space on a discrete time scale. Particles (beads) represent whole molecules, fragments of the polymer chain orfluid regions, rather than single atoms, and atomistic details are not considered. In other words, the internal degrees of freedom of particles are integrated out and instead of real interparticle interactions one introduces effective pairwise forces that depend on bead-bead distances. The total force acting on i -th bead from the interaction with all other j -th particles is given as a sum f i Æ X i, j h F C i j ÅF D i j ÅF R i j i , (45) where F C i j is a conservative force that identifies polymer beads, while the dissipative, FD i j , and random, F R i j , forces form a thermostat that keeps the temperature of the system constant. All the details of the 12601-19 S. Sokołowski, J. Ilnytskyi, O. Pizio method (e.g., justification of introduction of effective forces and explanation how to derive them) are given elsewhere [150–152]. A key property of the forces in DPD simulation is that they conserve the mo- mentum locally, so that hydrodynamic modes of the fluid emerge even for small particle numbers. The parametrization of conservative forces is density and temperature dependent and can be linked to the Flory-Huggins parameter for the case of mixtures (see, e.g., reference [154]). The gain of such a way of coarse-graining is twofold. Firstly, the number of moving entities is essentially reduced. Secondly, the effective forces in DPD aremuch “softer” than the atom-atom ones. This allows one for using much larger timesteps in the DPD simulation, in comparison with (semi-)atomistic MD discussed below. The conserva- tion of total momentum ensures a correct hydrodynamic behavior of the systems (i.e., the system satisfies the Langevin hydrodynamic equations). As a result, the method permits to access much longer time and length scales and ensures a faster micro-phase separation of non-miscible species, as compared with con- ventional MD simulations. Simulations of polymeric fluids in volumes up to 100 nm in linear dimension for tens of microseconds are now common [151]. The DPD modelling of confined fluids, or systems with walls, requires a special care. As it was dis- cussed by Visser et al. [155], the central problem is how to construct a solid wall for mesoscopic mod- elling. At this level of modelling, three boundary conditions must hold for a solid wall: (i) impenetrability (particles are not allowed to enter the wall), (ii) the wall should not artificially affect (i.e., by imposing an artificial structure) the fluid properties, and, in the case of flows, (iii) no-slip; the wall should impose velocities in accordance with Langevin hydrodynamics. Previous studies succeed to implement impen- etrability and no-slip boundary conditions. They usually used several layers of frozen particles residing on a surface [156–159] to ensure impenetrability. This, however, induces the ordering near a surface, which propagates into a substantial portion of the system. This effect has a physical justification for the atom-based simulation when one mimics the crystal-like structure of the wall, but should be interpreted as artificial for the meso-scale type of modelling. Quite recently, a new approach was introduced for modelling the walls in DPD studies [160]. It was assumed that, if the i -th bead with coordinates r i Æ {x i , y i ,z i } approaches the surface, e.g., the plane z Æ 0, it interacts with an “imaginary” bead that has the coordinates r j Æ {x i , y i ,0}. The total wall-bead interaction has all three contributions given by equation (45). Repulsion of the surface can be tuned by adjusting the parameters of conservative force, F C. The condition of totalmomentum conservation for the wall-bead interaction is applied by the following procedure. The force acquired by each wall due to the interaction with the approaching beads is accumulated during the evaluation of forces. Then, it is evenly distributed over all the beads inside the simulation box. The effective surface, therefore, is structureless and, thus, it does not introduce any artificial perturbation of the fluid structure near the wall. The details concerning DPD simulations can be found in references [151, 156, 160–164]. Let us now briefly summarize the most interesting DPD results for the systems involving tethered polymers. One of the first DPD studies of polymer-coated surfaces was carried out by Gibson et al. [165], who investigated the adsorption of colloidal particles on a surface modified with pinned polymers. In particular, it was found that the adsorption of colloidal particles was smaller when the size of the polymer relative to the colloidal particle and the density of the polymers increased. Moreover, the adsorption reduces when the polymer is well solvated. Malfreyt and Tildesley [163] simulated the brushes tethered at two walls of a slit in a good solvent. The density profiles of polymer segments across the pore were found to exhibit parabolic shape, and diffusion along the pore axis was significantly greater for the solvent particles in the middle of the slab than in the polymer region, since the solvent molecules were trapped within the entangled polymer. Later, Goujon, Malfreyt, and Tildesley [156] extended DPD method to the grand canonical ensemble and studied the compression of grafted polymer brushes in a good solvent. They also studied the compression of polymer brushes under shear and evaluated the friction coefficient as a function of compression, shear rate and properties of the solvent [164]. The constant and oscillatory flow of a fluid between two plates with tethered chains was simulated by Wijmans and Smit [166]. Pastorino et al. [167] investigated a short-chain melt between two brush-covered surfaces. The equi- librium situation, as well as the one under shear flow were explored. The polymers of both brush and melt were assumed identical. They studied the interdigitation of the melt and the brush and also consid- ered the orientation of bond vectors on different length scales, as well as radii of gyration, end-to-end vectors of free and grafted chains, and in the case of flow simulations, the velocity profiles. In the second 12601-20 Fluid-grafted chains interfaces work, Pastorino et al. [168] used the Gibbs criterion to localize the brush-melt interface and analyzed its equilibrium fluctuations in terms of a capillary wave Hamiltonian augmented by an elastic term that accounts for the deformability of the brush. However, one should stress that the simulation method used by Pastorino et al. [167, 168] was different from the standard DPD scheme [152]. Namely, they assumed that the interactions between the polymer segments and polymer segments with surfaces were of the Lennard-Jones type, and they performed molecular dynamics simulations with DPD thermostat. The effects of chain architecture on the structure of tethered rod-coil polymer layers was studied by Li et al. [169]. When immersed in a selective solvent for the coil blocks, rod blocks tend to form aggregates. Linear and Y-shaped polymers exhibited a similar aggregative behavior, but comb-like brushes were found to possess more diverse aggregative behavior compared to linear brushes. Surface structures with aggregates taking form of cones, cylinders, or layers of spheres were found. The behavior of grafted bi- nary polymer brushes with compatible components in the case of different chain lengths was investigated by Xue et al. [170]. They observed layered structures parallel to the surface indicating “phase separation”: short chains were suppressed in the layer adjacent to the surface, whereas longer chains were to much extent stretched. By slightly changing the solvent selectivity in preference of shorter chains, inversion of the layered structure was found. There have also been performed several DPD studies of the brushes tethered to non-uniform surfaces, see e.g., [138, 160, 161, 171, 172]. The systems involving brushes attached to some parts of the surface, the so-called structured brushes, are potentially very important in the manipulation of fluids at very short length scales. A prototype case consists of a solvent in contact with a surface, or confined between two identical plane-parallel substrates, decorated with stripes of tethered chains that alternate periodically in one direction (say, X ) and are infinite in the other transverse direction (Y ). The performed simulations indicate that heterogeneity of tethered layers has a great impact on the structure of confined fluid and on thermodynamic and dynamic properties of the systems. In particular, Ilnytskyi et al. [160] employed DPD to investigate the behavior of a binary mixture of A and B, exhibiting a demixing in a bulk phase, confined in slit-like pores with walls modified by the stripes of tethered brush of chains made of beads A. Themain issue was to determine possible morphologies that can be formed inside the pore depending on the geometrical parameters characterizing the system (the size of the pore and thewidth of the stripes), cf. figure 6. In order to describe the observedmorphologies, these authors calculated several characteristics, such as the density and local temperature profiles, the radii of gyration for the attached polymers, and the minimum polymer-polymer distances in the direction parallel and perpendicular to the pore walls. The summary of findings was presented as a sketch of the diagram of morphologies. The study of reference [160] has been extended to the case of different (so-called in- and out-of-phase) arrangements of the stripes at two surfaces and varying composition of the confined mixture [161]. For in-phase narrow stripes and narrow pores, the case is reduced to a quasi two-dimensional micro-phase separation within the mixture confined between the polymer-rich lamellae on each wall. On the other hand, for narrow pores and moderate stripe widths, the in-pillar separation is quasi one-dimensional. In the case of wider pores and very wide stripes, the pore interior splits into the regions showing quasi two- dimensional and bulk-like behavior. A number of solvent-mediated morphologies are observed when the stripes are weakly separated, both along the wall and across the pore. With an increase of the concentra- tion of A beads in the mixture, they bridge up the striped brush in either or in both directions. Certain effects of morphology switching are observed. On the application side, the concentration driven mor- phologies can be fixed permanently via crosslinking, if a mixture contains a short crosslinker compatible with the polymer. Fixed morphology could be used as a nanotemplate or as a microfluid channel [161]. An interesting work on the application of DPD simulation technique to describe brushes on nanopat- terned surfaces was published by Patra and Linse [138]. In themodel considered by them, polymer chains were randomly attached onto an infinite long stripe along the Y axis with finite width ¢ along the OX axis. The polymers were treated as jointed chains composed ofM spherical subunits (beads) connected by bonds. Each bead corresponded to a small part of a real polymer, approximately one Kuhn length [173], and was described by a soft repulsive potential. The simulations were carried out for a range of widths of the stripe, ¢, for several grafting densities ½C and for several numbers of beads,M . In figure 7 we show examples of the local densities of the beads as functions of X (along the surface) and z (in the direction perpendicular to the surface). The calculations were carried out for several values of ¢ (all distances are measured in the unit of the bead diameter). 12601-21 S. Sokołowski, J. Ilnytskyi, O. Pizio Figure 6. (Color online) Solvent-mediated morphologies formed by the polymer brush and the solvent mixture with an increase of the good solvent fraction (from left to right). Two cases are shown of the in- and out-of-phase arrangement of stripes, top and bottom rows of images, respectively. The plot is den- sity color-coded showing the polymer-rich domains (in blue) and good solvent-rich domains (in green). The “empty” space is occupied by a “bad” solvent. One should note solvent-mediated switching between lamellar and “meander”-like morphology in the bottom row of images. Reprinted with permission from reference [161]. Copyright CMP. As ¢ is increasing, the following trends in the system were observed: (i) an increase in the brush height h, (ii) an expansion of the region with high density close to the grafting surfaces, and (iii) an appearance of a considerable overshot of polymers outside the grafting region. The polymer density fell off rapidly when the perimeter of the brush was reached. Formally, the brush profile over a stripe can be evaluated by using AFM. However, the AFM tip, when approaching the surface, deforms the brush and, therefore, the height obtained from experiments differs from that for undeformed brush. Typical shapes of brush profiles from simulations and AFM measure- ments, as estimated by Patra and Linse [138], are displayed in figure 8. In contrast to the simulated profile, the experimental one does not display a clear plateau in the center of the stripe even though the brush is 10¥100 times as wide as it is high. Patra and Linse also present a strict derivation for the maximal decrease of the brush height outside the center. Since a polymer brush is a soft material, it bends under the force of the AFM tip. This behavior was also obtained by MD simulations of homogeneously grafted polymer brushes [174]. The simulations by Patra and Linse have been carried out in the domain M É ¢ É 5M . Their results suggest that this domain contains the transition from narrow polymer brushes to the regime of wide brushes. An improved understanding of this transition calls for a reassessment of the experimental con- ditions of the AFT technique, but even more exciting task would be to investigate the experimentally preferable orientation and conformational changes of chains near an edge by using, for example, labeled polymers. Many applications of nanopatterned polymer brushes would benefit from such a deeper un- derstanding of the molecular arrangement in narrow polymer brushes. Figure 7. The contours of the beads density for polymer lengthM Æ 100 and grafting density ½c Æ 0.11 at stripe width ¢Æ 15, 50, 200, and1 (from left to right). The contour curves mark the densities 0.01, 0.02, 0.05, 0.1, and 0.2. The location of stripes is indicated by the filled black areas. Note, that the scale of the abscissas differs among the panels. Reprinted with permission from reference [138]. Copyright ACS. 12601-22 Fluid-grafted chains interfaces Figure 8. Typical shapes of a nanopatterned brush as computed by simulation (solid line) and as mea- sured by AFM (broken line). Horizontal and vertical length scales have been indicated in the figure. The location of the stripe is indicated by thefilled black area. Reprinted with permission from reference [138]. Copyright ACS. 7.2. Atomistic and semi-atomistic modelling The strength of mesoscopic simulation methods, discussed in the previous subsection, is focused around their capability to cover nano-scale structure and the dynamics of polymer brush and the (mix- ture of) fluids contained within the pores. However, the soft interaction potentials, used in such studies, are typically of heuristic form and their relation to the real chemical substances is still an open ques- tion. Rigorous coarse-graining methods are still to be developed [175–178]. Therefore, the atom-based simulations have a twofold purpose: (i) to be used on smaller scale systems to serve as a feed for the esti- mates of coarse-grained interaction (so-called multiscale approach) and (ii) to be an independent tool for more “chemically detailed” simulation studies. An extensive literature is available on the applications of (semi-)atomistic models, reviewed earlier in references [143, 144, 179–181]. Here, we outline only some of the works that are of interest or have been omitted in previous reviews. One of the major tasks of computer simulations is to verify the analytic predictions from the SCF theory for grafted polymer layers. The MC simulation using the bond-fluctuation model was employed in [183] where the quantities describing the equilibrium structure of the brush are derived from the SCF theory and compared with the MC data without free parameter. In most cases, the results are in agreement with the SCF predictions. The causes of discrepancies are also discussed. The MD simulation of polymeric brushes in solvents of varying quality is presented in reference [91] where the resulting equilibrium structures are compared with the predictions of SCF approaches. The scaling of the brush height was found to agree with SCF predictions at temperatures T È T µ , where T µ is the µ temperature [2] of the solvent. For T Ç T µ , the agreement was poor, which was connected with the separation of the brush into solvent-rich and solvent-poor phases. This phase separation was found to disappear with an increasing grafting density and chain length. The effect of a weak attractive surface interaction on the structure of the brush was also studied. The entropy of polymer brushes in a good solvent limit was evaluated numerically in reference [184]. The cases of a free polymer brush, a compressed polymer brush and two polymer brushes facing each other were studied. The total number of configurations, entropy, and the force acting on the polymer were calculated for various chain lengths and grafting densities by means of MC simulations using an ef- ficient enrichment algorithm on a simple cubic lattice. It was demonstrated that the effect of an excluded volume plays an essential role and that the entropy of a free polymer brush drops as (spacing)¡2.4, while the entropy of a compressed polymer brush and two polymer brushes agrees very well with the analyt- ical SCF theory by Milner, Witten, and Cates [185]. The analysis has been also extended to the case of a polymer solvent. A MD simulation of polymeric brushes immersed in a melt of mobile polymer chains was presented [186]. The brush height and segment density profiles were evaluated for chains of length M grafted at one end to a solid surface, immersed in a melt of polymers of the same type. As the chain length P of free chains increases, there is a crossover from a wet to a dry brush in agreement with scaling and SCF theories. Since the interaction between all segments is identical, this crossover is driven purely by entropy. The simulation results were compared with earlier simulations. The exchange of chains between the brush and the bulk solution is of obvious relevance in under- standing the thermodynamics and the kinetics of the formation of a grafted layer. In reference [187], by including a finite head-group adsorption energy, the bond-fluctuation model has been used to simulate 12601-23 S. Sokołowski, J. Ilnytskyi, O. Pizio the polymer brush in a good solvent with the effect of chain exchange between the brush and the bulk. The self-adjusted surface coverage was measured for different values of chain lengths and head-group energies. It was found that the polymer chains in a grafted layer are replaced by introducing shorter chains of identical head groups, which supports some experimental studies. The irreversible adsorption of single chains grafted with one end to the surface was also studied by combining the scaling theory with computer simulations in reference [188]. Using scaling arguments, the authors derived relations between the overall time of adsorption the characteristic time of adsorption, and the chain length. To support the analysis, the MC simulations were performed using the bond fluctuation model. The sequence of adsorp- tion events is reproduced very well by simulations, and an analysis of various density profiles supports the theoretical model. To clarify the differences in the properties of grafted chains at various grafting densities, a num- ber of simulations were undertaken. In particular, the MC simulations were performed using bonded Lennard-Jones spheres. These simulations were aimed at observing the mushroom-to-brush transition at low grafting densities as well as the properties of a high density brush [189]. Through approximate scaling arguments, the location of the mushroom-to-brush transition was estimated and found to agree closely with the grafting density at which no more chains with mushroom dimensions can be inserted without overlap. Comparisons with analytical SCF theory for poor solvents were found to be much less favorable. In this case, the density profiles substantially differ from the SCF picture due to high segment densities and cluster formation. In reference [190], the static and dynamic properties of polymer brushes ofmoderate and high grafting densities were studied usingMD simulations and the Langevin-type model. The fluctuation dynamics in the direction lateral to the surface was found to be well described by Rouse scaling. The dependence of the dynamics on the grafting density was well described by the dynamic be- havior of thermal blobs for fluctuations in the direction both perpendicular and parallel to the surface. The pulling forces acting on the polymer bonds were also studied. The knowledge of the dependence of the bond forces on the grafting density is important to develop strategies for making brushes with high grafting densities. The MD method allowed the authors of ref- erence [190] to calculate the average vertical forces pulling the chains and the corresponding grafting energies. The binding potential between the nearest chains was the FENE potential, for details see refer- ence [190]. The obtained pulling force as a function of the average distance from the surface for chains of M Æ 64 segments is displayed in figure 9. The results were calculated at different grafting densities (shown in the plot). At the top of the brush, the stretching forces vanish, as long as the grafting density remains moderate. At high grafting densities, even the end bond displays some residual stretching as a result of strong interactions with neighboring chains. Moreover, at high and moderate grafting densities, one observes a plateau which indicates the brush regime. Figure 9. (Color online) Averaged bond forces along the z-direction for M Æ 64 as a function of the dis- tance from the surface at different grafting densities, given in the figure. Reprinted with permission from reference [190]. Copyright ACS. 12601-24 Fluid-grafted chains interfaces The wetting behavior of polymers on top of a brush consisting of identical polymers has been studied using the computational methodology developed for a direct measurement of a wetting transition and its order via the effective interface potential. The method also permits to estimate the contact angles in the non-wet state and to study the adsorption isotherms [182]. In the absence of long-range forces, the system shows a sequence of non-wet, wet, and non-wet states as the brush density is increased. Having included attractive long-range interactions, we can make the polymer wet the brush at all grafting densities, and thus both first- and second-order wetting transitions are observed. The study highlights a rich wetting behavior when competing adsorbent-substrate interactions of different scales are tuned over a broad range. Polyelectrolyte brushes have also been simulated using the MD method [191]. These were aimed at studying the fully charged polyelectrolyte brushes with salt, and the results are compared with SCF the- ory including finite stretching and volume effects [192]. The SCF approach was found to be capable of reproducing the brush heights at different grafting densities, salt concentrations and chain lengths at a semi-quantitative level. At high grafting densities, the density profiles obtained with both techniques ex- hibit a particularly close agreement, while at low densities, systematic deviations between their shapes are observed. The approximation of local electroneutrality, which the SCF approach is based on, was studied, and its implications were discussed. It is instructive to compare the local net charges of the vertical layers of the brush (figure 10). This quantity serves for validation of the local electroneutrality assumption made in the SCF approach. Not surprisingly, this assumption is satisfied rather well when the Bjerrum length is large (i.e., for strong electrostatics) and less within the weak charge regime. Throughout the osmotic regime (we recall that in the osmotic regime, the concentration of counterions inside the polyelectrolyte layer is higher than in the bulk solution [193, 194]), there exists a charge separation near the brush surface. The counterions form a charged layer just above the brush as a result of their osmotic pressure. It is this layer that generates the osmotic pull which contributes to the chain stretch. On the contrary, the local charge neutrality of the osmotic regime is rather well satisfied deep inside the brush. In the case of thick brush layers, when the surface contributes just a small fraction to the brush height, the approximation of local neutrality would turn even more accurate [191]. Figure 10. (Color online) Local net charge q(z) as a function of the vertical coordinate for different Bjer- rum lengths, given in the figure. Reprinted with permission from reference [191]. Copyright ACS. One of the topics extensively studied by using computer simulations is the adsorption of fluids on the surfaces modified by tethered chains [182, 184, 189] and the chromatographic process on columns filled with adsorbents modified by tethered oligomers [200–207]. A popular technique in the liquid chromatog- raphy is the Reverse Phase Liquid Chromatography (RPLC). The mobile phase in RPLC is a mixture of water and an organic solvent (usually methanol or acetonitrile). The stationary phase consists of a sup- port material (usually beads of silica of micrometer size) and a bonded phase. The bonded phase is most 12601-25 S. Sokołowski, J. Ilnytskyi, O. Pizio often built of C18 alkyl chains. The retention of solutes occurs primarily within or on top of the bonded alkyl chains of the stationary phase. The mechanism of retention in RPLC has been of interest for over three decades, but many of its as- pects are still highly debated [208, 209]. In general, solute retention in RPLC is driven by the free energy gradient on passing from the mobile phase into the stationary phase. One of the pioneering MD simula- tions of the transfer of a simple nonpolar solute from a water/methanol solvent mixture into a C18 sta- tionary phase at room temperature was carried out by Klatte and Beck [195]. These authors, in addition to a detailed examination of the local environment of the solute, computed the excess chemical potential profiles. The free energy change was consistent in magnitude with that expected for hydrophobic trans- fer from water rich to oil phases, but specific interfacial effects drew into question the bulk partitioning models. To elucidate the molecular-level structural features that control shape-selective separations, the MD was performed for chromatographic models with alkylsilane length and temperature analogous to ac- tual materials [196]. Alkyl chain order was found to increase both with an increase of the chain length and with a reduction of the temperature, as observed experimentally. Chromatography models of vari- ous chain lengths and over a temperature range that represents highly shape-selective RPLC stationary phases contained a series of well-defined rigid cavities. The size and depth of these “slots” increased for the C30 models, which may promote the enhanced separations of larger size shape-constrained solutes, such as carotenoids [196]. According to modern theories of RPLC, two main questions regarding the retention mechanism ap- pear: (i) is solute retention better described by partitioning into the stationary phase that resembles bulk liquid-liquid partitioning or by adsorption onto the stationary-phase surface? (ii) what is the driving force for the retention process, unfavorable interactions with the mobile phase or favorable interactions with the stationary phase? The answer of the above questions was considered in a series of papers by Siep- mann and co-workers [201–207], who employed an advanced MC simulation method. In particular, in an interesting recent work [206], the effects of stationary phase and the solute chain length, as well as the solute adsorption and separation were studied. The surface was modified by grafting dimethyl triacontyl (C30), dimethyl octadecyl (C18), dimethyl octyl (C8), and trimethyl (C1) to silane surface. Themodified sur- faces and the bare silane surface were in contact with a water/methanol (33% of methanol) mobile phase that contained normal alkanes and alcohols. The simulations were carried out using the Gibbs ensemble and the details are given in [206, 210, 211]. The computer simulations were carried out at T Æ 323 K and pressure of 10 atm, i.e., at the condi- tions typical of RPLC experiments. Figure 11 (panel a) gives typical snapshots of the configurations of the molecules in the selected systems. Simulating tetheredmolecules, Siempmann et al. [201–207] investigated how the structure of tethered layer changed with the length of the chains. In particular, they measured the dihedral angles, defined as the consecutive angles between carbons separated by two bonds and the normal to the silica substrate, as well as the so-called “gauche defects” [212] (defined as the fraction of dihedral angles in the entire chain deviating by more than 60± from the angle of the trans-conformer. They also studied the local composition and stationary phase thickness. For this purpose, density profiles for water, methanol, and stationary phase carbon atoms were evaluated, cf. figure 11 (panel b). The plots also show the Gibbs dividing surface (GDS, i.e., a plane defining the border between the mobile and stationary phases) and the width of the interfacial region. The latter was defined as the range of z-values where the total solvent density lies between 10% and 90% of its bulk value. However, due to very large oscillations in the total solvent density, it was not possible to determine a GDS and the interfacial region for the bare SiO 2 . One can see several similarities between the systems with the chains of different length. First, one observes a peak of the methanol local density just above the location of the GDS, i.e., there is a significant enrichment of the organic modifier in this region that extends deep into the interior of the tethered chains. The enrichment of methanol at the surface was similar for all lengths of chains. Second, there is a depletion of the total system density for all systems (see cyan lines), i.e., a dewetting effect at the surface. However, there are also some differences between the system with different chains. For the C1 sys- tem the interface is very sharp and the methanol enrichment is very pronounced. The solvent density shows an oscillating behavior while for other chains it decays smoothly in the interfacial region. For C8 12601-26 Fluid-grafted chains interfaces Figure 11. (Color online) (a): Examples of the simulation snapshots of the five systems obtained by Siep- mann et al. [206]. The stationary phase is shown as tubes with silicon in yellow, oxygen in orange, hy- drogen in white and carbon in black. The mobile phase is shown in the ball and stick representation with oxygen in red, hydrogen in white, and carbon in blue. The analyte molecules are depicted by large spheres with oxygen in red, hydrogen in white, and carbon in cyan. (b): System composition as a func- tion of distance from the silica surface. Stationary phase carbon density is shown in black, methanol in blue, and water in red. The total system and solvent densities are shown in cyan and purple, respectively. The GDS is shown by the orange dashed line and the interfacial region is shaded gray. Reprinted with permission from reference [206]. Copyright Elsevier. system, the interfacial width is larger than for the other systems, and one observes a more significant solvent penetration in the C8 system, as compared to the C18 and C30 stationary phases. This effect was interpreted by examining the snapshots in figure 11. The C8 phase is incapable of completely covering the silica surface, and thus the filaments of solvent molecules bridging the substrate and the mobile-phase region prevail. In the case of C18 and C30 phases, they form a continuous layer with respect to the solvent penetration, and solvent bridges are extremely rare for C18 phase. In other words, the C8 phase forms a more heterogeneous surface layer. It should be stressed that the water and methanol densities in the in- terior region of the C30 phase are about one order of magnitude smaller than the corresponding densities for the C18 phase. The simulations by Siepmann et al. indicated, that the retention of butane increases with an increas- ing chain length, which is in agreement with experimental data. The nonpolar alkane solute shows a mixed retention mechanism with contributions from adsorption in the interfacial region and partition- ing into the stationary phase. For the C1 phase, retention is dominated by adsorption. This phase is too thin to allow for partitioning. However, the adsorption of butane on a “hard” surface of C1 is different from adsorption on a more flexible, or “soft” surfaces involving longer chains. Siepmann et al. have dis- cussed not only the mechanism of retention of butane but also other molecules, and their calculations have clearly demonstrated that advanced computer simulations can be very helpful in explaining the mechanism of chromatographic separation. The absorption (confinement) of free linear chains into a polymer brush was studied with respect to chain size M and the Flory-Huggins type compatibility parameter Â with the brush by means of MC sim- ulations at various grafting densities using a bead-spring model. Different concentrations of free chains 12601-27 S. Sokołowski, J. Ilnytskyi, O. Pizio were examined. Contrary to the case of Â Æ 0, where all species are almost completely ejected by the polymer brush irrespective of their length M , it was found that for Â Ç 0, the absorbed amount under- goes a sharp crossover fromweak to strong absorption, discriminating between oligomers, 1ÉM É 8, and longer chains. For a moderately dense brush, the longer species, populate predominantly the deep inner part of the brush, whereas in a dense brush they penetrate into the “fluffy” tail of the dense brush only. Gyration radius Rg and end-to-end distance Re of absorbed chains thereby scale with length M as free polymers in the bulk. Three distinct regimes of penetration kinetics of free chains were found regarding the length M , in which the time of absorption grows with M at a different rate [197]. To summarize this part, one could just cite some simulation works that target technology-related as- pects of polymer brushes. For instance, the authors of reference [198] remark that “via computer simula- tions, we demonstrate how a densely grafted layer of polymers, a brush, could be turned into an efficient switch through chemical modification of some of its end-segments. In this way, a surface coating with reversibly switchable properties can be constructed. We analyze the fundamental physical principle be- hind its function, a recently discovered surface instability, and demonstrate that the combination of a high grafting density, an inflated end-group size and a high degree of monodispersity are conditions for an optimal functionality of the switch.” In this way, atomistic simulations of polymer brush involving ma- terials are expected to increase their predictive power and help to create a variety of smart, switchable, and multifunctional surfaces and thin films. Applications of these versatile and multifunctional brush coatings are envisioned in many areas including fluid control, microfluidics, and thin film sensors as well as some objects of interest in molecular biology [199]. 8. Summary To summarize, in this text we have reviewed the application of DF approaches and computer sim- ulation methods for a description of fluid-tethered polymer brushes interfaces. Several illustrations of important findings were given. The topic is very ample and, consequently, we were unable to describe it in full detail. Albeit the success in application of the theories outlined in this review, they will benefit of further improvements along different lines. In particular, a better description of the model of tethered polymer chains is necessary. Thermodynamic perturbation theory of the first order involved in the majority of present studies does not permit to capture structural peculiarities of several systems of interest. DF ap- proaches encounter difficulties in an accurate description of attractive interparticle and long-range elec- trostatic interactions. A comparison of theoretical predictions with computer simulation data is, there- fore, highly desirable to avoid incorrect interpretation of some theoretical results. As for technical issues necessary to be solved in the nearest future, it is of interest to include statistical mechanical analogue of mass action law into the available theories. This would permit to explore the effects of chemical associ- ation more in detail. We are convinced that theoretical developments in the area would permit to find novel phenomena faster and with less cost compared to other methods. Our review also presents selected applications of computer simulations in the studies of systems involving tethered chains. First, we briefly described some results obtained from DPD coarse-grained simulations and then we outlined the application of microscopic MC and MD simulations to explain the mechanism of separation in RPLC. With an increase in performance of computer facilities, the role of simulations will increase, and the simulation methods in many cases will replace not only theoretical approaches, but also some experiments. This is particularly true for the cases of complex systems and for out of equilibrium phenomena, such as fluid flows. One of possible directions toward the development of simulation techniques to describe brush-fluid interfaces would be the application of multiscale modelling using field-theoretic methodologies [213]. According to those methods, some fragments of polymers are treated on a fully atomistic (or even ab initio) level and the information is used for parametrization of coarse-grained potentials [214]. Such a technique has already been used to describe the formation of mi- celles and vesicles in amphiphilic systems [215]. However, when performing coarse-graining, one should always keep in mind that the coarse-grained systems are physically different from the initial, molecu- lar systems and that the coarse-graining methods cannot be used as a black box and require a thorough cross-checking. 12601-28 Fluid-grafted chains interfaces Acknowledgements This work was supported by the EU under IRSES Project STCSCMBS 268498. References 1. De Gennes P.-G., Macromolecules, 1980, 13, 1069; doi:10.1021/ma60077a009. 2. De Gennes P. G., Scaling Concepts in Polymer Physics, Cornell University Press, Ithaca, 1979. 3. Alexander S., J. Phys. (France), 1977, 38, 983; doi:10.1051/jphys:01977003808098300. 4. Binder K., Eur. J. Phys. E, 2002, 9, 293; doi:10.1140/epje/i2002-10076-2. 5. Schweizer K.S., Curro J.G., Adv. Chem. Phys., 1997, 98, 1; doi:10.1002/9780470141571. 6. Heine D.R., Grest G.S., Curro J.G., Adv. Polym. Sci., 2005, 173, 209; doi:10.1007/b99431. 7. Hall L.M., Jayaraman A., Schweizer K.S., Curr. Opin. Solid State Mat. Sci. 2010, 14, 38; doi:10.1016/j.cossms.2009.08.004. 8. Chandler D., In: The liquid state of matter: Fluids, Simple and Complex (Studies in statistical mechanics), Mon- troll E.W., Lebowitz J.L. (Eds.), North Holland, Amsterdam, 1982, 275. 9. Gast A.P., Leibler L., Macromolecules, 1986, 19, 686; doi:10.1021/ma00157a035. 10. Hasegawa R., Aoki Y., Doi M., Macromolecules, 1996, 29, 6656; doi:10.1021/ma960365x. 11. Semenov A.N., Macromolecules, 1993, 26, 2273; doi:10.1021/ma00061a021. 12. Jayaraman A., Schweizer K.S., Macromolecules, 2008, 41, 9430; doi:10.1021/ma801722m. 13. Jayaraman A., Schweizer K.S., J. Chem. Phys., 2008, 128, 164904; doi:10.1063/1.2907717. 14. Jayaraman A., Schweizer K.S., Langmuir, 2008, 24, 11119; doi:10.1021/la801432b. 15. Edwards S.F., Proc. Phys. Soc., 1965, 85, 613; doi:10.1088/0370-1328/85/4/301. 16. DiMarzio E.A., J. Chem. Phys., 1965, 42, 2101; doi:10.1063/1.1696251. 17. Vrij A., J. Polym. Sci. Part A2: Pol. Phys., 1968, 6, 1919; doi:10.1002/pol.1968.160061107. 18. Helfand E., Tagami Y., J. Chem. Phys., 1972, 56, 3592; doi:10.1063/1.1677735. 19. Helfand E., J. Chem. Phys., 1975, 62, 999; doi:10.1063/1.430517. 20. Schmid F., In: Computational Methods in Surface and Colloid Science, Borowko M. (Ed.), Marcel Dekker, New York, 2000, 631–685. 21. Matsen M.W., J. Phys.: Condens. Matter, 2002, 14, R21; doi:10.1088/0953-8984/14/2/201. 22. Fleer G.J., Cohen Stuart M.A., Scheutjens J.M.H.M., Cosgrove T., Vincent B., Polymers at Interfaces, Chapman and Hall, London, 1993. 23. Freed K., J. Chem. Phys., 1995, 103, 3230; doi:10.1063/1.470255. 24. Tang H., Freed K.F., J. Chem. Phys., 1991, 94, 1572; doi:10.1063/1.459962. 25. Lin E.K., Gast A.P., Macromolecules, 1996, 29, 390; doi:10.1021/ma9505282. 26. Lutsko J.F., Adv. Chem. Phys., 2010, 144, 1; doi:10.1002/9780470564318. 27. Chandler D., McCoy J.D., Singer S.J., J. Chem. Phys., 1986, 85, 5971; doi:10.1063/1.451510. 28. Chandler D., McCoy J.D., Singer S.J., J. Chem. Phys., 1986, 85, 5977; doi:10.1063/1.451511. 29. Ladanyi B.M., Chandler D., J. Chem. Phys., 1975, 62, 4308; doi:10.1063/1.431001. 30. Chandler D., Pratt L.R., J. Chem. Phys., 1976, 65, 2925; doi:10.1063/1.433529. 31. Woodward C.E., J. Chem. Phys., 1991, 94, 3183; doi:10.1063/1.459787. 32. Evans R., In: Fundamentals of Inhomogenous Fluids, Henderson D. (Ed.), Marcel Dekker, New York, 1992, 85– 175. 33. Woodward C.E., Forsman J., Phys. Rev. Lett., 2008, 100, 098301; doi:10.1103/PhysRevLett.100.098301. 34. Woodward C.E., Forsman J., Phys. Rev. E, 2006, 74, 010801(R); doi:10.1103/PhysRevE.74.010801. 35. Wertheim M.S., J. Stat. Phys., 1984, 35, 19; doi:10.1007/BF01017362. 36. Wertheim M.S., J. Stat. Phys., 1984, 35, 35; doi:10.1007/BF01017363. 37. Wertheim M.S., J. Stat. Phys., 1986, 42, 459; doi:10.1007/BF01127721. 38. Kierlik E., Rosinberg M.L., J. Chem. Phys., 1992, 97, 9222; doi:10.1063/1.463298. 39. Kierlik E., Rosinberg M.L., J. Chem. Phys., 1993, 99, 3950; doi:10.1063/1.466142. 40. Kierlik E., Rosinberg M.L., J. Chem. Phys., 1994, 100, 1716; doi:10.1063/1.466599. 41. Phan S., Kierlik E., Rosinberg M.L., Yu H., Stell G., J. Chem. Phys., 1993, 99, 5326; doi:10.1063/1.465976. 42. Yu Y.-X., Wu J., J. Chem. Phys., 2002, 117, 2368; doi:10.1063/1.1491240. 43. Bryk P., Sokołowski S., J. Chem. Phys., 2004, 121, 11314; doi:10.1063/1.1814075. 44. Bryk P., Pizio O., Sokołowski S., J. Chem. Phys., 2005, 122, 194904; doi:10.1063/1.1898484. 45. Li Z., Cao D., Wu J., J. Chem. Phys., 2005, 122, 224701; doi:10.1063/1.1924603. 46. Li Z., Cao D., Wu J., J. Chem. Phys., 2005, 122, 174708; doi:10.1063/1.1886685. 12601-29 http://dx.doi.org/10.1021/ma60077a009 http://dx.doi.org/10.1051/jphys:01977003808098300 http://dx.doi.org/10.1140/epje/i2002-10076-2 http://dx.doi.org/10.1002/9780470141571 http://dx.doi.org/10.1007/b99431 http://dx.doi.org/10.1016/j.cossms.2009.08.004 http://dx.doi.org/10.1021/ma00157a035 http://dx.doi.org/10.1021/ma960365x http://dx.doi.org/10.1021/ma00061a021 http://dx.doi.org/10.1021/ma801722m http://dx.doi.org/10.1063/1.2907717 http://dx.doi.org/10.1021/la801432b http://dx.doi.org/10.1088/0370-1328/85/4/301 http://dx.doi.org/10.1063/1.1696251 http://dx.doi.org/10.1002/pol.1968.160061107 http://dx.doi.org/10.1063/1.1677735 http://dx.doi.org/10.1063/1.430517 http://dx.doi.org/10.1088/0953-8984/14/2/201 http://dx.doi.org/10.1063/1.470255 http://dx.doi.org/10.1063/1.459962 http://dx.doi.org/10.1021/ma9505282 http://dx.doi.org/10.1002/9780470564318 http://dx.doi.org/10.1063/1.451510 http://dx.doi.org/10.1063/1.451511 http://dx.doi.org/10.1063/1.431001 http://dx.doi.org/10.1063/1.433529 http://dx.doi.org/10.1063/1.459787 http://dx.doi.org/10.1103/PhysRevLett.100.098301 http://dx.doi.org/10.1103/PhysRevE.74.010801 http://dx.doi.org/10.1007/BF01017362 http://dx.doi.org/10.1007/BF01017363 http://dx.doi.org/10.1007/BF01127721 http://dx.doi.org/10.1063/1.463298 http://dx.doi.org/10.1063/1.466142 http://dx.doi.org/10.1063/1.466599 http://dx.doi.org/10.1063/1.465976 http://dx.doi.org/10.1063/1.1491240 http://dx.doi.org/10.1063/1.1814075 http://dx.doi.org/10.1063/1.1898484 http://dx.doi.org/10.1063/1.1924603 http://dx.doi.org/10.1063/1.1886685 S. Sokołowski, J. Ilnytskyi, O. Pizio 47. Li Z., Wu J., J. Chem. Phys., 2007, 126, 144904; doi:10.1063/1.2715595. 48. MacDowell L.G., Bryk P., Phys. Rev. E, 2007, 75, 061609; doi:10.1103/PhysRevE.75.061609. 49. Malijevský Al., Bryk P., Sokołowski S., Phys. Rev. E, 2005, 72, 032801; doi:10.1103/PhysRevE.72.032801. 50. Xu X., Cao D., Zhang X., Wang W., Phys. Rev. E, 2009, 79, 021805; doi:10.1103/PhysRevE.79.021805. 51. Merlitz H., Ciu W., Wu C.-X., Sommer J.-U., Macromolecules, 2013, 46, 1248; doi:10.1021/ma302417j. 52. Bryk P., Roth R., Mecke K.R., Dietrich S., Phys. Rev. E, 2003, 68, 031602; doi:10.1103/PhysRevE.68.031602. 53. Cao D., Wu J., J. Chem. Phys., 2004, 121, 4210; doi:10.1063/1.1774983. 54. Cao D., Wu J., Langmuir, 2006, 22, 2712; doi:10.1021/la0527588. 55. Jiang T., Li Z., Wu J., Macromolecules, 2007, 40, 334; doi:10.1021/ma061939t. 56. Xu X., Cao D., Wu J., Soft Matter, 2010, 6, 4631; doi:10.1039/C0SM00034E. 57. Borówko M., Rżysko W., Sokołowski S., Staszewski T., J. Chem. Phys., 2007, 126, 214703; doi:10.1063/1.2743399. 58. Pizio O., Patrykiejew A., Sokołowski S., J. Phys. Chem. C, 2007, 111, 15743; doi:10.1021/jp0736847. 59. Patrykiejew A., Sokołowski S., Tscheliessnig R., Fischer J., Pizio O., J. Phys. Chem. B, 2008, 112, 4552; doi:10.1021/jp710978t. 60. Pizio O., Pusztai L., Sokołowska Z., Sokołowski S., J. Chem. Phys., 2009, 130, 134501; doi:10.1063/1.3103266. 61. Borówko M., Rżysko W., Sokołowski S., Staszewski T., J. Phys. Chem. B, 2009, 113, 4763; doi:10.1021/jp811143n. 62. Borówko M., Patrykiejew A., Sokołowski S., Staszewski T., Coll. Czech. Chem. Comm., 2010, 75, 221; doi:10.1135/cccc2009104. 63. Borówko M., Sokołowski S., Staszewski T., J. Chromat. A, 2011, 1218, 711; doi:10.1016/j.chroma.2010.12.029. 64. Pizio O., Borówko M., Rżysko W., Staszewski T., Sokołowski S., J. Chem. Phys., 2008, 128, 044702; doi:10.1063/1.2829247. 65. Borówko M., Patrykiejew A., Sokołowski S., J. Chem. Phys., 2011, 135, 054703; doi:10.1063/1.3617413. 66. Pizio O., Sokołowski S., Sokołowska Z., J. Chem. Phys., 2011, 134, 214702; doi:10.1063/1.3597773. 67. Emborsky C.P., Feng Z., Cox K.R., Chapman W.G., Fluid Phase Equilib., 2011, 306, 15; doi:10.1016/j.fluid.2011.02.007. 68. Marshall B.D., Chapman W.G., J. Phys. Chem. B, 2011, 115, 15036; doi:10.1021/jp2078677. 69. Jain S., Jog P., Weinhold J., Srivastava R., Chapman W.G., J. Chem. Phys., 2008, 128, 154910; doi:10.1063/1.2902976. 70. Yethiraj A., Woodward C.E., J. Chem. Phys., 1995, 102, 5499; doi:10.1063/1.469279. 71. Chandra N., Patra C.N., Yethiraj A., J. Chem. Phys., 2003, 118, 4702; doi:10.1063/1.1543141. 72. Müller M., MacDowell L.G., Yethiraj A., J. Chem. Phys., 2003, 118, 2929; doi:10.1063/1.1535893. 73. Müller M., MacDowell L.G., J. Phys.: Condens. Matter, 2003, 15, R609; doi:10.1088/0953-8984/15/19/201. 74. Daoulas K.Ch., Müller M., de Pablo J.J., Nealey P.F., Smith G.D., Soft Matter, 2006, 2, 573; doi:10.1039/B602610A. 75. Wang J., Müller M., Macromolecules, 2009, 42, 2251; doi:10.1021/ma8026047. 76. Yethiraj A., Fynewever H., Shew C.-Y., J. Chem. Phys., 2001, 114, 4323; doi:10.1063/1.1348031. 77. Cao D., Jiang T., Wu J., J. Chem. Phys., 2006, 124, 164904; doi:10.1063/1.2186323. 78. Chen H., Ye Z., Cai Y., Liu H., Hu Y., J. Phys. Chem. B, 2007, 111, 5927; doi:10.1021/jp068784b. 79. Bryk P., MacDowell L., J. Chem. Phys., 2009, 129, 104901; doi:10.1063/1.2976339. 80. Melenkevitz J., Curro J.G., Schweizer K.S., J. Chem. Phys., 1993, 99, 5571; doi:10.1063/1.465949. 81. Grayce C.J., Schweizer K.S., J. Chem. Phys., 1994, 100, 6846; doi:10.1063/1.467044. 82. Grayce C.J., Yethiraj A., Schweizer K.S., J. Chem. Phys., 1994, 100, 6857; doi:10.1063/1.467045. 83. Bryk P., MacDowell L.G., J. Chem. Phys., 2011, 135, 204901; doi:10.1063/1.3662139. 84. Israelachvili J.N., Intermolecular and surface forces, Academic Press, Boston, 1992. 85. Kreer T., Metzger S., Müller M., Binder K., Baschnagel J., J. Chem. Phys., 2004, 120, 4012; doi:10.1063/1.1642615. 86. Dimitrov D., Milchev A., Binder K., Macromol. Symp., 2007, 252, 47; doi:10.1002/masy.200750605. 87. Pastorino C., Binder K., Müller M., Macromolecules, 2009, 42, 401; doi:10.1021/ma8015757. 88. Shagolsem L., Sommer J.-U., Soft Matter, 2012, 8, 11328; doi:10.1039/C2SM26531A. 89. Jiang J., Xu X., Cao D., Phys. Rev. E, 2012, 86, 041805; doi:10.1103/PhysRevE.86.041805. 90. Milner S.T., Witten T.A., Cates M.E., Macromolecules, 1988, 21, 2610; doi:10.1021/ma00186a051. 91. Grest G.S., Murat M., Macromolecules, 1993, 26, 3108; doi:10.1021/ma00064a019. 92. Zhulina E., Borisov O.V., Pryamitsyn V.A., Birshtein T.M., Macromolecules, 1991, 24, 140; doi:10.1021/ma00001a023. 93. Eiser E., Klein J., Witten T.A., Fetters L.J., Phys. Rev. Lett., 1999, 82, 5076; doi:10.1103/PhysRevLett.82.5076. 94. Wijmans C.M., Leermakers F.A.M., Fleer G.J., J. Colloid Interface Sci., 1994, 167, 124; doi:10.1006/jcis.1994.1340. 95. Zhulina E., Singh C., Balasz A.C., Langmuir, 1999, 15, 3935; doi:10.1021/la981406g. 96. Pincus P., Macromolecules, 1991, 24, 2912; doi:10.1021/ma00010a043. 97. Balastre M., Schorr P., Yang J.C., Mays J.W., Tirrell M.V., Macromolecules, 2002, 35, 9480; doi:10.1021/ma011875g. 98. Borisov O.V., Birshtein T.M., Zhulina E.B., J. Phys. II France, 1991, 1, 521; doi:10.1051/jp2:1991186. 12601-30 http://dx.doi.org/10.1063/1.2715595 http://dx.doi.org/10.1103/PhysRevE.75.061609 http://dx.doi.org/10.1103/PhysRevE.72.032801 http://dx.doi.org/10.1103/PhysRevE.79.021805 http://dx.doi.org/10.1021/ma302417j http://dx.doi.org/10.1103/PhysRevE.68.031602 http://dx.doi.org/10.1063/1.1774983 http://dx.doi.org/10.1021/la0527588 http://dx.doi.org/10.1021/ma061939t http://dx.doi.org/10.1039/C0SM00034E http://dx.doi.org/10.1063/1.2743399 http://dx.doi.org/10.1021/jp0736847 http://dx.doi.org/10.1021/jp710978t http://dx.doi.org/10.1063/1.3103266 http://dx.doi.org/10.1021/jp811143n http://dx.doi.org/10.1135/cccc2009104 http://dx.doi.org/10.1016/j.chroma.2010.12.029 http://dx.doi.org/10.1063/1.2829247 http://dx.doi.org/10.1063/1.3617413 http://dx.doi.org/10.1063/1.3597773 http://dx.doi.org/10.1016/j.fluid.2011.02.007 http://dx.doi.org/10.1021/jp2078677 http://dx.doi.org/10.1063/1.2902976 http://dx.doi.org/10.1063/1.469279 http://dx.doi.org/10.1063/1.1543141 http://dx.doi.org/10.1063/1.1535893 http://dx.doi.org/10.1088/0953-8984/15/19/201 http://dx.doi.org/10.1039/B602610A http://dx.doi.org/10.1021/ma8026047 http://dx.doi.org/10.1063/1.1348031 http://dx.doi.org/10.1063/1.2186323 http://dx.doi.org/10.1021/jp068784b http://dx.doi.org/10.1063/1.2976339 http://dx.doi.org/10.1063/1.465949 http://dx.doi.org/10.1063/1.467044 http://dx.doi.org/10.1063/1.467045 http://dx.doi.org/10.1063/1.3662139 http://dx.doi.org/10.1063/1.1642615 http://dx.doi.org/10.1002/masy.200750605 http://dx.doi.org/10.1021/ma8015757 http://dx.doi.org/10.1039/C2SM26531A http://dx.doi.org/10.1103/PhysRevE.86.041805 http://dx.doi.org/10.1021/ma00186a051 http://dx.doi.org/10.1021/ma00064a019 http://dx.doi.org/10.1021/ma00001a023 http://dx.doi.org/10.1103/PhysRevLett.82.5076 http://dx.doi.org/10.1006/jcis.1994.1340 http://dx.doi.org/10.1021/la981406g http://dx.doi.org/10.1021/ma00010a043 http://dx.doi.org/10.1021/ma011875g http://dx.doi.org/10.1051/jp2:1991186 Fluid-grafted chains interfaces 99. Ross R.S., Pincus P., Macromolecules, 1992, 25, 2177; doi:10.1021/ma00034a018. 100. Zhulina E.B., Borisov O.V., Birshtein T.M., J. Phys. II France, 1992, 2, 63; doi:10.1051/jp2:1992113. 101. Borisov O.V., Zhulina E.B., Birshtein T.M., Macromolecules, 1994, 27, 4795; doi:10.1021/ma00095a021. 102. Zhulina E.B., Birshtein T.M., Borisov O.V., Macromolecules, 1995, 28, 1491; doi:10.1021/ma00109a021. 103. Wu T., Gong P., Szleifer I., Vlcek P., Subr V., Genzer J., Macromolecules, 2007, 40, 8756; doi:10.1021/ma0710176. 104. Currie E.P.K., Sieval A.B., Fleer G.J., Stuart M.A.C., Langmuir, 2000, 16, 8324; doi:10.1021/la991528o. 105. Witte K.N., Kim S., Won Y.-Y., J. Phys. Chem. B, 2009, 113, 11076; doi:10.1021/jp809814j. 106. Gong P., Wu T., Genzer J., Szleifer I., Macromolecules, 2007, 40, 8765; doi:10.1021/ma071018y. 107. Wu J., Li Z., Annu. Rev. Phys. Chem., 2007, 58, 85; doi:10.1146/annurev.physchem.58.032806.104650. 108. Li Z.D., Wu J.Z., J. Phys. Chem. B, 2006, 110, 7473; doi:10.1021/jp060127w. 109. Wang Z., Liu L., Neretnieks I., J. Phys.: Condens. Matter, 2011, 23, 175002; doi:10.1088/0953-8984/23/17/175002. 110. Wang Z., Liu L., Neretnieks I., J. Chem. Phys., 2011, 135, 244107; doi:10.1063/1.3672001. 111. Wang L., Liang H., Wu J.-Z., J. Chem. Phys., 2010, 133, 044906; doi:10.1063/1.3463426. 112. Jiang J., Feng J., Liu H., Hu Y., J. Chem. Phys., 2006, 124, 144908; doi:10.1063/1.2186316 113. Ballauff M., Borisov O., Curr. Opin. Colloid Interface Sci., 2006, 11, 316; doi:10.1016/j.cocis.2006.12.002. 114. Biesalski M., Johannsmann D., Ruhe J., J. Chem. Phys., 2004, 120, 8807; doi:10.1063/1.1690242. 115. Houbenov N., Minko S., Stamm M., Macromolecules 2003, 36, 5897; doi:10.1021/ma0341869. 116. Hehmeyer O.J., Stevens M.J., J. Chem. Phys., 2005, 122, 134909; doi:10.1063/1.1871937. 117. Tsori Y., Andelman D., Joanny J.-F., Europhys. Lett., 2008, 82, 46001; doi:10.1209/0295-5075/82/46001. 118. Adiga S.P., Brenner D.W., Nano Lett., 2002, 2, 567; doi:10.1021/nl025527j. 119. Adiga S.P., Brenner D.W., Nano Lett., 2005, 5, 2509; doi:10.1021/nl051843x. 120. Romet-Lemonne G., Daillant J., Guenoun P., Yang J., Mays J.W., Phys. Rev. Lett., 2004, 93, 148301; doi:10.1103/PhysRevLett.93.148301. 121. Ahrens H., Förster S., Helm C.A., Kumar N.A., Naji A., Netz R.R., Seidel C., J. Phys. Chem. B, 2004, 108, 16870; doi:10.1021/jp049553c. 122. Jiang D., Jin Z., Wu J.-Z., Nano Lett., 2011, 11, 5373; doi:10.1021/nl202952d. 123. Bazant M.Z., Storey B.D., Kornyshev A.A., Phys. Rev. Lett., 2011, 106, 046102; doi:10.1103/PhysRevLett.106.046102. 124. Chmiola J., Yushin G., Gogotsi Y., Portet C., Simon P., Taberna P.-L., Science, 2006, 313, 1760; doi:10.1126/science.1132195. 125. Largeot C., Portet C., Chmiola J., Taberna P.-L., Gogotsi Y., Simon P., J. Am. Chem. Soc., 2008, 130, 2730; doi:10.1021/ja7106178. 126. Feng G., Cummings P.T., J. Phys. Chem. Lett., 2011, 2, 2859; doi:10.1021/jz201312e. 127. Dobrynin A.V., Curr. Opin. Colloid Interface Sci., 2008, 13, 376; doi:10.1016/j.cocis.2008.03.006. 128. Pizio O., Sokołowski S., J. Chem. Phys., 2013, 138, 204715; doi:10.1063/1.4807777. 129. Pizio O., Sokołowski S., Sokołowska Z., J. Chem. Phys., 2012, 137, 234705; doi:10.1063/1.4771919. 130. Borówko M., Sokołowski S., Staszewski T., Sokołowska Z., Ilnytskyi J.M., J. Chem. Phys., 2012, 137 074707; doi:10.1063/1.4771919. 131. Fedorov M.V., Georgi N., Kornyshev A.A., Electrochem. Commun., 2010, 12, 296; doi:10.1016/j.elecom.2009.12.019. 132. Fedorov M.V., Kornyshev A.A., J. Phys. Chem. B, 2008, 112, 11868; doi:10.1021/jp803440q. 133. Wu J., Jiang T., Jiang D., Jin Z., Henderson D., Soft Matter, 2011, 7, 11222; doi:10.1039/C1SM06089A . 134. Henderson D., Wu J., J. Phys. Chem. B, 2012, 116, 2520; doi:10.1021/jp212082k. 135. Lockett V., Sedev R., Ralston J., HorneM., Rodopoulos T., J. Phys. Chem. C, 2008, 112, 7486; doi:10.1021/jp7100732. 136. Alam M.T., Islam M.M., Okajima T., Ohsaka T., J. Phys. Chem. C, 2009, 113, 6569; doi:10.1021/jp810865t. 137. Kondrat S., Perez C.R., Presser V., Gogotsi Y., Kornyshev A.A., Energy Environ. Sci., 2012, 5, 6474. doi:10.1039/c2ee03092f. 138. Patra M., Linse P., Nano Lett., 2006, 6, 133; doi:10.1021/nl051611y. 139. Wilder K., Quate C.F., Singh B., Alvis R., Arnold W.H., J. Vac. Sci. Technol. B, 1996, 14, 4004; doi:10.1116/1.588632. 140. Grest G.S., Murat M., In: Monte Carlo and Molecular Dynamics Simulations in Polymer Science, Binder K. (Ed.), Oxford Press, New York, 1995, 476–569. 141. Binder K., Milchev A., Baschnagel J., Annu. Rev. Mater. Sci., 1996, 26, 107; doi:10.1146/annurev.ms.26.080196.000543. 142. Slater G.W., Holm C., Chubynsky M.V., de Haan H.W., Dube A., Grass K., Hickey O.A., Kingsburry C., Sean D., Shendruk T.N., Zhan L., Electrophoresis, 2009, 30 792; doi:10.1002/elps.200800673. 143. Binder K., Milchev A., J. Polymer Sci. Part B: Polymer Physics, 2012, 50, 1515; doi:10.1002/polb.23168. 144. Binder K., Kreer T., Milchev A., Soft Matter, 2011, 7, 7159; doi:10.1039/C1SM05212H. 145. Neelov I.M., Binder K., Macromol. Theory Simul., 1995, 4, 119; doi:10.1002/mats.1995.040040108. 12601-31 http://dx.doi.org/10.1021/ma00034a018 http://dx.doi.org/10.1051/jp2:1992113 http://dx.doi.org/10.1021/ma00095a021 http://dx.doi.org/10.1021/ma00109a021 http://dx.doi.org/10.1021/ma0710176 http://dx.doi.org/10.1021/la991528o http://dx.doi.org/10.1021/jp809814j http://dx.doi.org/10.1021/ma071018y http://dx.doi.org/10.1146/annurev.physchem.58.032806.104650 http://dx.doi.org/10.1021/jp060127w http://dx.doi.org/10.1088/0953-8984/23/17/175002 http://dx.doi.org/10.1063/1.3672001 http://dx.doi.org/10.1063/1.3463426 http://dx.doi.org/10.1063/1.2186316 http://dx.doi.org/10.1016/j.cocis.2006.12.002 http://dx.doi.org/10.1063/1.1690242 http://dx.doi.org/10.1021/ma0341869 http://dx.doi.org/10.1063/1.1871937 http://dx.doi.org/10.1209/0295-5075/82/46001 http://dx.doi.org/10.1021/nl025527j http://dx.doi.org/10.1021/nl051843x http://dx.doi.org/10.1103/PhysRevLett.93.148301 http://dx.doi.org/10.1021/jp049553c http://dx.doi.org/10.1021/nl202952d http://dx.doi.org/10.1103/PhysRevLett.106.046102 http://dx.doi.org/10.1126/science.1132195 http://dx.doi.org/10.1021/ja7106178 http://dx.doi.org/10.1021/jz201312e http://dx.doi.org/10.1016/j.cocis.2008.03.006 http://dx.doi.org/10.1063/1.4807777 http://dx.doi.org/10.1063/1.4771919 http://dx.doi.org/10.1063/1.4771919 http://dx.doi.org/10.1016/j.elecom.2009.12.019 http://dx.doi.org/10.1021/jp803440q http://dx.doi.org/10.1039/C1SM06089A http://dx.doi.org/10.1021/jp212082k http://dx.doi.org/10.1021/jp7100732 http://dx.doi.org/10.1021/jp810865t http://dx.doi.org/10.1039/c2ee03092f http://dx.doi.org/10.1021/nl051611y http://dx.doi.org/10.1116/1.588632 http://dx.doi.org/10.1146/annurev.ms.26.080196.000543 http://dx.doi.org/10.1002/elps.200800673 http://dx.doi.org/10.1002/polb.23168 http://dx.doi.org/10.1039/C1SM05212H http://dx.doi.org/10.1002/mats.1995.040040108 S. Sokołowski, J. Ilnytskyi, O. Pizio 146. Saphiannikova M.G., Pryamitsyn V.A., Cosgrove T., Macromolecules, 1998, 31, 6662; doi:10.1021/ma971808f. 147. Tupitsyna A.I., Darinskii A.A., Birshtein T.M., Amoskov V.M., Emri I., Macromol. Theory Simul., 2004, 13, 771; doi:10.1002/mats.200400061. 148. Carlsson T., Kamerlin N., Arteca G.A., Elvingson C., Phys. Chem. Chem. Phys., 2011, 13, 16084; doi:10.1039/c1cp21433k. 149. Saphiannikova M.G., Pryamitsyn V.A., Birshtein T.M., Macromolecules, 2000, 33, 2740; doi:10.1021/ma9911400. 150. Español P., Warren P.B., Europhys. Lett., 1995, 30, 191; doi:10.1209/0295-5075/30/4/001. 151. Moeendarbary E., Ng T.Y., Zangeneh M., Int. J. Appl. Mech., 2009, 1, 737; doi:10.1142/S1758825109000381. 152. Frenkel D., Smit B., Understanding Molecular Simulation: From Algorithms to Applications (2nd Edn.), Aca- demic Press, London, 2002. 153. Pivkin I.V., Caswell B., Karniadakis G.E., In: Reviews in Computational Chemistry, Vol. 27, K. B. Lipkowitz (Ed.), John Wiley & Sons, Inc., Hoboken, NJ, 2011. 154. Groot R.D., Madden T.J., J. Chem. Phys., 1998, 108, 8713 doi:10.1063/1.476300. 155. Visser D.C., Hoefsloot H.C.J., Iedema P.D., J. Comput. Phys., 2005, 205, 626; doi:10.1016/j.jcp.2004.11.020. 156. Goujon F., Malfreyt P., Tildesley D.J., Chem. Phys. Phys. Chem., 2004, 5, 457; doi:10.1002/cphc.200300901. 157. Pivkin I.V., Karniadakis G.E., J. Comput. Phys., 2005, 207, 114; doi:10.1016/j.jcp.2005.01.006. 158. Pivkin I.V., Karniadakis G.E., J. Chem. Phys., 2006, 124, 184101; doi:10.1063/1.2191050. 159. Gibson J.B., Zhang K., Chen K., Chynoweth S., Manke C.W., Mol. Simul., 1999, 23, 1; doi:10.1080/08927029908022109. 160. Ilnytskyi J.M., Patsahan T., Sokołowski S., J. Chem. Phys., 2011, 134, 204903; doi:10.1063/1.3592562. 161. Ilnytskyi J.M., Sokołowski S., Patsahan T., Condens. Matter Phys., 2013, 16, No. 1, 13606; doi:DOI:10.5488/CMP.16.13606. 162. Groot R.D., Warren P.B., J. Chem. Phys., 1997, 107, 4423; doi:10.1063/1.474784. 163. Malfreyt P., Tildesley D.J., Langmuir, 2000, 16, 4732; doi:10.1021/la991396z. 164. Goujon F., Malfreyt P., Tildesley D.J., Mol. Phys., 2005, 103, 2675; doi:10.1080/00268970500134706. 165. Gibson J.B., Chen K., Chynoweth S., J. Colloid Interface Sci., 1998, 206, 464; doi:10.1006/jcis.1998.5734. 166. Wijmans C.M., Smit B., Macromolecules, 2002, 35 7138; doi:10.1021/ma020086b. 167. Pastorino C., Binder K., Kreer T., Müller M., J. Chem. Phys., 2006, 124, 064902; doi:10.1063/1.2162883. 168. Pastorino C., Binder K., Müller M., Macromolecules, 2009, 42, 401; doi:10.1021/ma8015757. 169. Li C.-S., Wu W.-C., Sheng Y.-S., Chen W.-C., J. Chem. Phys., 2008, 128, 154908; doi:10.1063/1.2904866. 170. Xue Y.-H., Liu H., Lu Z.-Y., Liang X.-Z., J. Chem. Phys., 2010, 132, 044903; doi:10.1063/1.3299730. 171. Petrus P., Lísal M., Brennan J.K., Langmuir, 2010, 26, 3695; doi:10.1021/la903200j. 172. Petrus P., Lísal M., Brennan J.K., Langmuir, 2010, 26, 14680; doi:10.1021/la102666g. 173. Grosberg A.Y., Khokhlov A.R., Statistical Physics of Macromolecules, AIP series in polymers and complex mate- rials, Woodbury, New York, 1996. 174. Murat M., Grest G.S., Macromolecules, 1996, 29, 8282; doi:10.1021/ma961267e. 175. Rühle V., Junghans C., Lukyanov A., Kremer K., Andrienko D., J. Chem. Theory Comput., 2009, 5, 3211; doi:10.1021/ct900369w. 176. Rühle V., Junghans C., Macromol. Theory Simul., 2011, 20, 472; doi:10.1002/mats.201100011. 177. Füchslin R.M., Fellermann H., Eriksson A., Ziock H.-J., J. Chem. Phys., 2009, 130, 214102; doi:10.1063/1.3143976. 178. Lyubartsev A.P., Karttunen M., Vattulainen I., Laaksonen A., Soft Mater., 2002, 1, 121; doi:10.1081/SMTS-120016746. 179. Milchev A., Rostiashvili V., Bhattacharya S., Vilgis T., Polymer Chain Adsorption on a Solid Surface: Scaling Arguments and Computer Simulations In: Nanophenomena at Surfaces, Michailov M., (Ed.), Springer Series in Surface Sciences, 2011, 47, 185; doi:10.1007/978-3-642-16510-8_8. 180. Isa L., Amstad E., Schwenke K., Del Gado E., Ilg P., Kröger M., Reimhult E., Soft Matter, 2011, 7, 7663; doi:10.1039/C1SM05407D. 181. Paturej J., Milchev A., Rostiashvili V.G., Vilgis T.A., Macromolecules, 2012, 45, 4371; doi:10.1021/ma202671n. 182. MacDowell L.G., Müller M., J. Chem. Phys., 2006, 124, 084907; doi:10.1063/1.2172597. 183. Lai P.-Y., Zhulina E.B., J. Phys. II France, 1992, 2, 547; doi:10.1051/jp2:1992148. 184. Ohno K., Sakamoto T., Mingawa T., Okabe Y., Macromolecules, 2007, 40, 723; doi:10.1021/ma0613234.30. 185. Milner S.T., Witten T.A., Cates M.E., Macromolecules, 1988, 22, 2610; doi:10.1021/ma00186a051. 186. Grest G.S., J. Chem. Phys., 1996, 105, 5532; doi:10.1063/1.472395. 187. Lai P.Y., J. Chem. Phys., 1993, 98, 669; doi:10.1063/1.464611. 188. Descas R., Sommer J.-U., Blumen A., J. Chem. Phys., 2006, 124, 94701; doi:10.1063/1.2159479. 189. Weinhold J.D., Kumar S.K., J. Chem. Phys., 1994, 101, 4312; doi:10.1063/1.467481. 190. He G.-L., Merlitz H., Sommer J.-U., Wu C.-X., Macromolecules, 2007, 40, 6721; doi:10.1021/ma070983l. 191. He S.-Z., Merlitz H., Chen L., Sommer J.-U., Wu C.-X., Macromolecules, 2010, 43, 7845; doi:10.1021/ma101230v. 12601-32 http://dx.doi.org/10.1021/ma971808f http://dx.doi.org/10.1002/mats.200400061 http://dx.doi.org/10.1039/c1cp21433k http://dx.doi.org/10.1021/ma9911400 http://dx.doi.org/10.1209/0295-5075/30/4/001 http://dx.doi.org/10.1142/S1758825109000381 http://dx.doi.org/10.1063/1.476300 http://dx.doi.org/10.1016/j.jcp.2004.11.020 http://dx.doi.org/10.1002/cphc.200300901 http://dx.doi.org/10.1016/j.jcp.2005.01.006 http://dx.doi.org/10.1063/1.2191050 http://dx.doi.org/10.1080/08927029908022109 http://dx.doi.org/10.1063/1.3592562 http://dx.doi.org/DOI:10.5488/CMP.16.13606 http://dx.doi.org/10.1063/1.474784 http://dx.doi.org/10.1021/la991396z http://dx.doi.org/10.1080/00268970500134706 http://dx.doi.org/10.1006/jcis.1998.5734 http://dx.doi.org/10.1021/ma020086b http://dx.doi.org/10.1063/1.2162883 http://dx.doi.org/10.1021/ma8015757 http://dx.doi.org/10.1063/1.2904866 http://dx.doi.org/10.1063/1.3299730 http://dx.doi.org/10.1021/la903200j http://dx.doi.org/10.1021/la102666g http://dx.doi.org/10.1021/ma961267e http://dx.doi.org/10.1021/ct900369w http://dx.doi.org/10.1002/mats.201100011 http://dx.doi.org/10.1063/1.3143976 http://dx.doi.org/10.1081/SMTS-120016746 http://dx.doi.org/10.1007/978-3-642-16510-8_8 http://dx.doi.org/10.1039/C1SM05407D http://dx.doi.org/10.1021/ma202671n http://dx.doi.org/10.1063/1.2172597 http://dx.doi.org/10.1051/jp2:1992148 http://dx.doi.org/10.1021/ma0613234.30 http://dx.doi.org/10.1021/ma00186a051 http://dx.doi.org/10.1063/1.472395 http://dx.doi.org/10.1063/1.464611 http://dx.doi.org/10.1063/1.2159479 http://dx.doi.org/10.1063/1.467481 http://dx.doi.org/10.1021/ma070983l http://dx.doi.org/10.1021/ma101230v Fluid-grafted chains interfaces 192. Biesheuvel P.M., de Vos W.M., Amoskov V.M., Macromolecules, 2008, 41, 6254; doi:10.1021/ma800828c. 193. Minko S., J. Macromol. Sci. C: Polymer Rev., 2006, 46, 397; doi:10.1080/15583720600945402. 194. Dedinaite A., In: Advanced Chemistry of Monolayers at Interfaces, Imae T. (Ed.), Elsevier, 2007, Ch. 2.; doi:10.1016/S1573-4285(06)14002-8. 195. Klatte S.J., Beck T.L., J. Chem. Phys., 1996, 100, 5931; doi:10.1021/jp953301h. 196. Lippa K.A., Sander L.C., Mountain R.D., Anal. Chem., 2005, 77, 7862; doi:10.1021/ac051085v. 197. Milchev A., Egorov S.A., Binder K., J. Chem. Phys., 2010, 132, 184905; doi:10.1063/1.3414996. 198. Merlitz H., He G.-L., Wu C.-X., Sommer J.-U., Phys. Rev. Lett., 2009, 102, 115702; doi:10.1103/PhysRevLett.102.115702. 199. Uhlmann P., Merlitz H., Sommer J.-U., StammM., Macromol. Rapid Commun., 2009, 30, 732; doi:10.1002/marc.200900113. 200. Zhang L., Sun Y., Biochem. Eng. J., 2010, 48, 408; doi:10.1016/j.bej.2009.12.003. 201. Zhang L., Sun L., Siepmann J.I., Schure M.R., J. Chromatogr. A, 2005, 1079, 127; doi:10.1016/j.chroma.2005.03.124. 202. Rafferty J.L., Siepmann J.L., Schure M.R., J. Chromatogr. A, 2011, 1218, 2203; doi:10.1016/j.chroma.2011.02.012. 203. Rafferty J.L., Siepmann J.L., Schure M.R., J. Chromatogr. A, 2008, 1204 11; doi:10.1016/j.chroma.2008.07.037. 204. Rafferty J.L., Siepmann J.L., Schure M.R., J. Chromatogr. A, 2009, 1216, 2320; doi:10.1016/j.chroma.2008.12.088. 205. Rafferty J.L., Siepmann J.L., Schure M.R., J. Chromatogr. A, 2011, 1218, 9183; doi:10.1016/j.chroma.2011.10.043. 206. Rafferty J.L., Siepmann J.L., Schure M.R., J. Chromatogr. A, 2012, 1223, 24; doi:10.1016/j.chroma.2011.11.039. 207. Rafferty J.L., Siepmann J.L., Schure M.R., Fluid Phase Equilib., 2010, 290, 25; doi:10.1016/j.fluid.2009.10.010. 208. Gritti F., Guiochon G., J. Chromatogr. A, 2006, 1128, 45; doi:10.1016/j.chroma.2006.06.038. 209. Guiochon G., Beaver L.A., J. Chromatogr. A, 2011, 1218, 8836; doi:10.1016/j.chroma.2011.09.008. 210. Panagiotopoulos A.Z., Quirke N., Stapleton M., Tildesley D.J., Mol. Phys., 1988, 63, 527; doi:10.1080/00268978800100361. 211. Panagiotopoulos A.Z., Fluid Phase Equil., 1992, 76, 97; doi:10.1016/0378-3812(92)85080-R. 212. Luedtke W.D., Landman U., J. Phys. Chem. B, 1998, 102, 6566; doi:10.1021/jp981745i. 213. Baeurle S.A., J. Math. Chem., 2009, 46, 363; doi:10.1007/s10910-008-9467-3. 214. Rühle V., Junghans C., Lukyanov A., Kremer K., Andrienko D., J. Chem. Theory Comput., 2009, 5, 3211; doi:10.1021/ct900369w. 215. Yamamoto S., Maruyama Y., Hyodo S., J. Chem. Phys., 2002, 116, 5842; doi:10.1063/1.1456031. 12601-33 http://dx.doi.org/10.1021/ma800828c http://dx.doi.org/10.1080/15583720600945402 http://dx.doi.org/10.1016/S1573-4285(06)14002-8 http://dx.doi.org/10.1021/jp953301h http://dx.doi.org/10.1021/ac051085v http://dx.doi.org/10.1063/1.3414996 http://dx.doi.org/10.1103/PhysRevLett.102.115702 http://dx.doi.org/10.1002/marc.200900113 http://dx.doi.org/10.1016/j.bej.2009.12.003 http://dx.doi.org/10.1016/j.chroma.2005.03.124 http://dx.doi.org/10.1016/j.chroma.2011.02.012 http://dx.doi.org/10.1016/j.chroma.2008.07.037 http://dx.doi.org/10.1016/j.chroma.2008.12.088 http://dx.doi.org/10.1016/j.chroma.2011.10.043 http://dx.doi.org/10.1016/j.chroma.2011.11.039 http://dx.doi.org/10.1016/j.fluid.2009.10.010 http://dx.doi.org/10.1016/j.chroma.2006.06.038 http://dx.doi.org/10.1016/j.chroma.2011.09.008 http://dx.doi.org/10.1080/00268978800100361 http://dx.doi.org/10.1016/0378-3812(92)85080-R http://dx.doi.org/10.1021/jp981745i http://dx.doi.org/10.1007/s10910-008-9467-3 http://dx.doi.org/10.1021/ct900369w http://dx.doi.org/10.1063/1.1456031 S. Sokołowski, J. Ilnytskyi, O. Pizio Опис iнтерфейсiв мiж флюїдом та приєднаними полiмерними ланцюжками: прогрес у застосуваннi теорiй функцiоналу густини та позаґраткових комп’ютерних симуляцiй С. Соколовський1, Я. Iльницький2, О. Пiзiо3 1 Вiддiл моделювання фiзико-хiмiчних процесiв, унiверситет Марiї Кюрi-Склодовської, 20031 Люблiн, Польща 2 Iнститут фiзики конденсованих систем НАН України, вул. Свєнцiцького, 1, 79011, Львiв, Україна 3 Iнститут хiмiї, Нацiональний автономний унiверситет м. Мехiко, Мехiко, Мексика Iнтерфейси мiж флюїдом та приєднаними полiмерними ланцюжками зустрiчаються у численних нано- скопiчних об’єктах. Дослiдження їх властивостей є цiкавим як з точки зору фундаментальних наукових знань, так i для низки застосувань, зокрема для проектування нанопристроїв iз спецiальним призначе- нням. В оглядi подано сучасний стан досягнень теоретичних пiдходiв в цiй областi (насамперед методiв, що базуються на функцiоналi густини), якi просунули наше розумiння мiкроскопiчних властивостей та- ких iнтерфейсiв. Цi теоретичнi дослiдження дозволяють описати залежностi адсорбцiйних, змочувальних та сольватацiйних сил та особливостi електро-iнтерфейсних явищ вiд термодинамiчного стану та хара- ктеристик приєднаних ланцюжкiв. Зроблено також огляд комп’ютерних симуляцiй таких систем, якi, як i теоретичнi дослiдження, спiвставлено з наявними експериментальними даними. Ключовi слова: приєднанi ланцюжки, полiмернi щiтки, флюїд, адсорбцiя, змочування, розчини електролiтiв, диференцiальна ємнiсть 12601-34 Introduction Overview of theoretical methods Scaling approaches. A few comments Integral equations Self-consistent field approaches Density functional approaches. Generalities Density functional theory Applications of the density functional approaches Phase transitions Solvation force Electrostatic density functionals Selected applications of the electrostatic density functional approach Computer simulations Mesoscopic modelling Atomistic and semi-atomistic modelling Summary

Description of fluid-tethered chains interfaces: advances in density functional theories and off-lattice computer simulations

Репозитарії

Схожі ресурси