Novel Trinuclear Copper(II) Complex: Crystal Structure at 100 K and Magnetic Properties of (R, S)-di[(6,7-dimethoxy-isoquinolin-1-yl)-(3,4-dimethoxy-phenyl)-methanolato]-tetra(2-hydroxybenzoato)-diaqua-tricopper dihydrate, [Cu3(C7H5O3)4(C20H20NO5)2(H2O)2]<middle dot>2(H2O)

Vol.:(0123456789) Journal of Chemical Crystallography (2024) 54:163–172 https://doi.org/10.1007/s10870-024-01009-2 ORIGINAL PAPER Novel Trinuclear Copper(II) Complex: Crystal Structure at100 K andMagnetic Properties of(R, S )‑ di[ (6, 7‑d ime tho xy‑ iso qui nol in‑ 1‑y l)‑ (3, 4‑d ime tho xy‑ phe nyl )‑m eth ano lat o]‑ tet ra( 2‑h ydr oxybenzoato)‑diaqua‑tricopper dihydrate, [Cu3(C7H5O3)4(C20H20NO5)2(H2O)2]·2(H2O) FedorValach1 · JánPavlik2 · IvanŠalitroš2,3 · MilanMelník2 · JozefKožíšek1 Received: 3 August 2023 / Accepted: 8 February 2024 / Published online: 28 March 2024 © The Author(s) 2024 Abstract The crystal structure of [Cu3(C7H5O3)4(C20H20NO5)2(H2O)2]·2(H2O) (1) and analysis of temperature and field dependence of magnetic susceptibility is reported in this work. The structure of 1 is composed of trinuclear complex units and water molecules. The middle copper atom occupies the center of symmetry. N, O-bonded (6,7-dimethoxy-isoquinolin-1yl)-(3,4-di methoxy-phenyl)-methanolato ligands, 2-hydroxybenzoates with bridging carboxylic groups, and oxo-bridged water molecules connect the middle Cu(II) atom with the terminal copper atoms. Two 2-hydroxybenzoates coordinate the terminal copper atoms via one carboxylic oxygen and an O atom of the hydroxyl group. The analysis of copper coordination by bond-valence sum approach and relevant structural correlation is consistent with hexacoordinated Cu(II) centers. Cu···Cu separation is 3.0269(3) Å. The magnetism of 1 shows a strong ferromagnetic interaction between the neighboring metallic centers accompanied by very weak antiferromagnetic intermolecular interactions. The complex units are mutually held by π···π stack interactions of 2-hydroxybenzoates and hydrogen bonds. * Fedor Valach v[email protected] 1 Faculty ofChemical, andFood Technology, Institute ofPhysical Chemistry andChemical Physics, Slovak University ofTechnology inBratislava, Radlinského 9, SK-81237Bratislava, SlovakRepublic 2 Present Address: Faculty ofChemical, andFood Technology, Institute ofInorganic Chemistry, Technology, andInorganic Materials, Slovak University ofTechnology inBratislava, Radlinského 9, SK-81237Bratislava, SlovakRepublic 3 Central European Institute ofTechnology, University ofTechnology inBrno, Purkyňova 123, 61200Brno, CzechRepublic 164 Journal of Chemical Crystallography (2024) 54:163–172 Graphical Abstract A new N,O bonded ligands, (R, S)-[(6,7-dimethoxy-isoquinolin-1-yl)-(3,4-dimethoxy-phenyl)-methanolate] coordinate the terminal atoms of the trinuclear copper(II) complex. Keywords Trinuclear copper(II) complex· Crystal structure· Carboxylate· Magnetic properties· Structure correlation· Inner coordination Introduction Structurally Closely Related Alkaloids toPapaverine (1-[(3,4-dimethoxyphenyl) methyl]-6,7dimethoxyisoquinoline) are important pharmacological drugs and biologically active compounds [1–3]. The molecule of papaverine (Scheme1) is unstable toward oxygen and can be oxidized easily to papaverinol and papaveraldine [4–7]. However, the crystal structure of papaverine hydrochloride [8] is composed of relatively stable nonoxidized papaverine molecules. The formation of copper complexes with the potential ligands contained in human organisms can cause many changes in cells [9]. To our knowledge, until now, with papaverine, only the crystal structure of the binuclear complex [Zn2(C6H5COO)4(pap)2] (pap = papaverine) was resolved [10]. Carboxylic copper(II) complexes may be mononuclear, binuclear, tetranuclear, or polynuclear [11–14]. Chemical coordination [15] defined by Werner [16] is an array of atoms surrounding the central atom because of chemical bonding. In transition element complexes, the atoms often interact with the central atom very weakly. In the quantum chemical definition, atomic valence is a measure of the degree of electron sharing of the given atom with the other atoms [17]. According to the residual valency of an atom, the atoms in a molecule may be subvalent, normal, or hypervalent [18]. Such classification deviates from the traditional valence of an atom in the crystal structure. It suggests using the atomic valence approximated by the bond valence sum as a parameter that is as close as possible to the integral value or “formal oxidation number“ [19, 20]. Scheme1 Molecule of papaverine 165Journal of Chemical Crystallography (2024) 54:163–172 The valence of a bond is the measure of interaction between electrons sharing the bonding of atoms. For the amorphous materials, the coordination spheres were limited [21, 22]. Beyond the distance, to the central atom around 6 Å, the refined parameters of dependence bond valance vs. interatomic distance do not change [23, 24]. Introducing chemical considerations, the concept of coordination requires estimating if an atom of a coordination array is bonded to the central atom or not. If such an atom significantly influences the interatomic distances of the remaining atoms of the array, it is primary bonded to the central atom. For the big family of copper complexes, the concept of coordination sphere plasticity [25] and equatorial-axial influence [26] reflect such influence via the central atom. Both effects are manifested in structural correlations of bond distances within the chromophores of transition element compounds. The structural correlations of coordination compounds with CuO6 and CuN6 chromophores approach the asymptotical values with the shortest equatorial bond length at around 1.9 Å [25]. Small changes in coordination polyhedra may cause significant changes in a hydrogen-bonding network [27]. In most cases, electronic distortions do not violate the valence sum rule [28, 29]. Importing dynamic features into molecular chemistry requires shifting from “static” to “dynamic” covalent bonds leading from molecular to constitutional dynamic chemistry [30]. As the magnetic properties of linearly arranged trinuclear hydroxy-bridged copper(II) complexes studied by the methods of density functional theory (DFT) showed, the angle Cu–O–Cu and out-of-plane displacement of the hydrogen atoms at the bridge are the key magnet factors that determine their magnetic behavior [31]. The magneto-structural study of binuclear copper(II) benzoates shows the dominance of antiferromagnetic spin interaction between metallic centers through bridging carboxylates over direct Cu···Cu interaction [32, 33]. These facts suggest the crucial structural influence of copper-ligand interactions within the coordination polyhedra of the central and terminal metallic centers. Here we report the crystal structure and magnetic properties of a novel compound consisting of trinuclear copper entities containing the alkaloid derived from papaverine. Experimental Section Preparation ofCrystals The light green product of the entitled compound was prepared by heating papaverine (0.01mol) with Cu(2HOC6H4COO)2·2H2O (0.01mol) in hot water and methanol. The mixture was stirred, filtered, and then left at room temperature. The solid product of the reaction was collected and recrystallized from hot methanol to give air-stable pale blue-turquoise crystals. The analytical data found/calculated are 12.42/12.54% (Cu), 53.40/53.74% (C), 4.40/4.6% (H), and 1.93/1.84% (N). Crystal Structure Solution andRefinement Single crystal data collection was performed on a STOE diffractometer at 100K. Crystal data and structure refinement details are summarized in Table1. The structure of1 was solved by the program SHELXS [34]. The structure was drawn by OLEX2 [35]. The refinements were carried out by SHELXL [36]. OMe group and water molecules disordered in two positions were refined under the constraint of their occupancy sums to unity. Positions of all hydrogen atoms were optimized under the constraint to ride on their parent atoms, with C–H bond length of 0.93 Å and O–H bond length of 0.82 Å. Magnetic Measurements Magnetic measurements were performed by a SQUID magnetometer MPMS-XL7 from Quantum Design in the RSO mode of detection. The temperature dependence of magnetic moment was recorded at 0.1 T as an external magnetic field and the temperature sweeping rate was set to 1K/min. The magnetic susceptibility of the polycrystalline sample was Table 1 The crystal data and details of crystal structure determination of 1 Chemical formula C68H64Cu3N2O24·2(H2O) Mr 1519.89 Crystal system, space group Triclinic, P¯1 a, b, c (Å) 10.3322(1), 13.4417(2), 13.6606(2) α, β, γ (°) 104.349(1), 111.855(1), 96.685(1) V (Å3) 1659.36(4) Z 1 Radiation type Cu Kα µ (mm−1) 1.857 Crystal size (mm) 0.48 × 0.16 × 0.15 Absorption correction Numerical Tmin, Tmax 0.469, 0.768 No. of measured, independent, and observed [I > 2σ(I)] reflections 64,315, 6353, 6328 Rint 0.022 (sin Θ/λ)max (Å−1) 0.654 R[F2 > 2σ(F2)], wR(F2), S 0.047, 0.131, 1.160 No. of reflections 6356 No. of parameters 510 No. of restraints 28 ∆ρmax, ∆ρmin (eÅ−3) 1.44, − 0.82 166 Journal of Chemical Crystallography (2024) 54:163–172 measured by the Faraday method in the temperature range of 77–300K. The mercury tetracyanato cobalt(II) was used as a calibrant [37]. Calculation Details Bond Valence Calculations ….For idealized copper hexacoordination of four equidistant equatorial in-plane Cu–N, O3 bonds and two equidistant outof-plane Cu–O′2 bonds (chromophore CuNO3O′2) the quantities rL, rS were defined. rS is the length of in-plane (equatorial) bonds and rL is the length of out-of-plane (axial) bonds. The bond valence sum (BVS) distribution was calculated as the sum of equatorial and axial contributions: where sCu–O and sCu–N are the bond valence functions for Cu–O and Cu–N bonds [38]: If BVS is kept constant, the contour of the function expressed by Eq.(1) with inserted Eq.(2) gives the dependence of rL vs. rS. For the idealized copper coordination of the chromophore Cun+NO3O′2 this dependence is where parameter n is the formal charge of a copper cation, Ro,Cu(n)–O is the bond valence parameter of a copper-oxygen bond, and sCu–N, sCu(n)–O are the functions [39]: The empirical bond-valence parameters Ro and b were used from the latest version of http:// www. iucr. org/ resou rces/ data/ dat as ets/ bondvalen ceparam eters. Equation(4) is appropriate for ionic, covalent, and intermediate ioniccovalent bonding, respectively [40]. The predicted axial Cu–O′(2×) bond length of the copper coordination of Cu2+O4O′2 chromophore was calculated as a numerical solution of the Eqs. [41, 42]: where BVSeq is the bond valance sum of equatorial Cu–O bonds calculated by Eq.(2). αi are the parameters of the function (2) for Cu–O bonds. (1) BVS( r S ,r L) =3s Cu−O( r S) +s Cu−N( r S) +2s Cu−O( r L) (2) s (r)= 𝛼 1 r + 𝛼 2 r 2+ 𝛼 3 r 3+ 𝛼 4 r 4+ 𝛼 5 r 5 (3) r L =R o,Cu(n)−O −bln [n−s Cu−N (r S ) −3 s Cu(n)−O( r S)] + b ln(2) (4) s (r)=exp [( R o −r ) ∕b ] (5) ( BVS eq −2 ) r 5 +2𝛼 1 r 4 +2𝛼 2 r 3 +2𝛼 3 r 2 +2𝛼 4 r+2𝛼 5 = 0 Magnetic Measurements The correction for the diamagnetism of the constituent atoms was calculated by Pascal’s constants [43]. The effective magnetic moments were calculated using the equation: To assess the values of magnetic exchange parameters the broken-symmetry DFT approach was employed. Major –O–CH3 (OMe) groups were used in calculations. The magnetic parameters were used in the form −JijŜiŜj. The final values of magnetic parameters were extracted using the adapted Eq.(7) of Yamaguchi and co-workers, for broken symmetry (BS) and high spin (HS) states [44, 46]: The fitting of the magnetic susceptibility and magnetization in 1 was performed with the program. PHI 3.1.1 [47]. Before the quantum-chemistry calculations, the positions of all hydrogen atoms were optimized using the method PBEh-3c [48]. All other atoms were kept in their positions as obtained from the crystal structure investigation. Calculation of magnetic exchange parameters was carried out within the program ORCA 4.2.0 [49, 50]. The second-order Douglass–Kroll–Hess correction (DKH2) [51] and the chainof-spheres approximation (RIJCOSX) [52] were set on. For all atoms, accounting for scalar relativistic effects Ahlrichs’s basis set DKH-def2-TZVP [53] was used with an auxiliary basis set def2/J [54]. The magnetic exchange parameters were obtained from DFT calculation utilizing the hybrid functional B3LYP [55–58]. In all calculations, the increased integration and fitting grids were used (grid6 and gridx6 in ORCA convention). Results andDiscussion Crystal Structure The crystal structure of compound 1 is composed of the tri copper units [Cu3(C7H5O3)4(C20H20NO5)2(H2O)2] (2) presented in Fig.1 and noncoordinated water molecules. Selected bond lengths and angles are collected in Table2. N, O -bonded(6,7-dimethoxyisoquinolin-1-yl)-(3,4dimethoxyphenyl)-methanolato ligand in this work is named papaverinyloxy ligand. The trinuclear unit shows chirality at C16. One OMe group attached to the phenyl ring is disordered in two positions of an occupancies ratio of 0.777(3)/0.223(3). The second OMe group is partially stabilized by the hydrogen bond. OW2–HW4‧‧‧O10, while its relative free methyl group is split in two positions of an occupancies ratio 0.500(3) /0.500(3). The noncoordinated water molecules of oxygens (6) 𝜇eff =2.83(𝜒 corr M ⋅T) 0.5 (7) J = ( E BS −E HS) ∕ ( <  S2> HS −<  S2> BS) 167Journal of Chemical Crystallography (2024) 54:163–172 OW2, and OW3 are disordered in two crystal sites with an occupancy ratio of 0.730(6) /0.270(6). The isoquinoline and phenyl rings of the papaverinyloxy ligand are significantly deviated from planarity. The angle between the least-squares planes of isoquinoline ring and the phenyl ring is 87.19(7)°. In papaverine hydrochloride, the corresponding angle is 80.13(8)°. In the structure of salicylic acid [60] the corresponding deviation is 0.0012(1) Å. The phenyl rings of 2-hydroxybenzoates do not deviate from planarity. But the atoms O6 of hydroxyl groups deviate from the least-squares planes of phenyl rings by 0.054(2) Å.The complex units 2 are connected via hydrogen bonds of OW2 water molecules connected into the infinite chains along [101] direction (Fig.2). The phenyl rings of adjacent symmetry-related 2-hydroxybenzoates assemble the complex units via π···π stack interactions [61, 62] into the chains along the [010] direction (Fig.3). The hydrogen bond lengths are presented in Table3. The inter-centroid distance and parallel shifts of phenyl rings are 3.7379(1) Å and 1.66901(5) Å. Slightly distorted trans-square coordination of Cu1 is created by two in-plane carboxylic O atoms and two oxygens of papaverinyloxy ligands. The five-membered ring containing Cu2 is slightly deviated from the planar geometry. The torsion angle N1–C1–C16–O7 is 18.3(3)°. The Cu2–O3 bond stabilizes the copper coordination of the 2-hydroxybenzoato ligand while the atom O2 remains free. The coordination polyhedron around Cu2 can be described as a distorted tetragonal bipyramid. The equatorial bonds of its inner coordination are Cu2–N1i, O1, O4, and O7i. The greatest deviation from the least-squares plane defined by the donor atoms is shown by atom N1i [0.079(2) Å]. Atom Cu2 is displaced from the basal plane toward the atom OW1 by 1.1406(1) Å. The hexacoordination of atom Cu2 is completed by the axial bonds Cu2–OW1 Fig. 1 The structure of complex 2 with the numbering of nonhydrogen atoms. The displacement ellipsoids are drawn at the 50% probability level. The equatorial mean planes of copper coordination polyhedra are shaded. The relative configuration at the C16 centers [59] is mentioned Table 2 Selected geometric parameters (Å, °) of the complex 2 Symmetry code: (i) − x + 1, − y + 1, − z + 2. Cu1-O5 1.940(2)Cu2–O4 1.945(2) Cu1–O7 1.959(2)Cu2–N1i1.970(2) Cu1–OW1 2.496(2)Cu2–OW1 2.442(2) Cu2–O1 1.925(2) Cu2–O3 2.873(2) Cu2–O7i1.930(2) O5–Cu1–O7 91.43(1) OW1–Cu1–O5 102.61(1) O1–Cu2–O4 91.10(1) N1–Cu2–O7 83.48(1) O1–Cu2–O3 64.43(1) N1–Cu2i–O3 92.65(1) Cu1–O7–Cu2i–O3i147.6(8) Cu1–O7–Cu2i–O1i− 74.8(1) Cu1–OW1i–Cu2i– O3i− 152.3(6) Cu1–OW1i–Cu2i– O1i160.0(6) Fig. 2 Schematical drawing of intermolecular hydrogen bonding in the crystal structure of 1 168 Journal of Chemical Crystallography (2024) 54:163–172 and Cu2–O3 of 2.445(2) Å and 2.872(2) Å, respectively. The angle formed by the Cu2–O3 bond and the basal leastsquares plane of the coordination polyhedron around Cu2 is 83.56(5)°. Bond Valence Analysis For copper atom valence 2+, the predicted lengths of Cu1–OW1 bonds calculated by Eq.(5) is 2.50(2) Å which does not deviate significantly from the experimental bond length of 2.497(2) Å. The coordination of Cu1 is completed by two out-of-plane water O atoms to an axially distorted tetragonal bipyramid of chromophore CuO4O′2. The bond valence sum of Cu1 calculated by formulae (2) is 2.05(2) v. u. Figure4 shows the contour graph of the copper atom valence (BVS) for idealized copper coordination of CuNO3O′2 chromophore. The numeric labels of isolines are the BVS values which are kept constant for relevant rL vs. rS dependence. The maximum bond lengths are rS(max) = 2.92(1) Å and rL(max) = 3.07(3) Å. For those distances, the contributions of equatorial and axial bonds to BVS in Eq.(1) yield zero values. In the regions of distances, the rS > rS(max) and rL > rL(max) these contributions are negative. The bond lengths rS(max) and rL(max) are the supreme limits of bonding distances within the inner copper coordination. If the bond valence sum of a copper atom is kept constant, any increase in the distance rS causes a decrease of rL and vice versa, i.e. such distortion of the inner coordination is plastic. 293 independent rS, rL values of the structures retrieved in CSD [64, 65] were approximated by the means of equatorial bond lengths (RS) and axial bond lengths (RL) [66]. The model structures were identified by the connectivity shown in Scheme2. The standard atomic radii and tolerances of CSD were used. Only the non-disordered structures with an R-factor less than 0.05 were accepted. The structures were inspected by the program ORTEP-3 for Windows [67]. REFCODES of model structures with RS and RL values are included in the Supplementary Information. The best approach to the contour line and the graph of Eq.(3) for the copper oxidation state 2.0 show points in the region rL ≥ rS. The weak bonds Cu2–O3 are included in the inner coordination of Cu2. The bond valences of these bonds calculated by Eq.(2) Fig. 3 Schematical drawing of interaromatic π···π interactions in the crystal structure of 1. Only the intramolecular hydrogen bonds and OMe groups of major occupancies are presented Table 3 Hydrogen-bond geometry (Å, º) of 1 Symmetry code: (i) − x + 1, − y + 1, − z + 2; (ii) − x, − y + 1, − z + 1 D–Ḣ ‧‧‧A D–H H‧‧‧A D‧‧‧A D–H‧‧‧A OW1–HW1‧‧‧O5i0.771(3) 2.810(4) 2.430(2) 111.8(2) OW1–HW2‧‧‧OW2 0.879(3) 1.827 (4) 2.698 (5) 170.8(2) OW2–HW3‧‧‧O2 0.849(3) 2.018(3) 2.807(4) 154.2(2) OW2–HW4‧‧‧O10ii 0.850(3) 2.077(2) 2.840(4) 159.8(2) Fig. 4 Diagram RLvs. RS (in Å) for the copper complexes of chromophore CuNO3O′2 compared with the calculated rLvs. rS dependences. Solid lines - contour graph of the function expressed by the Eqs.(1) and (2). Dashed lines – BVS functions expressed by the Eq.(3). ○ – model structures retrieved in CSD. ▲- the structure of 1. Spearman rank-order correlation coefficient [63] is − 0.844 Scheme2 The inner coordination of the model structures 169Journal of Chemical Crystallography (2024) 54:163–172 are 0.14(1) and 0.030(7) v.u. The bond valence sum of Cu2 is 2.08(1) v. u. Magnetic Properties In the trinuclear complex unit, three magnetic exchange pathways are to be considered: Pathways between the metallic centers Cu2‧‧‧Cu1, Cu1‧‧‧Cu2i and the pathway between more distant terminal copper atoms Cu2‧‧‧Cu2i (Fig.1). The resulting parameters of magnetic interaction are J12 = J23 = + 75.8 cm−1, and J13 = + 0.4 cm−1. The resulting Mullican spin densities are collected in Table4. As apparent, the DFT calculation predicts strong ferromagnetic interaction between adjacent centers in contrast with negligible interaction between the terminal copper centers. The resulting values were used as input for the following analysis of experimental magnetic functions of 1. The temperature and field dependence of magnetic functions are presented in Fig.5. The increase of susceptibility product at low temperatures indicates the presence of ferromagnetic exchange interaction, while the sudden decrease below 10K can be ascribed to the intermolecular antiferromagnetic interaction (molecular field). The magnetic functions were analyzed according to the following spin Hamiltonian [68]: where mB is the Bohr magneton, gi is the average value of the individual g-factors, Ŝiz is the projection of the spin operator of the i-th center, and zJ is the parameter of the molecular field. The two parameters of magnetic exchange interaction were fixed equal and the parameter between the terminal copper center Cu2 and Cu21 was completely omitted based on the DFT assessment (vide supra). The constraint g1 = g3 is due to the symmetry of complex 2. As can be seen, the fitted and calculated values of magnetic exchange interaction are in very good agreement. The g-factors of the magnetic centers gain typical values for copper(II) centers. Conclusions The trinuclear 2-hydroxybenzoato copper(II) complex with (6,7-dimethoxy-isoquinolin-1-yl)-(3,4-dimethoxy-phenyl)- methanol has been prepared and characterized. The copper atoms of the centrosymmetric trinuclear complex 2 are hexacoordinated. The novel N, O -bonded(6,7-dimethoxyisoquinolin-1-yl)-(3,4-dimethoxyphenyl)-methanolato (papaverinyloxy) ligands are bidentate. One-half of 2-hydroxybenzoates connect atoms Cu1 and Cu2 via bridging –COO groups. The remaining 2-hydroxybenzoates are bonded to the atoms Cu2 by the oxygens of nonbridging carboxylate groups. The positions of their phenyl rings are stabilized by the weak bonding of Cu2 atoms with the hydroxyl groups. A significant difference in Cu1 and Cu2 coordination geometries is consistent with the plasticity of the inner copper coordination of (8)  H =−J 12  S 1  S 2 −J 23  S 2  S 3 +m B g i  S iz −zJ  S i Table 4 Mullikan spin densities for individual broken symmetry spin states of 1 Symmetry code: (i) -x + 1, -y + 1, -z + 2 BS state ρ(Cu2) ρ(Cu1) ρ(Cu2i) βαα − 0.6550 0.6790 0.6570 αβα 0.6552 − 0.6786 0.6552 ααβ 0.6570 0.6790 − 0.6550 Fig. 5 The optimum fit of magnetic functions (solid lines) of compound 1. Experimental values are shown by circles 170 Journal of Chemical Crystallography (2024) 54:163–172 CuO4O’2 and CuNO3O′2 chromophores. The supreme limit for the equidistant equatorial bonds of the last chromophore is 2.92(1) Å. For the equidistant axial bonds, the supreme limit is 3.07(3) Å. If RL≥RS, the dependence RL vs. RS of copper inner coordination of Cun+NO3O′2 chromophore can be approximated by Eq.(3). The magnetism of complex 2 can be well explained by strong ferromagnetic interaction between the neighboring magnetic centers via the dominant pathways Cu1–O5–C35–O4–Cu2, Cu1–O7–Cu2i, and Cu1–OW1–Cu2. The plasticity of inner copper coordination and the papaverinyloxy ligand facilitate the formation of the trinuclear copper complex 2. By the copper activation of a papaverine methyl C–H bond, the hydrogen atom is replaced with bridging oxygen. Supplementary Information X-ray crystallographic data for the title compound (CCDC No.1999284) have been deposited with the Cambridge Crystallographic Data Center (http://www. ccdc.cam.ac.uk/).The datasets generated during and/or the current study are available from the corresponding author upon reasonable request. Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s1087002401009-2. Acknowledgements The article was created with the support of the Ministry of Education, Science, Research and Sport of the Slovak Republic within the Research and Development Operational Programme for the project “University Science Park of STU Bratislava”, ITMS 26240220084, co-funded by the European Regional Development Fund. We are also grateful to the HPC center at the Slovak University of Technology in Bratislava, which is a part of the Slovak Infrastructure of High-Performance Computing [SIVVP project, ITMS code 26230120002, funded by the European regional development funds (ERDF)]. We are grateful to the HPC center at the Slovak University of Technology in Bratislava, which is a part of the Slovak Infrastructure of High Performance Computing (SIVVP project, ITMS code 26230120002, funded by the European region development funds, ERDF), for the computational time and resources made available. This article was created with the support of the Ministry of Education, Science, Research and Sport of the Slovak Republic within the Research and Development Operational Programme. for the project “University Science Park of STU Bratislava”, ITMS 26240220084, for funding within the scheme Excellent research teams, co-funded by the European Regional Development Fund. The authors are grateful to Professor Tibor Gracza, and Dr. Daniel Vegh, for preliminary discussions. Author Contributions FV solved and refined the crystal structure and construcrted the structure correlations, JP and IŠ interpreted the data of magnetical measurements and quantum mechanical calculations, MM prepared the single crystal, JK refined the final structure model, FV and JP wrote the manuscript together. Funding Open access funding provided by The Ministry of Education, Science, Research and Sport of the Slovak Republic in cooperation with Centre for Scientific and Technical Information of the Slovak Republic. Funding for this research was provided by: Vedecká Grantová Agentúra MŠVV M SR a SAV (VEGA) (Grant No. 1/0175/23; Grant No. 1/0029/22), Agentúra na Podporu Vedy a Výskumu (Grant No. 20–0213, Grant No. 19–0087, Grant No. 22-0172, Grant No. DS-FR-22-0010). Data Availability The authors declare thet all other data supportig the findings of this study are available within the article and its supplementary information files. Declarations Competing Interests The authors declare no competing interests. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. References 1. Slavík J, Slavíková L (1994) Alkaloids from Papaver pinnatifidum Moris. Collect Czech Chem Commun 59:1879–1883 2. Utkan T, Sarioglu Y, Utkan NZ, Gönüllü NN, Yildirim MK (1996) Vascular smooth muscle reactivity and endothelium derived relaxing factor in experimental obstructive jaundice. Arch Physiol Biochem 104:30–35. https:// doi. org/ 10. 1076/ apab. 104.1. 30. 12864 3. Pavel IZ, Heller L, Sommerwerk S, Loesche A, Al-Harrasi A, Csuk R (2016) Drotaverine a concealed cytostatic! Arch Pharm 350:e1600289. https:// doi. org/ 10. 1002/ ardp. 20160 0289 4. Czyrski A, Girreser U, Hermann T (2013) 15N NMR chemical shifts in papaverine decomposition products. J Mol Struct 1036:111–114. https:// doi. org/ 10. 1016/j. molst ruc. 2012. 09. 072 5. Girreser U, Hermann TW, Piotrowska K (2003) Oxidation and degradation products of papaverine, part II[1]: investigations on the photochemical degradation of papaverine solutions. Arch Pharm Med Chem 336:401–405. https:// doi. org/ 10. 1002/ ardp. 20030 0751 6. Piotrowska K, Hermann TW, Augustyniak W (2006) Kinetics of photooxidation of papaverine hydrochloride and its major photooxidation products. J Pharm Biomed Anal 41:1391–1395. https:// doi. org/ 10. 1016/j. jpba. 2006. 02. 036 7. Hermann TW, Girreser U, Michalski P, Piotrowska K (2002) Oxidation and degradation products of papaverine. Part I: Gadamer and Schulemann’s papaverinol synthesis revisited. Archive Der Pharmazie 4:167–169. 8. Reynolds CD, Palmer RA, Gorinsky B (1973) Crystal and molecular structure of the alkaloid papaverine hydrochloride. J Cryst Mol Struct 4:213–218 9. Santini C, Pellei M, Gandin V, Porchia M, Tisato F, Marzano C (2014) Advances in copper complexes as anticancer agents. Am Chem Soc Chem Rev. https:// doi. org/ 10. 1021/ cr400 135x 10. Zeleňák V, Sabo M, Massa W, Černák J (2004) Tetra-µbenzoato-κ 8O:O′-bis(1-[(3,4-dimethoxyphenyl)methyl]-6, 7-dimethoxyisoquinoline-κNzinc(II)): the first crystal structure with papaverine as a ligand. Acta Crystallogr Sect C: Cryst Struct 171Journal of Chemical Crystallography (2024) 54:163–172 Commun 60:m85–m87. https:// doi. org/ 10. 1107/ S0108 27010 30288 3X 11. Catterick J, Thornton P (1977) Structures and physical properties of polynuclear carboxylates. Adv Inorg Chem Radiochem 20(08):291–362. https:// doi. org/ 10. 1016/ S00652792 12. Melník M (1981) Mono-, bi-, tetraand polynuclear copper(II) halogenocarboxylates. Coord Chem Rev 36:1–44. https:// doi. org/ 10. 1016/ S00108545(00) 80504-4 13. Melník M (1981) Binuclear caffeine adducts of Cu(II) acetate and cu(II) chloroacetates with unusually high antiferromagnetic interaction. J Inorg Nucl Chem 43:3035–3038. https:// doi. org/ 10. 1016/ 00221902(81) 80677X41 14. Melník M (1982) Structural isomerism of copper(II) compounds. Coord Chem Rev 47:239–261. https:// doi. org/ 10. 1016/ 00108545(82) 85033-9 15. Alig H, Trömel M (1992) Geometrische und chemische koordination. Z für Kristallographie 201:213–222. https:// doi. org/ 10. 1524/ zkri. 1992. 201.34. 213 16. Werner A (1893) Beitrag Zur Konstitution Anorganischer Verbindungen. Z für Anorganische Chemie 3:267–330. https:// doi. org/ 10. 1002/ zaac. 18930 030136 17. Gopinathan MS, Jug K (1983) Valency. I. A quantum chemical definition and properties. Theoret Chim Acta (Berl) 63:497–509. https:// doi. org/ 10. 1007/ BF023 94809 18. Gopinathan MS, Jug K (1983) Valency. II. Applications to molecules with first-row atoms. Theoret Chim Acta (Berl) 63:511–527. https:// doi. org/ 10. 1007/ PL000 20125 19. Jørgensen CK (1969) Stereochemistry and the choice between alternative descriptions of oxidation states. Chem Phys Lett 5:150–151. https:// doi. org/ 10. 1016/ 00092614(69) 80142-9 20. Jørgensen CK (1969) Oxidation numbers and oxidation states. Springer, Berlin, Heidelberg, New York 21. Adams S (2001) Relationship between bond valence and bond softness of alkali halides and chalcogenides. Acta Cryst B 57:278– 287. https:// doi. org/ 10. 1107/ S0108 76810 10030 68 22. Chen H, Adams S (2017) Bond softness sensitive bond-valence parameters for crystal structure plausibility tests. IUCrJ 4:614– 625. https:// doi. org/ 10. 1107/ S2052 25251 70102 11 23. Brown ID (2017) What is the best way to determine bond-valence parameters? IUCrJ 4:514–515. https:// doi. org/ 10. 1107/ s2052 25251 70117 82 24. Alcock NW (1972) Secondary bonding to nonmetallic elements. Advan Inorg Chem Radiochem 15(08):1–58. https:// doi. org/ 10. 1016/ S00652792 25. Gažo J, Bersuker IB, Garaj J, .Kabešová M, Kohout J, Langfelderová M, Melník M, Serátor M, Valach F (1976) Plasticity of the coordination sphere of copper(II) complexes, its manifestation and causes. Coord Chem Rev 19:253–298. https:// doi. org/ 10. 1016/ S00108545(00) 80317-3 26. Gažo J, Boča R, Jóna E, Kabešová M, Macášková L, Šima J, Pelikán P, Valach F (1982) Equatorial-axial interactions in complexes as a manifestation of mutual influence of ligands. Coord Chem Rev 43:87–131. https:// doi. org/ 10. 1016/ S00108545(00) 82093-7 27. Sundberg MR (2000) Effect of hydrogen bonding on coordinate bond. Semi-coordination in Werner-type copper(II) complexes revisited. Rev Inorg Chem 20:195–218. https:// doi. org/ 10. 1515/ REVIC. 2000. 20.3. 195 28. Brown ID (1992) Chemical and steric constraints in inorganic solids. Acta Cryst 48:553–572. https:// doi. org/ 10. 1107/ S0108 76819 20024 53 29. Brown ID (1981) The Bond-Valence Method: An Empirical Approach to Chemical Structure and Bonding. In: Structure and Bonding in Crystals. Edited by M. O’Keefee & A. Navrotsky, New York: Acadamic Press, Inc. pp. 1–30 30. Lehn JM (2012) Constitutional dynamic chemistry: Bridge from supramolecular chemistry to adaptive chemistry. Top Curr Chem 322:1–32. https:// doi. org/ 10. 1007/ 128_ 2011_ 256 31. Rodríguez-Fortea A, Ruiz E, Alemany P, Alvarez S (2003) Magneto-Structural correlations in trinuclear Cu(II) complexes: a density functional study. Monatshefte Chem 134:307–316. https:// doi. org/ 10. 1007/ s007060020523-6 32. Kawata T, Uekusa H, Ohba S, Furukawa T, Tokii T, Muto Y, Kato M (1992) Magneto-structural correlation in dimeric copper(II) benzoates. Acta Cryst B48:253–261. https:// doi. org/ 10. 1107/ S0108 76819 10126 97 33. Yamanaka M, Uekusa H, Ohba S, Saito Y, Iwata S, Kato M, Tokii T, Muto Y, Steward OW (1991) Correlation of electron density and spin-exchange interaction in dimeric copper(II) formates, acetates and silanecarboxylates. Acta Cryst B 47:344– 355. https:// doi. org/ 10. 1107/ S0108 76819 00124 59 34. Sheldrick GM (2008) A short history of SHELX. Acta Cryst A64:112–122. https:// doi. org/ 10. 1107/ S0108 76730 70439 30 35. Dolomanov OV, Bourhis RJ, Gildea RJ, Howard JAK, Puschmann H (2009) OLEX2: a complete structure solution, refinement and analysis program. J Appl Cryst 42:339–341. https:// doi. org/ 10. 1107/ S0021 88980 80427 26 36. Sheldrick GM (2015) Crystal structure refinement with SHELXL. Acta Crystallogr Sect C Struct Chem 71:3–8. https:// doi. org/ 10. 1107/ S2053 22961 40242 18 37. Figgis BN, Nyholm RS (1958) A convenient solid for calibration of the Gouy magnetic susceptibility apparatus. J Chem Soc 4:190–4191 38. Valach F (1999) A bond-valence approach to the semicoordination of copper-oxygen and copper-nitrogen complexes. Polyhedron 18:699–706. https:// doi. org/ 10. 1016/ S02775387(98) 00342-8 39. Brown ID, Altermatt D (1985) Bond-valence parameters obtained from a systematic analysis of the Inorganic Crystal structure database. Acta Cryst B41:244–247. https:// doi. org/ 10. 1107/ S0108 76818 50020 63 40. Urusov VS (1995) Semi-empirical groundwork of the bond‐ valence model. Acta Cryst B51:641–649. https:// doi. org/ 10. 1107/ S0108 76819 50034 17 41. Valach F, Tokarčík M, Maris T, Watkin DJ, Prout CK (2000) Bond-valence approach to the copper-copper and copper-oxygen bonding in binuclear copper(II) complexes: structure of tetrakis(2fluorobenzoato-O,O′)-bis(2-fluorobenzoate-O) dicopper(II). Z Kristallogr 215:56–60. https:// doi. org/ 10. 1524/ zkri. 2000. 215.1. 56 42. Valach F, Tokarčík M, Saunders A, Cowley A, Watkin DJ (2001) Structural isomerism of propionato copper(II) complex: Crystal and molecular structure of tetrakis(µ-o-propionato)bis(methyl 3-pyridyl-N-carbamate)dicopper(II) at 190 K. Polyhedron 20:1933–1937. https:// doi. org/ 10. 1016/ S02775387(01) 00784-7 43. König E (1966) Magnetic properties of coordination and organometallic transition metal compounds. Springer-, Berlin 44. Yamaguchi K, Takahara Y, Fueno T (1986) In: Applied Quantum Chemistry, V. H. Smith, V. Reidel (Ed.), Dordrecht, p. 155 45. Soda T, KitagawaY, Onishi T, Takano Y, Shigeta Y, Nagao H, Yoshika Y, Yamaguchi K (2000) Abinitio computations of effective exchange integrals for H-H, H-He-H and Mn2O2 complex: comparison of broken-symmetry approaches. Chem Phys Lett 319:223–230. https:// doi. org/ 10. 1016/ S00092614(00) 00166-4 46. Neese F (2009) Prediction of molecular properties and molecular spectroscopy with density functional theory: from fundamental theory to exchange-coupling. Coord Chem Rev 253:526–563. https:// doi. org/ 10. 1016/j. ccr. 2008. 05. 014 47. Chilton NF, Anderson RP, Turner LD, Soncini A, Murray KS (2013) PHI: a powerful new program for the analysis of anisotropic monomeric and exchange-coupled polynuclear dand