Next Article in Journal
Monitoring Test for Stability of Dependence Structure in Multivariate Data Based on Copula
Previous Article in Journal
Grey Coupled Prediction Model for Traffic Flow with Panel Data Characteristics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Possible Application of the Contribution of Aromaticity to Entropy: Thermal Switch

by
Romain Coustel
1,†,
Stéphane Carniato
2,† and
Gérard Boureau
2,*,†
1
Laboratoire de Chimie Physique et Microbiologie Pour l’Environnement, CNRS UMR 7564/Universitée de Lorraine, 405 Rue de Vandoeuvre, 54601 Villers-lès-Nancy, France
2
Laboratoire de Chimie Physique Matière et Rayonnement, Sorbonne Universités, UPMC Univ Paris 06, UMR 7614, 11 rue P. et M. Curie, 75005 Paris, France
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Entropy 2016, 18(12), 456; https://doi.org/10.3390/e18120456
Submission received: 2 March 2016 / Revised: 25 October 2016 / Accepted: 26 October 2016 / Published: 20 December 2016
(This article belongs to the Section Thermodynamics)

Abstract

:
It has been known for a long time that the loss of aromaticity of gaseous molecules leads to a large increase of the enthalpy and to a tiny increase of the entropy. Generally, the calculated transition temperature from an aromatic structure towards a non-aromatic structure at which these two contributions cancel is very high. The entropy associated to the loss of aromaticity of adsorbed molecules, such as pyridine on Si(100) and on Ge(100), is roughly the same while the associated enthalpy is much smaller, a consequence of which is a low transition temperature. This allows us to imagine monomolecular devices, such as thermal switches, based on the difference of electrical conductivity between aromatic and non-aromatic species adsorbed on Si(100) or on Ge(100).

1. Introduction

Since the pioneering work of Aviram and Ratner [1] discussing the use of a single organic molecule to construct a rectifier, molecular electronics has been a continuously expanding field [2]. A number of applications have been developed, such as the organic synthesis of conductive polymers [3]. In the present paper, we discuss a possible application of a curious property of aromaticity, the existence of a small negative contribution of aromaticity to entropy.
It is well known that aromatic compounds present much lower enthalpy of formation than the one expected for the corresponding Kekule structures. Similarly, more than sixty years ago, Franklin [4] predicted a weak and negative contribution of aromaticity to entropy. An example of this universal property is the adsorption of pyridine on Si(100): experimental results [5,6] showed that dative attachment which is the only way to preserve aromaticity is the favored form of adsorbed pyridine in dilute solutions at low temperatures. Recent theoretical analysis [7,8,9,10] confirmed that dative configuration is the more stable configuration. Two factors limit the domain of stability of dative attachment: repulsive interactions between dative bonds prevent them from being present in concentrated solutions while aromaticity contributes to a decrease in the entropy. When heated, aromatic grafted species (dative bonded) are converted to non-aromatic species which explains the vanishing of dative bonds at high temperatures even in dilute solutions as experimentally observed [5,6]. Since aromatic configuration is favored by lower enthalpy, this conversion appears driven by the entropic gain resulting from the vanishing of aromaticity [7,8] justifying here the deliberate use of possible negative contribution of aromaticity to entropy.
Recently, Casati et al. [11] show that when pressure is applied in aromatic compounds molecular stability which originates from the resonance between electronic configurations can be decreased in molecular crystals that gradually localizes one of the resonant configurations on compression.
This paper which is devoted to thermal switches used for circuit protection aims to go one step further by considering the combined effect of the temperature in addition to compression. Generally, thermal switches rely on a sharp drop in conductivity at the switch temperature, typically in the 300–500 K range. In the following, we shall discuss how the transition from an aromatic structure towards a non-aromatic structure can be leveraged to build molecular junctions for use as thermal switches. The negative contribution of aromaticity to entropy allows the existence of a transition temperature at which the entropy contribution compensates the energy contribution. Moreover, aromaticity is related to the cyclic delocalized distribution of π electrons [12]. Then, it may be expected that conductivity drops when aromaticity is lost at the transition temperature [13].

2. Thermal Switches

In order to discuss the relevance of this approach, let us use the results of earlier studies [7,8] for pyridine adsorbed on Si(100) and Ge(100), even if these materials are certainly not the easiest ones to handle. Of particular interest for our purpose is the spectroscopic study associated to DFT calculations of gold-molecule-gold junctions made by Yoshida et al. [14]. These authors have studied symmetric tolane-type molecules with pyridil as the anchoring group. Instead of gold, let us consider a tolane derived junction NC5H4–C≡C–C5H4N between two doped silicon samples. Two possible molecular devices are considered in the following: molecular junction with dative–dative (aromatic) attachment to Si(100) electrodes (see Figure 1) or dative–tetrapod (non-aromatic) attachment (see Figure 2).

3. Temperature Dependence of the Structure

It is of interest to plot the free energy (F = U − T·S) as a function of the distance between electrodes in both dative–dative and dative–tetrapod configurations.
The internal energy (U) has been evaluated as a function of the distance between electrodes (h). Electronic structure calculations have been completed with the GAMESS (US) software package. In the system considered here, the molecule is grafted between two semi-conducting silicon electrodes.
A Si21H20 cluster has been used to model the up electrode (with dative attachment). The top six silicon atoms compose the surface dimers while the remaining 15 Si atoms compose two subsurface layers which are hydrogen terminated to preserve the sp3 hybridization of the bulk diamond lattice. A Si15H18 cluster has been used to model the bottom electrode (with dative or tetrapod attachment). The top four Si atoms compose the surface dimers while the remaining eleven Si atoms compose two subsurface layers. We have considered different altitudes h (see Figure 3) between the two silicon electrodes in the range 19–22 Å.
Geometries of local minima on the potential energy surface were calculated at DFT (B3LYP) level of theory, with a 6–31 G* basis set for nitrogen, carbon and terminating hydrogens and a SBKJC pseudo-potential set for Si atoms. The geometry optimization consists to freeze the last silicon row of each electrode with full optimization of all the other atomic positions.
Our DFT simulations (see Figure 3) provide an estimate of the energy (the free energy at low temperature) for both configurations as function of h.
For large values of h, the dative–dative configuration is the most stable form. A minimum (EDD) of energy is found at h = 21.7 Å. At weaker altitudes the chain transforms to the dative–tetrapod configuration, with a minimum (EDT) of energy calculated at 20.3 Å.
It is interesting to note that both configuration present the same energy for h close to 21 Å.
The main point is: because of the dative–dative chain is longer than the dative–tetrapod one, an intermediate altitude at which both energy contributions are equal can be inferred. Our simplistic approach (zero point energy or extreme low-temperature conductivity effects are ignored here) gives an estimate of this altitude: knowledge of the exact value is not needed to further discuss their relative stability.
We assume the difference in entropy between aromatic (dative) and non-aromatic (tetrapod) configurations to be the one evaluated for pyridine in our earlier study [7,8] i.e., 17 cal·mol−1·K−1.
It is noteworthy that this small entropy difference is able to reverse the previous prediction derived from the energy difference alone. At moderate temperature (e.g., 300 K) the entropic term (−T·ΔS) is significant with respect to ΔU and, at fixed h = 21 Å (then fixed volume), the dative–tetrapod configuration is favored (see Figure 3).
This change of conformation of grafted molecule may have some consequences onto the local electronic conductivity across the molecular chain.
Analysis of the electronic density of the lowest unoccupied molecular orbital delocalized along the molecular chain is given in Figure 4. It appears that the conjugation within the molecular junction vanishes as the dative pattern is replaced by a tetrapod pattern and then, a change in conductivity is expected. At high temperature, the negative contribution of aromaticity to entropy favors non-aromatic structures; this property may be used to build thermal switches.

4. Summary and Conclusions

Let us briefly summarize our approach restricted to dative–dative and dative–tetrapod configurations of NC5H4–C≡C–C5H4N junction between two Si(100) electrodes. The key quantity is the difference in free energy, ΔF = ΔU − T·ΔS, between aromatic and non-aromatic forms. ΔU and ΔS were determined separately.
ΔU was evaluated from ab-initio (DFT) calculations. At large distances (h) between Si electrodes the dative–dative configuration is the most stable form. At weaker distances the dative–tetrapod becomes the most stable form. For a distance h between Si electrodes close to 21 Å the energy of the dative–dative configuration is found very close to the one of the dative–tetrapod configuration.
ΔS was assumed to be close to the contribution of aromaticity to entropy deduced for pyridine from thermodynamical measurements [7,8].
Then, the sign of ΔF was used to qualitatively predict the stabilization of dative–tetrapod configuration when temperature increases making possible to build thermal switches using the small negative contribution of aromaticity to entropy.

Acknowledgments

We gratefully acknowledge Y-Lan Boureau and Jean-Michel Mariot for a critical reading of the manuscript.

Author Contributions

All authors have read and approved the final manuscript; the authors contributed equally to this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aviram, A.; Ratner, M.A. Molecular rectifiers. Chem. Phys. Lett. 1974, 29, 277–283. [Google Scholar] [CrossRef]
  2. Ratner, M. A brief history of molecular electronics. Nat. Nanotechnol. 2013, 8, 378–381. [Google Scholar] [CrossRef] [PubMed]
  3. Akamatu, H.; Inokuchi, H.; Matsunaga, Y. Electrical Conductivity of the Perylene–Bromine Complex. Nature 1954, 173, 168–169. [Google Scholar] [CrossRef]
  4. Franklin, J.L. Calculation of Resonance Energies. J. Am. Chem. Soc. 1950, 72, 4278–4280. [Google Scholar] [CrossRef]
  5. Tao, F.; Qiao, M.H.; Wang, Z.H.; Xu, G.Q. Dative and Di−σ Binding States of Pyridine on Si(100) and Their Thermal Stability. J. Phys. Chem. B 2003, 107, 6384–6390. [Google Scholar] [CrossRef]
  6. Coustel, R.; Carniato, S.; Rochet, F.; Witkowski, N. Pyridine on Si(001)-(2×1): Density functional theory simulations compared with spectroscopic measurements. Phys. Rev. B 2012, 85, 035323. [Google Scholar] [CrossRef]
  7. Coustel, R.; Carniato, S.; Boureau, G. Thermodynamic factors limiting the preservation of aromaticity of adsorbed organic compounds on Si(100): Example of the pyridine. J. Chem. Phys. 2011, 134, 234708. [Google Scholar] [CrossRef] [PubMed]
  8. Coustel, R.; Carniato, S.; Boureau, G. Entropy, the Silent Killer of Aromaticity of Adsorbed Pyridine on Si(100) and Ge(100). J. Phys. Chem. C 2014, 118, 17505–17510. [Google Scholar] [CrossRef]
  9. Romeo, M.; Balducci, G.; Stener, M.; Fronzoni, G. N1s and C1s Near-Edge X-ray Absorption Fine Structure Spectra of Model Systems for Pyridine on Si(100): A DFT Simulation. J. Phys. Chem. C 2014, 118, 1049–1061. [Google Scholar] [CrossRef]
  10. Ng, W.K.H.; Liu, J.W.; Liu, Z.-F. Reaction Barriers and Cooperative Effects for the Adsorption of Pyridine on Si(100). J. Phys. Chem. C 2013, 117, 26644–26651. [Google Scholar] [CrossRef]
  11. Casati, N.; Kleppe, A.; Jephcoat, A.P.; Macchi, P. Putting pressure on aromaticity along with in situ experimental electron density of a molecular crystal. Nat. Commun. 2016, 7, 10901. [Google Scholar] [CrossRef] [PubMed]
  12. Poater, J.; Duran, M.; Solà, M.; Silvi, B. Theoretical Evaluation of Electron Delocalization in Aromatic Molecules by Means of Atoms in Molecules (AIM) and Electron Localization Function (ELF) Topological Approaches. Chem. Rev. 2005, 105, 3911–3947. [Google Scholar] [CrossRef] [PubMed]
  13. Wudl, F. Is there a relationship between aromaticity and conductivity? Pure Appl. Chem. 1982, 54, 1051–1058. [Google Scholar] [CrossRef]
  14. Yoshida, K.; Pobelov, I.V.; Manrique, D.Z.; Pope, T.; Mészáros, G.; Gulcur, M.; Bryce, M.R.; Lambert, C.J.; Wandlowski, T. Correlation of breaking forces, conductances and geometries of molecular junctions. Sci. Rep. 2015, 5, 9002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Dative–dative (aromatic) attachment of NC5H4–C≡C–C5H4N junction to Si(100).
Figure 1. Dative–dative (aromatic) attachment of NC5H4–C≡C–C5H4N junction to Si(100).
Entropy 18 00456 g001
Figure 2. Dative–tetrapod (non-aromatic) attachment of NC5H4–C≡C–C5H4N junction to Si(100).
Figure 2. Dative–tetrapod (non-aromatic) attachment of NC5H4–C≡C–C5H4N junction to Si(100).
Entropy 18 00456 g002
Figure 3. Free energy of dative–dative and dative–tetrapod configurations at 0 K as a function of the distance between electrodes. Free energy of dative–tetrapod configuration at 300 K was derived considering an additional entropic term −T·ΔS with ΔS = 17 cal·mol−1·K−1 [7,8] (dative–tetrapod configuration with h = 21.3 Å at 0 K has been arbitrarily chosen as reference).
Figure 3. Free energy of dative–dative and dative–tetrapod configurations at 0 K as a function of the distance between electrodes. Free energy of dative–tetrapod configuration at 300 K was derived considering an additional entropic term −T·ΔS with ΔS = 17 cal·mol−1·K−1 [7,8] (dative–tetrapod configuration with h = 21.3 Å at 0 K has been arbitrarily chosen as reference).
Entropy 18 00456 g003
Figure 4. Lowest unoccupied molecular orbital: dative–dative (a) and dative–tetrapod (b) configurations.
Figure 4. Lowest unoccupied molecular orbital: dative–dative (a) and dative–tetrapod (b) configurations.
Entropy 18 00456 g004

Share and Cite

MDPI and ACS Style

Coustel, R.; Carniato, S.; Boureau, G. A Possible Application of the Contribution of Aromaticity to Entropy: Thermal Switch. Entropy 2016, 18, 456. https://doi.org/10.3390/e18120456

AMA Style

Coustel R, Carniato S, Boureau G. A Possible Application of the Contribution of Aromaticity to Entropy: Thermal Switch. Entropy. 2016; 18(12):456. https://doi.org/10.3390/e18120456

Chicago/Turabian Style

Coustel, Romain, Stéphane Carniato, and Gérard Boureau. 2016. "A Possible Application of the Contribution of Aromaticity to Entropy: Thermal Switch" Entropy 18, no. 12: 456. https://doi.org/10.3390/e18120456

APA Style

Coustel, R., Carniato, S., & Boureau, G. (2016). A Possible Application of the Contribution of Aromaticity to Entropy: Thermal Switch. Entropy, 18(12), 456. https://doi.org/10.3390/e18120456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop