Next Article in Journal
Development of One-Step Non-Solvent Extraction and Sensitive UHPLC-MS/MS Method for Assessment of N-(n-Butyl) Thiophosphoric Triamide (NBPT) and N-(n-Butyl) Phosphoric Triamide (NBPTo) in Milk
Next Article in Special Issue
Computational Insights on the Potential of Some NSAIDs for Treating COVID-19: Priority Set and Lead Optimization
Previous Article in Journal
Chemical Composition and Antibacterial and Antioxidant Activity of a Citrus Essential Oil and Its Fractions
Previous Article in Special Issue
In Silico Studies of Some Isoflavonoids as Potential Candidates against COVID-19 Targeting Human ACE2 (hACE2) and Viral Main Protease (Mpro)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Deuteration and Temperature on Furan Ring Dynamics

1
Faculty of Chemistry, University of Wroclaw, Joliot–Curie 14, 50-383 Wroclaw, Poland
2
STC, Institute for Single Crystals, National Academy of Sciences of Ukraine, 60 Nauky Ave., 61001 Kharkiv, Ukraine
3
Faculty of Chemistry, V.N. Karazin Kharkiv National University, 4 Svobody Sq., 61077 Kharkiv, Ukraine
*
Authors to whom correspondence should be addressed.
Submission received: 5 April 2021 / Revised: 28 April 2021 / Accepted: 10 May 2021 / Published: 13 May 2021
(This article belongs to the Special Issue Molecular Modeling: Advancements and Applications)

Abstract

:
Despite significant progress in conformational analysis of cyclic molecules, the number of computational studies is still limited while most of that available in the literature data have been obtained long time ago with outdated methods. In present research, we have studied temperature driven conformational changes of the furan ring at three different temperatures. Additionally, the effect of deuteration on the ring dynamics is discussed; in addition, the aromaticity indices following the Bird and HOMA schemes are computed along all trajectories. Our ab initio molecular dynamic simulations revealed that deuteration has changed the furan ring dynamics and the obvious consequences; in addition, the shape and size of molecule are expected to be different.

Graphical Abstract

1. Introduction

Isotope effect and aromaticity are the key concepts in modern organic chemistry [1,2,3,4,5] with H/D isotope substitution being one of the fundamental aspects of NMR labeling, Neutron Scattering, and Mass Spectroscopy [6,7,8,9,10].
Furan is a cyclic aromatic ether found in cigarette smoke and used in the production of resins [11] and some agrochemicals and pharmaceuticals. Obtaining furan from “green” inedible lignocellulosic biomass is one of the promising sustainable, and industrially applicable alternatives to various petroleum-derived chemicals [12,13]. Recently, investigations on the presence of furan in coffee, juices, soups, and canned fruits or vegetables has gained widespread attention [14,15,16].
It is said that the presence of an aromatic conjugated system is responsible for bond equalization and most importantly the planarity of the ring. These statements seem to be in agreement with many experimental and theoretical studies [17,18,19] Therefore, for many decades, rings with aromatic character have been considered as flat and conformationally rigid molecules, until the recent theoretical studies [20,21,22,23] based on ab initio calculations, both at T = 0 K and with T > 0 K , that have shown the large amount of flexibility of aromatic molecules.
Our previous studies on isolated benzene molecule at 298 K exposed the fact that the population of this aromatic ring with planar geometry is very low (less than 7 % ). The molecule exists basically as a set of two families of flattened boat and half-boat conformations, which have average values of endocyclic dihedral angles of 7 , which allowed for explaining the results of experimental electronic spectra. Due to the high amount of symmetry of benzene molecules, electronic transitions should be forbidden However, out-of-plane deformations of benzene rings decrease the symmetry, allowing electronic transitions [24,25,26]. In light of these results, we asked the following questions: what is happening with the planarity of aromatic heterocyclic five–membered rings (i.e., furan), how does the temperature affect its dynamics, and how does deuteration of the ring influence its planarity?
Finally, we would like to point out that two competing mechanisms of odor detection—chemical and the spectral—are usually investigated by deuteration of odorants; see [27,28,29,30,31,32,33] for reference. However, our research exposed different perspectives on the influence of deuteration on the ring dynamics and thus also “volume/shape” of the molecule studied [34].

2. Results and Discussion

The full conformational space of five–membered ring includes pure 10 Envelope (E) and 10 Twist (T) conformations (and their intermediates) [35,36,37,38]. However, in our case, they can be reduced to six unique conformations due to the symmetry of the furan ring—see Supplementary Materials for the conformational circle.
Each single trajectory spanned around 100 ps, and, to asses changes in aromaticity, we have calculated two commonly used indices following the Bird and HOMA schemes [19,39,40] for each single step over all the trajectories. It is shown in Table 1.
It is seen that aromaticity decreases with growing temperature in both isotopomers, and the values of Bird and HOMA parameters are quite similar in two molecules at the same temperature. There is only one notable difference in the distribution (range) of values in both indices, it is wider for deuterofuran (FuD) at 50 K, and 300 K, but narrower at 500 K. It enables us to hypothesize that deuteration could be responsible for hampering the vibrations with maximal shortening or elongation of bonds.
Ring puckering parameters S defining degree of non–planarity and phase angle defining conformation type [41] were calculated for each single step of all trajectories. The frequencies of puckering amplitudes of the furan ring S obtained from our simulations at 50 K, 300 K and 500 K are shown in Figure 1. It is found that mean values of puckering amplitude S in deuterofuran are systematically lower as compared to furan in all cases. In addition, we observed a slight asymmetry in distribution over S for deuterofuran. At all temperatures, the median values are shifted to the area of smaller S as compared to furan, whilst the distribution width is almost the same. Thus, deuterofuran reveals long tails in the area of bigger S values. Therefore, the presence of D atoms makes the out-of-plane ring deformation less favorable, but the higher limit of flexibility is at the same level as in furan.
Simulations at the room temperature revealed that the population of planar conformation of the furan ring is around 39.5 % , which is much higher than that for benzene (cf. 6.6 % ), pyrimidine, or 1,2,4–triazine—cf. 30.1 % and 26.7 % , respectively [22]. However, this is to be expected since smaller cyclic molecules are geometrically more rigid. On the other hand, furan has much less aromatic character than benzene, and the population of planar geometry is much larger, which compels one to ask whether the concepts of aromatic organic chemistry still make sense in the light of the presented results.
Further analysis has found the distribution over the phase angle to be homogeneous in furan, and, surprisingly, it corresponds for deuterofuran to be more localized at three phase angles—0 , 180 and 360 , see Figure 2.
Thus, only deuterofuran tends to keep preferable conformations, even at higher temperatures, while furan does not demonstrate such localization and moves freely over all conformational space. The concurrent analysis of both puckering parameters did not reveal any significant difference between envelope and twist conformations in the areas of preferable conformations. All conformations that lead the oxygen atom out of the ring plane (first of all, Envelope O and Twist O–C α ) are more favorable than conformations that leave oxygen atoms in the plane.
The results of simulations show that the population of planar geometry for furan (deuterofuran) decreases significantly with the increase of the temperature: 99.85 (99.94)% versus 22.75 (32.27)% from 50 K to 500 K. The populations of six non–planar conformations are shown in Figure 3.

3. Materials and Methods

In order to investigate the conformational dynamics of furan and deuterated furan, we performed molecular dynamic calculations using ab initio molecular dynamics [42] using the efficient Car–Parrinello [43] propagation scheme as implemented in the CPMD program package [44]. The simulations were performed in the canonical ensemble at 50, 298, and 500 K. To control the temperature of the system, the NHC thermostat [45] was turned on and set at a frequency of 3000 (2500) cm 1 for furan (deuterated furan). A molecular dynamics time step of Δ t = 2 a.u. ( Δ t = 4 a.u.)—≈ 0.097 (≈ 0.048) fs—was used for the integration of the Car–Parrinello equations of motion using a fictitious mass parameter for the orbitals of 400 (700) a.u. together with the proper atomic masses for furan (deuterated furan). The supercell was a cubic box of 20.0 Å in length, and cluster boundary conditions were applied to properly treat the isolated system.
Conformational analysis was performed following the scheme developed by Zefirov, Palyulin, and Dashevskaya [41] that uses torsion angles instead of RMS displacements of atoms within the classical Cremer–Pople puckering formalism [35], in order to bypass some shortcomings of Cremer–Pople method related to non-equal bond lengths. There are a couple of different schemes for calculating puckering parameters; however, the difference between these methods is negligible for common heterocyclic rings like THF [46,47] and is likely negligible for furan. There are only two puckering parameters for a five-membered ring, one phase angle ψ that defines conformation, and puckering amplitude S that defines the degree of planarity of the ring. The phase angle ϕ that defines conformation and the puckering amplitude S that defines the degree of planarity of the ring are presented in the Supplementary Materials. For more details about the procedure used, please see [35,41] or the most recent paper by Chan [48].
For more details and equations, please see Supplementary Materials.

4. Conclusions

In this study, we investigated furan ring dynamics at three different temperatures also taking into account the effect of deuteration on ring puckering. It was found that furan is much more rigid than benzene and aromatic azines, while its mean degree of aromaticity is much lower. In addition, there is no strict relation: aromaticity vs. conformational dynamics of furan and deuterofuran at different temperatures. While degree of aromaticity is similar for both molecules, there are obvious differences in conformational behavior. It seems that, for molecules possessing an aromatic π -system, the particular degree of aromaticity does not simply predict the conformational dynamics and degree of non-planarity, especially at higher temperatures—other structural and dynamic conditions affect the activation of out-of-plane modes much more than aromatic π -conjugation. It is shown that a deuterated furan ring behaves differently in comparison to furan with puckering amplitude S being smaller, which strongly suggests that non-deuterated carbon atoms are ”volumetrically” larger. This is a very important conclusion that contradicts standard experimental assumptions that deuteration is leaving the size of molecules unchanged.

Supplementary Materials

The following are available online. Figure S1: The full conformational space of furan ring: pure 10 Envelopes (E) and 10 Twist (T) conformations and their intermediates. It is reduced to only 6 unique conformations due to the symmetry of furan ring see populations to all of them in Table S1, Figure S2: Distribution of I 5 indexes for furan and deuterofuran at three temperatures 50 K, 300 K and 500 K. The distributions at the same temperature for H-Furan and D-Furan are very similar and they simply overlap and to distinguish between them the y-scale has been changed that going up and down from ’0’ always positive values are shown and D-Furan is shown at bottom panel, Figure S3: Distribution of HOMA indexes for furan and deuterofuran at three temperatures 50 K, 300 K and 500 K. See description of Figure S2, Table S1: Population in % of all six conformations of the furan ring obtained from ab initio molecular dynamics at three temperatures 50 K, 300 K and 500 K. Note that the Envelope O and Twist C β –-C β conformations are doubly degenerated while all other are 4-fold degenerated, thus in order to compare the probabilities one must divide the percentage of each 4-fold degenerated conformations by 2.

Author Contributions

Conceptualization, P.D.; methodology, P.D. and I.V.O.; software, P.D. and I.V.O.; validation, Z.L.; writing—original draft preparation, P.D. and I.V.O.; writing—review and editing, P.D. and I.V.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Center Poland (2016/23/B/ST4/01099) (P.D.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was carried out with the support of the Interdisciplinary Center for Mathematical and Computational Modeling (ICM), University of Warsaw, under Grant No. G84-47 and at the Center of Informatics Tricity Academic Supercomputer & Network. The authors also gratefully acknowledge the Wroclaw Supercomputer Center (WCSS).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Von Ragué Schleyer, P. Introduction: Aromaticity. Chem. Rev. 2001, 101, 1115–1118. [Google Scholar] [CrossRef] [Green Version]
  2. Balaban, A.T.; Oniciu, D.C.; Katritzky, A.R. Aromaticity as a Cornerstone of Heterocyclic Chemistry. Chem. Rev. 2004, 104, 2777–2812. [Google Scholar] [CrossRef]
  3. Cyranski, M.K. Energetic Aspects of Cyclic Pi-Electron Delocalization: Evaluation of the Methods of Estimating Aromatic Stabilization Energies. Chem. Rev. 2005, 105, 3773–3811. [Google Scholar] [CrossRef] [PubMed]
  4. Krygowski, T.M.; Szatylowicz, H.; Stasyuk, O.A.; Dominikowska, J.; Palusiak, M. Aromaticity from the Viewpoint of Molecular Geometry: Application to Planar Systems. Chem. Rev. 2014, 114, 6383–6422. [Google Scholar] [CrossRef]
  5. Sola, M. Why Aromaticity Is a Suspicious Concept? Why? Front. Chem. 2017, 5, 22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Wiberg, K.B. The Deuterium Isotope Effect. Chem. Rev. 1955, 55, 713–743. [Google Scholar] [CrossRef]
  7. Kohen, A.; Limbach, H. Isotope Effects In Chemistry and Biology; CRC Press: Boca Raton, FL, USA, 2005. [Google Scholar]
  8. Bogle, X.S.; Singleton, D.A. Isotope-induced desymmetrization can mimic isotopic perturbation of equilibria. On the symmetry of bromonium ions and hydrogen bonds. J. Am. Chem. Soc. 2011, 133, 17172–17175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Swiderek, K.; Paneth, P. Binding Isotope Effects. Chem. Rev. 2013, 113, 7851–7879. [Google Scholar] [CrossRef]
  10. Perrin, C.L.; Shrinidhi, A.; Burke, K.D. Isotopic-Perturbation NMR Study of Hydrogen-Bond Symmetry in Solution: Temperature Dependence and Comparison of OHO and ODO Hydrogen Bonds. J. Am. Chem. Soc. 2019, 141, 17278–17286. [Google Scholar] [CrossRef]
  11. Hoydonckx, H.; Van Rhijn, W. Application of novel furan resins in composites. JEC Compos. Mag. 2008, 45, 34–35. [Google Scholar]
  12. Sanderson, K. Lignocellulose: A chewy problem. Nature 2011, 474, S12–S14. [Google Scholar] [CrossRef]
  13. Iroegbu, A.O.; Sadiku, E.R.; Ray, S.S.; Hamam, Y. Sustainable Chemicals: A Brief Survey of the Furans. Chem. Afr. 2020, 3, 481–496. [Google Scholar] [CrossRef] [Green Version]
  14. Moro, S.; Chipman, J.; Wegener, J.; Hamberger, C.; Dekant, W.; Mally, A. Furan in heat-treated foods: Formation, exposure, toxicity, and aspects of risk assessment. Mol. Nutr. Food Res. 2012, 56, 1197–1211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Waizenegger, J.; Winkler, G.; Kuballa, T.; Ruge, W.; Kersting, M.; Alexy, U.; Lachenmeier, D. Analysis and risk assessment of furan in coffee products targeted to adolescents. Food Addit. Contam. Part A 2012, 29, 19–28. [Google Scholar] [CrossRef] [PubMed]
  16. Alizadeh, M.; Barati, M.; Saleh-Ghadimi, S.; Roshanravan, N.; Zeinalian, R.; Jabbari, M. Industrial furan and its biological effects on the body systems. J. Food Biochem. 2018, 42, e12597. [Google Scholar] [CrossRef]
  17. Allcock, H.R. Recent advances in phosphazene (phosphonitrilic) chemistry. Chem. Rev. 1972, 72, 315–356. [Google Scholar] [CrossRef]
  18. Spackman, M.A. Molecular electric moments from x-ray diffraction data. Chem. Rev. 1992, 92, 1769–1797. [Google Scholar] [CrossRef]
  19. Krygowski, T.M.; Cyrański, M.K. Structural Aspects of Aromaticity. Chem. Rev. 2001, 101, 1385–1420. [Google Scholar] [CrossRef]
  20. Shishkin, O.V.; Pichugin, K.Y.; Gorb, L.; Leszczynski, J. Structural non-rigidity of six-membered aromatic rings. J. Mol. Struct. 2002, 616, 159–166. [Google Scholar] [CrossRef]
  21. Shishkin, O.V.; Dopieralski, P.; Omelchenko, I.V.; Gorb, L.; Latajka, Z.; Leszczynski, J. Dynamical Nonplanarity of Benzene. Evidences from the Car–Parrinello Molecular Dynamics Study. J. Phys. Chem. Lett. 2011, 2, 2881–2884. [Google Scholar] [CrossRef]
  22. Shishkin, O.V.; Dopieralski, P.; Omelchenko, I.V.; Gorb, L.; Latajka, Z.; Leszczynski, J. Entropy versus aromaticity in the conformational dynamics of aromatic rings. J. Mol. Model. 2013, 19, 4073–4077. [Google Scholar] [CrossRef]
  23. Omelchenko, I.V.; Shishkin, O.V.; Gorb, L.; Leszczynski, J. Aromaticity and conformational flexibility of five-membered monoheterocycles: Pyrrole-like and thiophene-like structures. Struct. Chem. 2016, 27, 101–109. [Google Scholar] [CrossRef]
  24. Berger, R.; Fischer, C.; Klessinger, M. Calculation of the Vibronic Fine Structure in Electronic Spectra at Higher Temperatures. 1. Benzene and Pyrazine. J. Phys. Chem. A 1998, 102, 7157. [Google Scholar] [CrossRef]
  25. Bernhardsson, A.; Forsberg, N.; Malmqvist, P.; Roos, B.O.; Serrano-Andres, L. A theoretical study of the 1B2u and 1B1u vibronic bands in benzene. J. Chem. Phys. 2000, 112, 2798–2809. [Google Scholar] [CrossRef] [Green Version]
  26. Li, J.; Lin, C.K.; Li, X.Y.; Zhu, C.Y.; Lin, S.H. Symmetry forbidden vibronic spectra and internal conversion in benzene. Phys. Chem. Chem. Phys. 2010, 12, 14967–14976. [Google Scholar] [CrossRef]
  27. Turin, L. A spectroscopic mechanism for primary olfactory reception. Chem. Senses 1996, 21, 773–791. [Google Scholar] [CrossRef] [Green Version]
  28. Franco, M.I.; Turin, L.; Mershin, A.; Skoulakis, E.M.C. Molecular vibration-sensing component in Drosophila melanogaster olfaction. Proc. Natl. Acad. Sci. USA 2011, 108, 3797–3802. [Google Scholar] [CrossRef] [Green Version]
  29. Block, E.; Jang, S.; Matsunami, H.; Sekharan, S.; Dethier, B.; Ertem, M.Z.; Gundala, S.; Pan, Y.; Li, S.; Li, Z.; et al. Implausibility of the vibrational theory of olfaction. Proc. Natl. Acad. Sci. USA 2015, 112, E2766–E2774. [Google Scholar] [CrossRef] [Green Version]
  30. Turin, L.; Gane, S.; Georganakis, D.; Maniati, K.; Skoulakis, E.M.C. Plausibility of the vibrational theory of olfaction. Proc. Natl. Acad. Sci. USA 2015, 112, E3154. [Google Scholar] [CrossRef] [Green Version]
  31. Block, E.; Jang, S.; Matsunami, H.; Batista, V.S.; Zhuang, H. Reply to Turin et al.: Vibrational theory of olfaction is implausible. Proc. Natl. Acad. Sci. USA 2015, 112, E3155. [Google Scholar] [CrossRef] [Green Version]
  32. Vosshall, L.B. Laying a controversial smell theory to rest. Proc. Natl. Acad. Sci. USA 2015, 112, 6525–6526. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Opitz, C.; Ahrné, E.; Goldie, K.N.; Schmidt, A.E.; Grzesiek, S. Deuterium induces a distinctive Escherichia coli proteome that correlates with the reduction in growth rate. J. Biol. Chem. 2019, 294, 2279–2292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Dunitz, J.D.; Ibberson, R.M. Is Deuterium Always Smaller than Protium? Angew. Chem. Int. Ed. 2008, 47, 4208–4210. [Google Scholar] [CrossRef]
  35. Cremer, D.; Pople, J.A. General definition of ring puckering coordinates. J. Am. Chem. Soc. 1975, 97, 1354–1358. [Google Scholar] [CrossRef]
  36. Cremer, D. Calculation of puckered rings with analytical gradients. J. Phys. Chem. 1990, 94, 5502–5509. [Google Scholar] [CrossRef]
  37. Ford, T.A. The evolution of the structural, vibrational and electronic properties of the cyclic ethers—On ring size. An ab initio study. J. Mol. Struct. 2014, 1073, 125–133. [Google Scholar] [CrossRef]
  38. Paoloni, L.; Rampino, S.; Barone, V. Potential-Energy Surfaces for Ring-Puckering Motions of Flexible Cyclic Molecules through Cremer–Pople Coordinates: Computation, Analysis, and Fitting. J. Chem. Theory Comput. 2019, 15, 4280–4294. [Google Scholar] [CrossRef]
  39. Bird, C.W. Heteroaromaticity, 5, a unified aromaticity index. Tetrahedron 1992, 48, 335–340. [Google Scholar] [CrossRef]
  40. Krygowski, T.M. Crystallographic studies of inter- and intramolecular interactions reflected in aromatic character of .pi.-electron systems. J. Chem. Inf. Comput. Sci. 1993, 33, 70–78. [Google Scholar] [CrossRef]
  41. Zefirov, N.S.; Palyulin, V.A.; Dashevskaya, E.E. Stereochemical studies. XXXIV. Quantitative description of ring puckering via torsional angles. The case of six-membered rings. J. Phys. Org. Chem. 1990, 3, 147–158. [Google Scholar] [CrossRef]
  42. Marx, D.; Hutter, J. Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  43. Car, R.; Parrinello, M. Unified Approach for Molecular Dynamics and Density-Functional Theory. Phys. Rev. Lett. 1985, 55, 2471–2474. [Google Scholar] [CrossRef] [Green Version]
  44. Copyright IBM Corp 1990-2019, Copyright MPI für Festkörperforschung Stuttgart 1997–2001. CPMD. Available online: http://www.cpmd.org (accessed on 23 March 2018).
  45. Martyna, G.J.; Klein, M.L.; Tuckerman, M. Nosé-Hoover chains: The canonical ensemble via continuous dynamics. J. Chem. Phys. 1992, 97, 2635–2643. [Google Scholar] [CrossRef]
  46. Altona, C.; Sundaralingam, M. Conformational analysis of the sugar ring in nucleosides and nucleotides. New description using the concept of pseudorotation. J. Am. Chem. Soc. 1972, 94, 8205–8212. [Google Scholar] [CrossRef]
  47. Han, S.J.; Kang, Y.K. Pseudorotation in heterocyclic five-membered rings: Tetrahydrofuran and pyrrolidine. J. Mol. Struct. Theochem 1996, 369, 157–165. [Google Scholar] [CrossRef]
  48. Chan, L.; Hutchison, G.R.; Morris, G.M. Understanding Ring Puckering in Small Molecules and Cyclic Peptides. J. Chem. Inf. Model. 2021, 61, 743–755. [Google Scholar] [CrossRef]
Figure 1. Relative frequency of puckering amplitude for furan (solid lines, furan-H) compared with deuterofuran (broken lines, furan-D) at three different temperatures: 50 K, 300 K, and 500 K derived from CPMD simulations. Inset represents studied molecules together with atoms labeling.
Figure 1. Relative frequency of puckering amplitude for furan (solid lines, furan-H) compared with deuterofuran (broken lines, furan-D) at three different temperatures: 50 K, 300 K, and 500 K derived from CPMD simulations. Inset represents studied molecules together with atoms labeling.
Molecules 26 02889 g001
Figure 2. Gradient maps over the puckering amplitude S and phase angle ψ , at three different temperatures: 50 K, 300 K, and 500 K derived from CPMD simulations for furan (furan–H) and deuterofuran (furan–D).
Figure 2. Gradient maps over the puckering amplitude S and phase angle ψ , at three different temperatures: 50 K, 300 K, and 500 K derived from CPMD simulations for furan (furan–H) and deuterofuran (furan–D).
Molecules 26 02889 g002
Figure 3. Furan and deuterofuran population of non-planar geometries decomposed into six different conformations: twist C β β (blue), envelope C β (brown), twist C α β (green), envelope C α , twist O–C α (red) and envelope O (dark blue) at three different temperatures: 50 K, 300 K, and 500 K derived from CPMD simulations. Non-planar populations at 50 K have been enlarged for visibility.
Figure 3. Furan and deuterofuran population of non-planar geometries decomposed into six different conformations: twist C β β (blue), envelope C β (brown), twist C α β (green), envelope C α , twist O–C α (red) and envelope O (dark blue) at three different temperatures: 50 K, 300 K, and 500 K derived from CPMD simulations. Non-planar populations at 50 K have been enlarged for visibility.
Molecules 26 02889 g003
Table 1. Aromaticity indices for furan (FuH) and deuterofuran (FuD) at three temperatures: 50 K, 300 K, and 500 K.
Table 1. Aromaticity indices for furan (FuH) and deuterofuran (FuD) at three temperatures: 50 K, 300 K, and 500 K.
I5HOMA
RangeMeanMedianRangeMeanMedian
FuH 50 K38.346.746.71.090.14 0.15
FuH 300 K104.741.941.93.42−0.10−0.06
FuH 500 K122.237.838.04.44−0.32−0.24
FuD 50 K42.646.746.71.150.14 0.15
FuD 300 K105.342.042.03.68−0.10−0.06
FuD 500 K113.637.938.04.36−0.32−0.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dopieralski, P.; Omelchenko, I.V.; Latajka, Z. Impact of Deuteration and Temperature on Furan Ring Dynamics. Molecules 2021, 26, 2889. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26102889

AMA Style

Dopieralski P, Omelchenko IV, Latajka Z. Impact of Deuteration and Temperature on Furan Ring Dynamics. Molecules. 2021; 26(10):2889. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26102889

Chicago/Turabian Style

Dopieralski, Przemyslaw, Iryna V. Omelchenko, and Zdzislaw Latajka. 2021. "Impact of Deuteration and Temperature on Furan Ring Dynamics" Molecules 26, no. 10: 2889. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26102889

Article Metrics

Back to TopTop