We study the conformational behavior of polarizable polymer chain under an external homogeneous electric field within the Flory type self-consistent field theory. We consider the influence of electric field on the polymer coil as well as on the polymer globule. We show that when the polymer chain conformation is a coil, application of external electric field leads to its additional swelling. However, when the polymer conformation is a globule, a sufficiently strong field can induce a globule-coil transition. We show that such “field-induced” globule-coil transition at the sufficiently small monomer polarizabilities goes quite smoothly. On the contrary, when the monomer polarizability exceeds a certain threshold value, the globule-coil transition occurs as a dramatic expansion in the regime of first-order phase transition. The developed theoretical model can be applied to predicting polymer globule density change under external electric field in order to provide more efficient processes of polymer functionalization, such as sorption, dyeing, and chemical modification

A coil-globule transition in dilute polymer solutions is one of the most famous phenomena in polymer physics. The importance of this phenomenon is explained by the fact that coil-globule transition plays a crucial role in recent technological advances from targeted delivery of drug molecules to their encapsulation in a polymer coil.1–3 On the other hand, understanding the conformational transitions of the single polymer chain in the solution is the basis of understanding the thermodynamic and structural properties of the so-called soft active materials, including dielectric elastomers,4 elastomeric gels,5 polyelectrolytes,6–17 copolymers,18,19 pH-sensitive hydrogels,20 and temperature-sensitive hydrogels.21 Today, the great efforts taken to develop a coil-globule transition theory have made a great contribution to understand this phenomenon. The coil-globule transition has been considered within a simple theory of self-consistent field as well as a sophisticated field-theoretic formalism.22–33 Most of these studies are based on the simple idea that the polymer coil shrinks lead to a collapse, when the solvent becomes poorer. It should be noted that the coil-globule transition usually occurs as first- or second-order phase transitions.34 

The possibilities to control the conformational transition by an external stimulus have attracted great attention in chemistry, biology, and material science. One of the most interesting external stimuli for different industrial applications enabling the control of the conformations of polymers, such as polyelectrolytes, dielectric elastomers, and liquid-crystalline polymers, is application of an external electric field. A lot of theoretical studies of the conformational behavior of single polyelectrolyte chain in the external homogeneous electric field have been published.35–41 In case of single polyelectrolyte chain, its 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 or a molecular polarizability of monomers that also make it possible to act on the polymer chain conformation by electric field application. Up to now, several theoretical works have been published where thermodynamic and structural properties of dielectric polymers in the bulk solution have been discussed.42–46 Nevertheless, to the best of our knowledge, until now there have been no the systematic theoretical studies of conformational behavior of a flexible dielectric polymer chain under the external electric field even at the level of mean-field approximation.

In this short communication, we undertake the first effort in this direction. We address the theoretical study of the polarizable polymer chain dissolved in a dielectric solvent under external homogeneous electric field within the Flory type self-consistent field theory.47 We study the influence of the external electric field on the polymer coil as well as on the polymer globule. We show that applying electric field leads to an additional swelling of the polymer coil. However, when the polymer chain conformation is a globule, we find out that a sufficiently large electric field can provoke globule-coil transition.

Let us consider a polarizable polymer chain dissolved in a solvent which we shall model as a continuous dielectric media with the dielectric permittivity εs. We assume that the polymer chain has a degree of polymerization N and that each monomer has a molecular polarizability γ. We also assume that the polymer chain is placed in a homogeneous electric field E. To study the conformations of the polymer chain in the external electric field, we shall use the simple Flory type47 self-consistent field theory, considering the radius of gyration Rg as a single order parameter. We also assume that the polymer chain occupies the volume which can be expressed through the gyration volume V g = 4 π R g 3 / 3 . Thus, a total free energy of the polymer chain can be written in the following form:

F ( R g ) = F conf ( R g ) + F vol ( R g ) + F el ( R g ) ,
(1)

where F conf ( R g ) is the free energy of the ideal polymer chain which can be calculated by the following interpolation formula:22,34,48

F conf ( R g ) = 9 4 k B T α 2 + 1 α 2 ,
(2)

where α = Rg/R0g is the expansion factor, R 0 g 2 = N b 2 / 6 is the mean-square radius of gyration of the ideal Gaussian polymer chain, b is the Kuhn length of the segment, kB is the Boltzmann constant, T is the absolute temperature; the contribution of monomer volume interactions to the total free energy takes the following form

F v o l ( R g ) = k B T N 2 B 2 V g + N 3 C 6 V g 2 ,
(3)

where B and C are the second and third virial coefficients of the monomer-monomer interactions, respectively; the electrostatic contribution can be written as

F el ( R g ) = V g ε p E 2 8 π ε s E 2 8 π ,
(4)

where εp is the effective dielectric permittivity inside the polymer volume, which can be determined in the mean-field approximation as49 

ε p = ε s + 4 π γ N V g ,
(5)

E = ε s E / ε p is the electric field inside the polymer volume, εs is the solvent dielectric permittivity, γ is the molecular polarizability of monomers, and E = |E| is the absolute value of the external electric field.

Using expression (5) and the expression for the local electric field E = ε s E / ε p inside the gyration volume, after some calculations, we obtain the following expression for the total free energy of the polymer chain in the solution under external electric field:

F k B T = 9 4 α 2 + 1 α 2 + N 2 B 2 V g + N 3 C 6 V g 2 N ε s γ E 2 2 k B T ε s + 4 π γ N V g .
(6)

Minimizing total free energy (6) with respect to the expansion factor α, after some algebra, we arrive at

α 5 α = 3 6 2 π b 3 B N + 27 C π 2 α 3 b 6 + 6 6 γ 2 ε s N E 2 k B T b 3 ε s + 18 6 γ α 3 b 3 N 2 .
(7)

As is seen from Eq. (7), the interaction of the induced dipoles of the polymer chain with the external electric field leads to the macromolecule swelling. The latter is well known in thermodynamics of dielectrics as an electrostriction effect.50 

Eq. (7) allows us to study the influence of the electric field on the conformation of a polarizable polymer chain. First, let us consider the case of good solvent when the polymer chain without the electric field is in coil conformation. In this case, B > 0 and the expansion factor α ≫ 1, so that the second term in the right-hand side of (7) (that is related to the triple monomer correlation) can be safely omitted. Then we get the following simplified equation for the expansion factor,

α 5 α = 3 6 2 π b 3 N B + 4 π γ 2 E 2 k B T ε s .
(8)

As one can see from (8), the presence of electric field inside the polymer coil only leads to an effective increase in the second virial coefficient and, consequently, its additional swelling. In case of strong electric field (E → ∞), we get the following evaluations for the expansion factor and the radius of gyration:

α 6 6 γ 2 k B T b 3 ε s 1 / 5 E 2 / 5 N 1 / 10 , R g b 2 γ 2 6 k B T ε s 1 / 5 E 2 / 5 N 3 / 5 .
(9)

It should be noted that in case of strong electric field, the size of the polymer coil is determined by the monomer polarizability, the temperature, the solvent dielectric permittivity, and the external electric field, but it does not depend on the parameters of monomer pair interactions.

Let us now consider a case of a poor solvent when the polymer chain without the electric field is in the globule conformation. In this regime, the second virial coefficient is a negative value (B < 0), while the expansion factor α ≪ 1. As is well known, in order to describe the conformational behavior of a polymer chain in the regime of poor solvent, we have to take into account the triple monomer correlations, leaving the contribution in (7) which is related to the third virial coefficient.34 In the case when E = 0, the expansion factor can be evaluated as α 3 6 C / π b 3 | B | 1 / 3 N 1 / 6 . Further, neglecting the terms in the left-hand side of (7) and introducing a notation α = sN−1/6, we obtain the following equation for parameter s:

B + 4 π γ 2 E 2 ε s k B T ε s + 18 6 γ s 3 b 3 2 + 3 6 C π s 3 b 3 = 0 .
(10)

At a sufficiently small electric field, equation (10) can be solved as

s 3 6 C π b 3 | B e f f ( E ) | 1 / 3 ,
(11)

where an effective second virial coefficient B eff ( E ) = B + 4 π γ 2 ε s E 2 / k B T ε s + 6 π γ | B | C 2 is introduced. As it follows from relation (11), as the electric field grows, the radius of gyration monotonically increases too. Solution (11) allows us to calculate the globule density,

ρ g = N V g 3 | B | 2 C 6 π γ 2 ε s E 2 C k B T ε s + 6 π γ | B | C 2 .
(12)

The first term in the right-hand side of Eq. (12) determines the globule density at zero electric field. The second term is the first correction on the electric field E. At a further electric field increase, the globule-coil transition may take place, which will be discussed detail in Subsection II C.

To proceed the numerical analysis of the above-mentioned in Sec. II B field-induced globule-coil transition, we introduce the following dimensionless variables: E ̃ = E ε s b 3 / k B T , B ̃ = B b 3 , C ̃ = C b 6 , and γ ̃ = γ b 3 / ε s . Thus, basic equation (7) can be rewritten as

α 5 α = 3 6 2 π N B ̃ + 4 π γ ̃ 2 E ̃ 2 1 + 18 6 γ ̃ α 3 N 2 + 27 C ̃ π 2 α 3 .
(13)

We choose parameters B ̃ and C ̃ so that polymer is in the globular conformation at zero electric field, namely, B ̃ = 0 . 25 and C ̃ = 0 . 4 .

Fig. 1 demonstrates the dependencies of expansion factor α on the external electric field E ̃ at different monomer polarizabilities γ ̃ , obtained by solving Eq. (13). As is seen, at the sufficiently small monomer polarizabilities, the electric field increase causes the expansion factor to grow smoothly. It should be noted that in the limit of the infinitely long polymer chain (N → ∞) at the sufficiently small monomer polarizability, the globule-coil transition occurs as a second-order phase transition. However, above some critical value γ ̃ c of the monomer polarizability, the expansion factor dependence on the electric field looks like an “S-shaped” curve. The latter means that the globule-coil transition at large monomer polarizabilities occurs as a first-order phase transition. The regions of metastable states and of unstable states are marked by the red and blue lines, respectively. It should be noted that in reality at certain electric field, an abrupt increase (black vertical lines) in the expansion factor should take place. We would also like to stress that the expansion factor values between which the abrupt increases occur were found to be the points of global minima of the total free energy. Fig. 1 also illustrates the binodal with the critical point (green line).

FIG. 1.

The phase diagram is presented in coordinates α vs E ̃ for the polarizable polymer chain. At the sufficiently small monomer polarizabilities, the electric field increase causes the expansion factor to grow smoothly. Above some critical value γ ̃ c of the monomer polarizability, the expansion factor dependence on the electric field looks like an “S-shaped” curve that indicates on the first-order phase transition. The regions of metastable states (blue lines) and of unstable states (red lines) are marked. The binodal is marked by the green line. The data are shown for B ̃ = 0 . 25 , C ̃ = 0 . 4 , N = 100.

FIG. 1.

The phase diagram is presented in coordinates α vs E ̃ for the polarizable polymer chain. At the sufficiently small monomer polarizabilities, the electric field increase causes the expansion factor to grow smoothly. Above some critical value γ ̃ c of the monomer polarizability, the expansion factor dependence on the electric field looks like an “S-shaped” curve that indicates on the first-order phase transition. The regions of metastable states (blue lines) and of unstable states (red lines) are marked. The binodal is marked by the green line. The data are shown for B ̃ = 0 . 25 , C ̃ = 0 . 4 , N = 100.

Close modal

In this communication, we have presented a simple self-consistent field theory of the conformational behavior of polarizable polymer chain in the dilute polymer solution under external homogeneous electric field. We have shown that in the regime of good solvent, when the polymer chain adopts a coil conformation, application of the electric field leads to additional swelling of the polymer coil. However, in the regime of poor solvent, when the polymer chain conformation is a globule, we have obtained that the electric field can induce a globule-coil transition in two qualitatively different ways. Namely, at the sufficiently small monomer polarizability, the globule-coil transition occurs quite smoothly. On the contrary, when the monomer polarizability exceeds a certain value, the globule-coil transition occurs as a first-order phase transition, i.e., as a discontinuous change of the radius of gyration.

We would like to stress that within the present model, we assume the isotropic dielectric response of the polymer chain. Such simplification is due to the fact that in this study, we would like to discuss only the principal possibility to change polarizable polymer chain conformation by electric field at simplest level of theory, leaving more rigorous consideration with an account of the response anisotropy for future research. Nevertheless, we believe that accounting for the dielectric response anisotropy could give some new interesting regimes.

In conclusion, we would like to speculate on the possible applications of the discussed phenomena. First, we believe that they might offer an additional possibility to control polymer conformation in the solutions for encapsulation of drug molecules into a polymer coil and their subsequent targeted delivery. Second, application of electric field, decreasing the polymer density, can help to provide more efficient processes of polymer functionalization, such as sorption, dyeing, and chemical modification.

This work was supported by grant from the President of the Russian Federation (No. MK-2823.2015.3).

1.
J.
Kost
and
R.
Langer
,
Adv. Drug Delivery Rev.
46
,
125
(
2001
).
2.
P.
Bawa
,
V.
Pillay
,
Y. E.
Choonara
, and
L. C.
du Toit
,
Biomed. Mater.
4
,
022001
(
2009
).
3.
S. D.
Fitzpatrick
,
E. F.
Lindsay
,
A.
Thakur
 et al,
Expert Rev. Med. Devices
9
(
4
),
339
(
2012
).
4.
Z.
Suo
,
Acta Mech. Solida Sin.
23
(
6
),
549
(
2010
).
5.
M.
Doi
,
J. Phys. Soc. Jpn.
78
,
052001
(
2009
).
6.
A. V.
Dobrynin
and
M.
Rubinstein
,
Prog. Polym. Sci.
30
,
1049
(
2005
).
7.
A. V.
Dobrynin
,
Curr. Opin. Colloid Interface Sci.
13
,
376
(
2008
).
8.
C.
Holm
,
J. F.
Joanny
,
K.
Kremer
,
R. R.
Netz
,
P.
Reineker
,
C.
Seidel
,
T. A.
Vilgis
, and
R. G.
Winkler
,
Adv. Polym. Sci.
66
,
67
(
2004
).
9.
V. Yu.
Boryu
and
I. Ya.
Erukhimovich
,
Macromolecules
21
(
11
),
3240
(
1988
).
10.
S. A.
Baeurle
,
M. G.
Kiselev
,
E. S.
Makarova
, and
E. A.
Nogovitsin
,
Polymer
50
(
7
),
1805
(
2009
).
11.
E. A.
Nogovitsyn
and
Yu. A.
Budkov
,
Russ. J. Phys. Chem. A
85
,
1363
1368
(
2011
).
12.
E. A.
Nogovitsin
and
Yu. A.
Budkov
,
Physica A
391
,
2507
(
2012
).
13.
Yu. A.
Budkov
,
E. A.
Nogovitstyn
, and
M. G.
Kiselev
,
Russ. J. Phys. Chem. A
87
,
638644
(
2013
).
14.
Yu. A.
Budkov
,
A. L.
Kolesnikov
,
E. A.
Nogovitsyn
, and
M. G.
Kiselev
,
Polym. Sci., Ser. A
56
(
5
),
697
(
2014
).
15.
Yu. A.
Budkov
,
A. L.
Kolesnikov
,
N.
Georgi
,
E. A.
Nogovitsyn
, and
M. G.
Kiselev
,
J. Chem. Phys.
142
,
174901
(
2015
).
16.
Yu. A.
Budkov
,
A. L.
Kolesnikov
,
N.
Georgi
,
E. A.
Nogovitsyn
, and
M. G.
Kiselev
,
J. Chem. Phys.
143
,
189903
(
2015
).
17.
A. L.
Kolesnikov
,
Yu. A.
Budkov
, and
E. A.
Nogovitsyn
,
J. Phys. Chem. B
118
,
13037
(
2014
).
18.
T.
Goldacker
,
V.
Abetz
,
R.
Stadler
,
I.
Erukhimovich
, and
L.
Leibler
,
Nature
398
,
137
(
1999
).
19.
L.
Wang
,
J.
Lin
, and
X.
Zhang
,
Polymer
54
,
3427
(
2013
).
20.
R.
Marcombe
,
S. Q.
Cai
,
W.
Hong
,
X. H.
Zhao
,
Y.
Lapusta
, and
Z. G.
Suo
,
Soft Matter
6
,
784
(
2010
).
21.
S. Q.
Cai
and
Z. G.
Suo
,
J. Mech. Phys. Solids
59
(
11
),
2259
(
2011
).
22.
A. Yu.
Grosberg
and
D. V.
Kuznetsov
,
Macromolecules
25
,
1970
(
1992
).
23.
T. M.
Birshtein
and
V. A.
Pryamitsyn
,
Macromolecules
24
,
1554
(
1991
).
24.
M. A.
Moore
,
J. Phys. A: Math. Gen.
10
(
2
),
305
(
1977
).
25.
P. G.
de Gennes
,
J. Phys. Lett.
36
(
2
),
55
(
1975
).
26.
I. C.
Sanchez
,
Macromolecules
12
(
5
),
980
(
1979
).
27.
28.
N. V.
Brilliantov
,
D. V.
Kuznetsov
, and
R.
Klein
,
Phys. Rev. Lett.
81
(
7
),
1433
(
1988
).
29.
A.
Halperin
and
E. B.
Zhulina
,
Europhys. Lett.
15
(
4
),
417
(
1991
).
30.
Yu. A.
Budkov
,
A. L.
Kolesnikov
,
N.
Georgi
, and
M. G.
Kiselev
,
J. Chem. Phys.
141
,
014902
(
2014
).
31.
Yu. A.
Budkov
,
I. I.
Vyalov
,
A. L.
Kolesnikov
 et al,
J. Chem. Phys.
141
,
204904
(
2014
).
32.
Yu. A.
Budkov
,
A. L.
Kolesnikov
,
N.
Georgi
, and
M. G.
Kiselev
,
Europhys. Lett.
109
,
36005
(
2015
).
33.
K.
Odagiri
and
K.
Seki
,
J. Chem. Phys.
143
,
134903
(
2015
).
34.
A. Yu.
Grosberg
and
A. R.
Khokhlov
,
Statistical Physics of Macromolecules
(
AIP
,
New York
,
1994
).
35.
M.
Muthukumar
,
J. Chem. Phys.
86
,
7230
(
1987
).
36.
R.
Podgornik
and
B.
Jonsson
,
Europhys. Lett.
24
,
501
(
1993
).
37.
O. V.
Borisov
,
E. B.
Zhulina
, and
T. M.
Birshtein
,
J. Phys. II
4
,
913
(
1994
).
38.
X.
Chatellier
and
J.-F.
Joanny
,
Phys. Rev. E
57
,
6923
(
1998
).
39.
R. R.
Netz
,
J. Phys. Chem. B
107
,
8208
(
2003
).
40.
N. V.
Brilliantov
and
C.
Seidel
,
Europhys. Lett.
97
,
28006
(
2012
).
41.
C.
Seidel
,
Yu. A.
Budkov
, and
N. V.
Brilliantov
,
Proc. Inst. Mech. Eng., Part N
227
,
142
(
2013
).
42.
M.
Muthukumar
,
J. Chem. Phys.
104
,
691
(
1996
).
43.
R.
Podgornik
,
Phys. Rev. E
70
,
031801
(
2004
).
44.
R.
Kumar
and
G. H.
Fredrickson
,
J. Chem. Phys.
131
,
104901
(
2009
).
45.
D. S.
Dean
and
R.
Podgornik
,
J. Chem. Phys.
136
,
154905
(
2012
).
46.
R.
Kumar
,
B. G.
Sumpter
, and
M.
Muthukumar
,
Macromolecules
47
,
6491
(
2014
).
47.
P.
Flory
,
Statistical Mechanics of Chain Molecules
(
Wiley-Interscience
,
New York
,
1969
).
48.
M.
Fixman
,
J. Chem. Phys.
36
(
2
),
306
(
1962
).
49.
Yu. A.
Budkov
,
A. L.
Kolesnikov
, and
M. G.
Kiselev
,
Europhys. Lett.
111
,
28002
(
2015
).
50.
L. D.
Landau
and
E. M.
Lifshitz
,
Electrodynamics of Continuous Media
,
A Course of Theoretical Physics
Vol.
8
(
Pergamon Press
,
Oxford
,
1960
).