Statistical theory of polarizable target compound impregnation into a polymer coil under the influence of an electric field

The paper presents a theoretical approach for describing the influence of an electric field on the conformation of an electrically neutral dielectric polymer chain dissolved in a dielectric solvent with an admixture of a target compound. Each monomer and each molecule of the target compound carries positive excess polarizability and the solvent is described as a continuous dielectric medium. The model is based on the Flory-type mean-field theory. We demonstrate non-monotonic dependences of the expansion factor and the concentration of the target compound on the strength of the electric field and molecular polarizability. Namely, the target compound concentration in the internal polymer volume as a function of electric field strength has pronounced maxima if the molecules are polarizable. In addition, the expansion factor of the non-polarizable polymer chain can be controlled by the electric field. The dependences of the expansion factor and target compound concentration on the monomer polarizability exhibit minima and intersection points. The intersection points correspond to the equality of dielectric permittivities in the bulk solution and in the internal polymer volume.


Introduction
It is well known that there are numerous external stimuli such as temperature, pressure, pH, composition, etc., which may affect the thermodynamic and structural properties of polymer solutions.These external stimuli are particularly important for the description of the so-called 'smart' or 'environmentally sensitive' polymers. 1,2These polymers are called smart because their properties can be changed by minor changes in the polymer surroundings.The abovementioned stimuli can produce changes in one or several properties such as shape, volume, solubility, surface characteristics, etc. 3 Moreover, the reversibility of these changes is the peculiarity of 'smart' materials.An electric field is one of the stimuli, despite its nature, which can be both artificial (a pair of electrodes) and natural (electric fields of biological membranes, for instance).It is also one of the possible stimuli for a controlled drug release system. 2,4The mechanism is based on the motion of electrically charged molecules under the influence of an electric field.The electric field affects both particles carrying a net electric charge and dielectric particles carrying permanent or induced dipole moments.The application of an electric field to dielectrics may lead to different effects, for instance, electrostriction and dielectrophoresis (electrostriction is the deformation of a dielectric medium under an electric field, and dielectrophoresis is the motion of particles in an inhomogeneous electric field 5,6 ).
The theoretical physics of dielectric polymers still remains a terra incognita.Indeed, there are several works devoted to the thermodynamic properties of dielectric polymer solutions and polymer brushes.Kumar et al. described the influence of monomer permanent dipole moments on the polymer mixture phase diagram shape. 7Dean and Podgornik formulated the theory of anisotropic van der Waals interactions of semi-flexible polymers and their influence on ordering. 8Podgornik also presented a theory describing the influence of electrostatic interactions of permanent dipoles on persistence length. 9Lu et al. analyzed static van der Waals interactions between two polarizable rigid polymers within the field-theoretic approach. 10n ref. 11, the authors formulated a statistical field theory of dielectric soft matter.In ref. 12, the authors developed a selfconsistent field theory of polar polymer brushes and analyzed its phase behavior depending on dipolar monomer interactions.
Theoretically, the electric field influence on various polymer properties has been studied mostly for polyelectrolyte solutions [13][14][15][16][17][18][19][20][21] and polyelectrolyte-based electrically responsive hydrogels as delivery systems (see, for instance, ref. 22-24).[25][26][27] Their response to the electric field is achieved by the presence of an electric charge on the polymer backbone.However, there is a wide class of neutral dielectric polymers with permanent dipoles and/or molecular polarizability of monomers that are exposed to an electric field.Nevertheless, dielectric polymeric materials under an electric field have theoretically been studied even much less extensively than the mentioned bulky dielectric polymers.Shkel and Klingenberg formulated generic models of solid material electrostriction. 28Suo suggested theoretical models of dielectric elastomers which are strongly deformable under an applied electric field. 29There are some works on the electrostriction of polymeric gels 30,31 and copolymer melts. 32,33ecent papers 34,35 have shown that the application of a homogeneous constant electric field leads to dielectric polymer expansion, regardless of the conformation (coil or globule).Moreover, these studies have investigated the influence of the many-body electrostatic dipole correlations of monomers on the conformational behavior of polymer chains. 35,36These papers have demonstrated that when the excess polarizability is positive, the application of an electric field to a polymer solution always leads to electrostriction swelling of the polymer chain.In recent papers, 37,38 the authors theoretically studied the influence of co-solvent polarizability on the differential capacitance of electric double layers.Using a local density functional approach, a modified Poisson-Boltzmann theory with an explicit account of polarizable co-solvent and excluded volume of the ions and the co-solvent was formulated.
The polymer coil swelling under the constant homogeneous electric field makes it reasonable to question: how will electrostriction affect the impregnation of some polarizable target compound (TC) inside a polarizable polymer coil?To the best of our knowledge, dilute solutions of dielectric polymers in the presence of the dielectric TC under an electric field have not been investigated systematically even at the level of the mean-field theory.A theoretical description of the TC concentration inside the polymer volume and polymer conformational behavior under an electric field is the main focus of the paper.
The paper is organized as follows: Section 2 presents the model and Section 3 the numerical results and discussion, followed by conclusions in Section 4.

Model
Now we will consider a simple model of a dilute polymer solution with an admixture of polarizable TC molecules.The solvent is treated as a dielectric medium with constant dielectric permittivity e s .As in our previous papers, [39][40][41][42] the whole volume of the solution is divided into two parts: the volume of the polymer chain (assessed by the polymer chain gyration volume) and the volume of the bulk solution.The impregnation process is described as a thermodynamic equilibrium between the bulk volume and the gyration volume of the polymer chain in terms of target compound concentration and polymer volume.The aim of this paper is to consider such a system under a static electric field taking into account the static positive excess polarizability of monomers and TC molecules.To reach our goal, we simplify the model, which assumes a linear and isotropic response of all molecules in the solution to the electric field.The latter assumption allows us to consider the grand thermodynamic potential for the gyration volume in the simplest form.Moreover, the effects discussed in this paper in general are independent of the polymer form and determined by the difference in dielectric permittivities in the gyration volume and bulk solution.We would like to stress that only our assumption of the isotropic dielectric response of the solution species makes it possible to consider the theory in terms of pressure instead of a stress tensor. 43

Model: bulk volume
We start from the free energy per particle of the TC in the bulk solution, which can be written as follows: where the concentration r of TC molecules and the temperature T are determined.The first contribution on the right hand side of ( 1) is the ideal gas free energy of the target compound, the second one is the interaction energy between the particles at the level of the second virial term and the third one is the electrostatic energy of the solution.Also, k B is the Boltzmann constant, L is the thermal de Broglie wavelength, w c is the second virial coefficient of the effective interactions between the TC molecules in the solvent medium, E is the strength of the electric field in the pure solvent, is the strength of the electric field in the solution (see below), e s is the dielectric permittivity of the solvent medium, and e is the dielectric permittivity in the bulk solution: 37,38 e = e s + 4pg c r, where the polarizability of the TC (g c ) is introduced.The third term on the right hand side of (1) describes the excess electrostatic free energy of the TC relative to the pure solvent.The electric field strength in the solution can be evaluated as follows: 43 Using these notations, the chemical potential and partial pressure of the TC molecules take the following form: and (5)

Model: internal polymer volume
Now let us discuss the model of a mixture in the gyration volume.We choose two independent order parameters: the gyration radius (R g ) of the polymer chain and the average number of the TC molecules c ) in the gyration volume.
0][41][42] The volume interactions between the molecules in the solution are described by the second virial coefficients (monomermonomer, TC-TC, monomer-TC interactions).Thus, bearing in mind all the above factors, one can construct the Gibbs thermodynamic potential as follows: Here is the ideal part of the free energy, where the first term is the conformation free energy of the Gaussian polymer chain within Fixman's approximation and the second term is the free energy of the ideal gas confined in the volume V g ; the contribution of the effective interactions between the species can be written as follows: 47 where N is the degree of polymerization, and w p , w c , and w pc are the second virial coefficients for monomer-monomer, TC-TC and monomer-TC interactions, respectively.It should be noted that in our consideration we assume that the degrees of freedom of the solvent molecules are integrated out, so the above mentioned second virial coefficients are renormalized with respect to the ''bare'' virial coefficients by the presence of the solvent.In other words, in accordance with the classical McMillan-Mayer theory, 48 all the considered second virial coefficients are related to the potentials of mean-force of the interaction between the species in the solvent medium.Thus, following the logic mentioned above (McMillan-Mayer theory), one can construct the solvation free energy of the polymer using partial pressure instead the total pressure.The total electrostatic free energy within the mean-field approximation consists of two terms: where F ex,p is the excess free energy of the mixture in the gyration volume relative to the bulk solution, and F ex,c is the excess free energy of the bulk solution in the gyration volume relative to the pure solvent.These contributions can be assessed within the mean-field approximation as follows: 49 and We would like to emphasize that the polymer is modeled as a polarizable sphere with the dielectric permittivity e p immersed in the dielectric medium with the dielectric permittivity e.Thus, the total electrostatic free energy is calculated as the excess electrostatic energy of the monomers together with the TC molecules contained in the gyration volume relative to the pure solvent occupying the same volume.m and P are, respectively, the chemical potential and partial pressure of the TC in the bulk solution satisfying expressions (4) and (5).Now, let us clarify the symbolic notations: R 0g 2 = Nb 2 /6 is the mean-square radius of gyration of the Gaussian polymer chain with the Kuhn length b, and E is the strength of the electric field.The gyration volume of the polymer chain is expressed through the gyration radius as: The dielectric permittivity of the solution in the gyration volume is expressed as: 37 where g c is the polarizability of one TC molecule, g p is the polarizability of monomers, r c is the TC concentration in the gyration volume, and r p is the concentration of monomers in the gyration volume.Now minimizing eqn (6) relative to the expansion factor of the polymer chain and the concentration of TC molecules in the gyration volume, one can obtain the following system of two coupled equations: where the expansion factor a = R g /R 0g is introduced.The latter system describes the behavior of the expansion factor a and concentration r c of the TC molecules near the polymer backbone depending on the electric field strength E, polarizabilities of the solution species g p and g c , second virial coefficients w p , w c and w pc , and the index of polymerization N. In eqn (15), the difference DP between the partial TC pressure inside the gyration volume which is expressed as: and the bulk pressure ( 5) is determined.The pressure difference DP consists of three parts: the first is the ideal contribution, the second is the contribution of the interactions between the target compound and the third is the electrostatic contribution.The exact form of the electrostatic contribution depends on the geometrical description of the polymer chain, while the first two terms hold the standard form of the virial expression. 47When the dielectric permittivity in the gyration volume is equal to the bulk permittivity e p = e, the electrostatic contributions in eqn ( 14) and ( 15) vanish.Let us now consider the limiting case of the expanded coil, i.e., when a c 1 and (r c À r)/r { 1.In the present limit, we obtain the following equation: where Aðr; EÞ e s 2 g p g c E 2 is the effective second virial coefficient of monomers in the solution under the electric field and the auxiliary function is also introduced.The results recently obtained by Budkov et al. 34,35 can be derived from eqn (17) by deleting the terms produced by the presence of TC molecules in the solution, i.e., when r = 0. Thus, in the absence of TC molecules in the solution, monomer polarizability is the only driving force for polymer chain swelling in the electric field.As we shall show below, an admixture of the TC molecules leads to non-monotonic dependences of both the expansion factor and TC concentration on the electric field and polarizability of the TC molecules.The TC concentration rc in the gyration volume and the expansion factor a of the polymer coil as functions of electric field strength are shown in Fig. 1 and 2. These calculations were done for the following set of TC polarizabilities gc = 0, gc = 0.2, and gc = 0.4 and fixed monomer polarizability gp = 0.8.As mentioned above, in the absence of an electric field the polymer coil exists in a collapsed conformation, whereas its application expands the coil.Such behavior is in agreement with our previous results. 34,35If the TC molecules are non-polarizable, then r c tends to the bulk value as the electric field strength increases.However, in the case of polarizable TC molecules, the picture changes.Namely, with an increase in the polarizability of TC molecules, their concentration in the gyration volume  monotonically decreases as the electric field strength increases.Moreover, at the non-zero polarizability of the TC molecules, the dependence of the TC concentration on the electric field strength has pronounced maxima in the vicinity of E E 2. These extrema are caused by introducing TC polarizability into the model.Now let us consider the behavior of r with the increase in the degree of polymerization N, analyzing eqn (14).It is evident that with N -N the second term in the exponent of eqn ( 14) vanishes.Thus, in the limit of infinite chain solution, eqn (14) results in r c = r.To illustrate that, Fig. 3 presents the dependence of TC concentration r c on the electric field for N = 100, 600, and 1000.As one can see, with the increase in the polymerization index the abovementioned nonmonotonic behavior of TC concentration disappears.

Regime II: cp o cc
Firstly, let us look at the dependence of a and rc on polymer polarizability at fixed electric field strength and other parameters.Fig. 4 and 5 show the concentration of TC molecules rc and expansion factor a as functions of gp monomer polarizability at different electric field strengths with fixed TC polarizability gc = 0.5.
Both functions have minima at non-zero monomeric polarizability.There is another interesting feature -these functions have intersection points that correspond to the equality of the dielectric permittivities in the bulk solution and the gyration volume.In order to prove this statement, we put e = e p .This substitution slightly changes the system of coupled equations.The only differences are: in eqn ( 14) the third term in the exponent is absent and in eqn (15) the third term is modified with DP = P c À P = k B T(r c À r + 1/2w c (r c 2 À r 2 )).Following these changes, let us write down: Using the new system, we obtain numerical values of gp E 0.09, a E 1.31 and rc E 0.12 which correspond to the intersections on the graphs.So, eqn (20) and ( 21) describe the case when neither the expansion factor nor the TC concentration depends on the electric field strength.Now let us discuss the case of a non-polarizable polymer chain, i.e. gp = 0, as opposed to TC molecules.Thus, it can be the case of a highly 'solvophobic' polymer relative to the TC molecules.The virial coefficients w ˜p = 0.5, w ˜c = 0.5, and w ˜pc = 3 are chosen in order to make this regime more pronounced.In Fig. 6, we present the behavior of the expansion factor with different values of TC polarizability, namely, gc = 0.2, gc = 0.4 and gc = 0.6.Also, in Fig. 7, the TC concentration behavior   inside the gyration volume is presented at the same TC polarizabilities.
Under low electric fields, the expansion factor remains nearly constant and then begins to increase monotonically.This regime can be described as follows.Due to the high 'solvophobicity' of the polymer chain the concentration of TC molecules inside the gyration volume is much lower compared to the value, which leads to polymeric coil compression.The electric field application tends to increase the TC concentration in the gyration volume, which, in turn, expands the polymer coil.It is worth noting that the TC molecule concentration in the gyration volume is always lower than the bulk value even at zero polymer solvophobicity (w pc = 0).That is due to the fact that the electric field in the gyration volume is weaker than its bulk value, so the TC molecules -carrying induced dipoles -prefer to be in the bulk solution rather than in the gyration volume.

Conclusions
In this work we have extended the theoretical description of internal and external stimuli for a polymer coil as a 'smart' polymer.We have described the impregnation of a polarizable solvophobic polymer coil with an additional amount of the polarizable target compound.
We tested two regimes: gp o gc and gp 4 gc .The first regime has shown (1) a monotonic increase in the expansion factor of the polymer coil with the increase of the electric field strength and (2) a maximum for the dependence of the target compound concentration in the gyration volume on the electric field.This maximum becomes more pronounced with the target compound polarizability growth.In the second regime, non-monotonic dependences of the expansion factor and the TC concentration on the polymer polarizability were demonstrated.
We also stress that the results described above were provided with the assumption of an isotropic dielectric response to the electric field, which allowed us to use the simple level of theoretical description.We assume that the effects related to the change in the coil shape should be important at a high strength of the electric field.In any case, our aim was to show the very possibility of controlling the polymer and the target compound concentration inside the polymer volume.We believe that all this information can contribute to a better understanding of the dielectric polymer behavior under an electric field.
We would like to stress that in this study we have neglected the contribution of the dipole correlations (for details see ref. 35) of the solution species.It can be shown that this contribution cannot qualitatively change the final outcomes of the present pure mean-field theory.The suggested theory also has an obvious limitation.Namely, it can be applied to polymer chains whose monomers are isotropically polarizable.In other words, in order to apply this theory, monomers and target compound molecules must have small non-diagonal components of polarizability tensors.The second limitation of the present theory is that it cannot describe the change in the polymer shape under an electric field.Indeed, we have introduced only one geometric parameter -radius of gyration -that makes it possible to control changes in the polymer volume as a whole, but it does not allow us to describe the polymer coil shape change.Nevertheless, in our opinion, this limitation is not crucial in the case, when small enough polymer coil deformations or electric field strengths (compared to the internal molecular electric field) are considered.The latter assumption is always satisfied under the experimentally accessible conditions.

For
numerical calculations the following set of dimensionless parameters has been introduced: w ˜p = w p /b 3 , w ˜c = w c /b 3 , w ˜pc = w pc /b 3 , gp = g p /(b 3 e s ), gc = g c /(b 3 e s ), and Ẽ ¼ E ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi e s b 3 =k B T p .The calculations have been made for the polymer chain with repulsive interactions between the monomers and TC molecules, that is to say w ˜pc 4 0. The numerical results are divided into two parts: gp 4 gc and gp o gc .3.1 Regime I: cp 4 cc The expansion factor a and the dimensionless concentration of TC molecules rc = r c b À3 in the gyration volume were plotted as functions of the electric field strength to demonstrate the influence of TC polarizability on them.The other parameters were fixed: N = 10 2 , w ˜p = 0.5, w ˜c = 0.5, w ˜pc = 1.5, and rc = 0.15.This combination of virial coefficients results in a collapsed conformation of the polymer coil at E -0 due to the polymer ''solvophobicity'' relative to the TC molecules.

Fig. 1
Fig. 1 Dependence of target compound concentration rc = r c b À3 on electric field strength Ẽ ¼ E ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi e s b 3 =k B T p in the gyration volume at constant polymer polarizability gp = g p /(b 3 e s ) = 0.8 and different co-solvent polarizabilities.

Fig. 2
Fig. 2 Dependence of expansion parameter a on electric field strength Ẽ ¼ E ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi e s b 3 =k B T p of the polymer coil at constant polymer polarizability gp = g p /(b 3 e s ) = 0.8 and different co-solvent polarizabilities.

Fig. 3
Fig. 3 Dependence of target compound concentration rc = rb À3 on the electric field strength Ẽ ¼ E ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi e s b 3 =k B T p in the gyration volume at constant polymer gp = g p /(b 3 e s ) = 0.8 and target compound gc = g c /(b 3 e s ) = 0.5 polarizabilities.

Fig. 4
Fig. 4 Expansion factor a as a function of monomer polarizability gp = g p / (b 3 e s ) for the solvophobic polymer at fixed values of electric field strength; target compound polarizability gc = g c /(b 3 e s ) = 0.5.

Fig. 5
Fig. 5 Concentration of target compound molecules in the gyration volume rc as a function of monomer polarizability gp = g p /(b 3 e s ) for the solvophobic polymer at fixed values of electric field strength; target compound polarizability gc = g c /(b 3 e s ) = 0.5.

Fig. 6
Fig.6Behavior of the solvophobic polymer with zero polarizability (g p = 0) in the presence of the polarizable target compound under an electric field.The polarizabilities of target compound molecules are taken as gc = g c /(b 3 e s ) = 0.2, 0.4, and 0.6.

Fig. 7
Fig. 7 Behavior of TC concentration inside the gyration volume of the solvophobic polymer with zero polarizability (g p = 0) under an electric field.The polarizabilities of target compound molecules are taken as gc = g c / (b 3 e s ) = 0.2, 0.4, and 0.6.