Next Article in Journal
Investigation of the Protective Effect for GcMAF by a Glycosidase Inhibitor and the Glycan Structure of Gc Protein
Next Article in Special Issue
The Gemstone Cyborg: How Diamond Films Are Creating New Platforms for Cell Regeneration and Biointerfacing
Previous Article in Journal
Indomethacin: Effect of Diffusionless Crystal Growth on Thermal Stability during Long-Term Storage
Previous Article in Special Issue
Synthesis of Triblock Polycarboxylate Superplasticizers with Well-Defined Structure and Its Dispersing Performance in β-Hemihydrate Gypsum
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 4
10.3390/molecules28041569
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structure and Vibrational Spectroscopy of C82 Fullerenol Valent Isomers: An Experimental and Theoretical Joint Study

by
Felix N. Tomilin
1,2,*,
Polina V. Artyushenko
2,3,4,
Irina A. Shchugoreva
3,4,
Anastasia V. Rogova
2,3,4,
Natalia G. Vnukova
1,2,
Grigory N. Churilov
1,2,
Nikolay P. Shestakov
1,2,
Olga N. Tchaikovskaya
5,
Sergei G. Ovchinnikov
1,2 and
Pavel V. Avramov
6,*
1
Kirensky Institute of Physics, Federal Research Center KSC SB RAS, Krasnoyarsk 660036, Russia
2
School of Non-Ferrous Metals and Materials Science, Siberian Federal University, Krasnoyarsk 660041, Russia
3
Laboratory for Digital Controlled Drugs and Theranostics, Federal Research Center Krasnoyarsk Scientific Center of the Siberian Branch of the RAS, Krasnoyarsk 660036, Russia
4
Laboratory for Biomolecular and Medical Technologies, Prof. V.F. Voino-Yasenetsky Krasnoyarsk State Medical University, Krasnoyarsk 660022, Russia
5
Department of Physics, Tomsk State University, Tomsk 634050, Russia
6
Department of Chemistry, Kyungpook National University, Daegu 41566, Republic of Korea
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(4), 1569; https://doi.org/10.3390/molecules28041569
Submission received: 27 December 2022 / Revised: 31 January 2023 / Accepted: 4 February 2023 / Published: 6 February 2023
(This article belongs to the Special Issue Feature Papers in Applied Chemistry)
Browse Figures
Review Reports Versions Notes

Abstract

:
Gd@C82OxHy endohedral complexes for advanced biomedical applications (computer tomography, cancer treatment, etc.) were synthesized using high-frequency arc plasma discharge through a mixture of graphite and Gd2O3 oxide. The Gd@C82 endohedral complex was isolated by high-efficiency liquid chromatography and consequently oxidized with the formation of a family of Gd endohedral fullerenols with gross formula Gd@C82O8(OH)20. Fourier-transformed infrared (FTIR) spectroscopy was used to study the structure and spectroscopic properties of the complexes in combination with the DFTB3 electronic structure calculations and infrared spectra simulations. It was shown that the main IR spectral features are formed by a fullerenole C82 cage that allows one to consider the force constants at the DFTB3 level of theory without consideration of gadolinium endohedral ions inside the carbon cage. Based on the comparison of experimental FTIR and theoretical DFTB3 IR spectra, it was found that oxidation of the C82 cage causes the formation of Gd@C82O28H20, with a breakdown of the integrity of the parent C82 cage with the formation of pores between neighboring carbonyl and carboxyl groups. The Gd@C82O6(OOH)2(OH)18 endohedral complex with epoxy, carbonyl and carboxyl groups was considered the most reliable fullerenole structural model.

Graphical Abstract

1. Introduction

Water-soluble fullerene derivatives attract considerable attention because of their potential for a variety of advanced biomedical applications. In particular, fullerenes reveal antiviral activity [1,2,3]. Shoji et. al. [2] found that C60 fullerene derivatives inhibit influenza A viral infection and can be considered as prospective candidates for the development of anti-influenza drugs. Water-soluble derivatives of fullerene C60 show a potential as biological antioxidants due to their chemically functionalized external surfaces by hydrophobic groups such as −OH and −COOH [4,5,6]. Fullerene derivatives can serve as drug delivery scaffolds with covalent or noncovalent linkages between the carbon cage and a bioactive moiety [7,8,9]. Highly iodinated fullerene derivatives [10] are of interest as potential imaging agents due to the high radioactive contrast of iodine ions. Gadolinium endohedral complexes of fullerenols are considered as prospective magnetic resonance imaging contrast agents [11,12,13,14,15]. Fullerene derivatives conjugated with single-stranded oligonucleotides called aptamers are considered a promising drug for cancer treatment [16] and several other diseases. The almost complete insolubility of pure fullerenes can be overcome by the attachment of various polar functional groups to the fullerene cages [17,18,19]. The presence and mutual arrangement of different functional groups affect the electronic structure, chemical activity and optical properties of fullerene derivatives [20,21,22].
Despite the desperate demand caused by a wide range of possible applications of such types of compounds in medicine, determination of the exact atomic structure of fullerene derivatives is still a challenging job. Density functional theory (DFT) calculations are a powerful tool to describe atomic and electronic structures and the stability of different derivatives of fullerenes [23,24,25,26], including endohedral metallofullerenes [27,28,29,30,31,32]. A joint experimental and theoretical approach proved to be a prospective method for the determination of substitution functional groups and their positions on fullerene cages, based on the comparison of experimental IR spectra of fullerene derivatives with theoretical ones. The dependence of the IR spectra upon the distribution of functional groups was successfully studied using DFT electronic structure calculations of hydroxylated [26,33] and oxidized [34] C60 derivatives. Since the vibration frequencies of the carbon cage and oxygen–containing groups lie in the range of 1000–4000 cm−1, this spectral energy region of fullerene derivatives is of special interest when determining the atomic structure of fullerenols by comparing the experimental and theoretical IR spectra [14,23]. For endohedral fullerenols, the guest metallic ions absorption bands do not fall into this region and do not need to be considered [14,35].
Among all types of fullerenes and their derivatives, with the number of carbon atoms in the parent carbon shells ranging from 16 to the hundreds [36], the C60 derivatives and complexes are the most studied [1,2,4,5,10,11,12,18] to date. However, of all endohedral fullerenes known so far, the C82 fullerene has become one of the most intensively studied cages [37,38,39]. Hydroxylation of a carbon cage is one of the common strategies used to solubilize the endohedral Gd@C82 complex [40]. Although considerable research has been devoted to the medical application of Gd@C82 and its derivatives [41,42,43], rather less attention has been paid to the determination of the atomic structure of functionalized Gd@C82 fullerene [22,44,45]. At the moment, the description of the atomic structure of fullerene is based on the prediction of the stability of various models by the method of density functional theory without comparisons with available experimental structural and spectroscopic data [22,44,45].
The Gd@C82OxHy complexes were successfully synthesized earlier [46], but their atomic structure has not been properly characterized yet. In a previous work [22], a number of atomistic models of endohedral fullerenols Gd@C82OxHy (x = 0 and 3; y = 8, 16, 24, 36 and 44) were proposed and studied using a combination of experimental mass spectroscopy techniques and density functional theory calculations. As a result of the combined experimental and theoretical investigation, the endohedral Gd@C82OxHy complex with 24 hydroxyl groups was considered the most prospective candidate for biomedical applications [22].
This study is devoted to the harsh–condition synthesis and joint experimental Fourier transform infrared (FTIR) absorption spectroscopy and the theoretical Density Functional Tight Binding (DFTB3) investigation of various Gd@C82OxHy complexes. The FTIR spectroscopy revealed an unexpected absorption band in the energy range from 1600 to 1800 cm−1 that could be assigned to double carbon–oxygen chemical bonds. For the sake of the proper interpretation of experimental FTIR data, a DFTB3 theoretical study of atomic and electronic structures and optical spectra of Gd@C82OxHy endohedral complexes was performed using the Density Functional Tight Binding approach. The structure of the proposed theoretical model of Gd@C82OxHy was found to be in good agreement with the experimental spectroscopic data. As a result, the influence of the key geometric parameters, the types of substitutional functional groups and their number in the fullerene cage were revealed by the DFTB3 method.

2. Results and Discussion

The harsh–condition synthesis of water–soluble fullerenols [46] is carried out under rather stringent conditions. This leads to a wide range of fullerenol molecules with different numbers and combinations of oxygen–containing functional groups. This, in turn, leads to difficulties in the determination of the atomic structure by means of X–ray electron diffraction. Therefore, various indirect methods should be and are used to verify the structures, which cannot give a direct and accurate answer about the atomic structure of the fullerenols.
It is reasonable to use infrared spectroscopy and X–ray photoelectron spectroscopy (XPS) to certify the water–soluble fullerenes. The FTIR spectroscopy method allows one to determine the functional groups, while XPS provides the opportunity to calculate the number of bond species of carbon and oxygen atoms. In the XPS fullerenol spectrum, the C1s bond energies extend to 5.3 eV and clearly deviate from the C1s peak of fullerene [47]. The large shift indicates a relatively high percentage of carbon with a higher degree of oxidation in fullerenol compared to fullerene. The wide range and spectral width also directly indicate that a significant number of isomers may be presented in the experimental mixture of the products. An analysis of the data provides the unique possibility to calculate the approximate atomic concentrations of different carbon and oxygen species with different bond types and oxidative states.
In the study, FTIR spectroscopy was used to determine the nature and coordination of the functional groups, since XPS does not provide detailed information about the diversity of functional groups. Certain bands in the experimental FTIR spectrum suggest that the C=O carbonyl, epoxy, −OH hydroxyl and −COOH carboxyl groups can be formed and coordinated to the surface of the fullerenol C82 carbon cage [47,48,49]. The atomistic theoretical models of possible C82 fullerenol isomers should take into account that the integrity of the parent fullerene cage could be broken by the appearance of C=O carbonyl and -OOH carboxyl groups coordinated to the surface of the molecule. Therefore, it is important to control that the fullerenol molecule does not break down, i.e., does not disintegrate into different species and fragments during geometry optimization in electronic structure calculations.
The approximate number of the functional groups was determined using XPS based on the proportion of carbon atoms chemically bonded to oxygen. Since the number of -OH hydroxyl groups attached to the fullerene must be even, the composition of a given fullerenol can be expressed as Gd@C82OxHy, where x:y ~ 1.5, so, in this work, the molecules C82O6(OOH)2(OH)18 and C82O8(OH)20 were constructed with an O:H ratio = 28:20. Since the main comparison of the theoretical models and experimental structures was done using infrared spectroscopy, the main goal in this work was to achieve the correct ratio of functional groups as it appears in the intensity of the corresponding spectral FTIR peaks. At the same time, gadolinium can be neglected in the calculation of IR spectra, since the concentration of one atom is very small compared to the rest, and, which could be even much more important, due to the huge Gadolinium atomic mass; the vibrations of gadolinium ion are in the long wavelength region of the spectrum, and its characteristic lines are poorly or not visible at all on the experimental spectrum. Therefore, the use of simpler methods such as the effective DFTB3 is more reasonable and faster.

2.1. FTIR Spectrum of Hydroxylated Gd Endohedral Complexes

The hydroxylated fullerene Gd endohedral complex Gd@C82 was synthesized and studied using FTIR spectroscopy. According to Reference [49], absorption bands from C–C (1554 cm−1), C=C (1621 cm–1), C=O (1703 cm−1) and C–O (1050–1150 cm−1) bonds were identified. The obtained experimental IR spectrum is presented in Figure 1 (black curve).
In order to reveal whether or not the previously studied structural model of Gd@C82O3(OH)24 [21] (Figure 2B) correctly represents the atomic structure of experimentally synthesized fullerenols, the corresponding theoretical DFTB3 (green curve) and experimental IR (black curve) spectra are compared in Figure 1. The experimental spectrum of the Gd@C82OxHy complex showed that the absorption bands provided by the C82 carbon cage lay off the bands of gadolinium ion (700–800 cm−1). Previously, it was shown that the presence of gadolinium ion slightly affected the carbon cage parameters [22]. Thus, it was assumed that excluding gadolinium ion from the theoretical model does not affect the part of the spectrum provided by the functional groups of the cage but can significantly simplify the calculations. Therefore, the gadolinium was removed from the model of Gd@C82O3(OH)24 (Figure 2C). Comparison of the FTIR experimental and DFTB3 theoretical IR spectra directly demonstrates the good agreement of the spectral lines in the spectral region 900–1400 cm−1, which reflects the epoxy and C–O groups’ vibrational modes, with a complete lack of intensity in the energy region (1400–1800 cm−1) assigned to the carboxylic and carbonyl groups and double and aromatic carbon–carbon bonds of the C82 carbon cage. Atomic coordinates of the C82O3(OH)24 structural model can be found in the Supplementary Materials Section.
Comparison of the experimental FTIR spectrum with the DFTB3 theoretical one shows that the structural model of C82O3(OH)24 does not satisfy the spectral features of hydroxylated fullerenes. For the C82O3(OH)24 spectral region for C–OH groups and carbon atoms of the cage is 1100–1350 cm−1, and for epoxy groups, the spectral region is around 850 cm−1 (Figure 1). In contrast to the theoretical spectrum, the experimental FTIR reveals a narrow very intense band of a complex shape in the spectral range from 1600 to 1800 cm−1. The spectral shoulder with the energy of 1650 cm−1 could be assigned to double C=C bond vibrations, whereas the main peak with the energy of 1750 cm−1 could be assigned to the C=O absorption band.
Fullerenes and their derivatives have polyhedral carbon cages in which sp2–carbons are directly bonded to three neighbors. Structurally, any fullerene cage cannot bear surface C=O groups without a breakdown of the carbon cage, since it violates the maximum carbon valence of four, and it could be assumed that the Gd@C82OxHy complex contained functional groups other than epoxy– and hydroxyls. The C=O bond’s presence in Gd@C82OxHy was confirmed by X–ray photoelectron spectroscopy as well [46]. Following the experimental FTIR spectrum, one can assume novel structural models of Gd@C82OxHy complexes with open–shell carbon cages caused by excess oxidation by oxygen during high–temperature synthesis.

2.2. Atomic Structure of C82O28H20 Valent Isomers

Following the experimental spectroscopic data of synthesized endohedral complexes and previous joint experimental and theoretical results [21], the Gd@C82 fullerene (Figure 2A) was considered the basis for a set of initial atomistic structural models. According to the experimental data (see the Experimental Methods Section), the ratio of hydroxyl groups to oxygen is 1:3, with the total number of oxygen atoms approximately equal to 28 [15]. Since the initial Gd@C82 fullerene belongs to the C2v group, 28 oxygen atoms were symmetrically distributed over the carbon cage. The oxygen atoms were added to the adjacent carbon atoms of six–membered rings (Figure 2E) and then transformed to hydroxyl, carbonyl and carboxyl functional groups. As a result of the modeling, two isomers of C82O28H20 containing C=O bonds were proposed to explain the peaks of experimental FTIR. The main difference between the isomers is how the C=O bond is coordinated to the C82 cage. The fundamental goal of the research is to study the main structural features of Gd@C82OxHy endohedral complexes and their spectroscopic characteristics.
The first C82O8(OH)20 valent isomer (I1) has 20 hydroxyl (−OH) and 8 carbonyl (C=O) groups (Figure 2F) distributed as eight carbonyls and hydroxyl pairs and six pairs of hydroxyls attached to the adjacent carbon atoms. Taking into account the limiting valence of a carbon atom equal to four sp3 C–C bonds, one could expect that the oxidation of a carbon atom to carbonyl must proceed with the breakdown of a chemical bond with one of its closest neighbors, with the consequent oxidation of the partner to at least hydroxyl and opening of the fullerene cage between the groups.
The second C82O6(OOH)2(OH)18 valent isomer (I2) has 18 hydroxyl (−OH), 6 carbonyl (C=O) and 2 carboxyl (−COOH) functional groups (Figure 2G). Two carboxyl groups were placed at opposite sides of the C82 cage. Carboxyl groups were formed by a carbon atom lying at the vertex of two neighboring hexagons and a pentagon. Thus, two initial C–C bonds were broken to form the −COOH carboxyl group, and two hydroxyl groups were attached to other ends of each broken C–C bond. Other functional groups were distributed in the same way as in I1. Similar to the C82O3(OH)24 model, gadolinium was removed from the new atomistic models. The atomic coordinates of both isomers can be found in the Supplementary Materials Section.
The atomic structures of both I1 and I2 isomers optimized at the DFTB3 level of theory are presented in Figure 2H,I. During the DFTB3 optimization, all carbonyl groups of the first isomer (Figure 2F) were transformed into epoxy groups keeping a closed parent C82 cage. Thus, the optimized structure of I1 has 20 hydroxyl (−OH) and 8 epoxy groups, in which the O atoms bridge the C–C bonds.
The DFTB3 equilibrium atomic structure of the I2 isomer is presented in Figure 2I. As a result of optimization, I2 has six pores in the cage, with two of them formed by a combination of two hexagons and one pentagon due to the formation of −COOH carboxyl groups. The mutual coordination of an I2 carboxyl and two hydroxyl groups prevented the formation of new C–O–C bonds between them with the formation of two epoxy groups. Each “carboxyl pore” has an adjacent pore emerging from the elongation of the distance between C=O and C–OH and combining two hexagons. Two other pores also formed by elongation of the distance between C=O and −OH and combining two hexagons on the C2v axis of the cage. Finally, the I2 isomer has an open C82 cage with 2 carboxyl, 2 epoxy, 4 carbonyl and 18 hydroxyl groups.

2.3. Theoretical IR Spectra of C82O28H20 Valent Isomers

The theoretical IR vibration spectra for both I1 and I2 isomers calculated at the DFTB3 level of theory are presented in Figure 1 with blue and red curves, respectively. Since both I1 and the previously proposed C82O3(OH)24 fullerenole are functionalized by hydroxyl and epoxy groups, one can see the similarity in their theoretical IR spectra, with some differences in the peak intensities caused by different amounts and different coordination of the hydroxyl and epoxy groups at the C82 cage. Comparison of the experimental FTIR spectrum (Figure 1) with the DFTB3 theoretical I1, one directly demonstrates a good agreement between the spectra in the 850–1350 cm−1 spectral region. In contrast, theoretical spectrum does not demonstrate any visible intensity in the carbonyl group spectral region (1750–1780 cm−1, Figure 1, blue curve), which is the most intensive region for the experimental one. One can conclude that the IR spectrum of the I1 isomer, which is characterized by the presence of eight epoxy groups, fails to describe the main spectral features of the Gd@C82OxHy complex.
The DFTB3 theoretical IR spectrum of the open−cage I2 isomer (Figure 1, red curve) demonstrates a much better agreement with the experimental FTIR one. The presence of the C=O carbonyl and −COOH carboxyl groups has a great impact on the IR spectrum of I2 in spectral region 1400–1800 cm−1. The DFTB3 calculated IR spectrum has characteristic bands at the 1600–1630 cm−1 (carbon cage), 1650–1670 cm−1 (vibrations of −COOH groups) and 1690–1720 cm−1 (C=O groups) (Figure 1, red curve) spectral regions. The good agreement of the DFTB3 theoretical I2 IR spectrum with the experimental FTIR one directly indicates the open−cage nature of the Gd@C82OxHy complex.

3. Methods

3.1. Experimental Methods

A mixture of Gd2O3 powder and graphite in a weighted ratio of 1:1 was processed in high–frequency arc plasma discharge by the sputtering of graphite electrodes with 3 mm axial holes [50]. A fullerene mixture was extracted from the carbon condensate by carbon disulfide in a Soxhlet apparatus. Using the experimental technique described in Reference [51], a mixture of Gd@C82 and higher fullerenes was isolated from the resulting solution, and the sample was dried and redissolved in toluene solvent. The Gd@C82 endohedral fullerene was isolated from the solution by high–efficiency liquid chromatography on an Agilent Technologies 1200 Series chromatograph. As it was shown [22] in the mass spectrum, only the main fraction of Gd@C82 is presented, well isolated from broadband noise. According to the method proposed in Reference [52], the −O and −OH groups were attached to the isolated endohedral metallofullerene. The number of these functional groups was calculated [46] using the X–Ray photoelectron Spectroscopy data from the fraction of carbon atoms chemically bonded to oxygen atoms. Taking into account that the number of −OH hydroxyl groups attached to the fullerene must be even [48], the composition of this product can be presented as Gd@C82Ox(OH)y (x = 10–12, y = 30–32, x + y = 40–42) [46]. The Fourier–transform infrared absorption spectrum (Figure 1, black curve) was recorded using the VERTEX 70 (Bruker Optik GMBH) spectrometer in the spectral region of 400–4000 cm−1 with a spectral resolution of 4 cm−1. To obtain the spectrum, round tablet samples of 0.5 mm thickness, 13 mm diameter and 0.140 g weight were prepared. The tablets were prepared as follows: Less than 0.00089 g sample weight of Gd@C82OxHy were thoroughly grounded with 0.14 g of KBr and subjected to cold pressing at 10,000 kg. The FTIR spectrometer was equipped by a global light source, wide band KBr beam splitter and RT–DLaTG detector (Bruker Optic GMBH).

3.2. Computational Details

The following computational procedure was used: (1) The fullerenol molecules were designed and constructed, (2) molecules were optimized to an equilibrium geometry and (3) the IR spectrum was calculated. The atomic structure and infrared spectra of free–standing C82O3(OH)24 hydroxylated fullerene and two valent isomers of C82O28H20, namely I1 (C82O8(OH)20) and I2 (C82O6(OOH)2(OH)18), were studied using the Density Functional Tight Binding [53] approach. To achieve high–quality vibrational spectra simulations and strict requirements for the structural optimization of low–dimensional clusters [54], the atomic structures of the Gd@C82O8(OH)20 and Gd@C82O6(OOH)2(OH)18 complexes (Figure 2) in the gas phase were optimized with the threshold of 10−4 Hartree/Bohr by DFTB3 using 3ob–3–1 parameters [53] using the GAMESS code [55]. Following the mandatory requirements of the Topology Conservation Theorem for low–dimensional lattices [54], the lowering or even the breakdown of high symmetry of low–dimensional atomic lattices due to any kind of internal mechanical stress must be taken into account. Due to this, all electronic structure calculations were performed without taking into account of any kind of point symmetry groups, and the C1 symmetry group was applied to all atomistic models of the C82 fullerene derivatives.
The DFTB3 is a fast and accurate parametrized quantum–mechanical method suitable to optimize extended molecules and clusters of dozens of atoms containing carbon, oxygen and hydrogen atoms [53] with high accuracy. The DFTB3 method is a very efficient approximation based on DFT parametrization, which provides a reasonably accurate atomic structure and force constants in comparison with the DFT and experimentally measured vibrational spectra [56] for a diverse range of molecules [57].

4. Conclusions

Chemically modified C82–based fullerenole endohedral complexes with Gd as a guest ion seem to be prospective agents for magnetic resonance imaging and cancer treatment due to the high contrast of Gd in magnetic resonance tomography, inert and protective carbon cages around Gd ions and chemical solubility in a biomedical environment. Gd@C82OxHy complexes were synthesized and investigated by joint experimental and theoretical techniques combining FTIR spectroscopy and theoretical DFTB3 electronic structure calculations. It was shown that consideration of the IR spectral mechanisms based on just the structure of an oxidized C82 cage to consider the force constants at the DFTB3 level of theory without consideration of gadolinium endohedral ion provides fundamental insights into the atomic structure and spectroscopic properties of Gd@C82OxHy complexes. A comparison of experimental FTIR and theoretical IR spectra of a number of C82 fullerenole isomers was performed for the sake of the determination of the atomic structure of a chemically modified Gd@C82OxHy cage. Based on the proposed theoretical structural models of C82OxHy fullerenols, theoretical DFTB3 IR spectra were calculated and compared with the experimental FTIR. It was shown that oxidation of the C82 cage with the formation of C82O28H20 should lead to the breakdown of some C–C bonds with the integrity breakdown of the parent C82 cage with the formation of pores between neighboring carbonyl and carboxyl groups. The Gd@C82O6(OOH)2(OH)18 endohedral complex with epoxy, carbonyl and carboxyl groups was considered as the most reliable fullerenole structural model.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28041569/s1: Atomic Cartesian coordinates for equilibrium structures of studied complexes of C82 fullerenols (in Å: X, Y, Z).

Author Contributions

Conceptualization, F.N.T., G.N.C., N.G.V. and P.V.A. (Pavel V. Avramov); methodology, O.N.T. and S.G.O.; software, P.V.A. (Polina V. Artyushenko), I.A.S. and A.V.R.; validation, G.N.C., N.P.S. and N.G.V.; formal analysis, N.G.V., N.P.S. and F.N.T.; investigation, F.N.T., P.V.A. (Polina V. Artyushenkoand), I.A.S., N.G.V. and A.V.R.; resources, S.G.O. and O.N.T.; data curation, F.N.T., P.V.A. (Polina V. Artyushenko) and G.N.C.; writing—original draft preparation, F.N.T., P.V.A. (Polina V. Artyushenko), N.G.V. and P.V.A. (Pavel V. Avramov); writing—review and editing, F.N.T., N.G.V., P.V.A. (Polina V. Artyushenko) and P.V.A. (Pavel V. Avramov); visualization, F.N.T.; supervision, F.N.T., S.G.O. and P.V.A. (Pavel V. Avramov); project administration, F.N.T. and S.G.O. and funding acquisition, O.N.T. and S.G.O. All authors have read and agreed to the published version of the manuscript.

Funding

Synthesis and spectroscopic study of the Gd@C82OxHy complexes were supported by the Ministry of Science and Higher Education of the Russian Federation under project FWES-2022-0005. Molecular design of the fullerene derivatives was supported by the National Research Foundation of the Republic of Korea, grant NRF 2021R1A2C1010455. DFTB3 electronic structure calculations were supported by Project FSWM-2020-0033 of the Russian Ministry of Science and Education.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available upon reasonable request from the corresponding authors.

Acknowledgments

The authors thank the JCSS Joint Super Computer Center of the Russian Academy of Sciences for providing supercomputers for computer simulations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dostalova, S.; Moulick, A.; Milosavljevic, V.; Guran, R.; Kominkova, M.; Cihalova, K.; Heger, Z.; Blazkova, L.; Kopel, P.; Hynek, D.; et al. Antiviral Activity of Fullerene C60 Nanocrystals Modified with Derivatives of Anionic Antimicrobial Peptide Maximin H5. Mon. Chem. Chem. Mon. 2016, 147, 905–918. [Google Scholar] [CrossRef]
  2. Shoji, M.; Takahashi, E.; Hatakeyama, D.; Iwai, Y.; Morita, Y.; Shirayama, R.; Echigo, N.; Kido, H.; Nakamura, S.; Mashino, T.; et al. Anti-Influenza Activity of C60 Fullerene Derivatives. PLoS ONE 2013, 8, e66337. [Google Scholar] [CrossRef]
  3. Marchesan, S.; Da Ros, T.; Spalluto, G.; Balzarini, J.; Prato, M. Anti-HIV Properties of Cationic Fullerene Derivatives. Bioorg. Med. Chem. Lett. 2005, 15, 3615–3618. [Google Scholar] [CrossRef] [PubMed]
  4. Galvan, Y.; Alperovich, I.; Zolotukhin, P.; Prazdnova, E.; Mazanko, M.; Belanova, A.; Chistyakov, V. Fullerenes as Anti-Aging Antioxidants. Curr. Aging Sci. 2017, 10, 56–67. [Google Scholar] [CrossRef]
  5. Sergeeva, V.; Kraevaya, O.; Ershova, E.; Kameneva, L.; Malinovskaya, E.; Dolgikh, O.; Konkova, M.; Voronov, I.; Zhilenkov, A.; Veiko, N.; et al. Antioxidant Properties of Fullerene Derivatives Depend on Their Chemical Structure: A Study of Two Fullerene Derivatives on HELFs. Oxidative Med. Cell. Longev. 2019, 2019, 4398695. [Google Scholar] [CrossRef]
  6. Roy, P.; Bag, S.; Chakraborty, D.; Dasgupta, S. Exploring the Inhibitory and Antioxidant Effects of Fullerene and Fullerenol on Ribonuclease A. ACS Omega 2018, 3, 12270–12283. [Google Scholar] [CrossRef]
  7. Sargazi, S.; ER, S.; Mobashar, A.; Gelen, S.S.; Rahdar, A.; Ebrahimi, N.; Hosseinikhah, S.M.; Bilal, M.; Kyzas, G.Z. Aptamer-Conjugated Carbon-Based Nanomaterials for Cancer and Bacteria Theranostics: A Review. Chem. Biol. Interact. 2022, 361, 109964. [Google Scholar] [CrossRef]
  8. Bolskar, R.D. Fullerenes for Drug Delivery BT—Encyclopedia of Nanotechnology; Bhushan, B., Ed.; Springer: Dordrecht, The Netherlands, 2016; pp. 1267–1281. [Google Scholar] [CrossRef]
  9. Al-Tamimi, B.H.; Farid, S.B.H. Fullerenes and Nanodiamonds for Medical Drug Delivery; Mallik, A.K., Ed.; IntechOpen: Rijeka, Croatia, 2021; Chapter 5. [Google Scholar] [CrossRef]
  10. Wharton, T.; Wilson, L.J. Highly-Iodinated Fullerene as a Contrast Agent For X-Ray Imaging. Bioorg. Med. Chem. 2002, 10, 3545–3554. [Google Scholar] [CrossRef]
  11. Grushko, Y.S.; Kozlov, V.S.; Sedov, V.P.; Kolesnik, S.G.; Lebedev, V.T.; Shilin, V.A.; Khodorkovsky, M.A.; Artamonova, T.O.; Shakhmin, A.L.; Shamanin, V.V.; et al. MRI-Contrasting System Based on Water-Soluble Fullerene/Gd-Metallofullerene Mixture. Fuller. Nanotub. Carbon Nanostruct. 2010, 18, 417–421. [Google Scholar] [CrossRef]
  12. Zhao, Z.; Zhen, M.; Zhou, C.; Li, L.; Jia, W.; Liu, S.; Li, X.; Liao, X.; Wang, C. A Gadofullerene Based Liver-Specific MRI Contrast Agent for an Early Diagnosis of Orthotopic Hepatocellular Carcinoma. J. Mater. Chem. B 2021, 9, 5722–5728. [Google Scholar] [CrossRef]
  13. Moreno-Vicente, A.; Mulet-Gas, M.; Dunk, P.W.; Poblet, J.M.; Rodríguez-Fortea, A. Probing the Formation of Halogenated Endohedral Metallofullerenes: Predictions Confirmed by Experiments. Carbon N. Y. 2018, 129, 750–757. [Google Scholar] [CrossRef]
  14. Zou, T.; Zhen, M.; Li, J.; Chen, D.; Feng, Y.; Li, R.; Guan, M.; Zhang, Y.; Shu, C.; Wang, C. The Effect of Hemiketals on the Relaxivity of Endohedral Gadofullerenols. RSC Adv. 2015, 5, 96253–96257. [Google Scholar] [CrossRef]
  15. Tóth, É.; Bolskar, R.D.; Borel, A.; González, G.; Helm, L.; Merbach, A.E.; Sitharaman, B.; Wilson, L.J. Water-Soluble Gadofullerenes:  Toward High-Relaxivity, PH-Responsive MRI Contrast Agents. J. Am. Chem. Soc. 2005, 127, 799–805. [Google Scholar] [CrossRef]
  16. Liu, Q.; Xu, L.; Zhang, X.; Li, N.; Zheng, J.; Guan, M.; Fang, X.; Wang, C.; Shu, C. Enhanced Photodynamic Efficiency of an Aptamer-Guided Fullerene Photosensitizer toward Tumor Cells. Chem.—Asian J. 2013, 8, 2370–2376. [Google Scholar] [CrossRef]
  17. Wakimoto, T.; Uchida, K.; Mimura, K.; Kanagawa, T.; Mehandjiev, T.R.; Aoshima, H.; Kokubo, K.; Mitsuda, N.; Yoshioka, Y.; Tsutsumi, Y.; et al. Hydroxylated Fullerene: A Potential Antiinflammatory and Antioxidant Agent for Preventing Mouse Preterm Birth. Am. J. Obstet. Gynecol. 2015, 213, 708.e1–708.e9. [Google Scholar] [CrossRef]
  18. Grebowski, J.; Konopko, A.; Krokosz, A.; DiLabio, G.A.; Litwinienko, G. Antioxidant Activity of Highly Hydroxylated Fullerene C60 and Its Interactions with the Analogue of α-Tocopherol. Free Radic. Biol. Med. 2020, 160, 734–744. [Google Scholar] [CrossRef]
  19. Semenov, K.N.; Andrusenko, E.V.; Charykov, N.A.; Litasova, E.V.; Panova, G.G.; Penkova, A.V.; Murin, I.V.; Piotrovskiy, L.B. Carboxylated Fullerenes: Physico-Chemical Properties and Potential Applications. Prog. Solid State Chem. 2017, 47–48, 19–36. [Google Scholar] [CrossRef]
  20. Yan, Q.-B.; Zheng, Q.-R.; Su, G. Theoretical Study on the Structures, Properties and Spectroscopies of Fullerene Derivatives C66X4 (X = H, F, Cl). Carbon N. Y. 2007, 45, 1821–1827. [Google Scholar] [CrossRef]
  21. Tang, C.; Chen, S.; Zhu, W.; Zhang, A.; Zhang, K.; Liu, M. Stabilities, Aromaticity, Infrared Spectra, and Optical Properties of Exohedral Fullerene Derivatives C76X18(X = H, F, Cl, and Br). Eur. Phys. J. D 2014, 68, 19. [Google Scholar] [CrossRef]
  22. Shakirova, A.A.; Tomilin, F.N.; Pomogaev, V.A.; Vnukova, N.G.; Churilov, G.N.; Kudryasheva, N.S.; Tchaikovskaya, O.N.; Ovchinnikov, S.G.; Avramov, P.V. Synthesis, Mass Spectroscopy Detection, and Density Functional Theory Investigations of the Gd Endohedral Complexes of C82 Fullerenols. Computation 2021, 9, 58. [Google Scholar] [CrossRef]
  23. Dawid, A.; Górny, K.; Gburski, Z. Water Solvent Effect on Infrared and Raman Spectra of C60(OH)24 Fullerenol Isomers: DFT Study. J. Phys. Chem. C 2017, 121, 2303–2315. [Google Scholar] [CrossRef]
  24. Wang, Z.; Gao, X.; Zhao, Y. Mechanisms of Antioxidant Activities of Fullerenols from First-Principles Calculation. J. Phys. Chem. A 2018, 122, 8183–8190. [Google Scholar] [CrossRef] [PubMed]
  25. Li, X.-J.; Yang, X.-H.; Song, L.-M.; Ren, H.-J.; Tao, T.-Z. A DFT Study on Structure, Stability, and Optical Property of Fullerenols. Struct. Chem. 2013, 24, 1185–1192. [Google Scholar] [CrossRef]
  26. Dawid, A.; Górny, K.; Gburski, Z. The Influence of Distribution of Hydroxyl Groups on Vibrational Spectra of Fullerenol C60(OH)24 Isomers: DFT Study. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2015, 136, 1993–1997. [Google Scholar] [CrossRef] [PubMed]
  27. Bai, H.; Gao, H.; Feng, W.; Zhao, Y.; Wu, Y. Interaction in Li@Fullerenes and Li+@Fullerenes: First Principle Insights to Li-Based Endohedral Fullerenes. Nanomaterials 2019, 9, 630. [Google Scholar] [CrossRef]
  28. Gao, H.; Sun, Y.; Zhang, J.; Wang, Q.; Wu, Y.; Bai, H. Understanding Endohedral Behaviors of Ten-Electron Atomic and Cluster System inside C60 from First-Principles. Phys. E Low-Dimens. Syst. Nanostruct. 2021, 127, 114532. [Google Scholar] [CrossRef]
  29. Dai, X.; Gao, Y.; Xin, M.; Wang, Z.; Zhou, R. The Ground State and Electronic Structure of Gd@ C82: A Systematic Theoretical Investigation of First Principle Density Functionals. J. Chem. Phys. 2014, 141, 244306. [Google Scholar] [CrossRef]
  30. Martínez-Flores, C.; Basiuk, V.A. Ln@C60 Endohedral Fullerenes: A DFT Analysis for the Complete Series from Lanthanum to Lutetium. Comput. Theor. Chem. 2022, 1217, 113878. [Google Scholar] [CrossRef]
  31. AlZahrani, A.Z. Cerium-Doped Endohedral Fullerene: A Density-Functional Theory Study. ISRN Condens. Matter Phys. 2012, 2012, 208234. [Google Scholar] [CrossRef]
  32. Senapati, L.; Schrier, J.; Whaley, K.B. Electronic Transport, Structure, and Energetics of Endohedral Gd@C82 Metallofullerenes. Nano Lett. 2004, 4, 2073–2078. [Google Scholar] [CrossRef]
  33. Velarde-Salcedo, M.V.; Gallo, M.; Guirado-López, R.A. Low Hydroxylated Fullerenes: Stability, Thermal Behavior, and Vibrational Properties. J. Phys. Chem. C 2018, 122, 13117–13129. [Google Scholar] [CrossRef]
  34. Palotás, J.; Martens, J.; Berden, G.; Oomens, J. Laboratory IR Spectra of the Ionic Oxidized Fullerenes C60O+ and C60OH+. J. Phys. Chem. A 2022, 126, 2928–2935. [Google Scholar] [CrossRef]
  35. Cui, R.; Li, J.; Huang, H.; Zhang, M.; Guo, X.; Chang, Y.; Li, M.; Dong, J.; Sun, B.; Xing, G. Novel Carbon Nanohybrids as Highly Efficient Magnetic Resonance Imaging Contrast Agents. Nano Res. 2015, 8, 1259–1268. [Google Scholar] [CrossRef]
  36. Shanbogh, P.P.; Sundaram, N.G. Fullerenes Revisited. Resonance 2015, 20, 123–135. [Google Scholar] [CrossRef]
  37. Wei, T.; Wang, S.; Liu, F.; Tan, Y.; Zhu, X.; Xie, S.; Yang, S. Capturing the Long-Sought Small-Bandgap Endohedral Fullerene Sc3N@C82 with Low Kinetic Stability. J. Am. Chem. Soc. 2015, 137, 3119–3123. [Google Scholar] [CrossRef]
  38. Roy, D.; Porwal, D.; Mukhopadhyay, K.; Ram, K. Theoretical Studies on C82 Fullerene with Encapsulated Germanium Metal. Synth. React. Inorg. Met. Nano-Metal Chem. 2006, 36, 43–46. [Google Scholar] [CrossRef]
  39. Shimada, T.; Ohno, Y.; Okazaki, T.; Sugai, T.; Suenaga, K.; Kishimoto, S.; Mizutani, T.; Inoue, T.; Taniguchi, R.; Fukui, N.; et al. Transport Properties of C78, C90 and Dy@C82 Fullerenes-Nanopeapods by Field Effect Transistors. Phys. E Low-Dimens. Syst. Nanostruct. 2004, 21, 1089–1092. [Google Scholar] [CrossRef]
  40. Mikawa, M.; Kato, H.; Okumura, M.; Narazaki, M.; Kanazawa, Y.; Miwa, N.; Shinohara, H. Paramagnetic Water-Soluble Metallofullerenes Having the Highest Relaxivity for MRI Contrast Agents. Bioconjug. Chem. 2001, 12, 510–514. [Google Scholar] [CrossRef]
  41. Meng, J.; Liang, X.; Chen, X.; Zhao, Y. Biological Characterizations of [Gd@C82(OH)22]n Nanoparticles as Fullerene Derivatives for Cancer Therapy. Integr. Biol. 2013, 5, 43–47. [Google Scholar] [CrossRef]
  42. Horiguchi, Y.; Kudo, S.; Nagasaki, Y. Gd@C82 Metallofullerenes for Neutron Capture Therapy—Fullerene Solubilization by Poly(Ethylene Glycol)-Block-Poly(2-(N, N-Diethylamino)Ethyl Methacrylate) and Resultant Efficacy in Vitro. Sci. Technol. Adv. Mater. 2011, 12, 44607. [Google Scholar] [CrossRef]
  43. Li, X.; Wang, C. The Potential Biomedical Platforms Based on the Functionalized Gd@C82 Nanomaterials. VIEW 2020, 1, e7. [Google Scholar] [CrossRef]
  44. Rodríguez-Zavala, J.G.; Tenorio, F.J.; Samaniego, C.; Méndez-Barrientos, C.I.; Peña-Lecona, F.G.; Muñoz-Maciel, J.; Flores-Moreno, R. Theoretical Study on the Sequential Hydroxylation of C82 Fullerene Based on Fukui Function. Mol. Phys. 2011, 109, 1771–1783. [Google Scholar] [CrossRef]
  45. Zhang, J.; Li, F.; Miao, X.; Zhao, J.; Jing, L.; Yang, G.; Jia, X. Most Stable Structures of Polyhydroxylated Endohedral Metallofullerene Gd@C82(OH)x (X = 1–24) from Density Function Theory. Chem. Phys. Lett. 2010, 492, 68–70. [Google Scholar] [CrossRef]
  46. Dudnik, A.I.; Vnukova, N.G.; Drokin, N.A.; Bondarev, V.S.; Shestakov, N.P.; Tomashevich, Y.V.; Churilov, G.N. Electrophysical Properties of Hydroxylated Endohedral Metallofullerene with Gadolinium. J. Phys. Chem. Solids 2019, 135, 109094. [Google Scholar] [CrossRef]
  47. Husebo, L.O.; Sitharaman, B.; Furukawa, K.; Kato, T.; Wilson, L.J. Fullerenols Revisited as Stable Radical Anions. J. Am. Chem. Soc. 2004, 126, 12055–12064. [Google Scholar] [CrossRef]
  48. Wang, B.-C.; Wang, H.-W.; Tso, H.-C.; Chen, T.-L.; Chou, Y.-M. Theoretical Studies of C70(OH)n (N = 14, 16, 18 and 20) Fullerenols. J. Mol. Struct. THEOCHEM 2002, 581, 177–186. [Google Scholar] [CrossRef]
  49. Li, J.; Wang, T.; Feng, Y.; Zhang, Y.; Zhen, M.; Shu, C.; Jiang, L.; Wang, Y.; Wang, C. A Water-Soluble Gadolinium Metallofullerenol: Facile Preparation, Magnetic Properties and Magnetic Resonance Imaging Application. Dalt. Trans. 2016, 45, 8696–8699. [Google Scholar] [CrossRef]
  50. Churilov, G.N.; Krätschmer, W.; Osipova, I.V.; Glushenko, G.A.; Vnukova, N.G.; Kolonenko, A.L.; Dudnik, A.I. Synthesis of Fullerenes in a High-Frequency Arc Plasma under Elevated Helium Pressure. Carbon N. Y. 2013, 62, 389–392. [Google Scholar] [CrossRef]
  51. Akiyama, K.; Hamano, T.; Nakanishi, Y.; Takeuchi, E.; Noda, S.; Wang, Z.; Kubuki, S.; Shinohara, H. Non-HPLC Rapid Separation of Metallofullerenes and Empty Cages with TiCl4 Lewis Acid. J. Am. Chem. Soc. 2012, 134, 9762–9767. [Google Scholar] [CrossRef]
  52. Chiang, L.Y.; Swirczewski, J.W.; Hsu, C.S.; Chowdhury, S.K.; Cameron, S.; Creegan, K. Multi-Hydroxy Additions onto C60 Fullerene Molecules. J. Chem. Soc. Chem. Commun. 1992, 24, 1791–1793. [Google Scholar] [CrossRef]
  53. Gaus, M.; Goez, A.; Elstner, M. Parametrization and Benchmark of DFTB3 for Organic Molecules. J. Chem. Theory Comput. 2013, 9, 338–354. [Google Scholar] [CrossRef]
  54. Avramov, P.V.; Kuklin, A.V. Topological and Quantum Stability of Low-Dimensional Crystalline Lattices with Multiple Nonequivalent Sublattices. New J. Phys. 2022, 24, 103015. [Google Scholar] [CrossRef]
  55. Barca, G.M.J.; Bertoni, C.; Carrington, L.; Datta, D.; De Silva, N.; Deustua, J.E.; Fedorov, D.G.; Gour, J.R.; Gunina, A.O.; Guidez, E.; et al. Recent Developments in the General Atomic and Molecular Electronic Structure System. J. Chem. Phys. 2020, 152, 154102. [Google Scholar] [CrossRef]
  56. Rüger, R.; Niehaus, T.; van Lenthe, E.; Heine, T.; Visscher, L. Vibrationally Resolved UV/Vis Spectroscopy with Time-Dependent Density Functional Based Tight Binding. J. Chem. Phys. 2016, 145, 184102. [Google Scholar] [CrossRef]
  57. Inakollu, V.S.S.; Yu, H. A Systematic Benchmarking of Computational Vibrational Spectroscopy with DFTB3: Normal Mode Analysis and Fast Fourier Transform Dipole Autocorrelation Function. J. Comput. Chem. 2018, 39, 2067–2078. [Google Scholar] [CrossRef]
Molecules 28 01569 g001 550
Figure 1. Experimental FTIR spectrum of Gd@C82OxHy and theoretical spectra of C82OxHy possible isomers. Top: Experimental FTIR spectrum of Gd@C82OxHy in the KBr matrix (black curve). Bottom: Theoretical IR spectra of free–standing C82O3(OH)24 (green curve). Middle: Theoretical IR spectrum of I1 isomer (C82O8(OH)20, blue curve), and theoretical IR spectrum of I2 isomer (C82O6(OOH)2(OH)18, red curve). Theoretical IR spectra calculated in the gas phase at the DFTB3 level of theory.
Figure 1. Experimental FTIR spectrum of Gd@C82OxHy and theoretical spectra of C82OxHy possible isomers. Top: Experimental FTIR spectrum of Gd@C82OxHy in the KBr matrix (black curve). Bottom: Theoretical IR spectra of free–standing C82O3(OH)24 (green curve). Middle: Theoretical IR spectrum of I1 isomer (C82O8(OH)20, blue curve), and theoretical IR spectrum of I2 isomer (C82O6(OOH)2(OH)18, red curve). Theoretical IR spectra calculated in the gas phase at the DFTB3 level of theory.
Molecules 28 01569 g001
Molecules 28 01569 g002 550
Figure 2. Some structural models of possible C82OxHy complexes. (A) Gd@C82; (B) Gd@C82O3(OH)24. (C,F,G) Initial structural models of C82O3(OH)24, where I1 is C82O8(OH)20 and I2 is C82O6(OOH)2(OH)18, valent isomers of the C82O28H20 complex, respectively. (D,H,I) Equilibrium geometries of C82O3(OH)24, I1 and I2 isomers of the C82O28H20 complex, respectively (atomic coordinates of I1 and I2 isomers can be found in the Supplementary Materials Section). Carbon, gadolinium, oxygen and hydrogen atoms are depicted in black, gray, red and white, respectively.
Figure 2. Some structural models of possible C82OxHy complexes. (A) Gd@C82; (B) Gd@C82O3(OH)24. (C,F,G) Initial structural models of C82O3(OH)24, where I1 is C82O8(OH)20 and I2 is C82O6(OOH)2(OH)18, valent isomers of the C82O28H20 complex, respectively. (D,H,I) Equilibrium geometries of C82O3(OH)24, I1 and I2 isomers of the C82O28H20 complex, respectively (atomic coordinates of I1 and I2 isomers can be found in the Supplementary Materials Section). Carbon, gadolinium, oxygen and hydrogen atoms are depicted in black, gray, red and white, respectively.
Molecules 28 01569 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tomilin, F.N.; Artyushenko, P.V.; Shchugoreva, I.A.; Rogova, A.V.; Vnukova, N.G.; Churilov, G.N.; Shestakov, N.P.; Tchaikovskaya, O.N.; Ovchinnikov, S.G.; Avramov, P.V. Structure and Vibrational Spectroscopy of C82 Fullerenol Valent Isomers: An Experimental and Theoretical Joint Study. Molecules 2023, 28, 1569. https://doi.org/10.3390/molecules28041569

AMA Style

Tomilin FN, Artyushenko PV, Shchugoreva IA, Rogova AV, Vnukova NG, Churilov GN, Shestakov NP, Tchaikovskaya ON, Ovchinnikov SG, Avramov PV. Structure and Vibrational Spectroscopy of C82 Fullerenol Valent Isomers: An Experimental and Theoretical Joint Study. Molecules. 2023; 28(4):1569. https://doi.org/10.3390/molecules28041569

Chicago/Turabian Style

Tomilin, Felix N., Polina V. Artyushenko, Irina A. Shchugoreva, Anastasia V. Rogova, Natalia G. Vnukova, Grigory N. Churilov, Nikolay P. Shestakov, Olga N. Tchaikovskaya, Sergei G. Ovchinnikov, and Pavel V. Avramov. 2023. "Structure and Vibrational Spectroscopy of C82 Fullerenol Valent Isomers: An Experimental and Theoretical Joint Study" Molecules 28, no. 4: 1569. https://doi.org/10.3390/molecules28041569

APA Style

Tomilin, F. N., Artyushenko, P. V., Shchugoreva, I. A., Rogova, A. V., Vnukova, N. G., Churilov, G. N., Shestakov, N. P., Tchaikovskaya, O. N., Ovchinnikov, S. G., & Avramov, P. V. (2023). Structure and Vibrational Spectroscopy of C82 Fullerenol Valent Isomers: An Experimental and Theoretical Joint Study. Molecules, 28(4), 1569. https://doi.org/10.3390/molecules28041569

Article Metrics

Back to TopTop