Insight into Markovnikov Reactions of Alkenes in Terms of Ab Initio and Molecular Face Theory

Electrophilic additions of hydrogen halides to alkenes in the gas phase are investigated with a high‐level ab initio method, MP2/6‐311+G(3df,2p). Based on this, the interesting features of these reactions along the IRC routes are characterized by the molecular face (MF) theory. For an alkene at the initial state, if the representative electron density (ED) encoded on the molecular face (MF) of the Markovnikov (M) carbon atom (the carbon with more hydrogen atoms) is larger than that of the anti‐Markovnikov (AM) carbon atom (the carbon with fewer hydrogen atoms), the electrophilic addition reaction is predicted to proceed along the Markovnikov addition route; in the reverse situtation, the anti‐Markovnikov addition route would be slightly preferred. It is then demonstrated that for a series of alkenes, the difference between activation energies of Markovnikov and anti‐Markovnikov addition routes [ΔE^#_(M?AM)] has a good linear correlation with sign(K_ED)K²_ED, where K_ED is characteristic of the electron density (ED) at the π region in the initial state of the alkenes. Interestingly, there is a good linear correlation between our sign (K_ED)K²_ED and the absolute values of difference in the core ionization energy between M and AM carbon atoms obtained by others (L. J. Sæthre, T. D. Thomas, S. Svensson J. Chem. Soc. P2 1997 , 2, 749.) in terms of the experimental study. In addition, the spatial dynamic changing features of the MF faces and interesting pictures of the electron transfer are clearly shown during the course of the electrophilic addition reactions. These results indicate that not only regioselectivity, but also activation energy and reactivity correlate with the π charge distribution in the initial state of the alkenes for electrophilic addition reactions.     

Valence state parameters of all transition metal atoms in metalloproteins—development of ABEEMσπ fluctuating charge force field

To promote accuracy of the atom‐bond electronegativity equalization method (ABEEMσπ) fluctuating charge polarizable force fields, and extend it to include all transition metal atoms, a new parameter, the reference charge is set up in the expression of the total energy potential function. We select over 700 model molecules most of which model metalloprotein molecules that come from Protein Data Bank. We set reference charges for different apparent valence states of transition metals and calibrate the parameters of reference charges, valence state electronegativities, and valence state hardnesses for ABEEMσπ through linear regression and least square method. These parameters can be used to calculate charge distributions of metalloproteins containing transition metal atoms (Sc‐Zn, Y‐Cd, and Lu‐Hg). Compared the results of ABEEMσπ charge distributions with those obtained by ab initio method, the quite good linear correlations of the two kinds of charge distributions are shown. The reason why the STO‐3G basis set in Mulliken population analysis for the parameter calibration is specially explained in detail. Furthermore, ABEEMσπ method can also quickly and quite accurately calculate dipole moments of molecules. Molecular dynamics optimizations of five metalloproteins as the examples show that their structures obtained by ABEEMσπ fluctuating charge polarizable force field are very close to the structures optimized by the ab initio MP2/6–311G method. This means that the ABEEMσπ/MM can now be applied to molecular dynamics simulations of systems that contain metalloproteins with good accuracy. © 2014 Wiley Periodicals, Inc.     

Gold(I)/Gold(III) Catalysis that Merges Oxidative Addition and π‐Alkene Activation

Heteroarylation of alkenes with aryl iodides was efficiently achieved with a (MeDalphos)AuCl complex through Au^I/Au^III catalysis. The possibility to combine oxidative addition of aryl iodides and π‐activation of alkenes at gold is demonstrated for the first time. The reaction is robust and general (>30 examples including internal alkenes, 5‐, 6‐, and 7‐membered rings). It is regioselective and leads exclusively to trans addition products. The (P,N) gold complex is most efficient with electron‐rich aryl substrates, which are troublesome with alternative photoredox/oxidative approaches. In addition, it provides a very unusual switch in regioselectivity from 5‐exo to 6‐endo cyclization between the Z and E isomers of internal alkenols.     

Preparation of σ- and π-Benzyl Type Derivatives of Molybdenum

The reaction of -LiM₀(CO)₂Cp with α-chloromethylnaphthalene gives a yellow, crystalline σ-benzyl type derivative of molybdenum (1). Irridation of (1) in n-hexane gives a red π-benzyl type derivative (2). Both complexes (1) and (2) are stable in air.     

Enolate Stabilization by Anion–π Interactions: Deuterium Exchange in Malonate Dilactones on π‐Acidic Surfaces

Of central importance in chemistry and biology, enolate chemistry is an attractive topic to elaborate on possible contributions of anion–π interactions to catalysis. To demonstrate the existence of such contributions, experimental evidence for the stabilization of not only anions but also anionic intermediates and transition states on π‐acidic aromatic surfaces is decisive. To tackle this challenge for enolate chemistry with maximal precision and minimal uncertainty, malonate dilactones are covalently positioned on the π‐acidic surface of naphthalenediimides (NDIs). Their presence is directly visible in the upfield shifts of the α‐protons in the ¹H NMR spectra. The reactivity of these protons on π‐acidic surfaces is measured by hydrogen–deuterium (H–D) exchange for 11 different examples, excluding controls. The velocity of H–D exchange increases with π acidity (NDI core substituents: SO₂R>SOR>H>OR>OR/NR₂>SR>NR₂). The H–D exchange kinetics vary with the structure of the enolate (malonates>methylmalonates, dilactones>dithiolactones). Moreover, they depend on the distance to the π surface (bridge length: 11–13 atoms). Most importantly, H–D exchange depends strongly on the chirality of the π surface (chiral sulfoxides as core substituents; the crystal structure of the enantiopure (R,R,P)‐macrocycle is reported). For maximal π acidity, transition‐state stabilizations up to ?18.8 kJ mol^?1 are obtained for H–D exchange. The Brønsted acidity of the enols increases strongly with π acidity of the aromatic surface, the lowest measured pK_a=10.9 calculates to a ΔpK_a=?5.5. Corresponding to the deprotonation of arginine residues in neutral water, considered as “impossible” in biology, the found enolate–π interactions are very important. The strong dependence of enolate stabilization on the unprecedented seven‐component π‐acidity gradient over almost 1 eV demonstrates quantitatively that such important anion–π activities can be expected only from strong enough π acids.     

Intramolecular π–π Interactions in Flexibly Linked Partially Fluorinated Bisarenes in the Gas Phase

Three compounds with phenyl and pentafluorophenyl rings bridged by (CH₂)₃ and (CH₂)₂SiMe₂ units were synthesized by hydrosilylation and C−C coupling reactions. Their solid‐state structures are dominated by intermolecular π stacking interactions, primarily leading to dimeric or chain‐type aggregates. Analysis of free molecules in the gas phase by electron diffraction revealed the most abundant conformer to be significantly stabilized by intramolecular π–π interactions. For the silicon compounds, structures characterized by σ–π interactions between methyl and pentafluorophenyl groups are second lowest in energy and cannot be excluded completely by the gas electron diffraction experiments. C₆H₅(CH₂)₃C₆F₅, in contrast, is present as a single conformer. The gas‐phase structures served as a reference for the evaluation of a series of (dispersion‐corrected) quantum‐chemical calculations.     

Fourfold alkyl wrapping of a copper(II) porphyrin thwarts macrocycle π–π stacking in a compact supramolecular package

Porphyrins are valuable constituents in optoelectronic, catalytic, and other applications, yet control of intermolecular π–π stacking is invariably essential to attain the desired properties. Superstructures built onto the porphyrin, often via meso‐aryl groups, can afford facial encumbrance that suppresses π–π stacking, although some molecular designs have provided insufficient facial coverage and many have entailed cumbersome syntheses. In this study, a copper(II) porphyrin bearing four meso substituents, namely, {10,20‐bis[2,6‐bis(octyloxy)phenyl]‐5,15‐dibromoporphinato}copper(II), [Cu(C₆₄H₈₂Br₂N₄O₄)], was prepared by metalation of the corresponding free‐base porphyrin and was characterized by single‐crystal X‐ray diffraction. The crystal structure reveals a dihedral angle of 111.1 (2)° for the plane of the meso‐aryl group relative to the plane of the porphyrin, with both aryl groups tilted in the same direction. Each of the four octyloxy groups exhibits a gauche conformation for the –OCH₂CH₂– unit but is extended with four or five anti (–CH₂CH₂–/H) conformations thereafter, causing each octyl group to span the dimension of the macrocycle. In a global frame of reference where the two Br atoms define the north/south poles and the two aryl groups are at antipodes on the equator, two octyl groups of one aryl unit project over the northern hemisphere (covering pyrroles A and B), whereas those of the other aryl unit project over the southern hemisphere (covering pyrroles C and D). Together, the four octyl groups ensheath the two faces of the porphyrin in a self‐wrapped assembly. The closest approach of the Cu atom to an octyl methylene C atom (position 6) is 3.5817 (18) Å, the mean separations of neighboring porphyrin planes are 8.059 (4) and 4.693 (8) Å along the a and c axes, respectively, and the center‐to‐center distances between the Cu atoms of neighboring porphyrins are 10.2725 (4), 12.2540 (6), and 12.7472 (6) Å along the a, b, and c axes, respectively. The Hirshfeld surface analysis and two‐dimensional (2D) fingerprint plots provide information concerning contact interactions in the supramolecular assembly of the solid crystal.     

Enantio‐ and Site‐Selective α‐Fluorination of N‐Acyl 3,5‐Dimethylpyrazoles Catalyzed by Chiral π–CuII Complexes

Catalytic enantioselective α‐fluorination reactions of carbonyl compounds are among the most powerful and efficient synthetic methods for constructing optically active α‐fluorinated carbonyl compounds. Nevertheless, α‐fluorination of α‐nonbranched carboxylic acid derivatives is still a big challenge because of relatively high pK_a values of their α‐hydrogen atoms and difficulty of subsequent synthetic transformation without epimerization. Herein we show that chiral copper(II) complexes of 3‐(2‐naphthyl)‐l ‐alanine‐derived amides are highly effective catalysts for the enantio‐ and site‐selective α‐fluorination of N‐(α‐arylacetyl) and N‐(α‐alkylacetyl) 3,5‐dimethylpyrazoles. The substrate scope of the transformation is very broad (25 examples including a quaternary α‐fluorinated α‐amino acid derivative). α‐Fluorinated products were converted into the corresponding esters, secondary amides, tertiary amides, ketones, and alcohols with almost no epimerization in high yield.     

Ab initio study of the π–π interactions between CO2 and benzene,pyridine, and pyrrole

The π–π interactions between CO₂ and three aromatic molecules, namely benzene (C₆H₆), pyridine (C₅H₅N), and pyrrole (C₄H₅N), which represent common functional groups in metal‐organic/zeoliticimidazolate framework materials, were characterized using high‐level ab initio methods. The coupled‐cluster with single and double excitations and perturbative treatment of triple excitations (CCSD(T)) method with a complete basis set (CBS) was used to calibrate Hartree–Fock, density functional theory, and second‐order M?ller–Plesset (MP2) with resolution of the identity approximation calculations. Results at the MP2/def2‐QZVPP level showed the smallest deviations (only about 1 kJ/mol) compared with those at the CCSD(T)/CBS level of theory. The strength of π–π binding energies (BEs) followed the order C₄H₅N > C₆H₆ ～ C₅H₅N and was roughly correlated with the aromaticity and the charge transfer between CO₂ and aromatic molecule in clusters. Compared with hydrogen‐bond or electron donor–acceptor interactions observed during BE calculations at the MP2/def2‐QZVPP level of theory, π–π interactions significantly contribute to the total interactions between CO₂ and aromatic molecules. © 2013 Wiley Periodicals, Inc.     

Comparative studies on group III σ‐hole and π‐hole interactions

The σ‐hole of M₂H₆ (M = Al, Ga, In) and π‐hole of MH₃ (M = Al, Ga, In) were discovered and analyzed, the bimolecular complexes M₂H₆···NH₃ and MH₃···N₂P₂F₄ (M = Al, Ga, In) were constructed to carry out comparative studies on the group III σ‐hole interactions and π‐hole interactions. The two types of interactions are all partial‐covalent interactions; the π‐hole interactions are stronger than σ‐hole interactions. The electrostatic energy is the largest contribution for forming the σ‐hole and π‐hole interaction, the polarization energy is also an important factor to form the M···N interaction. The electrostatic energy contributions to the interaction energy of the σ‐hole interactions are somewhat greater than those of the π‐hole interactions. However, the polarization contributions for the π‐hole interactions are somewhat greater than those for the σ‐hole interactions. © 2016 Wiley Periodicals, Inc.     

Gold(III) π‐Allyl Complexes

Gold(III) π‐complexes have been authenticated recently with alkenes, alkynes, and arenes. The key importance of Pd^II π‐allyl complexes in organometallic chemistry (Tsuji–Trost reaction) prompted us to explore gold(III) π‐allyl complexes, which have remained elusive so far. The (P,C)Au^III(allyl) and (methallyl) complexes 3 and 3′ were readily prepared and isolated as thermally and air‐stable solids. Spectroscopic and crystallographic analyses combined with detailed DFT calculations support tight quasi‐symmetric η³‐coordination of the allyl moiety. The π‐allyl gold(III) complexes are activated towards nucleophilic additions, as substantiated with β‐diketo enolates.     

Structure determination of three furan‐substituted benzimidazoles and calculation of π–π and C—H...π interaction energies

The synthesis and structural characterization of 2‐(furan‐2‐yl)‐1‐(furan‐2‐ylmethyl)‐1H‐benzimidazole [C₁₆H₁₂N₂O₂, (I)], 2‐(furan‐2‐yl)‐1‐(furan‐2‐ylmethyl)‐1H‐benzimidazol‐3‐ium chloride monohydrate [C₁₆H₁₃N₂O₂⁺·Cl⁻·H₂O, (II)] and the hydrobromide salt 5,6‐dimethyl‐2‐(furan‐2‐yl)‐1‐(furan‐2‐ylmethyl)‐1H‐benzimidazol‐3‐ium bromide [C₁₈H₁₇N₂O₂⁺·Br⁻, (III)] are described. Benzimidazole (I) displays two sets of aromatic interactions, each of which involves pairs of molecules in a head‐to‐tail arrangement. The first, denoted set (Ia), exhibits both intermolecular C—H...π interactions between the 2‐(furan‐2‐yl) (abbreviated as Fn) and 1‐(furan‐2‐ylmethyl) (abbreviated as MeFn) substituents, and π–π interactions involving the Fn substituents between inversion‐center‐related molecules. The second, denoted set (Ib), involves π–π interactions involving both the benzene ring (Bz) and the imidazole ring (Im) of benzimidazole. Hydrated salt (II) exhibits N—H...OH₂...Cl hydrogen bonding that results in chains of molecules parallel to the a axis. There is also a head‐to‐head aromatic stacking of the protonated benzimidazole cations in which the Bz and Im rings of one molecule interact with the Im and Fn rings of adjacent molecules in the chain. Salt (III) displays N—H...Br hydrogen bonding and π–π interactions involving inversion‐center‐related benzimidazole rings in a head‐to‐tail arrangement. In all of the π–π interactions observed, the interacting moieties are shifted with respect to each other along the major molecular axis. Basis set superposition energy‐corrected (counterpoise method) interaction energies were calculated for each interaction [DFT, M06‐2X/6‐31+G(d)] employing atomic coordinates obtained in the crystallographic analyses for heavy atoms and optimized H‐atom coordinates. The calculated interaction energies are −43.0, −39.8, −48.5, and −55.0 kJ mol⁻¹ for (Ia), (Ib), (II), and (III), respectively. For (Ia), the analysis was used to partition the interaction energies into the C—H...π and π–π components, which are 9.4 and 24.1 kJ mol⁻¹, respectively. Energy‐minimized structures were used to determine the optimal interplanar spacing, the slip distance along the major molecular axis, and the slip distance along the minor molecular axis for 2‐(furan‐2‐yl)‐1H‐benzimidazole.     

A comparative experimental and theoretical investigation of hydrogen‐bond,halogen‐bond and π–π interactions in the solid‐state supramolecular assembly of 2‐ and 4‐formylphenyl arylsulfonates

To explore the operational role of noncovalent interactions in supramolecular architectures with designed topologies, a series of solid‐state structures of 2‐ and 4‐formylphenyl 4‐substituted benzenesulfonates was investigated. The compounds are 2‐formylphenyl 4‐methylbenzenesulfonate, C₁₄H₁₂O₄S, 3a , 2‐formylphenyl 4‐chlorobenzenesulfonate, C₁₃H₉ClO₄S, 3b , 2‐formylphenyl 4‐bromobenzenesulfonate, C₁₃H₉BrO₄S, 3c , 4‐formylphenyl 4‐methylbenzenesulfonate, C₁₄H₁₂O₄S, 4a , 4‐formylphenyl 4‐chlorobenzenesulfonate, 4b , C₁₃H₉ClO₄S, and 4‐formylphenyl 4‐bromobenzenesulfonate, C₁₃H₉BrO₄S, 4c . The title compounds were synthesized under basic conditions from salicylaldehyde/4‐hydroxybenzaldehydes and various aryl sulfonyl chlorides. Remarkably, halogen‐bonding interactions are found to be important to rationalize the solid‐state crystal structures. In particular, the formation of O…X (X = Cl and Br) and type I X…X halogen‐bonding interactions have been analyzed by means of density functional theory (DFT) calculations and characterized using Bader's theory of `atoms in molecules' and molecular electrostatic potential (MEP) surfaces, confirming the relevance and stabilizing nature of these interactions. They have been compared to antiparallel π‐stacking interactions that are formed between the arylsulfonates.     

“Assessing the RAFT Equilibrium Constant via Model Systems: An EPR Study”—Response to a Comment

We have presented an EPR‐based approach for deducing the RAFT equilibrium constant, K_eq, of a dithiobenzoate‐mediated system [Meiser, W. and Buback M. Macromol. Rapid Commun. 2011 , 32, 1490]. Our value is by four orders of magnitude below K_eq from ab initio calculations for the identical monomer‐free system. Junkers et al. [Macromol. Rapid Commun. 2011 , 32, 1891] claim that our EPR approach would be model dependent and our data could be equally well fitted by assuming slow addition of radicals to the RAFT agent and slow fragmentation of the so‐obtained intermediate radical as well as high cross‐termination rate. By identification of all side products, our EPR‐based method is shown to be model independent and to provide reliable K_eq values, which demonstrate the validity of the intermediate radical termination model.     

The tropolone–isobutylamine complex: a hydrogen‐bonded troponoid without dominant π–π interactions

Tropolone long has served as a model system for unraveling the ubiquitous phenomena of proton transfer and hydrogen bonding. This molecule, which juxtaposes ketonic, hydroxylic, and aromatic functionalities in a framework of minimal complexity, also has provided a versatile platform for investigating the synergism among competing intermolecular forces, including those generated by hydrogen bonding and aryl coupling. Small members of the troponoid family typically produce crystals that are stabilized strongly by pervasive π–π, C—H…π, or ion–π interactions. The organic salt (TrOH·iBA) formed by a facile proton‐transfer reaction between tropolone (TrOH) and isobutylamine (iBA), namely isobutylammonium 7‐oxocyclohepta‐1,3,5‐trien‐1‐olate, C₄H₁₂N⁺·C₇H₅O₂⁻, has been investigated by X‐ray crystallography, with complementary quantum‐chemical and statistical‐database analyses serving to elucidate the nature of attendant intermolecular interactions and their synergistic effects upon lattice‐packing phenomena. The crystal structure deduced from low‐temperature diffraction measurements displays extensive hydrogen‐bonding networks, yet shows little evidence of the aryl forces (viz. π–π, C—H…π, and ion–π interactions) that typically dominate this class of compounds. Density functional calculations performed with and without the imposition of periodic boundary conditions (the latter entailing isolated subunits) documented the specificity and directionality of noncovalent interactions occurring between the proton‐donating and proton‐accepting sites of TrOH and iBA, as well as the absence of aromatic coupling mediated by the seven‐membered ring of TrOH. A statistical comparison of the structural parameters extracted for key hydrogen‐bond linkages to those reported for 44 previously known crystals that support similar binding motifs revealed TrOH·iBA to possess the shortest donor–acceptor distances of any troponoid‐based complex, combined with unambiguous signatures of enhanced proton‐delocalization processes that putatively stabilize the corresponding crystalline lattice and facilitate its surprisingly rapid formation under ambient conditions.     

Transition‐Metal π‐Ligation of a Tetrahalodiborane

The reaction of tetraiododiborane (B₂I₄) with trans‐[Pt(BI₂)I(PCy₃)₂] gives rise to the diplatinum(II) complex [{(Cy₃P)(I₂B)Pt}₂(μ₂:η³:η³‐B₂I₄)], which is supported by a bridging diboranyl dianion ligand [B₂I₄]^2?. This complex is the first transition‐metal complex of a diboranyl dianion, as well as the first example of intact coordination of a B₂X₄ (X=halide) unit of any type to a metal center.     

The peptide NCbz‐Val‐Tyr‐OMe and aromatic π–π interactions

The peptide N‐benzyloxycarbonyl‐L‐valyl‐L‐tyrosine methyl ester or NCbz‐Val‐Tyr‐OMe (where NCbz is N‐benzyloxycarbonyl and OMe indicates the methyl ester), C₂₃H₂₈N₂O₆, has an extended backbone conformation. The aromatic rings of the Tyr residue and the NCbz group are involved in various attractive intra‐ and intermolecular aromatic π–π interactions which stabilize the conformation and packing in the crystal structure, in addition to N—H...O and O—H...O hydrogen bonds. The aromatic π–π interactions include parallel‐displaced, perpendicular T‐shaped, perpendicular L‐shaped and inclined orientations.     

Interplay between Metal⋅⋅⋅π Interactions and Hydrogen Bonds: Some Unusual Synergetic Effects of Coinage Metals and Substituents

The ternary systems of C₂H₄ (C₂H₂ or C₆H₆)‐MCN‐HF (M=Cu, Ag, Au) and the respective binary systems were investigated to study the interplay between metal???π interactions and hydrogen bonds. The metal???π interactions in C₂H₄‐MCN become stronger with the irregular order Ag<Cu<Au, while the hydrogen bonds in MCN‐HF become weaker following the same order. The metal???π interactions are weakened as the H atoms in the π system are replaced with electron‐withdrawing groups and enhanced by electron‐donating groups. Type 1 of these ternary systems, in which MCN acts as Lewis base and acid simultaneously, is more stable than type 2, in which C₂H₄ acts as a double Lewis base. Negative cooperativity is present in type 2 ternary systems with a weakening of the metal???π interactions and the hydrogen bonds. Positive cooperativity is found in type 1 ternary systems with an enhancement of the metal???π interactions and the hydrogen bonds, except for C₂(CN)₄‐AuCN‐HF‐1. The weaker metal???π interaction in C₆H₆‐AuCN has a greater enhancing effect on the hydrogen bond in AuCN‐HF than those in C₂H₄‐AuCN and C₂H₂‐AuCN. These synergetic effects were analyzed with the natural bond orbital and energy decomposition.     

A Stable Cyclic (R2SnAu)3 Anion Having In‐Plane σ‐Möbius Aromaticity

A cyclic (R₂SnAu)₃ anion ( 3^? , R₂Sn=2,2,5,5‐tetrakis(trimethylsilyl)‐1‐stannacyclopentane‐1,1‐diyl) has been synthesized as a stable blue salt with K⁺(THF)₆ through the reaction of stable dialkylstannylene 1 with R′₃PAuCl (R′=Et, Ph) followed by the reduction with KC₈. Crystallographic and NMR analysis shows that the six‐membered (SnAu)₃ ring of 3^? is planar and highly symmetric with an equal distance of six Au?Sn bonds. A UV/Vis spectrum of 3^? in hexane reveals an intense absorption maximum at 598 nm. While cyclic Au₃^? with four valence electrons is known as unstable anti‐aromatic anion, 3^? with three divalent tin ligands is stable σ aromatic anion with an unprecedented Möbius orbital array as predicted by the perturbation MO and CCSD analysis of 3^? .     

Synthesis and Crystal Structure of a [Mo8O26]4– Cluster Derivative with 4‐MePyH+. First β‐Octamolybdate Derivative with π–π Stacking

Dioxobis(pyridine‐2‐thiolate‐N, S)molybdenum(VI) (MoO₂(Py‐S)₂), reacts with of 4‐methylpyridine (4‐MePy) in acetonitrile, by slow diffusion, to afford the title compound. This has been characterized by elemental analysis, IR and ¹H NMR spectroscopy. The X‐ray single crystal structure of the complex is described. Structural studies reveal that the molecular structure consists of a β‐Mo₈O₂₆ polyanion with eight MoO₆ distorted edge‐shared octahedra with short terminal Mo–O bonds (1.692–1.714 Å), bonds of intermediate length (1.887–1.999 Å) and long bonds (2.150–2.473 Å). Two different types of hydrogen bonds have been found: N–H···O (2.800–3.075 Å) and C–H···O (3.095–3.316 Å). The presence of π–π stacking interactions and strong hydrogen bonds are presumably responsible for the special disposition of the pyridinic rings around the polyanion cluster.