Open AccessArticle

Polar Vector Property of the Stationary State of Condensed Molecular Matter

Jürg Hulliger

^*,

Luigi Cannavacciuolo

and

Mathias Rech

Department of Chemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland

Author to whom correspondence should be addressed.

Symmetry 2014, 6(4), 844-850; https://doi.org/10.3390/sym6040844

Submission received: 3 June 2014 / Revised: 15 September 2014 / Accepted: 26 September 2014 / Published: 13 October 2014

(This article belongs to the Special Issue Crystal Symmetry and Structure)

Download

Browse Figures

Versions Notes

Abstract

Crystalline phases undergoing 180° orientational disorder of dipolar entities in the seed or at growing (hkl) faces will show a polar vector property described by ∞/mm symmetry. Seeds and crystals develop a bi-polar state (∞/mm), where domains related by a mirror plane m allow for a ∞m symmetry in each domain. The polarity of domains is due to energetic favorable interactions at the object-to-nutrient interface. Such interactions are well reproduced by an Ising Hamiltonian. Two-dimensional Monte Carlo simulations performed for real molecules with full long-range interactions allow us to calculate the spatial distribution of the electrical polarization P_el. The investigation has been extended to liquid droplets made of dipolar entities by molecular dynamics simulations. We demonstrate the development of an

m \bar{\infty}

quasi bi-polar state leading to a charged surface.

Keywords:

molecular crystals; polarity; parity; continuous groups; orientational disorder

1. Introduction

Recent experimental and theoretical work [1] lead us to pose a fundamental question: Can mono-domain polar molecular crystals exist [2,3]? Looking at data of the CSD, there seems to be no doubt of a positive answer: roughly 17% of the known structures of molecular compounds [4] (no polymers and proteins) were refined to a polar space group. If we exclude ionic systems (salts) to only consider neutral small molecules, there remain about 9% of the polar structures. Although twinning and disorder are currently addressed, we notice that only a small number of structural investigations have provided some knowledge on: (i) the absolute orientation of the polar axis in relation to the morphology [5–12]; and (ii) just in a few cases, the spatial distribution of the electrical polarization P_el has been analyzed [13–15]. As a matter of fact, we do not know much about the polar domain states of (non-ferroelectric) molecular crystals. By analyzing both the 180° orientational disorder of dipolar entities and the forms of twinning related to the spatial distribution of the polarization, we can, however, answer the aforementioned key question: for crystals expressing a polar morphology by their (hkl) faces, a mono-domain, i.e., a homogeneous polarization distribution, may be observed experimentally [5–12]. This is a kinetic effect due to the critical slowing down occurring in real growth that hinders the reversal transition into a bi-polar state, as predicted by statistical mechanics calculations [2,3] (Figure 1a). This very form of twinning process can, in fact, transform an initial mono-domain polar seed into a final object featuring macroscopic domains of opposite polarities. In this article, we focus on understanding the bi-polar state of condensed matter constituted by polar molecules by symmetry and statistical mechanics arguments. Because of the close analogy, the analysis will also include liquid droplets (Figure 1b) and crystals made thereof.

2. Symmetry and Conservation

The laws of physics are found invariant under the consecutive action of C (charge conjugation), P (parity operation, space inversion, reflection) and T (time reversal) symmetry operations [16]. Molecules and their condensed states are described by quantum and statistical mechanics. The Hamiltonian H_mol of a molecular system commutes with each of these operations. Similarly, this holds for molecular crystals made stable by inter-molecular potentials. In our context, parity P is of particular relevance: both the polar vector operators μ_el (dipole moment of molecules) and P_el (polarization of molecular crystals) change their direction (sign) under the action of P (so, C, not T) [17]. When attempting to calculate observable properties, such as μ_el or P_el, we encounter the issue of integrating a state equation over time and space, because for every experiment, a certain time-space domain is committed. In quantum mechanics, the dipole moment μ_el is defined by

\int Ψ * e r Ψ d V

(e: charge of the electron; r: distance vector). Because of:

\int Ψ * e r Ψ d V = e \int Ψ * r Ψ d V = e \int r {| Ψ |}^{2} d V

(1)

where the last integrand is the product of a function, r, with parity odd (−) and Ψ^*Ψ = |Ψ|² with parity even (+), and since the integral is calculated over the whole space, the expectation value is zero. This is a short hand proof of a statement we can find in textbooks on quantum mechanics [18]: “Quantum-mechanical systems in the stationary state do not have permanent electric dipole moments”. For most molecules, Born–Oppenheimer-based calculations provide a description that can attribute a rigid nuclear frame. This may be extended to a vibronic state; however, the amplitudes of nuclear motion will only cover a small part of the full space of coordinates. Consequently, such a time averaged single object of fixed nuclear conformation (transforming under a polar point group) shows a non-zero dipole moment μ_el. Assuming quasi-rigid objects carrying a dipole moment, thermal energy, i.e., temperature, is the variable that will lead to averaging. Consequently, a vector property showing P = −1 for single objects will not be observed for a stationary state of an ensemble. However, droplets adopt a bi-polar type of state (Figure 1b), when projecting single dipole moments onto diameters for a summation to yield 〈P_el〉_zone. Preliminary molecular dynamics (MD) calculations for carbon monoxide in the 2D liquid state show that the oxygen atom of CO molecules is preferentially pointing out of the boundary, whereas towards the center, averaging of 〈P_el〉 takes place.

When a liquid-to-solid transition occurs, packing may affect a unique polar axis. Along this axis, physical measurements will find polar vector and higher order uneven rank tensor properties, irrespective of what is giving rise to the polarity: isolated molecules or molecular field-modified entities [19]. Here, we just address the symmetry of the projection of 〈μ_el〉 onto this unique axis. The axis itself represents ∞/m(C_∞h) symmetry [20]. Because of the finite size of crystals, the mirror plane m gets localized in the middle of the object. This means that the crystal splits into two domains with m perpendicular to the unique axis. A mono-domain polar state is thus not allowed: a polar packing, i.e., the average polarization, is expected to develop a bi-polar state, described by ∞/mm(D_∞h) symmetry [20] (Figure 1a). In each domain, a polar vector type property of ∞m symmetry is allowed. For a detailed discussion of the eigensymmetry and the composite symmetry of twinned crystals, see [21] and Figure 3.3.10.2. therein. Considering real crystal structures, e.g., nucleation processes leading to n/m (n = 2, 3, 4, 6) point groups and growth along the axis n, the stochastic process transforms pro-chiral faces into chiral sectors. Because of the mirror plane m passing through a seed, sectors involving the +n or the −n direction show an enantiomorphic relation (for a discussion of other point groups and low index faces, see [22]).

So far, the bi-polar state for molecular crystals was experimentally demonstrated by scanning pyroelectric and phase-sensitive second harmonic microscopy [1]. Here, the measurement of the circular dichroism of each sector could add a third proof for a stochastic mechanism of symmetry lowering.

What may be a mechanism by which a molecular packing can attain a bi-polar state? By the 180° orientational disorder of dipolar building blocks: Monte Carlo (MC) simulations [2,3] have suggested a growth-induced stochastic reversal transition for building blocks, when starting from a polar seed, ending up in a bi-polar state. At this point, we turn to the seed itself. What will be the polar state of dipolar molecules aggregating in a nucleus leading to further crystallization? MC simulations based on a simple interaction scheme [2,3] have shown that a 3D system adopts a bi-polar state. Most recently, MC simulations with full long range interactions confirmed [23] this behavior for a 2D array of acceptor-donor-substituted benzene molecules (Figure 2).

For both MC-type simulations, many random initial orientational states for the entities are employed for the average of the quantities of interest. From there, the system evolves into a bi-polar state. Analogously, when starting from an ordered, but centric seed, MC growth also results in a bi-polar state [2,3], producing a certain degree of polar alignment in each of the domains. As presented so far, the bi-polar state was demonstrated by simulations. An Ising [24] model describing a system of interacting adjacent chains allows us to understand that opposite average polarities can exist in each of the domains. The spin Hamiltonian H (Equation (2)) has been derived in two dimensions, under very general assumptions: (i) for any given state configuration S_i = ±1, i = 0,…, N, Equation (2) is the energy of that configuration; (ii) periodic boundary conditions on the transverse direction and a free boundary condition on the longitudinal one.

H = ϵ + Δ E_{f} \sum_{j} (S_{j} (0) - S_{j} (N)) + \sum_{i j} [S_{j + 1} (i) S_{j} (i) Δ E_{⊥} / 2 - S_{i + 1} (j) S_{i} (j) (Δ E_{A} + Δ E_{D}) / 4]

(2)

where ∆E_⊥ = E_par − E_antip, where E_par (resp. E_antip) is the energy of two next-neighbor spins parallel (resp. anti-parallel), due to the transverse interaction. The sums are taken over all lattice points, specifically, i and j run over the longitudinal and transverse direction, respectively; ∊ is a constant term. Whereas the bilinear terms (longitudinal and transverse interactions) can induce polar alignment of vectors, they do not provide a definite direction of polarity 〈P_el〉. Absolute polarity (directionality) is introduced by the linear effective particle operators S₀ and S_N in H. These terms are formally equivalent to the action of an electrical field E, representing and, thus, introducing ∞m [20,25], i.e., polar vector symmetry. The energy parameter ∆E_f [2,3] associated with these two linear terms in H reflect the polar symmetry of interacting molecules: in a simple description of the interactions, ∆E_f is equal to E_AA − E_DD, i.e., the difference in the energy for defect formation in a single chain: E_AA : −A ⋯ A−; E_DD : −D ⋯ D− (A: acceptor; D: donor). Therefore, at each border of the unique axis, an effective field (∞m) [25] is reducing the symmetry of the axis to yield a polar vector property (∞m; also chiral ∞) in each domain. It can be shown that the introduction of the ∆E_f energy parameter is not restricted to nearest neighbor interactions, but can be extended to a lattice sum.

3. Summary and Conclusions

One main point emerges from this study: in the case of a crystal lattice, (i) the early (nucleation) or (ii) the later state (growth) can accommodate 180° orientational defects of dipolar entities (Figure 1b); both the seed and the final crystal will show a polar vector property described by ∞/mm (D_∞h) [20] symmetry, i.e., a bi-polar state (Figure 1a), where domains related by a mirror plane m show a polar vector property (∞m, 〈P_el〉) induced by the linear terms ∆E_f S_i(i = 0, N) in H (Equation (2)). This general behavior predicted by H was reproduced by different types of MC and MD simulations, assuming either a simple nearest neighbor interaction scheme [2,3] or long-range interactions of real molecules (Figure 2).

Did we properly answer our initial question? There is one important point left out so far: in case molecules by kinetic control, i.e., nucleation at high supersaturation/supercooling, form a mono-domain polar seed (instead of a well-thermalized bi-polar object), this state may persist up to a macroscopic size, because the probability to undergo the reversal transition [2,3] may be too small to occur up to the final size that a crystal can attain. It is, therefore, likely to find crystals representing a polar morphology and a mono-domain state. Nevertheless, our final conclusion is the mono-domain polar molecular crystals, for reasons of symmetry and orientational defects (configurational entropy), should not exist.

Acknowledgments

This work was supported by the SNF, project No. 200021_129472/1. We thank Helmut Klapper and Theo Hahn for a helpful discussion on continuous groups and Jürg Hauser for data searching.

Author Contributions

Jürg Hulliger: general theoretical frame; Luigi Cannavacciulolo: MC simulations, Hamiltonian; Mathias Rech: Hamiltonian and basic physics.

Conflicts of Interest

The authors declare no conflict of interest.

References

Hulliger, J.; Wüst, T.; Brahimi, K.; Burgener, M.; Aboulfadl, H. A Stochastic principle behind polar properties of condensed molecular matter. New J. Chem. 2013, 37, 2229–2235. [Google Scholar]
Hulliger, J.; Wüst, T.; Brahimi, K.; Martinez Garcia, J.C. Can mono domain polar molecular crystals exist? Cryst. Growt. Des. 2012, 12, 5211–5218. [Google Scholar]
Hulliger, J.; Wüst, T.; Rech, M. Symmetry and the polar state of condensed molecular matter. Z. Kristallogr. 2013, 228, 607–610. [Google Scholar]
Allen, F.H. The Cambridge Structural Database: A quarter of a million crystal structures and rising. Acta Crystallogr. Sect. B. 2002, 58, 380–388. [Google Scholar]
Curtin, D.Y.; Paul, I. Chemical consequences of the polar axis in organic chemistry. Chem. Rev. 1981, 81, 525–541. [Google Scholar]
Davey, R.J.; Milisavljevic, B.; Bourne, J.R. Solvent Interactions at Crystal Surfaces: The Kinetic Story of α-Resorcinol. J. Phys. Chem. 1988, 92, 2032–2036. [Google Scholar]
Chen, B.D.; Garside, J.; Davey, R.J.; Maginn, S.J.; Matsuoka, M. Growth of m-Chloronitrobenzene Crystals in the Presence of Tailor-Made Additives: Assignment of the Polar Axes from Morphological Calculations. J. Phys. Chem. 1994, 98, 3215–3221. [Google Scholar]
Buguo, W.; Weizhuo, Z.; Cheng, Y.; Erwei, S.; Zhiwen, Y. Analysis of growth habits of polar organic crystal N-4-nitrophenyl-(L)-prolinol (NPP) based on the incorporation of growth units. J. Cryst. Growth. 1996, 160, 375–381. [Google Scholar]
Ryu, G.; Yoon, C.S. Seeded supercooled melt growth and polar morphology of organic nonlinear optical crystal, meta-nitroaniline (mNA). J. Cryst. Growth. 1998, 191, 190–198. [Google Scholar]
Hong, H.-K.; Park, J.W.; Lee, K.-S.; Yoon, C.S. Growth of highly nonlinear optical organic crystal, 3-methyl-4-methoxy-4′-nitrostilbene (MMONS). J. Cryst. Growth. 2005, 277, 509–517. [Google Scholar]
Cuppen, H.M.; Beurskens, S.; Kozuka, S.; Tsukamoto, K.; Smits, J.M.M.; de Gelder, R.; Grimbergen, R.F.P.; Meekes, H. Crystal Structure and Growth Behavior of Aspartame Form I-A. Cryst. Growt. Des. 2005, 5, 917–923. [Google Scholar]
Srinivasan, K.; Sherwood, J.N. Asymmetric Growth of α-Resorcinol Crystals: Comparison of Growth from the Vapor Phase and from Aqueous Solution. Cryst. Growt. Des. 2005, 5, 1359–1370. [Google Scholar]
Batagiannis, A.; Wübbenhorst, M.; Hulliger, J. Piezo- and pyroelectric microscopy. Curr. Op. Solid State Mater. Sci. 2010, 14, 107–115. [Google Scholar] [Green Version]
Behrnd, N.-R.; Couderc, G.; Wübbenhorst, M.; Hulliger, J. Scanning pyroelectric microscopy revealing the spatial polarity distribution in topologically centric crystals of trans-4-chloro-4′-nitrostilbene. Phys. Chem. Chem. Phys. 2006, 8, 4132–4137. [Google Scholar]
Labat, G.; Behrnd, N.-R.; Couderc, G.; Bonin, M.; Tsuwi, J.; Batagannis, A.; Berger, R.; Bertoni, M.; Prodi-Schwab, A.; Hulliger, J. Polymorphism polar morphology and absolute structure determination of 4-iodo-4′-nitrobiphenyl (INBP). Cryst. Eng. Comm. 2010, 12, 1252–1262. [Google Scholar]
Feynman, R.P.; Leighton, R.B.; Sands, M. Lectures o. Physics., 17th ed; Addison-Wesley: Reading, MA, USA, 1972; Volume I, III. [Google Scholar]
Wagnière, G.H. On Chirality and the Universa. Asymmetry.; Wiley: Zürich, Switzerland, 2007. [Google Scholar]
Blatt, J.M.; Weisskopf, V. Theoretical Nuclea. Physics.; Wiley: New York, NY, USA, 1952; Chapter 7; pp. 23–25. [Google Scholar]
Spackmann, M.A.; Munshi, P.; Dittrich, B. Dipole moment enhancement in molecular crystals from X-ray diffraction data. Chem. Phys. Chem. 2007, 8, 2051–2063. [Google Scholar]
Hahn, T. International Tables for Crystallography, Vol. A Space-Group Symmetry, Tab. 10.1 A Space-Group Symmetry, Tab. 10.1.4.2.; Kluwer Acad. Publ.: Dordrecht, NL, 2002; In previous work [2,3] we were using the notation by Shubnikov et al. (∞/∞m), being replaced here by the international symbol. [Google Scholar]
Authier, A. International Tables for Crystallography Volume D: Physical Properties o. Crystals.; Kluwer Academic Publishers: Dordrecht, the Netherlands, 2003; Chapters 3.2–3.4. [Google Scholar]
Gervais, C.; Hulliger, J. Impact of Surface Symmetry: On the prediction of Growth-Induced Properties. Cryst. Growt. Des. 2007, 7, 1925–1935. [Google Scholar]
Cannavacciuolo, L.; Hulliger, J. University of Bern, Bern, Switzerland. . Monte Carlo simulations demonstrating the bi-polar state for a 2D nucleus made of real molecules undergoing long rang. interactions. 2014; Unpublished work. [Google Scholar]
Greiner, W.; Neise, L.; Stöcker, H. Thermodynamik und Statistisch. Mechanik., 2nd ed; Harri Deutsch: Frankfurt, Germany, 1993; Volume 9. [Google Scholar]
Paufler, P. Physikalisch. Kristallographie.; Akad. Verlag: Berlin, Germany, 1986; p. 29. [Google Scholar]

Figure 1. (a) Bi-polar state ∞/mm showing domains (∞m). The polarization may be oriented “down”(as shown) or “up”; (b) Quasi bi-polar state of a stationary spherical object

(m \bar{\infty}, K_{h})

Figure 1. (a) Bi-polar state ∞/mm showing domains (∞m). The polarization may be oriented “down”(as shown) or “up”; (b) Quasi bi-polar state of a stationary spherical object

(m \bar{\infty}, K_{h})

Figure 2. Average polarization vs. lattice site for an MC 2D-simulation, showing the bi-polar state formation. The full line is the fit with an exponential function. The system is composed of 21 × 5 molecules (1-bromo-4′-cyano-tetrafluoro-benzene) arranged in five parallel chains.

© 2014 by the authors; licensee MDPI, Basel, Switzerland This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Hulliger, J.; Cannavacciuolo, L.; Rech, M. Polar Vector Property of the Stationary State of Condensed Molecular Matter. Symmetry 2014, 6, 844-850. https://doi.org/10.3390/sym6040844

AMA Style

Hulliger J, Cannavacciuolo L, Rech M. Polar Vector Property of the Stationary State of Condensed Molecular Matter. Symmetry. 2014; 6(4):844-850. https://doi.org/10.3390/sym6040844

Chicago/Turabian Style

Hulliger, Jürg, Luigi Cannavacciuolo, and Mathias Rech. 2014. "Polar Vector Property of the Stationary State of Condensed Molecular Matter" Symmetry 6, no. 4: 844-850. https://doi.org/10.3390/sym6040844

APA Style

Hulliger, J., Cannavacciuolo, L., & Rech, M. (2014). Polar Vector Property of the Stationary State of Condensed Molecular Matter. Symmetry, 6(4), 844-850. https://doi.org/10.3390/sym6040844

Article Menu

Polar Vector Property of the Stationary State of Condensed Molecular Matter

Abstract

1. Introduction

2. Symmetry and Conservation

3. Summary and Conclusions

Acknowledgments

Author Contributions

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI