TY - JOUR

T1 - Longitudinal polarizability of long polymeric chains

T2 - Quasi-one-dimensional electrostatics as the origin of slow convergence

AU - Kudin, Konstantin N.

AU - Car, Roberto

AU - Resta, Raffaele

PY - 2005

Y1 - 2005

N2 - The longitudinal linear polarizability α (N) of a stereoregular oligomer of size N is proportional to N in the large- N limit, provided the system is nonconducting in that limit. It has long been known that the convergence of α (N) N to the asymptotic α∞ value is slow. We show that the leading term in the difference between α (N) N and α∞ is of the order of 1N. The difference [α (N) -α (N-1)], as well as αcenter (N) (when computationally accessible), also converge to α∞, but faster, the leading term being of the order of 1 N2. We also present evidence that in these cases the power law convergence behavior is due to quasi-one-dimensional electrostatics, with one exception. Specifically, in molecular systems the difference between α (N) N and α∞ has not just one but two sources of the O (1N) term, with one being due to the aforementioned Coulomb interactions, and the second due to the short ranged exponentially decaying perturbations on chain ends. The major role of electrostatics in the convergence of the remainders is demonstrated by means of a Clausius-Mossotti-type classical model. The conclusions derived from the model are also shown to be applicable in molecular systems, by means of test-case ab initio calculations on linear stacks of H2 molecules, and on polyacetylene chains. The implications of the modern theory of polarization for extended systems are also discussed.

AB - The longitudinal linear polarizability α (N) of a stereoregular oligomer of size N is proportional to N in the large- N limit, provided the system is nonconducting in that limit. It has long been known that the convergence of α (N) N to the asymptotic α∞ value is slow. We show that the leading term in the difference between α (N) N and α∞ is of the order of 1N. The difference [α (N) -α (N-1)], as well as αcenter (N) (when computationally accessible), also converge to α∞, but faster, the leading term being of the order of 1 N2. We also present evidence that in these cases the power law convergence behavior is due to quasi-one-dimensional electrostatics, with one exception. Specifically, in molecular systems the difference between α (N) N and α∞ has not just one but two sources of the O (1N) term, with one being due to the aforementioned Coulomb interactions, and the second due to the short ranged exponentially decaying perturbations on chain ends. The major role of electrostatics in the convergence of the remainders is demonstrated by means of a Clausius-Mossotti-type classical model. The conclusions derived from the model are also shown to be applicable in molecular systems, by means of test-case ab initio calculations on linear stacks of H2 molecules, and on polyacetylene chains. The implications of the modern theory of polarization for extended systems are also discussed.

UR - http://www.scopus.com/inward/record.url?scp=24144465702&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=24144465702&partnerID=8YFLogxK

U2 - 10.1063/1.1871934

DO - 10.1063/1.1871934

M3 - Article

C2 - 15847502

AN - SCOPUS:24144465702

SN - 0021-9606

VL - 122

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

IS - 13

M1 - 134907

ER -