From Organic Fragments to Photoswitchable Catalysts: The OFF-ON Structural Repository for Transferable Kernel-Based Potentials.

Structurally and conformationally diverse databases are needed to train accurate neural networks or kernel-based potentials capable of exploring the complex free energy landscape of flexible functional organic molecules. Curating such databases for species beyond "simple" drug-like compounds or molecules composed of well-defined building blocks (e.g., peptides) is challenging as it requires thorough chemical space mapping and evaluation of both chemical and conformational diversities. Here, we introduce the OFF-ON (organic fragments from organocatalysts that are non-modular) database, a repository of 7869 equilibrium and 67,457 nonequilibrium geometries of organic compounds and dimers aimed at describing conformationally flexible functional organic molecules, with an emphasis on photoswitchable organocatalysts. The relevance of this database is then demonstrated by training a local kernel regression model on a low-cost semiempirical baseline and comparing it with a PBE0-D3 reference for several known catalysts, notably the free energy surfaces of exemplary photoswitchable organocatalysts. Our results demonstrate that the OFF-ON data set offers reliable predictions for simulating the conformational behavior of virtually any (photoswitchable) organocatalyst or organic compound composed of H, C, N, O, F, and S atoms, thereby opening a computationally feasible route to explore complex free energy surfaces in order to rationalize and predict catalytic behavior.

Tuning the Thermal Stability and Photoisomerization of Azoheteroarenes through Macrocycle Strain*.

Azobenzene and its derivatives are one of the most widespread molecular scaffolds used in a range of modern applications, as well as in fundamental research. After photoexcitation, azo-based photoswitches revert back to the most stable isomer on a timescale ( t 1 / 2 ) that determines the range of potential applications. Attempts to bring t 1 / 2 to extreme values prompted the development of azobenzene and azoheteroarene derivatives that either rebalance the E- and Z-isomer stabilities, or exploit unconventional thermal isomerization mechanisms. In the former case, one successful strategy has been the creation of macrocycle strain, which tends to impact the E/Z stability asymmetrically, and thus significantly modify t 1 / 2 . On the bright side, bridged derivatives have shown an improved optical switching owing to the higher quantum yields and absence of degradation. However, in most (if not all) cases, bridged derivatives display a reversed thermal stability (more stable Z-isomer), and smaller t 1 / 2 than the acyclic counterparts, which restricts their potential interest to applications requiring a fast forward and backwards switch. In this paper, the impact of alkyl bridges on the thermal stability of phenyl-azoheteroarenes is investigated by using computational methods, and it is revealed that it is indeed possible to combine such improved photoswitching characteristics while preserving the regular thermal stability (more stable E-isomer), and increased t 1 / 2 values under the appropriate connectivity and bridge length.

The Photoisomerization Pathway(s) of Push-Pull Phenylazoheteroarenes*.

Azoheteroarenes are the most recent derivatives targeted to further improve the properties of azo-based photoswitches. Their light-induced mechanism for trans-cis isomerization is assumed to be very similar to that of the parent azobenzene. As such, they inherited the controversy about the dominant isomerization pathway (rotation vs. inversion) depending on the excited state (nπ* vs. ππ*). Although the controversy seems settled in azobenzene, the extent to which the same conclusions apply to the more structurally diverse family of azoheteroarenes is unclear. Here, by means of non-adiabatic molecular dynamics, the photoisomerization mechanism of three prototypical phenyl-azoheteroarenes with increasing push-pull character is unraveled. The evolution of the rotational and inversion conical intersection energies, the preferred pathway, and the associated kinetics upon both nπ* and ππ* excitations can be linked directly with the push-pull substitution effects. Overall, the working conditions of this family of azo-dyes is clarified and a possibility to exploit push-pull substituents to tune their photoisomerization mechanism is identified, with potential impact on their quantum yield.

Thermal spin crossover in Fe(ii) and Fe(iii). Accurate spin state energetics at the solid state.

The thermal spin crossover (SCO) phenomenon refers to an entropy-driven spin transition in some materials based on d6-d9 transition metal complexes. While its molecular origin is well known, intricate SCO behaviours are increasingly common, in which the spin transition occurs concomitantly to e.g. phase transformations, solvent absorption/desorption, or order-disorder processes. The computational modelling of such cases is challenging, as it requires accurate spin state energies in the solid state. Density Functional Theory (DFT) is the best framework, but most DFT functionals are unable to balance the spin state energies. While a few hybrid functionals perform better, they are still too expensive for solid-state minima searches in moderate-size systems. The best alternative is to dress cheap local (LDA) or semi-local (GGA) DFT functionals with a Hubbard-type correction (DFT+U). However, the parametrization of U is not straightforward due to the lack of reference values, and because ab initio parametrization methods perform poorly. Moreover, SCO complexes undergo notable structural changes upon transition, so intra- and inter-molecular interactions might play an important role in stabilizing either spin state. As a consequence, the U parameter depends strongly on the dispersion correction scheme that is used. In this paper, we parametrize U for nine reported SCO compounds (five based on FeII, 1-5 and four based on FeIII, 6-9) when using the D3 and D3-BJ dispersion corrections. We analyze the impact of the dispersion correction treatments on the SCO energetics, structure, and the unit cell dimensions. The average U values are different for each type of metal ion (FeIIvs. FeIII), and dispersion correction scheme (D3 vs. D3-BJ) but they all show excellent transferability, with mean absolute errors (MAE) below chemical accuracy (i.e. MAE <4 kJ mol-1). This enables a better description of SCO processes and, more generally, of spin state energetics, in materials containing FeII and FeIII ions.

Exploring chemical space in the search for improved azoheteroarene-based photoswitches.

In quest of improved photo switches, azoheteroarenes have emerged as a potential alternative to azobenzenes. However, to date the number and types of these species that have been subjected to study are insufficient to provide an in-depth understanding of the photochemical effects brought about by different substituents. Here, we computationally screen the optical properties and thermal stabilities of 512 azoheteroarenes that consist of eight different N-containing heteroarenes combined with 64 substitution patterns. The most promising compounds are identified and their properties rationalized based on the nature of the azoheteroarene core and the location and type of substitution patterns.

Revising the common understanding of metamagnetism in the molecule-based bisdithiazolyl BDTMe compound.

The BDTMe molecule-based material is the first example of a thiazyl radical to exhibit metamagnetic behavior. Contrary to the common idea that metamagnetism occurs in low-dimensional systems, it is found that BDTMe magnetic topology consists of a complex 3D network of almost isotropic ferromagnetic spin-ladders that are coupled ferromagnetically and further connected by some weaker antiferromagnetic interactions. Calculated magnetic susceptibility χT(T) data is in agreement with experiment. Calculated M(H) data clearly show the typical sigmoidal shape of a metamagnet at temperatures below 2 K. The calculated critical field becomes more apparent in the dM/dH(H) plot, being in very good agreement with experiment. Our computational study concludes that the magnetic topology of BDTMe is preserved throughout the entire experimental range of temperatures (0-100 K). Accordingly, the ground state is the same irrespective of the temperature at which we study the BDTMe crystal. Revising the commonly accepted understanding of a metamagnet explained as ground state changing from antiferromagnetic to ferromagnetic, the Boltzmann population of the different states is here suggested to be the key concept: at 2 K the ground singlet state has more weight (24%) than at 10 K (1.5%), where excited states have an important role. Changes in the antiferromagnetic interactions that couple the ferromagnetic skeleton of BDTMe will directly affect the population of the distinct states that belong to a given magnetic topology and thus its magnetic response. Accordingly, this strategy could be valid for a wide range of bisdithiazolyl BDT-compounds whose magnetism can be tuned by means of weak antiferromagnetic interactions.

Spin State Chemistry: Modulation of Ligand p Ka by Spin State Switching in a [2×2] Iron(II) Grid-Type Complex.

The iron(II) [2×2] grid complex Fe-8H has been synthesized and characterized. It undergoes spin-crossover (SCO) upon deprotonation of the hydrazine-based terpyridine-like ligand. The deprotonation patterns have been determined by X-ray crystallography and 1H NMR spectroscopy and discussed in relation to the spin state of the iron(II) centers, which influences greatly the p Ka of the ligand. The synthesis of the magnetically silent zinc(II) analogue is also reported, and its (de)protonation behavior has been characterized to serve as a reference for the study of the FeII grid complexes. DFT computations have also been performed in order to investigate how the successive deprotonation of the bridging ligands affects the SCO behavior within the grid.

Absorption Spectroscopy and Photophysics of a ReI -dppz Probe for DNA-Mediated Charge Transport.

Optical properties of [Re(CO)3 (dppz)(py)]+ (dppz=dipyrido[3,2-a:2',3'-c]phenazine; py=pyridine) in acetonitrile, water and DNA have been investigated based on DFT, time-dependent-DFT (TD-DFT)/ conductor-like screening model, with and without explicit solvent molecules, and molecular dynamics. Whereas implicit solvent model is not appropriate to model optical properties of dppz-substituted metal complexes, adding explicit solvent molecules in interaction with dppz stabilizes the metal-to-ligand-charge-transfer (MLCT) transitions. Classical molecular dynamics simulations point to an important conformational flexibility, as evidenced by the coexistence of two conformers A and B. When considering the conformational sampling, the lowest band of the absorption spectrum is red-shifted and broadened up to 500 nm in agreement with the experimental spectra supporting important dynamical effects. The absorption spectra of [Re(CO)3 (dppz)(py-R)]+/ GC-DNA and [Re(CO)3 (dppz)(py-R)]+ /AT-DNA (R=CH2 -CH2 -COO- ) intercalated in both major or minor grooves exhibit a lowest energy charge separated (CS) band at about 600 nm and 500 nm, respectively, corresponding mainly to excitations from guanine and adenine to dppz. These states may play a central role into DNA-mediated charge transport processes. The over stabilization of the lowest 3 ILdppz state of [Re(CO)3 (dppz)(py)]+ in water as compared to acetonitrile could be responsible for the quenching of emission in water.

Cooperativity in Spin Crossover Systems. An Atomistic Perspective on the Devil's Staircase.

Cooperativity is key in defining the shape (i.e., gradual, abrupt, or hysteretic) of thermally driven spin transitions in magnetic switches. Despite its importance, there is very little information on its atomistic origin, which hinders the rational design of materials displaying a bistability region (i.e., hysteresis). In this paper, we investigate the spin transition of two solvatomorphs of [Fe(2-pic)3]Cl2, an Fe(II)-complex displaying thermal spin crossover (SCO) from a low-spin (LS) to a high-spin (HS) state with either gradual or abrupt two-step character. To do it, we apply a novel computational protocol to study the cooperativity of SCO compounds from DFT calculations, which does not rely on a priori assumptions on the studied system. The distinct shape of the spin transition is successfully captured, and the atomistic origin of cooperativity is traced back to geometrical distortions of the Fe-N6 core in case of the solvatomorph exhibiting an abrupt transition. According to our calculations, HS and LS molecules contribute differently to cooperativity, which results in a complex energetic evolution of the spin transition that cannot be described by the common Slichter-Drickamer model. The present work opens new avenues for the study of cooperativity of SCO systems having a chemically oriented perspective and demonstrates that quantum chemistry calculations can discriminate the shape of a spin transition, while providing insight into the atomistic underlying factors that contribute to its cooperative behavior.

A Spin-Crossover Molecular Material Describing Four Distinct Thermal Pathways.

Spin-crossover (SCO) molecular solids are valued switchable materials for their common abrupt and reversible thermal transitions, large thermal hysteresis, or guest-dependent effects. These properties usually involve crystallographic transitions coupled to the SCO events. These phenomena are of great value for the understanding of solid-state transformations and also for exploiting them. We present here a lattice of the complex [FeL(bbp)](ClO4)2 (1; L and bbp are tris-imine ligands) featuring an unprecedented rich succession of SCO and crystallographic phase transformations. Magnetometry measurements unveil a thermally irreversible sequence of spin conversions that delineate four different thermal pathways. All of these are single-crystal-to-single-crystal processes and can thus be monitored by single crystal X-ray diffraction using one unique specimen. Fresh crystals of 1 contain one molecule of acetone per Fe center (1·ac) that abandons the lattice upon warming at the same time that a SCO from an ordered mixed spin state (1:1 high spin/low spin; HS/LS) to a fully HS state, 1α, occurs. This crystallographic phase, accessed through a template effect by the solvent, converts into another one, 1ß, upon cooling, as triggered by a HS to LS SCO. Warming of 1ß induces a new SCO (LS to ordered HS/LS) coupled to another crystallographic phase transition, 1ß â†’ 1γ. The fully HS state of 1γ can not be reached before decomposition of the compound. Instead, this phase cycles between the HS/LS and the LS states through superimposable pathways, different from that of the prerequired 1ß â†’ 1γ phase change. Analysis of the thermal variation of the free energy, G, through density functional theory methods provides trends in agreement with the observation of these transformations and clarifies the possible metastable nature of the various phases identified. This unique behavior allows the access to four different magnetic responses depending on the thermal history of the sample, within a given range of temperatures near the ambient conditions.

Pairing-up viologen cations and dications: a microscopic investigation of van der Waals interactions.

The microscopic origin of van der Waals- and magnetic-interactions in 4,4' methyl viologen cation-based units (MV+˙ and MV2+) was inspected using wave function (variational DDCI and perturbative MP2, CASPT2) and density functional theory (DFT) calculations. The analysis deepens the comprehension of the magnetic behavior of experimental bis-viologen cyclophanes ([CYC]2(+˙)), in which the MV+˙ units are connected through alkyl linkers of different lengths. The formation of the so-called long-multicenter bonds in such radical dimers, responsible for the quenching of the magnetic response, was analyzed in [MV2]2(+˙). Dynamical correlation effects, accessible from second-order perturbation corrections, were decisive in observing a bonding regime characterized by an equilibrium distance of 3.3 Šand a 45 kJ mol-1 dissociation energy. At larger intermolecular distances, our calculations on [MV2]2(+˙) indicate that the singlet and triplet states are energetically competing (i.e. weak exchange interactions, JAB). Despite the absence of any clear bonding regime at the MP2 level, the puzzling association of two di-cations into [MV2]4+ is anticipated at 3.3 Šusing weakly screened point charges (ε = 1.5) to account for the Coulomb interactions between the solvated subunits. The main conclusion is that both dispersion interactions and environment effects are required to overcome the Coulomb repulsion associated with doubly-charged species. All these data provide some complementary insights into the nature and amplitude of interactions between cation and dication units, and their relevance in various experimental manifestations.

Luminescent Dinuclear Copper(I) Complexes as Potential Thermally Activated Delayed Fluorescence (TADF) Emitters: A Theoretical Study.

The excited state properties of a series of binuclear NHetPHOS-Cu(I) complexes (NHetPHOS) have been investigated by means of density functional theory (DFT) and time-dependent DFT (TD-DFT). It is shown that experimental trends observed in powder, generally explored via S1 and T1 excited state energetics and S1 ⇔ T1 intersystem crossing (ISC) efficiency, are hardly analyzed on the basis of excited state properties calculated in solution. Indeed, several local minima corresponding to various structural deformations are evident on the lowest excited state potential energy surfaces (PES) when solvent correction is applied, leading to a four-state thermally activated delayed fluorescence (TADF) mechanism. In contrast, preliminary simulations performed in the solid point to the reduction of nuclear flexibility and consequently to a rather simple two-state model.

Bistability in Organic Magnetic Materials: A Comparative Study of the Key Differences between Hysteretic and Non-hysteretic Spin Transitions in Dithiazolyl Radicals.

Dithiazolyl (DTA)-based radicals have furnished many examples of organic spin-transition materials, some of them occurring with hysteresis and some others without. Herein, we present a combined computational and experimental study aimed at deciphering the factors controlling the existence or absence of hysteresis by comparing the phase transitions of 4-cyanobenzo-1,3,2-dithiazolyl and 1,3,5-trithia-2,4,6-triazapentalenyl radicals, which are prototypical examples of non-bistable and bistable spin transitions, respectively. Both materials present low-temperature diamagnetic and high-temperature paramagnetic structures, characterized by dimerized (⋅⋅⋅A-A⋅⋅⋅A-A⋅⋅⋅)n and regular (⋅⋅⋅A⋅⋅⋅A⋅⋅⋅A⋅⋅⋅A⋅⋅⋅)n π-stacks of radicals, respectively. We show that the regular π-stacks are not potential energy minima but average structures arising from a dynamic inter-conversion between two degenerate dimerized configurations: (⋅⋅⋅A-A⋅⋅⋅A-A⋅⋅⋅)n ↔(-A⋅⋅⋅A-A⋅⋅⋅A-)n . The emergence of this intra-stack dynamics upon heating gives rise to a second-order phase transition that is responsible for the change in the dominant magnetic interactions of the system. This suggests that the promotion of a (⋅⋅⋅A-A⋅⋅⋅A-A⋅⋅⋅)n ↔(-A⋅⋅⋅A-A⋅⋅⋅A-)n dynamics is a general mechanism for triggering spin transitions in DTA-based materials. Yet, this intra-stack dynamics does not suffice to generate bistability, which also requires a rearrangement of the intermolecular bonds between the π-stacks via a first-order phase transition.

Lattice-Solvent Effects in the Spin-Crossover of an Fe(II)-Based Material. The Key Role of Intermolecular Interactions between Solvent Molecules.

The spin transition of Fe(II) complexes is the subject of intensive synthetic and computational efforts. In this manuscript, we analyze the spin crossover (SCO) of [Fe(E-dpsp)2]2+ (1), which features a spin transition depending on the cocrystallizing solvent molecules. Whereas the use of acetone results in a hysteretic spin transition at ∼170 K, the use of propylene carbonate (PC) results in a permanent diamagnetic signal up to 300 K. By means of DFT+U+D2 calculations in the solid state of the material, we unravel the reasons for such different behavior. Our results allow us to ascribe the relatively low transition temperature of 1(BF4)2·acetone to the distorted arrangement of the SCO molecules in the low-spin state of the material. In turn, intermolecular interactions play the primary role in the case of 1(BF4)2·2PC. In particular, we found that solvent-solvent interactions actively promote the stability of the low-spin state due to the formation of PC dimers. These dimers would appear at larger distances in the high-spin phase, with the subsequent loss of phase stability. This is yet another proof of how subtle is the spin transition phenomenon in Fe(II)-based architectures.

Electron transport through a spin crossover junction. Perspectives from a wavefunction-based approach.

The present paper reports the application of a computational framework, based on the quantum master equation, the Fermi's golden Rule, and conventional wavefunction-based methods, to describe electron transport through a spin crossover molecular junction (Fe(bapbpy) (NCS)2, 1, bapbpy = N-(6-(6-(Pyridin-2-ylamino)pyridin-2-yl)pyridin-2-yl)-pyridin-2-amine). This scheme is an alternative to the standard approaches based on the relative position and nature of the frontier orbitals, as it evaluates the junction's Green's function by means of accurate state energies and wavefunctions. In the present work, those elements are calculated for the relevant states of the high- and low-spin species of 1, and they are used to evaluate the output conductance within a given range of bias- and gate-voltages. The contribution of the ground and low-lying excited states to the current is analyzed, and inspected in terms of their 2S + 1 Ms-states. In doing so, it is shown the relevance of treating not only the ground state in its maximum-Ms projection, as usually done in most computational-chemistry packages, but the whole spectrum of low-energy states of the molecule. Such improved representation of the junction has a notable impact on the total conductivity and, more importantly, it restores the equivalence between alpha and beta transport, which means that no spin polarization is observed in the absence of Zeeman splitting. Finally, this work inspects the strong- and weak-points of the suggested theoretical framework to understand electron transport through molecular switchable materials, identifies a pathway for future improvement, and offers a new insight into concepts that play a key role in spintronics.

On the zeroth-order hamiltonian for CASPT2 calculations of spin crossover compounds.

Complete active space self-consistent field theory (CASSCF) calculations and subsequent second-order perturbation theory treatment (CASPT2) are discussed in the evaluation of the spin-states energy difference (ΔH(elec)) of a series of seven spin crossover (SCO) compounds. The reference values have been extracted from a combination of experimental measurements and DFT + U calculations, as discussed in a recent article (Vela et al., Phys Chem Chem Phys 2015, 17, 16306). It is definitely proven that the critical IPEA parameter used in CASPT2 calculations of ΔH(elec), a key parameter in the design of SCO compounds, should be modified with respect to its default value of 0.25 a.u. and increased up to 0.50 a.u. The satisfactory agreement observed previously in the literature might result from an error cancellation originated in the default IPEA, which overestimates the stability of the HS state, and the erroneous atomic orbital basis set contraction of carbon atoms, which stabilizes the LS states.

Understanding the Influence of the Electronic Structure on the Crystal Structure of a TTF-PTM Radical Dyad.

The understanding of the crystal structure of organic compounds, and its relationship to their physical properties, have become essential to design new advanced molecular materials. In this context, we present a computational study devoted to rationalize the different crystal packing displayed by two closely related organic systems based on the TTF-PTM dyad (TTF = tetrathiafulvalene, PTM = polychlorotriphenylmethane) with almost the same molecular structure but a different electronic one. The radical species (1), with an enhanced electronic donor-acceptor character, exhibits a herringbone packing, whereas the nonradical protonated analogue (2) is organized forming dimers. The stability of the possible polymorphs is analyzed in terms of the cohesion energy of the unit cell, intermolecular interactions between pairs, and molecular flexibility of the dyad molecules. It is observed that the higher electron delocalization in radical compound 1 has a direct influence on the geometry of the molecule, which seems to dictate its preferential crystal structure.

Towards an accurate and computationally-efficient modelling of Fe(II)-based spin crossover materials.

The DFT + U methodology is regarded as one of the most-promising strategies to treat the solid state of molecular materials, as it may provide good energetic accuracy at a moderate computational cost. However, a careful parametrization of the U-term is mandatory since the results may be dramatically affected by the selected value. Herein, we benchmarked the Hubbard-like U-term for seven Fe(ii)N6-based pseudo-octahedral spin crossover (SCO) compounds, using as a reference an estimation of the electronic enthalpy difference (ΔHelec) extracted from experimental data (T1/2, ΔS and ΔH). The parametrized U-value obtained for each of those seven compounds ranges from 2.37 eV to 2.97 eV, with an average value of U = 2.65 eV. Interestingly, we have found that this average value can be taken as a good starting point since it leads to an unprecedented mean absolute error (MAE) of only 4.3 kJ mol(-1) in the evaluation of ΔHelec for the studied compounds. Moreover, by comparing our results on the solid state and the gas phase of the materials, we quantify the influence of the intermolecular interactions on the relative stability of the HS and LS states, with an average effect of ca. 5 kJ mol(-1), whose sign cannot be generalized. Overall, the findings reported in this manuscript pave the way for future studies devoted to understand the crystalline phase of SCO compounds, or the adsorption of individual molecules on organic or metallic surfaces, in which the rational incorporation of the U-term within DFT + U yields the required energetic accuracy that is dramatically missing when using bare-DFT functionals.

Elucidating the 2D magnetic topology of the 'metal-radical' TTTA⋠Cu(hfac)2 system.

The TTTA⋅Cu(hfac)2 polymer (1; in which TTTA = 1,3,5-trithia-2,4,6-triazapentalenyl, and hfac = (1,1,1,5,5,5)-hexafluoroacetylacetonate) is one of the most prominent examples of the rational use of the 'metal-radical' synthetic approach to achieve ferromagnetic interactions. Experimentally, the magnetic topology of 1 could not be fully deciphered. Herein, the first-principles bottom-up procedure was applied to elucidate the nature and strength of the magnetic JAB exchange interactions present in 1. The computed JAB values give rise to a 2D magnetic topology of ferromagnetic dimers (+11.9 cm(-1)) coupled through weaker antiferromagnetic interactions (-3.0 and -3.2 cm(-1)) in two different spatial directions. The hitherto unknown origin of the antiferromagnetic interdimer interactions is thus unveiled. By using the 2D magnetic topology, the agreement between calculated and experimental χT(T) data is extraordinary. In the metal-radical TTTA⋅Cu(hfac)2 compound, the computational model transcends the local dimer cluster model owing to strong interactions between metal centers and organic radicals, thereby creating a de facto biradical. In addition, it is shown that the magnetic topology cannot be inferred from the polymeric [TTTA⋅⋅⋅Cu(hfac)2]n crystal motif, that is, from its chemical coordination pattern. Instead, one should think in terms of magnetic building blocks, namely, the de facto biradicals.

Linear or cyclic clusters of Cu(II) with a hierarchical relationship.

The polydentate ligand 2,6-bis(5-(2-hydroxyphenyl)-pyrazol-3-yl)-pyridine, H4L, exhibits a series of coordination pockets favoring the establishment of metal sequences with predetermined motifs, together with a degree of flexibility for the formation of clusters with various overall topologies. With Cu(II) under strong basic conditions it has a marked tendency to stabilize a cyclic [Cu16L8] cluster. The sequential formation of this compound via [Cu7L8](2-) intermediates, recognized in its structure, is suggested by crystallographic evidence, which shows the persistent formation of the complex salt (NBu4)2[Cu7L8] in the presence of the organic cation. Also, the crystallographic identification of the related cluster [Cu11L5(OH)2(py)12] from similar reaction conditions underscores the rich multiplicity of species attainable from this simple reaction system.