Stability trends in carbocation intermediates stemming from germacrene A and hedycaryol

• Naziha Tarannam,
• Prashant Kumar Gupta,
• Shani Zev and
• Dan Thomas Major

Beilstein J. Org. Chem. 2024, 20, 1189–1197, doi:10.3762/bjoc.20.101

• reaction trajectory towards this pathway and make sure that it stays along this pathway and does not rearrange to a similar, more stable carbocation. Here, using density functional theory (DFT) calculations, we provide a rationale for the relative stability of the intermediate carbocations formed from

• and hedycaryol were prepared from their corresponding SMILES string using the RDKit library [30]. Electronic structure calculations on these carbocations were performed using the hybrid DFT functional M06-2X [31] and with the range-separated hybrid meta-GGA functional ωB97M-V [32] with the 6-31+G(d,p

• ) basis set [33]. Additionally, to check the reliability of the energy calculations single point calculations were also performed with the larger def2-TZVPP [34] basis set for both functionals and to account for intramolecular non-covalent interactions, D3 dispersion corrections were added to the M06-2X

Two-fold addition reaction of silylene to C₆₀: structural and electronic properties of a bis-adduct

• Masahiro Kako,
• Masato Kai,
• Masanori Yasui,
• Michio Yamada,
• Yutaka Maeda and
• Takeshi Akasaka

Beilstein J. Org. Chem. 2024, 20, 1179–1188, doi:10.3762/bjoc.20.100

• silirane (silacyclopropane) derivative of fullerenes. Electrochemical measurements confirmed that the redox potentials of 3 are shifted cathodically compared to those of the parent mono-adduct 2. Density functional theory (DFT) calculations provided the basis for the electronic properties of compound 3

• -silylene adduct 3 of C60 based on spectroscopic measurements, X-ray crystallography, electrochemical analyses, and theoretical calculations. Results and Discussion Synthesis of bis-silylene adduct 3 The synthesis of the silylene adduct was conducted using a modified literature procedure [15]. A degassed

• . These results confirm that the two addends in 3 are both silirane structures, which represents the first example of a crystallographic analysis of silirane derivatives of fullerenes. Theoretical calculations Hirsch and co-workers reported the bis-functionalization of C60 using the Bingel–Hirsch, the

Manganese-catalyzed C–C and C–N bond formation with alcohols via borrowing hydrogen or hydrogen auto-transfer

• Mohd Farhan Ansari,
• Atul Kumar Maurya,
• Abhishek Kumar and
• Saravanakumar Elangovan

Beilstein J. Org. Chem. 2024, 20, 1111–1166, doi:10.3762/bjoc.20.98

• polyethyleneimine derivatives by coupling ethylene glycol and ethylenediamine using manganese-pincer catalyst Mn1 (1 mol %), t-BuOK (10 mol %) in toluene at 150 °C (Scheme 21) [48]. The mechanistic investigation based on the experimental and DFT calculations suggested a BH pathway. First, dehydrogenation of the

• amount of t-BuOK (20 mol %) in t-AmOH at 135 °C for 24 h allowed for the investigation of several alkyl and aryl ketones with amines and vicinal diols yielding good to excellent yields of the desired pyrroles (Scheme 70) with water and hydrogen gas being the only byproducts. DFT calculations suggested

Kinetically stabilized 1,3-diarylisobenzofurans and the possibility of preparing large, persistent isoacenofurans with unusually small HOMO–LUMO gaps

• Qian Liu and
• Glen P. Miller

Beilstein J. Org. Chem. 2024, 20, 1099–1110, doi:10.3762/bjoc.20.97

• Qian Liu Glen P. Miller Department of Chemistry, University of New Hampshire, 23 Academic Way, Durham, New Hampshire 03864-3598, USA 10.3762/bjoc.20.97 Abstract DFT calculations demonstrate that an isoacenofuran of any size possesses a smaller HOMO–LUMO gap than the corresponding acene bearing an

• . However, calculations indicate that the center ring of pentacene and other acenes is the most aromatic [2] and therefore resonance form P3 of Figure 1 is most significant. By contrast, isotetracenofuran and other isoacenofurans can only be drawn in one neutral, closed-shell resonance form in which only

• isopentacenofuran series (compounds 17–20). They possess HOMO–LUMO gaps between 1.25 and 1.36 eV, well below the 1.77 eV value calculated for hexacene. These Spartan ’20 [8] calculations confirm our expectation of a more highly delocalized π-system in any isoacenofuran (Figure 2) compared to the corresponding acene

Novel route to enhance the thermo-optical performance of bicyclic diene photoswitches for solar thermal batteries

• Akanksha Ashok Sangolkar,
• Rama Krishna Kadiyam and
• Ravinder Pawar

Beilstein J. Org. Chem. 2024, 20, 1053–1068, doi:10.3762/bjoc.20.93

• barriers were accounted using the DLPNO-CCSD(T) and (8,8)-CASPT2 methods, respectively. The photophysical properties including the absorption onset, excitation wavelengths, and the absorption intensities were extensively investigated with the time-dependent calculations. The result provides information on

• is associated with dissociation of the photoproduct. To validate the interconnection of the TS with the metastable photoisomer and the expected degradation product, IRC calculations were performed. The IRC calculations confirm that the obtained TS interlinks the metastable photoproduct with an

• undesired byproduct (Figure 3). Calculations reveal that the barrier for thermal dissociation of the photoproduct type-Id (209.43 kJ/mol) is significantly larger than the TBR barrier accounted for all the studied photoswitches. Nevertheless, the results reflect the propensity of the type-Id system to yield

Structure–property relationships in dicyanopyrazinoquinoxalines and their hydrogen-bonding-capable dihydropyrazinoquinoxalinedione derivatives

• Tural N. Akhmedov,
• Ajeet Kumar,
• Daken J. Starkenburg,
• Kyle J. Chesney,
• Khalil A. Abboud,
• Novruz G. Akhmedov,
• Jiangeng Xue and
• Ronald K. Castellano

Beilstein J. Org. Chem. 2024, 20, 1037–1052, doi:10.3762/bjoc.20.92

• illustrated in Figures S37–S41 (Supporting Information File 1). This serves as good evidence that the synthesized compounds exhibit no signs of aggregation in DMSO within the assessed concentration range (i.e., 2.5–45 μM). Electronic structure calculations Gas-phase computational investigations were conducted

• to estimate energies of the frontier molecular orbitals (FMOs), specifically the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), for DCPQs 1a–7a as well as 1b–7b. The density functional theory (DFT) calculations were performed using the B3LYP/6-31+G

• ). Cyclic voltammetric measurements of DCPQs, specifically 4a and 6a, and DPQDs 4b and 6b, support the overall trends observed through the electronic structure calculations. X-ray crystallography Single crystal X-ray diffraction data could be obtained for 2b, 5b, and 6b. Figure 6 illustrates the single

Evaluation of the enantioselectivity of new chiral ligands based on imidazolidin-4-one derivatives

• Jan Bartáček,
• Karel Chlumský,
• Jan Mrkvička,
• Lucie Paloušová,
• Miloš Sedlák and
• Pavel Drabina

Beilstein J. Org. Chem. 2024, 20, 684–691, doi:10.3762/bjoc.20.62

• tridentate chiral ligands: PyBOX [15] and PyBidine [16]. Additionally, we explored further modifications by substituting the pyridine moiety with an imidazole ring in ligand III (Figure 1), motivated by DFT calculations (see Supporting Information File 1, part S4) which confirmed that imidazole has a similar

• catalysed by copper(II) complex of ligand IV. Asymmetric aldol reactions of various aldehydes with ketones catalysed by compound IV. Supporting Information Supporting Information File 62: General information and experimental data of prepared compounds, copies of 1H and 13C NMR spectra and DFT calculations

Enhanced reactivity of Li⁺@C₆₀ toward thermal [2 + 2] cycloaddition by encapsulated Li⁺ Lewis acid

• Hiroshi Ueno,
• Yu Yamazaki,
• Hiroshi Okada,
• Fuminori Misaizu,
• Ken Kokubo and
• Hidehiro Sakurai

Beilstein J. Org. Chem. 2024, 20, 653–660, doi:10.3762/bjoc.20.58

• of this straightforward and selective reaction for Li+@C60. Results and Discussion We began by performing density functional theory (DFT) calculations to screen the substrates with suitable HOMO levels for the thermal [2 + 2] cycloaddition with Li+@C60. The structures of several kinds of possible

• a HOMO level of −6.08 eV or higher, such as styrene derivatives 1, 2, and 3. Based on the results of theoretical calculations, styrene derivatives 1, 2, and 3 were selected as the reactants for the thermal [2 + 2] cycloaddition with Li+@C60. For comparison, we also investigated the reaction of

HPW-Catalyzed environmentally benign approach to imidazo[1,2-a]pyridines

• Luan A. Martinho and
• Carlos Kleber Z. Andrade

Beilstein J. Org. Chem. 2024, 20, 628–637, doi:10.3762/bjoc.20.55

• protocol showed a lower E-factor by far, which makes it a cleaner and environmentally friendly process. Through the mass intensity parameter, our procedure also showed a lower mass loss compared to the total raw material used to produce a given mass of product. The calculations of these parameters can be

• conditions and Green Chemistry metrics for the GBB reactions under μw and conventional heating conditions.a Supporting Information Supporting Information File 44: Typical experimental procedures, FTIR, NMR and mass spectra of all compounds and GC metrics calculations. Funding The authors thank the

A laterally-fused N-heterocyclic carbene framework from polysubstituted aminoimidazo[5,1-b]oxazol-6-ium salts

• Andrew D. Gillie,
• Matthew G. Wakeling,
• Bethan L. Greene,
• Louise Male and
• Paul W. Davies

Beilstein J. Org. Chem. 2024, 20, 621–627, doi:10.3762/bjoc.20.54

• radius 3.5 Å, mesh spacing 0.10, H atoms removed for calculations. Colour coding represents positioning of steric bulk relative to the centre of the sphere, scale in Å. Use of AImOxAuCl 13 in catalysis. aYields are calculated from the 1H NMR spectra against an internal standard unless otherwise stated

Chemical and biosynthetic potential of Penicillium shentong XL-F41

• Ran Zou,
• Xin Li,
• Xiaochen Chen,
• Yue-Wei Guo and
• Baofu Xu

Beilstein J. Org. Chem. 2024, 20, 597–606, doi:10.3762/bjoc.20.52

• inapplicability of the NOESY experiment. To establish the absolute configuration of compound 1, electron circular dichroism (ECD) calculations were conducted using the time-dependent density functional theory (TDDFT) approach at a B3LYP/6-311G (d,p) (IEFPCM) level (Figure 4). Considering the uncertainty of the

Possible bi-stable structures of pyrenebutanoic acid-linked protein molecules adsorbed on graphene: theoretical study

• Yasuhiro Oishi,
• Motoharu Kitatani and
• Koichi Kusakabe

Beilstein J. Org. Chem. 2024, 20, 570–577, doi:10.3762/bjoc.20.49

• . The results of the nudged elastic band calculations show that the energy difference between the conformations is about 0.35 eV and the activation barrier on the path is about 0.55 eV (Figure 2) [9]. We will refer to the pathway on the potential energy surface with such an activation barrier as an

• . A substitution effect that does not effectively change the potential energy surface is isotopic substitution. In fact, the adiabatic potential surface defined for a hypothetical motion at the absolute zero is not changed by a change in nuclear mass. Even in density functional theory calculations

• , 2, and the conformation at the activation barrier top, respectively. Acknowledgements We are grateful to Prof. H. Ogi, Prof. M. Ohtani, Prof. N. Nakamura, Prof. K. Tanigaki, Prof. S. Hagiwara, and Prof. K. Nakajima for fruitful discussions. The calculations were done at the computer centers of

Synthesis of photo- and ionochromic N-acylated 2-(aminomethylene)benzo[b]thiophene-3(2Н)-ones with a terminal phenanthroline group

• Vladimir P. Rybalkin,
• Sofiya Yu. Zmeeva,
• Lidiya L. Popova,
• Irina V. Dubonosova,
• Olga Yu. Karlutova,
• Oleg P. Demidov,
• Alexander D. Dubonosov and
• Vladimir A. Bren

Beilstein J. Org. Chem. 2024, 20, 552–560, doi:10.3762/bjoc.20.47

• )methyl)-N-phenylacetamide was used as an actinometer for the quantum yield calculations (φfl = 0.60 ± 0.005) [14][33]. For preparative purposes, a Sweko IP65 led emitter (SUL-S1-20W-230-4000K-WH) was used. Spectral-fluorescent experiments were performed using solutions in acetonitrile or DMSO in quartz

• g/cm3, 17777 reflections measured (9.61° ≤ 2Θ ≤ 152.768°), 4053 unique (Rint = 0.0201, Rσ = 0.0153) that were used in all calculations. The final R1 was 0.0307 (I ≥ 2σ(I)) and wR2 was 0.0813 (all data). Reflections were recorded and unit cell parameters were determined and refined using the

Pseudallenes A and B, new sulfur-containing ovalicin sesquiterpenoid derivatives with antimicrobial activity from the deep-sea cold seep sediment-derived fungus Pseudallescheria boydii CS-793

• Zhen Ying,
• Xiao-Ming Li,
• Sui-Qun Yang,
• Hong-Lei Li,
• Xin Li,
• Bin-Gui Wang and
• Ling-Hong Meng

Beilstein J. Org. Chem. 2024, 20, 470–478, doi:10.3762/bjoc.20.42

• structures were solved by direct methods using the SHELXTL software package [18][19]. All non-hydrogen atoms were refined anisotropically. The H atoms were located by geometrical calculations, and their positions and thermal parameters were fixed during the structure refinement. The absolute structures were

Enhanced host–guest interaction between [10]cycloparaphenylene ([10]CPP) and [5]CPP by cationic charges

• Eiichi Kayahara,
• Yoshiyuki Mizuhata and
• Shigeru Yamago

Beilstein J. Org. Chem. 2024, 20, 436–444, doi:10.3762/bjoc.20.38

• , [10]CPP⊃[5]CPP2+. While the same host–guest complex consisted of neutral CPPs, [10]CPP⊃[5]CPP, was already reported, the cationic complex showed an about 20 times higher association constant in (CDCl2)2 at 25 °C (103 mol L−1). Electrochemical and photophysical analyses and theoretical calculations

• photophysical and electrochemical analyses and DFT calculations. Therefore, this CT is most likely the origin of the increased Ka value. We also discuss the charged double-layer structure, as determined by X-ray crystallographic analysis. Results and Discussion The size-complementary interaction between CPP2

• absorption bands in the visible region and characteristic broad absorption bands extending to the NIR region (≈2000 nm) (Figure 4, black line). According to time-dependent density functional theory (DFT) calculations at the ωB97-XD/6-31G* level of theory, the NIR absorption of the most thermodynamically

Facile approach to N,O,S-heteropentacycles via condensation of sterically crowded 3H-phenoxazin-3-one with ortho-substituted anilines

• Eugeny Ivakhnenko,
• Vasily Malay,
• Pavel Knyazev,
• Nikita Merezhko,
• Nadezhda Makarova,
• Oleg Demidov,
• Gennady Borodkin,
• Andrey Starikov and
• Vladimir Minkin

Beilstein J. Org. Chem. 2024, 20, 336–345, doi:10.3762/bjoc.20.34

• –H proton. Therefore, three tautomeric forms are possible for 5 (Scheme 4), one of which, the 7H-tautomer 7b, inevitably adopts a bipolar or biradical structure. According to the data from the DFT calculations performed at the B3LYP/6-311++G(d,p) approximation (Figure S6, Supporting Information File

• reaction products as 12H-quinoxaline[2,3-b]phenoxazines was confirmed through DFT calculations, X-ray crystallography, and NMR spectroscopy. Electronic absorption spectra (Table 2 and Figures 5–7) and electrochemical properties (Table 3) of the heteropentacyclic compounds 4a–h, 5a–c, 6a,b, and 10c revealed

• ). The DFT calculations [32] were performed using the Gaussian 16 program package [33] with the B3LYP functional [34] and the 6-311++G(d,p) basis set. The structures corresponding to minima on the potential energy surface and states with broken symmetry [20] were found through complete optimization of

Unveiling the regioselectivity of rhodium(I)-catalyzed [2 + 2 + 2] cycloaddition reactions for open-cage C₇₀ production

• Cristina Castanyer,
• Anna Pla-Quintana,
• Anna Roglans,
• Albert Artigas and
• Miquel Solà

Beilstein J. Org. Chem. 2024, 20, 272–279, doi:10.3762/bjoc.20.28

• opening. Keywords: cycloadditions; DFT calculations; [70]fullerene; open-cage fullerenes; rhodium; Introduction The discovery of C60 (buckminsterfullerene) in 1985 [1] initiated the search for possible technological applications of fullerenes. Nowadays, applications for these carbon-based molecules have

• group was substituted by a mesyl substituent and BIPHEP was used as a model phosphine ligand instead of Tol-BINAP to reduce the computational cost. The calculations, conducted at the B3LYP-D3/cc-pVTZ-PP(SMD=o-DCB)//B3LYP-D3/cc-pVDZ-PP level (see full computational details in Supporting Information File

• in both cases. This preference was rationalized by means of DFT calculations. Moreover, the photooxidation of one of the C–C double bonds of the new bis(fulleroids) afford open-cage C70 derivatives having a 12-membered ring opening. It is noteworthy to mention that examples of open-cage C70

Photochromic derivatives of indigo: historical overview of development, challenges and applications

• Gökhan Kaplan,
• Zeynel Seferoğlu and
• Daria V. Berdnikova

Beilstein J. Org. Chem. 2024, 20, 228–242, doi:10.3762/bjoc.20.23

• '-dinitrobenzoylindigo 18b were provided by de Melo and co-workers in 2020 using modern femtosecond spectroscopic techniques as well as TD-DFT calculations [60]. In 1989, Takahashi and co-workers reported an important study of geometry and conformational stability of indigo and its N,N'-disubstituted derivatives [61][62

• along with theoretical calculations [43]. It was found that N-aryl substituents in combination with the polar solvents result in the rapid non-radiative deactivation upon excitation of the initial E-isomer of indigo. At the same time, N-aryl substituents significantly increased the barrier for the

Perspectives on push–pull chromophores derived from click-type [2 + 2] cycloaddition–retroelectrocyclization reactions of electron-rich alkynes and electron-deficient alkenes

• Michio Yamada

Beilstein J. Org. Chem. 2024, 20, 125–154, doi:10.3762/bjoc.20.13

• also performed theoretical calculations, unveiling that the ring-opening step presented a higher energy barrier compared with the ring-formation step, particularly when the DCV compounds lacked aryl substituents. Subsequently, they successfully isolated the corresponding cyclobutene intermediate 13

• on the fullerene sphere. This was followed by a nucleophilic attack of the fullerene carbanion on one of the cyano groups, ultimately yielding 33. Computational calculations also suggested that the rearrangement was facilitated by the destabilization of the reactant due to steric interactions

• . The absolute configurations of the axially chiral TCBD and DCNQ derivatives were ascertained by a comparative analysis of the experimental CD spectra against the spectra derived from time-dependent density functional theory (TD-DFT) calculations. The axial chirality was stabilized by the steric

Electron-beam-promoted fullerene dimerization in nanotubes: insights from DFT computations

• Laura Abella,
• Gerard Novell-Leruth,
• Josep M. Ricart,
• Josep M. Poblet and
• Antonio Rodríguez-Fortea

Beilstein J. Org. Chem. 2024, 20, 92–100, doi:10.3762/bjoc.20.10

• pressures and high temperatures [13][14][15]. We aim to shed light in these reaction mechanisms and energy profiles by using complementary methodologies as standard density functional theory (DFT) calculations and first-principles Car–Parrinello molecular dynamics (CPMD) simulations. Firstly, we have

• (periodic calculations using VASP code) by comparing with other more standard procedures (non-periodic calculations using ADF code; see Computational methods for details) when computing the energy profile for the radical cation C60 + C60+. Besides, we computed the profile for the neutral dimerization C60

• -TEM images and movies will provide us a deeper understanding of fullerene coalescence processes in peapods. Computational Methods The Amsterdam Density Functional (ADF) code [21][22] was used for the electronic structure calculations and to optimize reactants, products, intermediates and transition

Optimizing reaction conditions for the light-driven hydrogen evolution in a loop photoreactor

• Pengcheng Li,
• Daniel Kowalczyk,
• Johannes Liessem,
• Mohamed M. Elnagar,
• Dariusz Mitoraj,
• Radim Beranek and
• Dirk Ziegenbalg

Beilstein J. Org. Chem. 2024, 20, 74–91, doi:10.3762/bjoc.20.9

• an exhaust gas line. To align the results using photocatalysts from different batches, an adjustment factor of 2.212 was calculated. The calculations are shown in Supporting Information File 1, Figure S8. Design of experiments A DOE analysis was used to evaluate the effect of four operating

Cycloaddition reactions of heterocyclic azides with 2-cyanoacetamidines as a new route to C,N-diheteroarylcarbamidines

• Pavel S. Silaichev,
• Tetyana V. Beryozkina,
• Vsevolod V. Melekhin,
• Valeriy O. Filimonov,
• Andrey N. Maslivets,
• Vladimir G. Ilkin,
• Wim Dehaen and
• Vasiliy A. Bakulev

Beilstein J. Org. Chem. 2024, 20, 17–24, doi:10.3762/bjoc.20.3

• (sigma) = 0.0978) which were used in all calculations. The final R1 was 0.0693 (I > 2σ(I)) and wR2 was 0.1849 (all data). The refined twin ratio was 0.6748(18):0.3252(18). The experiment was accomplished on the automated X-ray diffractometer «Xcalibur Ruby» with CCD detector following standard procedures

Studying specificity in protein–glycosaminoglycan recognition with umbrella sampling

• Mateusz Marcisz,
• Sebastian Anila,
• Margrethe Gaardløs,
• Martin Zacharias and
• Sergey A. Samsonov

Beilstein J. Org. Chem. 2023, 19, 1933–1946, doi:10.3762/bjoc.19.144

• analysis of the contacts between protein and GAG molecules established in the course of the simulation. Binding free energy calculations MM/GBSA (molecular mechanics generalized born surface area) model igb = 2 [48] from AMBER20 was used for free energy calculations on the trajectories obtained from RS

• manner more sound for the applied protocol, is to estimate PMF along the reaction coordinate. Such calculations also support the conclusion that the applied protocol did not succeed in returning the system to the original binding pose (Supporting Information File 1, Figures S9–S11). In comparison to the

Beyond n-dopants for organic semiconductors: use of bibenzo[d]imidazoles in UV-promoted dehalogenation reactions of organic halides

• Kan Tang,
• Megan R. Brown,
• Chad Risko,
• Melissa K. Gish,
• Garry Rumbles,
• Phuc H. Pham,
• Oana R. Luca,
• Stephen Barlow and
• Seth R. Marder

Beilstein J. Org. Chem. 2023, 19, 1912–1922, doi:10.3762/bjoc.19.142

• potentials of the halides that can be reduced in this way, quantum-chemical calculations, and steady-state and transient absorption spectroscopy suggest that UV irradiation accelerates the reactions via cleavage of the dimers to the corresponding radical monomers. Keywords: dehalogenation; n-dopant

• weak shoulder at ca. 400 nm (ε400 ≈ 150 M−1 cm−1), whereas (Cyc-DMBI)2 exhibits only a strong feature with a maximum at 327 nm (εmax = 13000 M−1 cm−1). TD-DFT calculations (M06/6-31G(d,p), isolated molecules) qualitatively reproduce the different behavior of the two dimers: for (Cyc-DMBI)2 the S0→S1

• longer-time spectra with Cyc-DMBI• having a blue-shifted absorption relative to N-DMBI•, consistent with the observed features arising from the monomeric radicals. Intersystem crossing (ISC) from the singlet to triplet excited state is also a possibility and TD-DFT calculations suggest that the dimer T1

Controlling the reactivity of La@C₈₂ by reduction: reaction of the La@C₈₂ anion with alkyl halide with high regioselectivity

• Yutaka Maeda,
• Saeka Akita,
• Mitsuaki Suzuki,
• Michio Yamada,
• Takeshi Akasaka,
• Kaoru Kobayashi and
• Shigeru Nagase

Beilstein J. Org. Chem. 2023, 19, 1858–1866, doi:10.3762/bjoc.19.138

• analysis of the single-bonded adduct 3a revealed the addition site of the p-methoxybenzyl group on La@C2v-C82. Theoretical calculations indicated that the addition site carbons in neutral La@C2v-C82 have high spin density, whereas those in the La@C2v-C82 anion do not have high charge densities. Thus, the

• @Ih-C80]2− = −1.62 V vs Fc+/Fc; C602− = −1.50 V vs Fc+/Fc)). The different reactivity of [M3N@Ih-C80]2− was explained by theoretical calculations. The charge density of the highest occupied molecular orbital (HOMO) was more highly localized on the fullerene cage for [Lu3N@Ih-C80]2−, whereas it was

• 2299232 (for 3a) contains the supplementary crystallographic data for this study. This data was provided free of charge by the joint Cambridge Crystallographic Data Centre and Fachinformationszentrum Karlsruhe Access Structures Service. Theoretical calculations POAV (θσπ−90º) and charge densities values