Next Article in Journal
Understanding Avobenzone Crystallization in Sunscreen Formulations: Role of Metal Oxide-Driven Nucleation and Stabilization Strategies
Next Article in Special Issue
Structural and Electronic Properties of U5M+ and T5M+ (U = Uracil, T = Thymine, M = Ag and Au) Cluster Cations
Previous Article in Journal
Photo-Induced Degradation of Priority Air Pollutants on TiO2-Based Coatings in Indoor and Outdoor Environments—A Mechanistic View of the Processes at the Air/Catalyst Interface
Previous Article in Special Issue
Synthesis and Crystal Structure Analysis of Some Aromatic Imines of Syringaldehyde
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two Modifications of Nitrilotris(methylenephenylphosphinic) Acid: A Polymeric Network with Intermolecular (O=P–O–H)3 vs. Monomeric Molecules with Intramolecular (O=P–O–H)3 Hydrogen Bond Cyclotrimers

1
Institut für Anorganische Chemie, TU Bergakademie Freiberg, D-09596 Freiberg, Germany
2
Institut für Analytische Chemie, TU Bergakademie Freiberg, D-09596 Freiberg, Germany
3
Institute of Resource Ecology, Helmholtz-Zentrum Dresden-Rossendorf eV, D-01328 Dresden, Germany
4
Zentrum für Effiziente Hochtemperaturstoffwandlung (ZeHS), TU Bergakademie Freiberg, D-09596 Freiberg, Germany
*
Author to whom correspondence should be addressed.
Crystals 2024, 14(7), 662; https://doi.org/10.3390/cryst14070662
Submission received: 28 June 2024 / Revised: 9 July 2024 / Accepted: 16 July 2024 / Published: 19 July 2024
(This article belongs to the Special Issue Different Kinds of Hydrogen Bonds in Crystal Structures)

Abstract

:
Nitrilotris(methylenephenylphosphinic) acid (NTPAH3) was silylated using hexamethyldisilazane to produce the tris(trimethylsilyl) derivative NTPA(SiMe3)3. From the latter, upon alcoholysis in chloroform, NTPAH3 could be recovered. Thus, a new modification of that phosphinic acid formed. Meanwhile, NTPAH3 synthesized in aqueous hydrochloric acid crystallized in the space group P3c1 with the formation of O-H⋅⋅⋅O H-bonded networks (NTPAH3P), in chloroform crystals in the space group R3c formed (NTPAH3M), the constituents of which are individual molecules with exclusively intramolecular O-H⋅⋅⋅O hydrogen bonds. Both solids, NTPAH3P and NTPAH3M, were characterized by single-crystal X-ray diffraction, multi-nuclear (1H, 13C, 31P) solid-state NMR spectroscopy, and IR spectroscopy as well as quantum chemical calculations (both of their individual constituents as isolated molecules as well as in the periodic crystal environment). In spite of the different stabilities of their constituting molecular conformers, the different crystal packing interactions rendered the modifications of NTPAH3P and NTPAH3M similarly stable. In both solids, the protons of the acid are engaged in cyclic (O=P–O–H)3 H-bond trimers. Thus, the trialkylamine N atom of this compound is not protonated. IR and 1H NMR spectroscopy of these solids indicated stronger H-bonds in the (O=P–O–H)3 H-bond trimers of NTPAH3M over those in NTPAH3P.

Graphical Abstract

1. Introduction

Nitrilotris(methylenephosphinic) acids are essentially related to nitrilotris(acetic) acid I (Figure 1), which has been shown to crystallize as a zwitterionic compound [1,2]. The protonation of the α-amine N atom by one of the acidic groups is a common feature, which is very well known from biogenic α-amino acids [3]. Upon replacing carboxylic with phosphinic (compound II, [4]) or phosphonic acid groups (compound III [5]), the slightly more acidic P-containing acid groups (for example, PhMePOOH [6] and arylphosphinic acids of the type Aryl(H)POOH [7] were shown to be more acidic than benzoic acid) serve as protonators, and this is also true for nitrilotris(methylenephosphonic) acid IV [8] and compound V [9]. In contrast to the double-zwitterionic nature of V, cyclohexane-1,2-diamine derivative VI [10] has mono-zwitterionic features. The lone pair of its second N atom is involved in an N–H⋅⋅⋅N hydrogen bond and thus less susceptible to protonation.
In previous studies [11,12] nitrilotris(methylenephenylphosphinic) acid (NTPAH3) was shown to be an interesting tripodal ligand. In the course of our ongoing investigations of Si-O-P compounds [13,14,15], phosphinic acids revealed interesting coordination features at silicon [16]. In this context, we prepared both NTPAH3 and its trimethylsilyl derivative NTPA(SiMe3)3 as starting materials for further investigations. Thereby, we encountered highly interesting features of the molecular structure of NTPAH3 itself in the solid state, which distinguish this acid from related aminoalkyl functionalized phosphinic and phosphonic (and carboxylic) acids. Both the preference of formation of O–H⋅⋅⋅O hydrogen bonds over the protonation of the intramolecular base (i.e., the alkylamine N atom) as well as the capability of solvent-dependent formation of intramolecular vs. intermolecular (O=P–O–H)3 H-bond trimers render crystalline NTPAH3 an interesting hydrogen bonding system. While the former stabilizes monomeric NTPAH3, the latter renders this acid an interesting candidate as a building block for hydrogen-bonded organic frameworks (HOFs) [17,18,19]. As hydrogen bonding itself is of current interest (e.g., with respect to spectroscopic properties [20] and as a bond type that may complement other bond types in H-bond cross-linked materials [21]), we decided to have a closer look at this particular H-bond system of NTPAH3.

2. Materials and Methods

2.1. General Considerations

Starting materials NH4Cl (ORG-Laborchemie, Bunde, Germany, 99.7%), phenylphosphinic acid (Sigma-Aldrich, Steinheim, Germany, 99%), paraformaldehyde (thermo scientific, Kandel, Germany, 96%), hydrochloric acid (VWR chemicals, Fontenay-sous-bois, France, 37%), hexamethyldisilazane (abcr, Karlsruhe, Germany, 98.5%), chloroform, stabilized with amylenes (Fisher Scientific, Loughborough, UK, 99.8%), (CD3)2SO (deutero, Kastellaun, Germany, 99.8%) and CDCl3 (deutero, Kastellaun, Germany, 99.8%) were obtained from commercial suppliers. The following solvents were dried prior to use: chloroform and (CD3)2SO were stored over 3 Å molecular sieves (Carl Roth, Karlsruhe, Germany); CDCl3 was distilled from CaH2 (Sigma-Aldrich, Steinheim, Germany, 95%) and kept under argon atmosphere.
Synthesis and characterization of NTPA(SiMe3)3 were carried out under an atmosphere of dry argon utilizing standard Schlenk techniques. Solution NMR spectra (1H, 13C, 29Si, 31P) (cf. Figures S1–S7 in the supporting information) were recorded on a Nanobay 400 MHz spectrometer (Bruker Biospin, Ettlingen, Germany). 1H, 13C and 29Si chemical shifts are reported relative to Me4Si (0 ppm), either as internal reference or in accordance with the shift of the solvent signal ((CD3)2SO at δ 39.5 ppm in case of the 13C NMR spectrum of NTPAH3), 31P chemical shifts are reported relative to 85% H3PO4 (from an external calibration) set at 0 ppm. Solid-state NMR experiments (for spectra cf. Figures S8–S13 in the supporting information) were performed on an Avance HD 400 MHz WB spectrometer (Bruker Biospin, Ettlingen, Germany) with ZrO2 rotors. The 13C and 31P spectra were recorded using a 4 mm triple resonance CP MAS DVT probe with a spinning rate of 10 kHz. The 1H spectra were measured on a 2.5 mm CP MAS VTN probe with a spinning rate of 30 kHz. The 31P NMR chemical shift was referenced to 85% H3PO4 using NH4H2PO4 as external secondary standard. The 13C chemical shift was referenced to TMS using adamantane as external secondary standard. Cross-polarization NMR experiments were carried out with 1 ms contact time for 31P and 2 ms for 13C. An 80% ramp was used for 31P and a 70% ramp for 13C. The recycle delay of all 13C and 31P CP experiments was 30 s and for 1H 10 s. IR spectra were measured at room temperature using a Nicolet 380 instrument (Thermo Fisher, Waltham, MA, USA). The solid material was ground with dry KBr in a mortar and pressed into a pellet. For each measurement, 32 scans were collected using wavenumbers from 400 cm−1 to 4000 cm−1. The spectra were corrected for background (prior to each measurement the background was collected using a freshly prepared KBr pellet).
For single-crystal X-ray diffraction analyses, crystals were selected under an inert oil and mounted on a glass capillary (which was coated with silicone grease). Diffraction data were collected on a Stoe IPDS-2/2T diffractometer (STOE, Darmstadt, Germany) using Mo Kα-radiation. Data integration was performed with the STOE software XArea version 2.3. The structures were solved using SHELXT and refined with the full-matrix least-squares methods of F2 against all reflections with SHELXL-2019/3 [22,23,24]. All non-hydrogen atoms were anisotropically refined, C-bound hydrogen atoms were isotropically refined in idealized position (riding model), and hydrogen atoms of OH groups were located as residual electron density peaks and were refined without restraints. For details of data collection and refinement, see Appendix A, Table A1. Graphics of molecular structures were generated with ORTEP-3 [25,26] and POV-Ray 3.7 [27]. CCDC 2365392 (NTPAH3P), 2365391 (NTPAH3M) and 2365393 (NTPA(SiMe3)3) contain the supplementary crystal data for this article. These data can be obtained free of charge from the Cambridge Crystallographic Data Centre via https://www.ccdc.cam.ac.uk/structures/ (accessed on 25 June 2024).
The geometry optimizations of isolated molecules were carried out with ORCA 5.0.3 [28] using the restricted PBE0 functional with a relativistically recontracted Karlsruhe basis set ZORA-def2-TZVPP [29,30] for all atoms, scalar relativistic ZORA Hamiltonian [31,32], atom-pairwise dispersion correction and with the Becke–Johnson damping scheme (D3BJ) [33,34]. Calculations were started from the molecular structures obtained by single-crystal X-ray diffraction analysis. Numerical frequency calculations were performed to prove convergence at the local minimum after geometry optimization and to obtain the Gibbs free energy (293.15 K). The crystal structures of NTPAH3P, and NTPAH3M were optimized with TURBOMOLE rev. V7-7-1 [35] using the restricted PBE functional with a basis set pob-TZVP and RI-J auxiliary basis set for all atoms. A grid size of m3 and atom-pairwise dispersion correction with a Becke–Johnson damping scheme (D3BJ) [33,34] were used. During the periodic boundary optimization, the redundant coordinates were optimized and cell parameters were held constant. At the final geometry, a single-point calculation using PBE0 functional was performed. QTAIM analyses were performed with Multiwfn 3.8 [36] starting from the results obtained with the PBE0 functional. Graphics were generated using Chemcraft version 1.8 (Build 164) [37] or TURBOMOLE rev. V7-7-1 [35]. The R-factor for the match of 3-fold symmetry of the individual optimized molecular structures was evaluated with Chemcraft [37].

2.2. Syntheses and Characterization

Compound NTPAH3P (C21H24NO6P3). This compound was synthesized by following a published protocol [38] using NH4Cl (1.411 g, 26.37 mmol), paraformaldehyde (4.752 g, 158.2 mmol), phenylphosphinic acid (11.24 g, 4.232 mmol), 60 mL of water and 60 mL of hydrochloric acid. Yield: 5.98 g (12.5 mmol 47%) (cf. lit. 51.5% [38]) M.p. 238 °C (cf. lit. 233-236 °C [38]). As the original report [38] did not list NMR and IR spectroscopic data, we list the data obtained from our product. NTPAH3 in solution: 1H NMR (400 MHz, (CD3)2SO): δ (ppm) 13.07 (s, 3H, OH), 7.67-7.57 (m, 9H, Ph-o,p), 7.48 (m, 6H, Ph-m), 3.52 (d, 2JPH = 4.4 Hz, 6H, CH2); 13C NMR (101 MHz, (CD3)2SO): δ (ppm) 132.2 (Ph-p), 131.4 (d, 2JPC = 10.6 Hz, Ph-o), 130.9 (d, 1JPC = 131.5 Hz, Ph-i), 128.4 (d, 3JPC = 12.8 Hz, Ph-m), 58.9 (dt, 1JPC = 108.1 Hz, 3JPC = 5.2 Hz, CH2); 31P{1H} NMR (162 MHz, (CD3)2SO): δ (ppm) 33.9. Solid NTPAH3P: IR: υ ¯ (cm−1) 419 (m), 476 (s), 537 (s), 694 (s), 752 (s), 763 (w), 796 (m), 886 (s), 955 (s), 987 (m), 1029 (w), 1071 (m), 1099 (m), 1099 (m), 1126 (m), 1157 (m), 1182 (s), 1246 (s), 1313 (w), 1326 (w), 1344 (m), 1394 (w), 1424 (s), 1440 (s), 1488 (w), 1592 (m), 1684 (b), 2126 (b), 2266 (b), 2560 (b), 2802 (m), 2846 (w), 2968 (m), 3056 (w).
Compound NTPA(SiMe3)3 (C30H48NO6P3Si3). In a 100 mL Schlenk flask with magnetic stirrer and reflux condenser, 4.793 g (10.0 mmol) of NTPAH3P were suspended in 5 mL (ca. 3.85 g, 23.9 mmol) of hexamethyldisilazane and heated with stirring. The reaction mixture was refluxed under inert atmosphere until the evolution of ammonia had ceased and a clear colorless solution was obtained. This solution was allowed to cool to room temperature within few minutes to afford a white precipitate, whereupon the supernatant was removed with a syringe and the white solid was dried in vacuo. Yield: 6.80 g (9.77 mmol, 98%).
The P atoms of this compound are centers of chirality. In CDCl3 solution, this compound forms a set of diastereomeric pairs of enantiomers, the RRR/SSS pair, which was encountered in the crystal structure, and the RRS/SSR pair. We observed an RRR/SSS: RRS/SSR ratio of 85:15. Because of the chemical identity of the three phosphinic acid moieties, the former pair gives rise to a single set of signals. The latter gives rise to a double set of signals in relative 2:1 intensity ratio (these additional signals are shown in the Supplementary Material in Figures S4–S7 and listed in Table S1). RRR/SSS pair: 1H NMR (400 MHz, CDCl3): δ (ppm) = 7.43 (appt. tq, 3JHH = 7.5 Hz, 4JHH + 5JPH = 1.3 Hz, 3 H, Ph-p), 7.16 (appt. td, 4JPH = 3.6 Hz, 6 H, Ph-m), 7.05 (ddd, 3JPH = 12.2 Hz, 3JHH = 8.1 Hz, 4JHH = 1.4 Hz, 6 H, Ph-o), 4.41 (tt, 2JHH + 2JPH = 14.7, 4JPH = 2.6 Hz, 1H, CH2), 2.60 (dd, 2JHH = 14.8, 2JPH = 11.4 Hz, 1H, CH2), 0.17 (s, 9H, Si(CH3)3); 13C{1H} NMR (100 MHz, CDCl3): δ (ppm) 132.0 (d, 2JPC = 10.5 Hz, Ph-o), 131.6 (d, 1JPC = 132.3 Hz, Ph-i), 131.0 (d, 4JPC = 2.8 Hz, Ph-p), 127.8 (d, 3JPC = 13.1 Hz, Ph-m), 57.2 (dt, 1JPC = 122.6 Hz, 3JPC = 9.8 Hz, CH2), 1.27 (d, 3JPC = 1.2 Hz, Si(CH3)3); 29Si{1H} NMR (79 MHz, CDCl3): δ (ppm) = 22.2 (d, 3JPSi = 9.3 Hz); 31P{1H} NMR (162 MHz, CDCl3): δ (ppm) = 29.4.
A coarse crystalline product (used for single-crystal X-ray diffraction analyses) was obtained according to a related protocol but using more hexamethyldisilazane as solvent, i.e., 2.08 g (4.34 mmol) of NTPAH3P in 20 mL (ca. 15.4 g, 95.6 mmol) of hexamethyldisilazane. The greater amount of solvent caused a lower yield of isolated crystalline product (2.24 g, 3.22 mmol, 74%).
Compound NTPAH3M (C21H24NO6P3). In a 10 mL Schlenk flask, NTPA(SiMe3)3 (0.569 g, 0.818 mmol) was dissolved in chloroform (5 mL) and anhydrous EtOH (1 mL) was added carefully (chloroform phase layered with ethanol), whereupon a small amount of white precipitate formed upon contact between the two phases within the first minute. Eventually, long crystalline needles of the product formed upon standing overnight. Thereafter, the solvent was decanted, and the white solid was dried in vacuo. Yield: 0.361 g (0.753 mmol, 92%). M.p. 238 °C. IR: υ ¯ (cm−1) 481 (m), 559 (s), 690(s), 712 (m), 722 (m), 734 (m), 762 (m), 957 (s), 1024 (w), 1035 (w), 1068 (w), 1097 (m), 1126 (s), 1154 (m), 1269 (w), 1310 (m), 1410 (m) 1439 (m), 1538 (b), 1595 (m), 2883 (m), 2937 (w), 2992 (w), 3051 (m).

3. Results and Discussion

3.1. Syntheses of NTPAH3 and Its Trimethylsilyl Ester NTPA(SiMe3)3

The nitrilotriphosphinic acid NTPAH3 in our study was prepared according to a previously published protocol [38] (Scheme 1), and the product obtained in this procedure (referred to as modification NTPAH3P; the superscript index P indicates its H-bonded polymeric nature in the solid state) was used for our study. Silylation of NTPAH3 in excess hexamethyldisilazane (which served both as the silylating reagent and as the solvent) afforded the tris(trimethylsilyl) derivative NTPA(SiMe3)3, which crystallized from the reaction solution upon cooling to room temperature. This compound is highly soluble in chloroform and, in this solution, highly sensitive toward protolysis (e.g., alcoholysis). Thus, deliberate alcoholysis of NTPA(SiMe3)3 in chloroform afforded long needles of NTPAH3 in a new modification (referred to as modification NTPAH3M, the superscript index M indicates its monomeric nature caused by exclusively intramolecular H-bonds).

3.2. Crystallographic Analysis of the Molecular Structures of NTPA(SiMe3)3, NTPAH3P and NTPAH3M

Compound NTPA(SiMe3)3 crystallized in the trigonal space group type R 3 ¯ with two independent thirds of molecules in the asymmetric unit (Figure 2, Appendix A Table A1). Their N atoms are located on crystallographic 3-fold rotation axes (Figure 2b,d). The Si–O bond lengths (1.6756(14) and 1.6760(14) Å) are similar to those encountered with trimethylsilyl phosphate O=P(OSiMe3)3 (1.674(3) Å, [39]) and a derivative of a cyclic phosphinic acid, (Me3SiO)(O=)P[-C(SiMe3)2–CH2–CH=C(SiMe3)-] (1.672(2) Å, [40]). The latter was the only hit in a CSD search [41] for crystallographically characterized triorganosilyl derivatives of phosphinic acids. Thus, the crystal structure of NTPA(SiMe3)3 bears some novelty in this regard. Further metric data of these molecules (Table 1) will be discussed in the context of the further analysis of the acid NTPAH3.
Crystals of both modifications NTPAH3P and NTPAH3M were analyzed by single-crystal X-ray diffraction (Figure 3 and Figure 4, Table 1 and Table A1). Both modifications are trigonal (space group types P3c1 and R3c, respectively), and their asymmetric units consist of one third of a molecule of NTPAH3 with its nitrogen atom located on a crystallographically imposed 3-fold rotation axis. Even though the three P-atoms within one molecule exhibit the same chirality, both crystal structures accommodate the respective (R,R,R)/(S,S,S)-pair of enantiomers because of the c glide in their space groups types. In both cases, the N atom is not protonated; i.e., the phosphinic acid moieties are retained. Hence, the molecules are not zwitterionic. This was an unexpected observation, as in both cases the trialkylamine N atoms’ lone pairs are essentially vacant (they are not involved in any other hydrogen bonds, the structures do not contain any other H-bond donors such as ammonium ions, water molecules etc., which could interfere with the N-located lone pair). This can be attributed to the involvement of all of the P(O)(OH) groups’ hydrogen atoms in cyclic hydrogen bond patterns, i.e., R ( 12 ) 3 3 motifs [42] of the type (O=P–O–H)3 (vide infra). The difference between the molecular conformations of NTPAH3 contained in its two modifications is essentially based on an approximate 110° torsion about the CH2–P bond (N–C–P–C torsion angles are 66.7(2)° in NTPAH3P and 176.3(1)° in NTPAH3M), as shown in Figure 3. This conformational difference gives rise to entirely different inter- and intramolecular interactions. While in modification NTPAH3P, the three phenyl groups form a calix shape and the P(O)(OH) moieties point away from the central N atom (Figure 4a,c), in modification NTPAH3M, the phenyl groups represent the periphery, and the three P(O)(OH) moieties are oriented toward the central 3-fold rotation axis (Figure 4b,d). In the latter case, this orientation results in the formation of a cyclic hydrogen bond system (indicated by the thin blue lines in Figure 4d). In modification NTPAH3P, a related R ( 12 ) 3 3 hydrogen bond motif is formed, also about a crystallographically imposed 3-fold rotation axis, but in an intermolecular manner (cf. Figure 4e). The molecular conformation of NTPAH3 in NTPAH3P resembles the conformation of NTPA(SiMe3)3.
For a comparison of the individual molecules, a set of selected corresponding bond lengths and angles of compounds NTPA(SiMe3)3, NTPAH3P and NTPAH3M is listed in Table 1. The N-C-P angles in NTPA(SiMe3)3 and NTPAH3P are noticeably wider than the tetrahedral angle, and we interpret this widening as originating from the steric congestion about the CH2 groups. This angle is much closer to the tetrahedral angle in the cage structure of NTPAH3M, which lacks that kind of steric repulsion of the N- and P-bound substituents. In the compounds studied here, the coordination spheres about the P atoms are distorted tetrahedrally, which is to be expected from the VSEPR effect of the P=O bond. Thus, the widest angles are associated with X-P=O (X = C, O) moieties. Interestingly, the O-P=O angles are slightly wider than the C-P=O angles, and two reasons can be considered for this phenomenon: The C(Ph)-P=O angle is particularly small (smaller than the C(CH2)-P=O angle) because of C–H⋅⋅⋅O=P attraction, which involves the P=O moiety and an ortho-H atom of the phenyl group. The C(CH2)-P=O angle is still smaller than the O-P=O angle. The P–O bond may feature partial multi-bond character (in case of the P–O–Si moiety because of the electropositive silyl substituent, which renders the P–O-bound oxygen electron-rich, and in the case of the P–O–H moiety, because of the hydrogen bridge, which renders the P–O-bound oxygen electron-rich and also lowers the P=O double bond character). That would explain the enhanced O–P=O repulsion vs. relatively lowered C–P=O repulsion.
The inter- vs. intramolecular O–H⋅⋅⋅O hydrogen bonding in NTPAH3P and NTPAH3M, respectively, results in further differences between the two crystal structures: In NTPAH3P the C-N-C angles are 109.6(2)°. Thus, they are similar to the C-N-C angles in NTPA(SiMe3)3 (110.4(1) and 111.1(1)°). In NTPA3M, the corresponding angles are significantly wider (116.0(1)°). We attribute this angle widening to the strain associated with the cage formation by the intramolecular H-bonds. The metrics of the R ( 12 ) 3 3 motifs are different, too. With O-O-distances of 2.487(3) Å in the O–H⋅⋅⋅O motifs of NTPAH3P its intermolecular hydrogen bonds appear to be weaker than the intramolecular H-bonds in NTPAH3M with O-O-distances of 2.450(2) Å. The P=O and P–O bonds seem to exhibit notable differences, too (1.494(2) and 1.547(2) Å, respectively, in NTPAH3P vs. 1.510(2) and 1.533(2) Å, respectively, in NTPAH3M). As these features may suffer from unresolved disorder effects (simultaneous cyclic disordering of O–H⋅⋅⋅O vs. O⋅⋅⋅H–O), we employed computational analyses for closer inspection of these hydrogen bond systems (cf. Section 3.3.).
The different molecular conformations of NTPAH3 in its modifications result in different intra- and intermolecular interactions in the solid state. In particular, in NTPAH3M the O–H⋅⋅⋅O hydrogen bonds are an intramolecular feature, whereas in NTPAH3P they contribute to the intermolecular packing. To account for related kinds of interactions in a comparable manner, a Hirshfeld surface analysis was performed for the asymmetric units (which are one third molecule in each case), and therefore the O–H⋅⋅⋅O hydrogen bonds were accounted for in the same manner (Selected dnorm isosurface plots of the Hirshfeld surface analyses are contained in the Supplementary Material Figures S14 and S15). The dnorm ranges (−1.092…+2.229 for NTPAH3P, −1.094…+1.167 for NTPAH3M) vary noticeably, and the pronounced upper value of +2.229 in NTPAH3P indicates less efficient packing in this modification (well in accordance with the volume per molecule, which is much larger for NTPAH3P than for NTPAH3M with 596 vs. 549 Å3, respectively). Figure 5 shows the fingerprint plots (de vs. di) of these analyses. Note: Because of the analysis of the asymmetric unit, the N–C bonds (and adjacent intramolecular N-H interactions) are also shown in the fingerprint plots, even though they are intramolecular features and determine the lowest dnorm values in both analyses. The O-H interactions are dominated by the O–H⋅⋅⋅O hydrogen bonds. Therefore, these contributions are representative of the (intermolecular and intramolecular, respectively) O–H⋅⋅⋅O hydrogen bond motifs in NTPAH3P and NTPAH3M, respectively. (C–H⋅⋅⋅O contacts are also present in both modifications; they contribute to the longer distance O-H interactions. For comparison of fingerprint plots (de vs. di) of all O-H interactions see Figure S16). Whereas these parts of the fingerprint plots are similar for the two modifications, great differences are evident for the long-range interactions part of the maps, where C-C, C-H and H-H interactions can be found, and C-C interactions are responsible for the greatest difference. Figure 6 shows that the long-range C-C interactions in NTPAH3P are caused by the calix shape created by the three phenyl groups about the 3-fold rotation axis. The faces of the phenyl rings inside the calix are thus protected from closer packing. In contrast, in NTPAH3M, the phenyl groups are more accessible for intermolecular packing from all sides. In addition to the closer C-C interactions, a larger fraction of C-C contacts on the packing of NTPAH3M results therefrom (5.8% in NTPAH3M, 2.8% in NTPAH3P) vs. a larger fraction of C-H contacts in NTPAH3P (16.8% in NTPAH3P, 12.4% in NTPAH3M), see diagram in Figure 7. In contrast, the O-H interactions (and thus the O–H⋅⋅⋅O hydrogen bonds in particular) contribute to surprisingly similar extent (26.2% in NTPAH3P, 27.8% in NTPAH3M).
Note: As the O-bound hydrogen atoms in the crystal structures of NTPAH3P and NTPAH3M were refined without restraints, the O–H bond lengths obtained from crystal structure refinement varied noticeably (0.62(4) Å in NTPAH3P, 0.93(5) Å in NTPAH3M). It is well known that there are problems associated with locating H atoms from X-ray diffraction data. In order to rule out the possibility that this deviation of different refined O–H bond lengths has greater influence on our Hirshfeld surface analysis, we performed the same analysis with structure models that were obtained by a refinement with idealized OH groups (riding model, H atom attached with the ShelX-code AFIX 147). In principle, it reflects the same findings (same appearance of fingerprint maps, same dnorm ranges). A representation of the contributions of contacts (corresponding to Figure 7, including the models with idealized OH groups) is included in the Supplementary Material, Figure S17.

3.3. Computational Analyses of NTPAH3P and NTPAH3M

For energetic comparison of the two modifications of NTPAH3 and for closer inspection of the H-bond situation and its influence on bond lengths, we have optimized the molecular structures of NTPAH3P and NTPAH3M both as isolated molecules (NTPAH3POPT-1 and NTPAH3MOPT-1) and in their periodic crystal environment (NTPAH3POPT-CRYST and NTPAH3MOPT-CRYST). For details of the computational analyses, please see Section 2.1. and the Supplementary Material (Figures S18–S21, Tables S2–S5).
For both approaches of optimization of the molecular structures, the atomic coordinates obtained from crystallographic analyses were used as starting point, and optimizations were performed without constraints of symmetry operations. With respect to the latter, in all cases the 3-fold molecular symmetry was essentially retained in the local minima molecular conformations (with R-factors for the 3-fold symmetry match of 0.082 (NTPAH3MOPT-1), 0.008 (NTPAH3POPT-1), 0.028 (NTPAH3MOPT-CRYST) and 0.013 (NTPAH3POPT-CRYST), with an R-factor of 0.000 representing a perfect C3 symmetry).
Comparison of the total energies revealed that the two modifications exhibit very similar stability with NTPAH3MOPT-CRYST being only marginally less stable than NTPAH3POPT-CRYST by 0.34 kcal⋅mol−1 (PBE) or 1.66 kcal·mol−1 (single-point calculation with PBE0). As an isolated molecule, NTPAH3MOPT-1 is noticeably more stable than NTPAH3POPT-1 by −20.6 kcal·mol−1 (PBE0). This difference is mainly attributed to the presence vs. absence of the O–H···O hydrogen bonds in NTPAH3MOPT-1 and NTPAH3POPT-1, respectively. Nonetheless, the energy involved in formation of the R ( 12 ) 3 3 H-bond trimer can be expected to exceed that value. Dimerization energies of phosphinic acids (formation of R ( 8 ) 2 2 motifs by two O–H⋅⋅⋅O hydrogen bonds) were experimentally determined at 21 ± 6 to 60 ± 10 kcal·mol−1 [44,45], which would imply H-bond energies around 10.5 to 30 kcal·mol−1 per H-bond and thus energies around 31.5 to 90 kcal·mol−1 for related H-bond cyclo-trimers. However, the energy of the molecular conformation of NTPAH3POPT-1 also gains contributions from intramolecular interactions, which are associated with the formation of the (P-Ph)3 calix and which are absent in NTPAH3MOPT-1. These contributions serve as some compensation to the expected energy difference. Therefore, we analyzed the wave function features of the H-bonds to evaluate their energetic properties. Using quantum theory of atoms in molecules (QTAIM) analyses, we observed bond-critical points at the O–H···O hydrogen bonds. Emamian et al. [46] reported that via the electron density at the bond critical point, the hydrogen bond strength can be estimated. For NTPAH3MOPT-1, an average H-bond energy of 17 kcal·mol−1 per H-bond is observed, and for each of the sets, NTPAH3MOPT-CRYST and NTPAH3POPT-CRYST, it is ca. 20 kcal·mol−1. These values are in the above-mentioned expected range of 10.5 to 30 kcal·mol−1 per H-bond. The energetic difference between the H-bonds of 0.4 kcal·mol−1 (20.3 kcal·mol−1 for NTPAH3MOPT-CRYST and 19.9 kcal·mol−1 for NTPAH3POPT-CRYST) vaguely indicates stronger H-bonds in NTPAH3M but is too marginal for a reliable conclusion.
Closer inspection of the bonds involved in the H-bond systems in NTPAH3P and NTPAH3M (Table 2) indicates differences between the strength of the R ( 12 ) 3 3 H-bond trimers in these two modifications. While in the silyl derivative NTPA(SiMe3)3 the P–O bond is ca. 0.10 Å longer than the P=O bond, this bond length difference is less pronounced in the crystal structures of the acids NTPAH3P and NTPAH3M (with differences of 0.05 and 0.02 Å, respectively). The enhanced averaging of P–O and P=O bond lengths in NTPAH3M could be an indication of a stronger H-bond system with respect to NTPAH3P, and the shorter O···O distance (by 0.04 Å) in NTPAH3M provides further support. These observations may (at least in part) also originate from unresolved O–H···O vs. O⋅⋅⋅H–O disorder in the crystal structures, and therefore further support is essential to confirm this hypothesis. Comparison of NTPAH3POPT-1 (single molecule devoid of H-bonds) and NTPAH3MOPT-1 (single molecule with intramolecular H-bond system) already shows that P=O and P–O bonds are elongated and shortened, respectively, upon H-bond formation, and elongation of the O–H bond is observed as well. However, the additional interactions in the crystal packing cause further changes, which is evident from the comparison of NTPAH3MOPT-1 and NTPAH3MOPT-CRYST. Longer P=O and P–O bonds (with respect to the related bonds in the isolated molecule) are observed, and features of stronger H-bond systems are revealed. In NTPAH3MOPT-CRYST, a longer O–H bond, shorter H⋅⋅⋅O distance and an overall shorter O⋅⋅⋅O distance is observed. Comparison of the two optimized molecular structures in the crystal environment (NTPAH3POPT-CRYST and NTPAH3MOPT-CRYST) reveals that these modifications exhibit very similar H-bond features. Whereas in NTPAH3M the slightly longer O–H bond may indicate a stronger H-bond system, the slightly longer O⋅⋅⋅O distance is indicative of the opposite. Hence, the structural differences are less pronounced than expected from the crystallographic analysis. Therefore, further support will be delivered by solid-state 1H NMR and IR spectroscopy (cf. Section 3.4). Last but not least, the different C-N-C bond angles in NTPAH3P and NTPAH3M were also confirmed by the computational analyses. Hence, the wide C-N-C angles are a characteristic feature of the cage structure of the molecule in NTPAH3M.

3.4. IR and Solid-State NMR Spectroscopic Analyses of NTPAH3P and NTPAH3M

The IR spectra of NTPAH3P and NTPAH3M (Figure 8) are characteristic fingerprints, which clearly distinguish between the two modifications. The most prominent difference is associated with the OH groups of these compounds. As reported earlier [47], the vibrational modes of phosphinic acids´ OH groups in cyclic hydrogen bond situations may give rise to an “A-B-C”-pattern of broad bands of high intensity. These bands are located around 2560, 2200 and 1680 cm−1 for NTPAH3P and around 2630, 2190 and 1540 cm−1 for NTPAH3M. Notably, the “C” band is shifted to lower wavenumbers and gains in relative intensity over “A” and “B” bands for NTPAH3M. In the literature [47] it is interpreted as a sign of the increasing strength of the set of hydrogen bonds involved. Furthermore, a splitting of the “B” band is particularly obvious for NTPAH3P, and in the literature that splitting has been mentioned in the context of weaker H-bonds.
1H MAS NMR spectroscopy of these solids also supports the presence of stronger H-bonds in NTPAH3M. As shown in Figure 9, the signal of the OH protons (at δiso 13.6 ppm for NTPAH3P and at δiso 16.6 ppm for NTPAH3M) is sufficiently well resolved from the shift range of the aryl protons (the broad signal in the range 4–10 ppm) and thus clearly shows that the OH signal of NTPAH3M is shifted by Δδ ca. +3 ppm. It is well known that 1H NMR signals of protons in H-bond situations are shifted to higher frequencies and that this shift is more pronounced for nuclei in stronger H-bonds [48,49]. This has been underlined by correlations of decreasing OH groups’ 1H NMR shielding with shortening in O⋅⋅⋅O distances in series of carboxylic acids and phosphonic acids [48]. Moreover, this shift difference of the OH signal represents the greatest difference in the 1H MAS NMR spectra of these two compounds.
The 13C MAS NMR (Figure 10) and 31P MAS NMR spectra (Figure 11) reveal further spectroscopic differences between NTPAH3P and NTPAH3M. In the former, the notable shift difference Δδ of ca. 10 ppm between the signals of the CH2 carbon atoms (at δiso 54 ppm for NTPAH3P and at δiso 64 ppm for NTPAH3M) indicates the involvement of these groups in different bonding situations in the two modifications (i.e., as part of a cage motif in NTPAH3M vs. absence of this motif in NTPAH3P). In contrast, the shift differences of the aromatic 13C signals (in the range 125–135 ppm) are less pronounced. The 31P MAS NMR signals (at δiso 34.2 ppm for NTPAH3P and at δiso 37.7 ppm for NTPAH3M) also allow us to assign the modification to the respective solid. Both fingerprint features of the two modifications, the 13C signal of the CH2 groups and the 31P signals, are sufficiently well separated to rule out the presence of significant amounts of the other modification in the respective sample.

4. Conclusions

In general, nitrilotris(methylenephenylphosphinic) acid (NTPAH3) can be silylated using hexamethyldisilazane to afford the tris(trimethylsilyl) derivative, NTPA(SiMe3)3. The latter is sensitive toward protolysis, and alcoholysis in chloroform liberates the acid NTPAH3. The different solvents used for synthesis of NTPAH3 (aqueous hydrochloric acid vs. chloroform) give rise to the crystallization of two different modifications of this acid, i.e., NTPAH3P and NTPAH3M, respectively. As both modifications exhibit similar stability (shown by computational analysis), the formation of these modifications can be attributed to the different crystallization kinetics in these different solvents, which may depend on the predominant conformer in solution. Whereas aqueous environment may foster H-bond interactions of NTPAH3 with the solvent, thus supporting the conformation encountered in NTPAH3P crystallized therefrom, chloroform cannot be involved into intermolecular H-bonds to similar extent, thus fostering the formation of intramolecular H-bonds in NTPAH3 and supporting the conformation encountered in NTPAH3M crystallized therefrom. The absence of zwitterionic species in the crystals of NTPAH3 can be attributed to the special feature of cyclic (O=P–O–H)3 H-bond trimers (i.e., R ( 12 ) 3 3 motifs), which engage all H-bond donor and acceptor sites of the phosphinic acid moieties, thus leaving no acid for protonation of the trialkylamine N atom. Therefore, we conclude that the 3-fold symmetry of the NTPAH3 molecule, which supports crystallization in trigonal space groups and thus also supports formation of (O=P–O–H)3 H-bond trimers (driven by crystal symmetry) represents a key to this rather unexpected appearance of a non-zwitterionic aminomethylphosphinic acid in its solid-state modifications. The crystal structure of NTPA(SiMe3)3, which is devoid of related H-bond systems, underlines this trend of crystallization of the NTPA motif in a trigonal crystal symmetry.
Molecular networks with a high density of hydrogen bonds (e.g., in layered silicon hydrogen phosphates [14] and various other H-bonded networks [50,51,52,53,54,55]) may exhibit good proton conductivity. Thus, further modification of the NTPA motif toward a covalently or coordinatively bound network with retention of O–H⋅⋅⋅O H-bond motifs and the additional trialkylamine N atoms may give rise to H+-conductive materials. Further investigations in this direction are currently under way.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cryst14070662/s1: NMR spectra (1H, 13C{1H}, 31P{1H} and, where applicable, 29Si{1H}) of solutions of compounds NTPAH3 (Figures S1–S3) and NTPA(SiMe3)3 (Figures S4–S7, Table S1), of solid NTPAH3P and NTPAH3M (Figures S8–S13), Hirshfeld surface analysis dnorm isosurface views of the asymmetric units of NTPAH3P and NTPAH3M (Figures S14 and S15), comparison of the Hirshfeld surface analysis fingerprint maps for O-H-contacts of the asymmetric units of NTPAH3P and NTPAH3M (Figure S16), comparison of the dependence of percentage of interatomic contacts on the refinement model: free vs. idealized H-refinement of OH groups (Figure S17), atomic coordinates, total energies and graphical representations of the optimized molecular structures of NTPAH3P and NTPAH3M as isolated molecules (Figures S18 and S19, Tables S2 and S3) and in the periodic crystal environment (Figures S20 and S21, Tables S4 and S5).

Author Contributions

Conceptualization, J.W.; funding acquisition, E.B., R.G. and E.K.; investigation, S.K., E.B., R.G. and J.W.; writing—original draft preparation, J.W.; writing—review and editing, S.K., E.B., R.G., E.K. and J.W.; visualization, S.K. and J.W.; supervision, E.B., E.K. and J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the German Science Foundation (DFG) under Project 509900844 (Neue Si-O-P-N Verbindungen und Materialien: Synthesen, Strukturen und Materialeigenschaften) and by the German Federal Ministry of Environment, Nature Conservation, Nuclear Safety and Consumer Protection (BMUV) under Project 1501667 (Am-BALL).

Data Availability Statement

CCDC 2365392 (NTPAH3P), 2365391 (NTPAH3M) and 2365393 (NTPA(SiMe3)3) contain the supplementary crystal data for this article. These data can be obtained free of charge from the Cambridge Crystallographic Data Centre via https://www.ccdc.cam.ac.uk/structures/ (accessed on 25 June 2024).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Crystallographic data from data collection and refinement for the nitrilotriphosphinic acid N(CH2P(Ph,O,OH))3 (NTPAH3) in its modifications NTPAH3P and NTPAH3M, and for its tris(trimethylsilyl) ester N(CH2P(Ph,O,OSiMe3))3 NTPA(SiMe3)3.
Table A1. Crystallographic data from data collection and refinement for the nitrilotriphosphinic acid N(CH2P(Ph,O,OH))3 (NTPAH3) in its modifications NTPAH3P and NTPAH3M, and for its tris(trimethylsilyl) ester N(CH2P(Ph,O,OSiMe3))3 NTPA(SiMe3)3.
ParameterNTPAH3P 1NTPAH3MNTPA(SiMe3)3 2
FormulaC21H24NO6P3C21H24NO6P3C30H48NO6P3Si3
Mr479.32479.32695.87
T(K)180(2)180(2)160(2)
λ(Å)0.710730.710730.71073
Crystal systemtrigonaltrigonaltrigonal
Space groupP3c1R3c R 3 ¯
a = b(Å)9.3792(3)18.2554(5)12.1818(3)
c(Å)15.6549(6)11.4177(3)92.431(3)
V3)1192.65(9)3295.3(2)11,878.8(7)
Z2612
ρcalc(g·cm−1)1.341.451.17
μMoKα (mm−1)0.30.30.3
F(000)50015004440
θmax(°), Rint28.0, 0.066228.0, 0.035025.0, 0.0492
Completeness100%100%99.9%
Reflns collected14,28118,163255,15
Reflns unique193217484686
Restraints110
Parameters9998265
GoF1.0361.0551.106
χFlack0.05(6)−0.01(3)n/a
R1, wR2 [I > 2σ(I)]0.0317, 0.07110.0224, 0.06020.0376, 0.1017
R1, wR2 (all data)0.0425, 0.07420.0243, 0.06100.0550, 0.1066
Largest peak/hole (e·Å−3)0.21, −0.210.21, −0.150.30, −0.24
1 The batch of crystals contained long trigonal prisms as well as pseudo-hexagonal plates. The data set was collected from a prism. The plates represent the same modification, but they are twins according to the twin law (1 1 0, 0–1 0, 0 0–1). 2 The batch of crystals consisted of both thick plates and truncated bipyramids, both of which show trigonal symmetry. The data set was collected from a small cap of a truncated bipyramidal crystal. The central parts of the truncated bipyramids exhibited layering, which indicated stacking faults. Data collection from the larger trigonal plates revealed essentially the same unit cell parameters, but superstructure reflections along l indicated a primitive c axis of the same length. Those data sets allowed for structure solution in both space group types R 3 ¯ and P3, but both models refined to unsatisfactory R-values (above 0.20) and resulted in noticeable residual electron density peaks. We therefore interpret these crystals as two-phase twins, which result from large contributions of stacking faults along c and give rise to different local space group symmetries (i.e., space group types R 3 ¯ and P3). We are aware that these kinds of twins can be refined, for example it has been done for a 1:1 twin of space group types P212121 and Pmnb [56], but such an analysis was beyond the scope of our study as crystal picking and collection of a data set in space group R 3 ¯ already provided satisfactory information of the molecular structure of compound NTPA(SiMe3)3.

References

  1. Stanford, R.H., Jr. The Crystal Structure of Nitrilotriacetic Acid. Acta Crystallogr. 1967, 23, 825–832. [Google Scholar] [CrossRef]
  2. Skrzypczak-Jankun, E.; Smith, D.A. Nitrilotriacetic Acid, C6H9NO6. Acta Crystallogr. C 1994, 50, 1097–1099. [Google Scholar] [CrossRef]
  3. Hoppe, B.; Martens, J. Aminosäuren—Bausteine des Lebens. Chem. Unserer Zeit 1983, 17, 4153. [Google Scholar] [CrossRef]
  4. Urbanovský, P.; Kotek, J.; C’ísařová, I.; Hermann, P. Selective and clean synthesis of aminoalkyl-H-phosphinic acids from hypophosphorous acid by phospha-Mannich reaction. RSC Adv. 2020, 10, 21329–21349. [Google Scholar] [CrossRef] [PubMed]
  5. Kasser, J.; Nazarov, A.A.; Hartinger, C.G.; Wdziekonski, B.; Dani, C.; Kuznetsov, M.L.; Arion, V.B.; Keppler, B.K. A one step/one pot synthesis of N,N-bis(phosphonomethyl)amino acids and their effects on adipogenic and osteogenic differentiation of human mesenchymal stem cells. Bioorg. Med. Chem. 2009, 17, 3388–3393. [Google Scholar] [CrossRef] [PubMed]
  6. Baldwin, R.A.; Cheng, M.T. Organometalloid azides. IV. Preparation and reactions of N,N′-[p-arylenebis(diphenylphosphoranylidyne)]bis(P-phenyl phosphonamidic azides. J. Org. Chem. 1967, 32, 2636–2639. [Google Scholar] [CrossRef]
  7. Quin, L.D.; Dysart, M.R. Arylphosphinic Acids: Dissociation Constants and Reaction with Diazomethane. J. Org. Chem. 1962, 27, 1012–1014. [Google Scholar] [CrossRef]
  8. Daly, J.J.; Wheatley, P.J. The Crystal and Molecular Structure of Nitrilotrimethylene Triphosphonic Acid. J. Chem. Soc. A 1967, 1967, 212–221. [Google Scholar] [CrossRef]
  9. Moschona, A.; Plesu, N.; Mezei, G.; Thomas, A.G.; Demadis, K.D. Corrosion protection of carbon steel by tetraphosphonates of systematically different molecular size. Corros. Sci. 2018, 145, 135–150. [Google Scholar] [CrossRef]
  10. Bailly, T.; Burgada, R.; Lecouvey, M.; Neuman, A.; Prangé, T. Trans 1,2 diaminocyclohexane as a template in the synthesis of ligands for transition metal and actinide in vivo detoxification. ARKIVOC 2003, 9, 140–149. [Google Scholar] [CrossRef]
  11. Cecconi, F.; Ghilardi, C.A.; Gili, P.; Midollini, S.; Lorenzo Luis, P.A.; Lozano-Gorrìn, A.D.; Orlandini, A. Complexation of nickel(II) and lead(II) cations with the tripodal nitrilo-tris(methylenephenylphosphinic) acid (H3L). X-ray crystal structure of the dimer [Ni(HL)(DMSO)]2·2DMSO. Inorg. Chim. Acta 2001, 319, 67–74. [Google Scholar] [CrossRef]
  12. Cecconi, F.; Ghilardi, C.A.; Lorenzo Luis, P.A.; Midollini, S.; Orlandini, A.; Dakternieks, D.; Duthie, A.; Dominguez, S.; Berti, E.; Vacca, A. Complexes of the tripodal nitrilotrimethylenetrisphosphonic (H6L) and P,P′,P″-triphenylnitrilotrimethylenetrisphosphinic (H3L°) acids with the copper(II) ion. Synthesis and characterization of [Hpy][Cu(H3L)(H2O)] and [Cu(HL°)(py)]2⋅2Me2CO. J. Chem. Soc. Dalton Trans. 2001, 30, 211–217. [Google Scholar] [CrossRef]
  13. Jähnigen, S.; Brendler, E.; Böhme, U.; Kroke, E. Synthesis of silicophosphates containing SiO6-octahedra under ambient conditions—Reactions of anhydrous H3PO4 with alkoxysilanes. Chem. Commun. 2012, 48, 7675–7677. [Google Scholar] [CrossRef]
  14. Kowalke, J.; Arnold, C.; Ponomarev, I.; Jäger, C.; Kroll, P.; Brendler, E.; Kroke, E. Structural Insight into Layered Silicon Hydrogen Phosphates Containing [SiO6] Octahedra Prepared by Different Reaction Routes. Eur. J. Inorg. Chem. 2019, 2019, 828–836. [Google Scholar] [CrossRef]
  15. Viehweger, C.; Kowalke, J.; Brendler, E.; Schwarzer, S.; Vogt, C.; Kroke, E. Five- and six-fold coordinated silicon in silicodiphosphonates: Short range order investigation by solid-state NMR spectroscopy. New J. Chem. 2020, 44, 4613–4620. [Google Scholar] [CrossRef]
  16. Kowalke, J.; Wagler, J.; Viehweger, C.; Brendler, E.; Kroke, E. Ionic Dissociation of SiCl4: Formation of [SiL6]Cl4 with L = Dimethylphosphinic Acid. Chem. Eur. J. 2020, 26, 8003–8006. [Google Scholar] [CrossRef] [PubMed]
  17. Mohan, B.; Singh, G.; Gupta, R.K.; Sharma, P.K.; Solovev, A.A.; Pombeiro, A.J.L.; Ren, P. Hydrogen-bonded organic frameworks (HOFs): Multifunctional material on analytical monitoring. Trends Anal. Chem. 2024, 170, 117436. [Google Scholar] [CrossRef]
  18. Lin, R.-B.; Chen, B. Hydrogen-bonded organic frameworks: Chemistry and functions. Chem 2022, 8, 2114–2135. [Google Scholar] [CrossRef]
  19. Chen, L.; Zhang, B.; Chen, L.; Liu, H.; Hu, Y.; Qiao, S. Hydrogen-bonded organic frameworks: Design, applications, and prospects. Mater. Adv. 2022, 3, 3680–3708. [Google Scholar] [CrossRef]
  20. Samuel, H.S.; Nweke-Maraizu, U.; Etim, E.E. Understanding Intermolecular and Intramolecular Hydrogen Bonds: Spectroscopic and Computational Approaches. J. Chem. Rev. 2023, 5, 439–465. [Google Scholar] [CrossRef]
  21. Song, P.; Wang, H. High-Performance Polymeric Materials through Hydrogen-Bond Cross-Linking. Adv. Mater. 2020, 32, 1901244. [Google Scholar] [CrossRef] [PubMed]
  22. Sheldrick, G.M. SHELXT—Integrated space-group and crystal-structure determination. Acta Crystallogr. A 2015, 71, 3–8. [Google Scholar] [CrossRef] [PubMed]
  23. Sheldrick, G.M. Program for the Refinement of Crystal Structures; SHELXL-2019/3; University of Göttingen: Göttingen, Germany, 2019. [Google Scholar]
  24. Sheldrick, G.M. A short history of SHELX. Acta Crystallogr. A 2008, 64, 112–122. [Google Scholar] [CrossRef] [PubMed]
  25. Farrugia, L.J. ORTEP-3 for windows—A version of ORTEP-III with a graphical user interface (GUI). J. Appl. Crystallogr. 1997, 30, 565. [Google Scholar] [CrossRef]
  26. Farrugia, L.J. WinGX and ORTEP for Windows: An update. J. Appl. Crystallogr. 2012, 45, 849–854. [Google Scholar] [CrossRef]
  27. POV-RAY (Version 3.7), Trademark of Persistence of Vision Raytracer Pty. Ltd., Williamstown, Victoria (Australia). Copyright Hallam Oaks Pty. Ltd., 1994–2004. Available online: http://www.povray.org/download/ (accessed on 28 June 2021).
  28. Neese, F. Software update: The ORCA program system—Version 5.0. WIREs Comput. Mol. Sci. 2022, 12, e1606. [Google Scholar] [CrossRef]
  29. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef] [PubMed]
  30. Pantazis, D.A.; Neese, F. All-electron basis sets for heavy elements. WIREs Comput. Mol. Sci. 2014, 4, 363–374. [Google Scholar] [CrossRef]
  31. van Lenthe, E.; Baerends, E.J.; Snijders, J.G. Relativistic regular two-component Hamiltonians. J. Chem. Phys. 1993, 99, 4597–4610. [Google Scholar] [CrossRef]
  32. van Wüllen, C. Molecular density functional calculations in the regular relativistic approximation: Method, application to coinage metal diatomics, hydrides, fluorides and chlorides, and comparison with first-order relativistic calculations. J. Chem. Phys. 1998, 109, 392–399. [Google Scholar] [CrossRef]
  33. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456. [Google Scholar] [CrossRef] [PubMed]
  34. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104–154119. [Google Scholar] [CrossRef] [PubMed]
  35. Balasubramani, S.G.; Chen, G.P.; Coriani, S.; Diedenhofen, M.; Frank, M.S.; Franzke, Y.J.; Furche, F.; Grotjahn, R.; Harding, M.E.; Hättig, C.; et al. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J. Chem. Phys. 2020, 152, 184107. [Google Scholar] [CrossRef] [PubMed]
  36. Lu, T.; Chen, F. Multiwfn: A Multifunctional Wavefunction Analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  37. Chemcraft, Version 1.8 (Build 164). 2016. Available online: https://www.chemcraftprog.com/ (accessed on 19 September 2015).
  38. Plaza, A.I.; Grim, S.O.; Motekaitis, R.J.; Martell, A.E. ETHYLENEBIS(NITRILODIMETHYLENE)TETRAKIS(PHENYLPHOSPHINIC ACID). In Inorganic Syntheses; Basolo, F., Ed.; McGraw-Hill, Inc.: New York, NY, USA, 1976; Volume 16, pp. 199–202. [Google Scholar] [CrossRef]
  39. Nickelsen, J.; van Gerven, D.; Wickleder, M.S. In Situ single crystal growth of tris(trimethylsilyl)phosphate, (Me3SiO)3PO. Z. Anorg. Allg. Chem. 2022, 648, e202200184. [Google Scholar] [CrossRef]
  40. Ishida, S.; Hirakawa, F.; Iwamoto, T. Reactions of a Stable Phosphinyl Radical with Stable Aminoxyl Radicals. Chem. Lett. 2015, 44, 94–96. [Google Scholar] [CrossRef]
  41. This Refers to a Search in the Cambridge Structure Database Using ConQuest Version 2023.3.0 (Build 392113). Available online: https://www.ccdc.cam.ac.uk/solutions/software/conquest (accessed on 10 June 2024).
  42. Bernstein, J.; Davis, R.E.; Shimoni, L.; Chang, N.-L. Patterns in Hydrogen Bonding: Functionality and Graph Set Analysis in Crystals. Angew. Chem. Int. Ed. 1995, 34, 1555–1573. [Google Scholar] [CrossRef]
  43. Spackman, P.R.; Turner, M.J.; McKinnon, J.J.; Wolff, S.K.; Grimwood, D.J.; Jayatilaka, D.; Spackman, M.A. CrystalExplorer: A program for Hirshfeld surface analysis, visualization and quantitative analysis of molecular crystals. J. Appl. Cryst. 2021, 54, 1006–1011. [Google Scholar] [CrossRef]
  44. Asfin, R.E.; Denisov, G.S.; Poplevchenkov, D.N.; Tokhadze, K.G.; Velikanova, T.V. IR υ(OH) Band and Dimerization of Phosphorus Acids in the Gas Phase and Solid State. Pol. J. Chem. 2002, 76, 1223–1231. [Google Scholar]
  45. Asfin, R.E.; Denisov, G.S.; Tokhadze, K.G. The infrared spectra and enthalpies of strongly bound dimers of phosphinic acids in the gas phase. (CH2Cl)2POOH and (C6H5)2POOH. J. Mol. Struct. 2002, 608, 161–168. [Google Scholar] [CrossRef]
  46. Emamian, S.; Lu, T.; Kruse, H.; Emamian, H. Exploring Nature and Predicting Strength of Hydrogen Bonds: A Correlation Analysis Between Atoms-in-Molecules Descriptors, Binding Energies, and Energy Components of Symmetry-Adapted Perturbation Theory. J. Comput. Chem. 2019, 40, 2868–2881. [Google Scholar] [CrossRef] [PubMed]
  47. Hadzi, D. Infrared spectra of strongly hydrogen-bonded systems. Pure Appl. Chem. 1965, 11, 435–454. [Google Scholar] [CrossRef]
  48. Harris, R.K.; Jackson, P.; Merwin, L.H.; Say, B.J.; Hägele, G. Perspectives in high-resolution solid-state nuclear magnetic resonance, with emphasis on combined rotation and multiple-pulse spectroscopy. J. Chem. Soc. Faraday Trans. 1 1988, 84, 3649–3672. [Google Scholar] [CrossRef]
  49. Berglund, B.; Vaughan, R.W. Correlations between proton chemical shift tensors, deuterium quadrupole couplings, and bond distances for hydrogen bonds in solids. J. Chem. Phys. 1980, 73, 2037–2043. [Google Scholar] [CrossRef]
  50. Chand, S.; Elahi, S.M.; Pal, A.; Das, M.C. Metal–Organic Frameworks and Other Crystalline Materials for Ultrahigh Superprotonic Conductivities of 10–2 s cm−1 or Higher. Chem. Eur. J. 2019, 25, 6259–6269. [Google Scholar] [CrossRef] [PubMed]
  51. Hao, B.-B.; Wang, X.-X.; Zhang, C.-X.; Wang, Q. Two Hydrogen-Bonded Organic Frameworks with Imidazole Encapsulation: Synthesis and Proton Conductivity. Cryst. Growth. Des. 2021, 21, 3908–3915. [Google Scholar] [CrossRef]
  52. Li, S.; Liu, Y.; Li, L.; Liu, C.; Li, J.; Ashraf, S.; Li, P.; Wang, B. Enhanced Proton Conductivity of Imidazole-Doped Thiophene-Based Covalent Organic Frameworks via Subtle Hydrogen Bonding Modulation. ACS Appl. Mater. Interfaces 2020, 12, 22910–22916. [Google Scholar] [CrossRef] [PubMed]
  53. Staiger, A.; Paren, B.A.; Zunker, R.; Hoang, S.; Häußler, M.; Winey, K.I.; Mecking, S. Anhydrous Proton Transport within Phosphonic Acid Layers in Monodisperse Telechelic Polyethylenes. J. Am. Chem. Soc. 2021, 143, 16725–16733. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, K.; Li, T.; Zeng, H.; Zou, G.; Zhang, Q.; Lin, Z. Ionothermal Synthesis of Open-Framework Metal Phosphates Using a Multifunctional Ionic Liquid. Inorg. Chem. 2018, 57, 8726–8729. [Google Scholar] [CrossRef]
  55. Zhang, K.-M.; Lou, Y.-L.; He, F.-Y.; Duanm, H.-B.; Huang, X.-Q.; Fan, Y.; Zhao, H.-R. The water-mediated proton conductivity of a 1D open framework inorganic-organic hybrid iron phosphate and its composite membranes. Inorg. Chem. Commun. 2021, 134, 109032. [Google Scholar] [CrossRef]
  56. Caldwell, L.M.; Hill, A.F.; Rae, A.D.; Willis, A.C. Alkynylselenolatoalkylidynes: [Mo(≡CSeC≡CR)(CO)2{HB(pzMe2)3}] (R = CMe3, SiMe3; pzMe2 = 3,5-Dimethylpyrazol-1-yl). Organometallics 2008, 27, 341–345. [Google Scholar] [CrossRef]
Figure 1. Selected multi-functional α-amino acids of the carboxylic, phosphinic and phosphonic acid types (IVI) and the nitrilotris(methylenephenylphosphinic) acid (NTPAH3) studied herein. The current study will reveal the location of the acidic protons in the latter compound.
Figure 1. Selected multi-functional α-amino acids of the carboxylic, phosphinic and phosphonic acid types (IVI) and the nitrilotris(methylenephenylphosphinic) acid (NTPAH3) studied herein. The current study will reveal the location of the acidic protons in the latter compound.
Crystals 14 00662 g001
Scheme 1. Syntheses of compounds NTPAH3 (in its modifications NTPAH3P and NTPAH3M) and NTPA(SiMe3)3.
Scheme 1. Syntheses of compounds NTPAH3 (in its modifications NTPAH3P and NTPAH3M) and NTPA(SiMe3)3.
Crystals 14 00662 sch001
Figure 2. Molecular structure of NTPA(SiMe3)3 with thermal displacement ellipsoids at the 50% probability level (for clarity, hydrogen atoms are not depicted). Sections (a,c) show a selected view of molecule 1 and 2, respectively, with the atoms of the asymmetric unit labelled, (b,d) show the molecules viewed along the crystallographic c axis. The two independent nitrogen atoms (N1, N2) are located on crystallographic 3-fold rotation axes; the whole molecules are generated by the symmetry operations (−x+y, −x+1, z) and (–y+1, x–y+1, z) for molecule 1, (–y+2, x–y+1, z) and (−x+y+1, −x+2, z) for molecule 2.
Figure 2. Molecular structure of NTPA(SiMe3)3 with thermal displacement ellipsoids at the 50% probability level (for clarity, hydrogen atoms are not depicted). Sections (a,c) show a selected view of molecule 1 and 2, respectively, with the atoms of the asymmetric unit labelled, (b,d) show the molecules viewed along the crystallographic c axis. The two independent nitrogen atoms (N1, N2) are located on crystallographic 3-fold rotation axes; the whole molecules are generated by the symmetry operations (−x+y, −x+1, z) and (–y+1, x–y+1, z) for molecule 1, (–y+2, x–y+1, z) and (−x+y+1, −x+2, z) for molecule 2.
Crystals 14 00662 g002
Figure 3. Thermal ellipsoid plots (50% probability) of the asymmetric units of NTPAH3 in its modifications NTPAH3P (a) and NTPAH3M (b), view along the P–CH2 bond.
Figure 3. Thermal ellipsoid plots (50% probability) of the asymmetric units of NTPAH3 in its modifications NTPAH3P (a) and NTPAH3M (b), view along the P–CH2 bond.
Crystals 14 00662 g003
Figure 4. Molecular structures of NTPAH3 in its modifications NTPAH3P (in sections (a,c,e,f)) and NTPAH3M (in (b,d)) with thermal displacement ellipsoids at the 50% probability level (for clarity, C-bound hydrogen atoms are not depicted). Sections (a,b) show a selected view of each molecule with the atoms of the asymmetric unit labelled, (c,d) show the molecules viewed along the crystallographic c axis. Their nitrogen atoms (N1) are located on crystallographic 3-fold rotation axes, in both cases the whole molecules are generated by the symmetry operations (–y+1, x–y+1, z) and (−x+y, −x+1, z). For NTPAH3M these symmetry operations create a trimeric intramolecular hydrogen bond system (shown in (d), labels with * and ** indicate symmetry equivalents). In NTPAH3P intermolecular hydrogen bonds in a polymeric network (shown in (f)) give rise to a related (O=P–O–H)3 motif (magnified section of interest from (f) shown in (e), these symmetry equivalents * and ** are generated by symmetry operations –y+1, x–y, z and −x+y+1, −x+1, z, respectively).
Figure 4. Molecular structures of NTPAH3 in its modifications NTPAH3P (in sections (a,c,e,f)) and NTPAH3M (in (b,d)) with thermal displacement ellipsoids at the 50% probability level (for clarity, C-bound hydrogen atoms are not depicted). Sections (a,b) show a selected view of each molecule with the atoms of the asymmetric unit labelled, (c,d) show the molecules viewed along the crystallographic c axis. Their nitrogen atoms (N1) are located on crystallographic 3-fold rotation axes, in both cases the whole molecules are generated by the symmetry operations (–y+1, x–y+1, z) and (−x+y, −x+1, z). For NTPAH3M these symmetry operations create a trimeric intramolecular hydrogen bond system (shown in (d), labels with * and ** indicate symmetry equivalents). In NTPAH3P intermolecular hydrogen bonds in a polymeric network (shown in (f)) give rise to a related (O=P–O–H)3 motif (magnified section of interest from (f) shown in (e), these symmetry equivalents * and ** are generated by symmetry operations –y+1, x–y, z and −x+y+1, −x+1, z, respectively).
Crystals 14 00662 g004
Figure 5. Fingerprint plots (de vs. di) of the Hirshfeld surface analyses (performed with CrystalExplorer version 21.5, revision 608bb32 [43]) of (a) the asymmetric unit of the crystal structure of NTPAH3P and (b) the asymmetric unit of the crystal structure of NTPAH3M. Note: The N–C interactions correspond to the N–C bonds, which connect the asymmetric unit with the two further symmetry equivalent parts of the molecule.
Figure 5. Fingerprint plots (de vs. di) of the Hirshfeld surface analyses (performed with CrystalExplorer version 21.5, revision 608bb32 [43]) of (a) the asymmetric unit of the crystal structure of NTPAH3P and (b) the asymmetric unit of the crystal structure of NTPAH3M. Note: The N–C interactions correspond to the N–C bonds, which connect the asymmetric unit with the two further symmetry equivalent parts of the molecule.
Crystals 14 00662 g005
Figure 6. Selected section of the Hirshfeld surface analysis (performed with CrystalExplorer version 21.5, revision 608bb32 [43]) of the asymmetric unit of the crystal structure of NTPAH3P showing (a) the area of C–C interactions on the isosurface (dnorm) and (b) the corresponding section of the fingerprint plot (de vs. di).
Figure 6. Selected section of the Hirshfeld surface analysis (performed with CrystalExplorer version 21.5, revision 608bb32 [43]) of the asymmetric unit of the crystal structure of NTPAH3P showing (a) the area of C–C interactions on the isosurface (dnorm) and (b) the corresponding section of the fingerprint plot (de vs. di).
Crystals 14 00662 g006
Figure 7. Percentage distribution (with respect to the area on the Hirshfeld surface) of the contacts of the asymmetric units of the crystal structures of NTPAH3P and NTPAH3M. Note: The N–C and N–H interactions (shown in the black box) correspond to the N–C bonds and interactions of N with their CH2 H atoms. Thus, they are intramolecular interactions in both cases.
Figure 7. Percentage distribution (with respect to the area on the Hirshfeld surface) of the contacts of the asymmetric units of the crystal structures of NTPAH3P and NTPAH3M. Note: The N–C and N–H interactions (shown in the black box) correspond to the N–C bonds and interactions of N with their CH2 H atoms. Thus, they are intramolecular interactions in both cases.
Crystals 14 00662 g007
Figure 8. IR spectra of (a) NTPAH3P and (b) NTPAH3M (in KBr pellets, recorded at room temperature).
Figure 8. IR spectra of (a) NTPAH3P and (b) NTPAH3M (in KBr pellets, recorded at room temperature).
Crystals 14 00662 g008
Figure 9. 1H MAS NMR spectra of (a) NTPAH3P and (b) NTPAH3M recorded at MAS frequencies of 30 kHz. (The probe background has not been subtracted).
Figure 9. 1H MAS NMR spectra of (a) NTPAH3P and (b) NTPAH3M recorded at MAS frequencies of 30 kHz. (The probe background has not been subtracted).
Crystals 14 00662 g009
Figure 10. 13C CP MAS NMR spectra of (a) NTPAH3P and (b) NTPAH3M recorded at MAS frequencies of 10 kHz. The asterisks (*) indicate spinning sidebands of the aryl-C signals.
Figure 10. 13C CP MAS NMR spectra of (a) NTPAH3P and (b) NTPAH3M recorded at MAS frequencies of 10 kHz. The asterisks (*) indicate spinning sidebands of the aryl-C signals.
Crystals 14 00662 g010
Figure 11. 31P MAS NMR spectra of (a) NTPAH3P and (b) NTPAH3M recorded at MAS frequencies of 10 kHz. The asterisks (*) indicate spinning sidebands.
Figure 11. 31P MAS NMR spectra of (a) NTPAH3P and (b) NTPAH3M recorded at MAS frequencies of 10 kHz. The asterisks (*) indicate spinning sidebands.
Crystals 14 00662 g011
Table 1. Selected corresponding bond lengths (Å) and angles (deg) in the crystal structures of NTPA(SiMe3)3, NTPAH3P and NTPAH3M.
Table 1. Selected corresponding bond lengths (Å) and angles (deg) in the crystal structures of NTPA(SiMe3)3, NTPAH3P and NTPAH3M.
NTPA(SiMe3)3 1NTPAH3PNTPAH3M
N–C1.473(2)/1.469(2)1.469(3)1.467(2)
P–C(CH2)1.801(2)/1.802(2)1.806(3)1.817(2)
P–C(Ph)1.788(2)/1.797(2)1.796(3)1.786(2)
P–O1.575(2)/1.571(2)1.547(2)1.533(2)
P=O1.469(2)/1.471(2)1.494(2)1.510(2)
C-N-C110.37(12)/111.11(12)109.56(17)116.04(9)
N-C-P113.21(14)/113.08(14)115.74(18)109.29(15)
C-P-C107.12(10)/107.92(9)109.81(13)109.32(9)
O-P-O115.57(9)/116.07(9)115.79(15)114.40(11)
C(CH2)-P-O(H/Si)100.15(8)/99.73(8)98.48(13)107.29(1)
C(CH2)-P=O114.93(9)/114.58(9)113.47(12)109.62(10)
C(Ph)-P-O(H/Si)104.79(10)/104.74(9)109.33(13)106.67(10)
C(Ph)-P=O112.98(11)/112.57(10)109.46(14)109.45(9)
ΣX-P=O 2343.5/343.2338.7333.5
1 Two individual sets of parameters are given in accordance with the asymmetric unit. 2 Sum of angles about P=O (X = O, C).
Table 2. Selected interatomic distances (Å) and angles (deg) in NTPA(SiMe3)3 (from XRD) as well as in NTPAH3P and NTPAH3M (data from XRD and optimized molecular structures).
Table 2. Selected interatomic distances (Å) and angles (deg) in NTPA(SiMe3)3 (from XRD) as well as in NTPAH3P and NTPAH3M (data from XRD and optimized molecular structures).
CompoundP=OP–OO–HH⋅⋅⋅OO⋅⋅⋅OC-N-C
NTPA(SiMe3)3 (XRD)1.469(2)1.574(2)---110.4(1)
1.471(2)1.571(2)---111.1(1)
NTPAH3P (XRD)1.494(2)1.547(2)- 1- 12.487(3)109.6(2)
NTPAH3M (XRD)1.510(2)1.533(2)- 1- 12.450(2)116.0(1)
NTPAH3POPT-1 21.4761.6080.96--112.5
NTPAH3MOPT-1 21.4961.5531.031.472.495116.8
NTPAH3POPT-CRYST 21.5591.6071.061.412.469109.3
NTPAH3MOPT-CRYST 21.5581.6101.071.412.480116.9
1 The O–H bond lengths from XRD analyses are not taken into consideration for comparison because of the inherent problem of uncertainty of determination of H-atom positions by X-ray diffraction. 2 For molecules or structures with multiple independent bonds or atom distances of the same kind, the average value is reported here.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Knerr, S.; Brendler, E.; Gericke, R.; Kroke, E.; Wagler, J. Two Modifications of Nitrilotris(methylenephenylphosphinic) Acid: A Polymeric Network with Intermolecular (O=P–O–H)3 vs. Monomeric Molecules with Intramolecular (O=P–O–H)3 Hydrogen Bond Cyclotrimers. Crystals 2024, 14, 662. https://doi.org/10.3390/cryst14070662

AMA Style

Knerr S, Brendler E, Gericke R, Kroke E, Wagler J. Two Modifications of Nitrilotris(methylenephenylphosphinic) Acid: A Polymeric Network with Intermolecular (O=P–O–H)3 vs. Monomeric Molecules with Intramolecular (O=P–O–H)3 Hydrogen Bond Cyclotrimers. Crystals. 2024; 14(7):662. https://doi.org/10.3390/cryst14070662

Chicago/Turabian Style

Knerr, Steven, Erica Brendler, Robert Gericke, Edwin Kroke, and Jörg Wagler. 2024. "Two Modifications of Nitrilotris(methylenephenylphosphinic) Acid: A Polymeric Network with Intermolecular (O=P–O–H)3 vs. Monomeric Molecules with Intramolecular (O=P–O–H)3 Hydrogen Bond Cyclotrimers" Crystals 14, no. 7: 662. https://doi.org/10.3390/cryst14070662

APA Style

Knerr, S., Brendler, E., Gericke, R., Kroke, E., & Wagler, J. (2024). Two Modifications of Nitrilotris(methylenephenylphosphinic) Acid: A Polymeric Network with Intermolecular (O=P–O–H)3 vs. Monomeric Molecules with Intramolecular (O=P–O–H)3 Hydrogen Bond Cyclotrimers. Crystals, 14(7), 662. https://doi.org/10.3390/cryst14070662

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop