Elastomers are used as dielectric layers contained between the parallel conductive plates of capacitors. The introduction of filler particles into an elastomer changes its permittivity ε. When particle organization in a composite is intentionally varied, this alters its capacitance. Using numerical simulations, we examine how conductive particle chains introduced into polydimethylsiloxane (PDMS) alter ε. The effects of filler volume fraction ψ, interparticle d and interchain spacing a, zigzag angle θ between adjacent particles and overall chain orientation, particle size r, and clearance h between particles and the conductive plates are characterized. When filler particles are organized into chainlike structures rather than being just randomly distributed in the elastomer matrix, ε increases by as much as 85%. When particles are organized into chainlike forms, ε increases with increasing ψ and a, but decreases with increasing d and θ. A composite containing smaller particles has a higher ε when ψ<9% while larger particles provide greater enhancement when ψ is larger than that value. To enhance ε, adjacent particles must be interconnected and the overall chain direction should be oriented perpendicular to the conductive plates. These results are useful for additive manufacturing on electrical applications of elastomeric composites.

Due to their low cost, lighter weight, and superior hydrophobicity, electrically insulating rubber-like materials, such as silicone rubber, are often used as the dielectric medium in capacitors. A dielectric medium is characterized by its permittivity, which is measured using a capacitor. At steady state, the relative permittivity, i.e., the ratio of the capacitance for a dielectric medium to that for a vacuum,1 is also called the static relative permittivity.2 The properties of its individual constituents influence the bulk properties of a composite material. The addition of conductive filler particles into an elastomer matrix improves the bulk relative permittivity of the composite, predictions of which can be based on mixing rules or effective medium theory.3–5 However, since predictive methods only consider the volume fractions of the constituents, the geometrical arrangement of the filler particles, which also influences the bulk properties of the composite,6,7 is not appropriately represented.

There are several means to organize the filler, e.g., using optical tweezers,8 electron beams,9 or field-assisted self-assemblies.10–12 The former two methods require expensive tools, whereas the latter method uses simpler tools,13 such as application of an external magnetic field on a liquid polymer that is infused with suspended particles in ferrofluid.10,11,14–17 This creates a basis for a ‘material printer’. For instance, starting with a liquid polymer-particle mixture, desired bulk material attributes can be imparted to a material by a priori specifying the required particle arrangements. Enhancing the electric permittivity by organizing the filler particles in an elastomeric composite is expected to be beneficial for flexible electronics,18 3D printed circuit boards19,20 and sensors,21,22 data storage,23 optics,24 and artificial muscles.25 

The simplest and most common particle organization using field-assisted control is a chainlike structure.12,16,26 Upon application of an external field, filler particles within the liquid prepolymer align along the field direction to form straight12 or zigzag chains.12,27–29 Therefore, we arrange conductive particles in chainlike structures that lie perpendicular to the conductive plates within polydimethylsiloxane (PDMS), as shown in Figure 1. Numerical simulations are employed to investigate the effects of different particle arrangements and sizes on the static relative permittivity of the composite using the COMSOL Multiphysics® package along with its electrostatics physics module. Our simulations assume uniform spherical particles, a steady electrostatic condition, and the absence of defect, fatigue and hysteresis. Despite these assumptions, the results are useful for designing and understanding how the electric permittivity of polymer-particle composites can be improved.

FIG. 1.

Model of a parallel capacitor consisting of two circular plates of radius R placed at a distance L apart. The capacitor is filled with a dielectric composite consisting of polydimethylsiloxane (PDMS) and spherical conductive particles. For simplicity, particle organization is limited to a chainlike structure characterized by interparticle spacing d, interchain spacing a, particle radius r, clearance h between the outermost particle of a chain and its adjacent plate, zigzag angle θ, and volume fraction of the particles θ.

FIG. 1.

Model of a parallel capacitor consisting of two circular plates of radius R placed at a distance L apart. The capacitor is filled with a dielectric composite consisting of polydimethylsiloxane (PDMS) and spherical conductive particles. For simplicity, particle organization is limited to a chainlike structure characterized by interparticle spacing d, interchain spacing a, particle radius r, clearance h between the outermost particle of a chain and its adjacent plate, zigzag angle θ, and volume fraction of the particles θ.

Close modal

The relative permittivity of a composite (ε) is predicted by simulating a parallel capacitor at steady state without any current flow,

(1)

where V denotes the applied voltage, υρυdυ total charge, C overall capacitance, ε0 absolute permittivity of vacuum, A projected surface area of the plates, and L distance between the plates. The PDMS matrix has relative permittivity of 2.69.30 Since the relative permittivity of metal is infinite,31,32 we assume a large value, 2.8×105, for the conductive particles as an appropriate idealization. Our models are verified for pure PDMS and for PDMS-particle composites with randomly distributed particles using various approaches, as explained in the supplementary material. We also ensure that the applied voltages do not lead to electric field strengths inside the composites that exceed the breakdown strength for the PDMS matrix.33 When particle chains are oriented parallel to the plates, ε lies within 3% of that obtained with randomly distributed particles. Therefore, only the results for chains oriented perpendicular to the plates are discussed.

We examine the influence of the following parameters on ε: (1) filler particle volume fraction ψ, (2) interparticle spacing d, (3) interchain spacing a, (4) clearance h between the outermost particle and its adjacent conductive plate, (5) zigzag angle θ (the acute angle between the center-to-center direction of adjacent particles and the overall chain direction), and (6) size of the filler particles r, all of which are illustrated in the schematic of Figure 1.

Organizing filler particles into a chainlike structure within a dielectric composite enhances ε. For the same ψ, this ε enhancement diminishes rapidly as the interparticle spacing d/2r increases and ε approaches the corresponding value for a composite that contains only randomly distributed particles, as shown in Figure 2a. This result is explained through the interparticle long-range dipolar interaction that decays proportionally to m3, where m denotes the average distance between neighboring particles.34 A composite with higher ψ has larger ε, which agrees with the mixing rule concept.5 The ε enhancement is most pronounced for d/2r=0 when the particles form continuous chains and it is negligible when d/2r>0, i.e., when discontinuous chains are formed. Thus, when the composite contains continuous filler particle chains, ε enhancement is more significant than when discontinuous chains are included in it. Hence, results for composites that include continuous particle chains are discussed below.

FIG. 2.

(a) Organizing the filler particles into a chainlike structure increases the relative permittivity ε. When the interparticle spacing d/2r is increased, ε decreases and approaches the value for the composite that includes only randomly distributed particles. When d/2r=0, the particles form continuous particle chains that enhance ε more than discontinuous chains do, which occurs when d/2r>0. The larger the volume fraction ψ, the higher ε is. The green squares, blue diamonds, and red circles correspond to ψ= 1.9%, 8.7%, and 13%, respectively. (b) The clearance h has negligible influence on ε. The entire particle organization can be placed off-center perpendicular to the plate without influencing ε. The clearance for particle organization that is centered between the two plates is hc. The farther the ratio h/hc is from unity, the more off-center is the placement of the particle organization. Considering symmetry, results for h/hc1 are presented for ψ= 8.7%.

FIG. 2.

(a) Organizing the filler particles into a chainlike structure increases the relative permittivity ε. When the interparticle spacing d/2r is increased, ε decreases and approaches the value for the composite that includes only randomly distributed particles. When d/2r=0, the particles form continuous particle chains that enhance ε more than discontinuous chains do, which occurs when d/2r>0. The larger the volume fraction ψ, the higher ε is. The green squares, blue diamonds, and red circles correspond to ψ= 1.9%, 8.7%, and 13%, respectively. (b) The clearance h has negligible influence on ε. The entire particle organization can be placed off-center perpendicular to the plate without influencing ε. The clearance for particle organization that is centered between the two plates is hc. The farther the ratio h/hc is from unity, the more off-center is the placement of the particle organization. Considering symmetry, results for h/hc1 are presented for ψ= 8.7%.

Close modal

Consider that the two clearances between each plate of a capacitor and their nearest particles form two smaller capacitors that are connected in series by a continuous conductive particle chain, i.e.,

(2)

where V1 and V2 denote the electric potentials across each of the two capacitors. No matter how V1 and V2 change, the electric potential V across the two plates is constant. Therefore, when the entire filler organization is placed off-center, changing h does not significantly influence ε, as shown in Figure 2b. Thus, for both continuous and discontinuous chains, ε enhancement is independent of h.

Reducing a diminshes ε enhancement, as shown in Figure 3a. This result is counterintuitive with the general understanding that the capacitance, hence ε, should increase with decreasing distance between two conductive objects. However, the statement is only true for two conductive objects with equal and opposite charges.35 Since each conductive chain herein has zero net charge, this leads to the different observed behavior.

FIG. 3.

(a) Permittivity enhancement ε/εz increases with increasing interchain spacing a. When aax, which is the maximum allowable value for a chainlike arrangement at any volume fraction, permittivity enhancement is largest. The distribution of the surface charge density in (b)-(d) shows that as a decreases, like charges at the ends of the particle chains are repelled towards the middle of the chains, which decreases the induced charges at chain ends, as illustrated by the smaller and non-symmetric blue and purple regions in (d). The smaller the induced charges, the lower the capacitance is, hence the smaller the relative permittivity.

FIG. 3.

(a) Permittivity enhancement ε/εz increases with increasing interchain spacing a. When aax, which is the maximum allowable value for a chainlike arrangement at any volume fraction, permittivity enhancement is largest. The distribution of the surface charge density in (b)-(d) shows that as a decreases, like charges at the ends of the particle chains are repelled towards the middle of the chains, which decreases the induced charges at chain ends, as illustrated by the smaller and non-symmetric blue and purple regions in (d). The smaller the induced charges, the lower the capacitance is, hence the smaller the relative permittivity.

Close modal

Essentially, for the same applied voltage, as a decreases, charge interactions between adjacent particle chains reduce the electric charges at the chain ends. This results in a lower overall capacitance and therefore smaller ε for the composite. Consider two continuous particle chains in a dielectric composite that are placed far enough to avoid long-range dipolar interactions.34 The electric charges induced by an applied voltage V at the two ends of each chain are ±Q0, which can be visualized by examining the distribution of surface charge density (Figure 3bd). When the chains are moved closer, i.e., a is decreased, like charges located at the ends of the chains repel each other and migrate towards the middle of these chains. Here, charges from opposite ends of the chain intersperse and neutralize their polarity (green region in Figure 3c and d). The remaining charges at the two ends of each chain are now ±Q1 where Q1<Q0, shown as the smaller and asymmetrical blue and purple (higher surface charge densities) regions in Figure 3d. For a capacitor of specific overall geometry and applied voltage, there is a proportional relationship between its electric charges vρvdv=Q and ε, cf. Equation (1). Therefore, with decreasing a, charge interactions within the chains diminish the overall charge at their ends, diminishing ε enhancement, as shown in Figure 3a. Since a has the maximum value ax for any ψ (see the supplementary material for details), Figure 3a also implies that when the particle chains are placed evenly with a/ax=1, ε enhancement is largest.

Increasing ψ increases ε. The ε enhancement is more significant for a composite containing continuous particle chains than one that includes randomly distributed particles (εz), as shown in Figure 4a. When ψ< 9%, a composite containing smaller filler particles has higher ε than one containing larger particles, whereas when ψ> 9% the behavior is opposite because chains formed with smaller particles have a greater overall surface area.36 In addition, the surface morphologies of the interfaces between the matrix and smaller filler particles are more intricate, leading to a higher interfacial polarization,37,38 a size advantage that diminishes as ψ increases. Due to the competing influences of increasing ψ and charge interactions discussed above, a composite containing smaller particles has an inherently smaller ax when ψ> 9%, resulting in more vigorous charge interactions that decrease the ε enhancement.

FIG. 4.

(a) When particles are arranged into continuous chains (solid lines) or distributed randomly (dashed lines), the relative permittivity ε increases with volume fraction ψ for all three particle sizes, r=0.4μm (green squares), 0.5 μm (blue diamonds), and 0.68 μm (red circles). ε is higher when particles form continuous chains, i.e., d/2r=0. When the particles are randomly distributed, ε is smaller. For smaller particles (green squares), ε is higher when ψ< 9%, but this size advantage diminishes as ψ increases. This is because, when ψ> 9%, charge interactions between adjacent chains that contain smaller particles are stronger than for chains with larger particles, which diminishes ε enhancement. (b) The ε enhancement for composites that have filler particle organizations is compared with counterparts containing randomly distributed particles εz. Maximum enhancement occurs at ψ= 7%, 9%, and 11% for particle radii r=0.4μm, 0.5 μm, and 0.68 μm, respectively. This figure demonstrates that organized particle arrangements are more beneficial than randomly distributed filler particles for permittivity enhancement.

FIG. 4.

(a) When particles are arranged into continuous chains (solid lines) or distributed randomly (dashed lines), the relative permittivity ε increases with volume fraction ψ for all three particle sizes, r=0.4μm (green squares), 0.5 μm (blue diamonds), and 0.68 μm (red circles). ε is higher when particles form continuous chains, i.e., d/2r=0. When the particles are randomly distributed, ε is smaller. For smaller particles (green squares), ε is higher when ψ< 9%, but this size advantage diminishes as ψ increases. This is because, when ψ> 9%, charge interactions between adjacent chains that contain smaller particles are stronger than for chains with larger particles, which diminishes ε enhancement. (b) The ε enhancement for composites that have filler particle organizations is compared with counterparts containing randomly distributed particles εz. Maximum enhancement occurs at ψ= 7%, 9%, and 11% for particle radii r=0.4μm, 0.5 μm, and 0.68 μm, respectively. This figure demonstrates that organized particle arrangements are more beneficial than randomly distributed filler particles for permittivity enhancement.

Close modal

For smaller ψ, ax is larger so that interchain charge interactions are negligible. These interactions become stronger with increasing ψ for continuous particle chains regardless of particle size. Therefore, for organized particles, ε increases with ψ and ax but the higher ψ is, the smaller ax becomes whereas εz depends on ψ only. Consequently, there exists an optimal value of ψ for ε enhancement that maximizes the benefit of using organized particles as compared to simply distributing them randomly. The permittivity reaches a maximum value ε/εz1.85 for any particle size at an optimal ψ value, as shown in Figure 4b. For instance, when r=0.68μm, the maximum value of ε/εz occurs at ψ= 11% while the maxima for r=0.4μm and 0.5 μm lie at ψ= 7% and 9%, respectively.

The particle chains discussed above are straight and aligned (θ=0°) with overall chain directions perpendicular to the conductive plates. In Figure 5, we demonstrate the influence of θ=0° to 60° and ψ on ε for r=0.5μm. For continuous particle chains (Figure 5a), ε decreases with increasing θ and even more rapidly with increasing ψ. Without varying the number of particles in a chain, the larger θ is, the closer are the particles within a chain. Since the electric potential over the entire surface of the chain is the same, for a chain that is oriented perpendicular to the equipotential surfaces, the closer the particles are, the lower is the distortion of the equipotential surfaces. This leads to smaller electric field strengths and therefore smaller ε for the composite. For discontinuous chains, ε is more or less independent of θ, as shown in Figure 5b, the largest difference in ε being 0.265. Hence, through various combinations of ψ and θ for the same matrix and filler material, it is possible to design elastomer-particle composites with ε in the range from 3 to 8 (Figure 5c).

FIG. 5.

Relative permittivity ε as a function of zigzag angle ε and volume fraction ε for (a) continuous chains (d/2r=0) and (b) discontinuous chains (d/2r=0.2), both composed of particles with radius r=0.5μm. For continuous chains, ε decreases with increasing θ. The larger ψ is, the greater is the decrease in ε. For d/2r=0.2, as shown in (b), θ has negligible effect on ε with the largest difference in ε being 0.265. (c) Surface plot of ε, ψ, and θ for r=0.5μm. The top surface corresponds to the results for d/2r=0 whereas that on the bottom is for d/2r=0.2. The same overall behavior is observed, irrespective of the constituent particle size. (See supplementary material for the corresponding figures for r=0.4μm and 0.68 μm.)

FIG. 5.

Relative permittivity ε as a function of zigzag angle ε and volume fraction ε for (a) continuous chains (d/2r=0) and (b) discontinuous chains (d/2r=0.2), both composed of particles with radius r=0.5μm. For continuous chains, ε decreases with increasing θ. The larger ψ is, the greater is the decrease in ε. For d/2r=0.2, as shown in (b), θ has negligible effect on ε with the largest difference in ε being 0.265. (c) Surface plot of ε, ψ, and θ for r=0.5μm. The top surface corresponds to the results for d/2r=0 whereas that on the bottom is for d/2r=0.2. The same overall behavior is observed, irrespective of the constituent particle size. (See supplementary material for the corresponding figures for r=0.4μm and 0.68 μm.)

Close modal

In summary, the manner in which conductive filler particles contained within a PDMS elastomeric matrix are arranged influences the electric permittivity of the composite material. Organizing the particles into chainlike structures enhances the bulk permittivity by 85% in comparison with that for a composite containing the same filler, but which is randomly distributed. The value of ψ that is most beneficial when organized particles are used is 9% α 2%. The bulk permittivity increases with increasing ψ and a(x), but decreases with increasing d and ψ. When ψ< 9%, a composite containing smaller particles has higher ε, while for ψ> 9%, introduction of larger particles increases ε. The enhancement in ε is more gradual beyond ψ 9% due to the competition between effects induced by increasing ψ, which increases ε, and decreasing ax, which decreases ε. When adjacent particles are interconnected, ε enhancement is more pronounced because continuous particle chains that are oriented perpendicular to the equipotential surfaces increase the induced charge of the composite, hence its bulk permittivity. These results are relevant for the design of soft dielectric materials with different permittivities when identical material compositions are used but their filler particles morphologies are varied, such as for additive manufacturing flexible electronics and 3D printed circuits.

See supplementary material for the model verifications, the determination of ax, and results for zigzag configurations with varying particle sizes.

This work was supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (NSERC) under grant number RGPIN-2014-04066, and the Mitacs Globalink program.

1.
R. J.
Young
and
P. A.
Lovell
,
Introduction to polymers
(
CRC press
,
2011
).
2.
M.
Alger
,
Polymer Science Dictionary
(
Springer Science & Business Media
,
1996
).
3.
D.
Ponnamma
,
C. W.
Kishor Kumar Sadasivuni
,
S.
Thomas
, and
M. A.-A.
AlMa’adeed
,
Flexible and Stretchable Electronic Composites
(
Springer International Publishing
,
Cham
,
2016
).
4.
A. J.
Pickles
and
M. B.
Steer
,
IEEE Access
1
,
523
(
2013
).
5.
P.
Barber
,
S.
Balasubramanian
,
Y.
Anguchamy
,
S.
Gong
,
A.
Wibowo
,
H.
Gao
,
H. J.
Ploehn
, and
H.-C.
zur Loye
,
Materials
2
,
1697
(
2009
).
6.
D.
Fragouli
,
A.
Das
,
C.
Innocenti
,
Y.
Guttikonda
,
S.
Rahman
, and
L.
Liu
,
ACS Appl. Mater. Interfaces
6
,
4535
(
2014
).
7.
P. J.
Tsai
,
S.
Ghosh
,
P.
Wu
, and
I. K.
Puri
,
ACS Appl. Mater. Interfaces
(
2016
).
8.
M. J.
Guffey
and
N. F.
Scherer
,
Nano Lett.
10
(
11
),
4302
(
2010
).
9.
A. A.
Tseng
,
K.
Chen
,
C. D.
Chen
, and
K. J.
Ma
, “
Electronics packaging manufacturing
,”
IEEE Transactions on
26
(
2
),
141
(
2003
).
10.
A. R.
Abdel Fattah
,
S.
Ghosh
, and
I. K.
Puri
,
J. Magn. Magn. Mater.
401
,
1054
(
2016
).
11.
A. R.
Abdel Fattah
,
S.
Ghosh
, and
I. K.
Puri
,
ACS Appl. Mater. Interfaces
8
(
17
),
11018
(
2016
).
12.
R. R.
Collino
,
T. R.
Ray
,
R. C.
Fleming
,
C. H.
Sasaki
,
H.
Haj-Hariri
, and
M. R.
Begley
,
Extreme Mechanics Letters
5
,
37
(
2015
).
13.
R. L.
Truby
and
J. A.
Lewis
,
Nature
540
(
7633
),
371
(
2016
).
14.
S.
Ghosh
,
M.
Tehrani
,
M. S.
Al-Haik
, and
I. K.
Puri
,
Materials
8
(
2
),
474
(
2015
).
15.
S.
Ghosh
and
I. K.
Puri
,
Faraday Discuss.
181
,
1
(
2014
).
16.
S.
Ghosh
and
I. K.
Puri
,
Soft Matter
9
,
2024
(
2013
).
17.
R.
Ganguly
and
I. K.
Puri
,
Adv. Appl. Mech.
41
,
293
(
2007
).
18.
A. C.
Siegel
,
S. T.
Phillips
,
M. D.
Dickey
,
N.
Lu
,
Z.
Suo
, and
G. M.
Whitesides
,
Adv. Funct. Mater.
20
,
28
(
2010
).
19.
F. J.
Juskey
,
A. B.
Suppelsa
, and
D. W.
Dorinski
, U.S. patent 5,264,061 (
1993
).
20.
H. F.
Schmidt
,
S.
Rauchmaul
, and
J.
Bednarz
, U.S. patent 5008496 A (
1991
).
21.
A.
Frutiger
,
J. T.
Muth
,
D. M.
Vogt
,
Y.
Mengüç
,
A.
Campo
,
A. D.
Valentine
,
C. J.
Walsh
, and
J. A.
Lewis
,
Adv. Mater.
27
,
2440
(
2015
).
22.
J. T.
Muth
,
D. M.
Vogt
,
R. L.
Truby
,
Y.
Mengüç
,
D. B.
Kolesky
,
R. J.
Wood
, and
J. A.
Lewis
,
Adv. Mater.
26
(
36
),
6202
(
2014
).
23.
Z. Z.
Sun
,
A.
López
, and
J.
Schliemann
,
J. Appl. Phys.
109
,
104303
(
2011
).
24.
J.
Ge
,
L.
He
,
Y.
Hu
, and
Y.
Yin
,
Nanoscale
3
,
177
(
2011
).
25.
H.
Stoyanov
,
M.
Kollosche
,
S.
Risse
,
D. N.
McCarthy
, and
G.
Kofod
,
Soft Matter
7
,
194
(
2011
).
26.
A.
Boudenne
,
Y.
Mamunya
,
V.
Levchenko
,
B.
Garnier
, and
E.
Lebedev
,
Eur. Polym. J.
63
,
11
(
2015
).
27.
Y.
Han
,
W.
Hong
, and
L. E.
Faidley
,
Int. J. Solids Struct.
50
,
2281
(
2013
).
28.
T.
Borbáth
,
S.
Günther
,
D.
Yu Borin
,
T.
Gundermann
, and
S.
Odenbach
,
Smart Mater. Struct.
21
(
10
),
105018
(
2012
).
29.
E.
Coquelle
and
G.
Bossis
,
Journal of Advanced Science
17
,
132
(
2005
).
30.
J.
Trajkovikj
,
J.-F.
Zurcher
, and
A. K.
Skrivervik
, in
Antennas & Propagation Conference (LAPC), 2012 Loughborough Loughborough
,
2012
, pp.
1
4
.
31.
S. J.
van Enk
,
Journal of Modern Optics
42
,
321
(
1995
).
32.
G. L.
Klimchitskaya
,
U.
Mohideen
, and
V. M.
Mostepanenko
,
Rev. Mod. Phys.
81
,
1827
(
2009
).
33.
P.
Brochu
and
Q.
Pei
,
Macromol. Rapid Commun.
31
,
10
(
2010
).
34.
S. P.
Gubin
,
G. B.
Khomutov
, and
Y. A.
Koksharov
, in
Magnetic Nanoparticles
, edited by
S. P.
Gubin
(
Wiley-VCH Verlag GmbH & Co. KGaA
,
Weinheim, Germany
,
2009
), pp.
117
195
.
35.
H. A.
Haus
and
J. R.
Melcher
, in
Electromagnetic Fields and Energy
(
Prentice-Hall, Englewood Cliffs
,
NJ
,
1989
).
36.
Prateek
,
V. K.
Thakur
, and
R. K.
Gupta
,
Chem. Rev.
116
,
4260
(
2016
).
37.
N.
Xu
,
L.
Hu
,
Q.
Zhang
,
X.
Xiao
,
H.
Yang
, and
E.
Yu
,
ACS Appl. Mater. Interfaces
7
,
27373
(
2015
).
38.
J.
Yuan
,
S.
Yao
,
W.
Li
,
A.
Sylvestre
, and
J.
Bai
,
The Journal of Physical Chemistry C
118
,
22975
(
2014
).

Supplementary Material