Inorganic Salts of N-phenylbiguanidium(1+)—Novel Family with Promising Representatives for Nonlinear Optics

Seven inorganic salts containing N-phenylbiguanide as a prospective organic molecular carrier of nonlinear optical properties were prepared and studied within our research of novel hydrogen-bonded materials for nonlinear optics (NLO). All seven salts, namely N-phenylbiguanidium(1+) nitrate (C2/c), N-phenylbiguanidium(1+) perchlorate (P-1), N-phenylbiguanidium(1+) hydrogen carbonate (P21/c), bis(N-phenylbiguanidium(1+)) sulfate (C2), bis(N-phenylbiguanidium(1+)) hydrogen phosphate sesquihydrate (P-1), bis(N-phenylbiguanidium(1+)) phosphite (P21), and bis(N-phenylbiguanidium(1+)) phosphite dihydrate (P21/n), were characterised by X-ray diffraction (powder and single-crystal X-ray diffraction) and by vibrational spectroscopy (FTIR and Raman). Two salts with non-centrosymmetric crystal structures—bis(N-phenylbiguanidium(1+)) sulfate and bis(N-phenylbiguanidium(1+)) phosphite—were further studied to examine their linear and nonlinear optical properties using experimental and computational methods. As a highly SHG-efficient and phase-matchable material transparent down to 320 nm and thermally stable to 483 K, bis(N-phenylbiguanidium(1+)) sulfate is a promising novel candidate for NLO.

The most efficient way of designing and preparing these molecular materials consists of using a two-step approach, as described previously in [12]: (a) molecular engineering, i.e., synthesis (selection) of a proper molecular carrier with NLO properties and (b) crystal structure engineering, i.e., using suitable intermolecular interactions and/or co-crystallisation partners to grow crystals with the desired symmetry and NLO properties.
The symmetry of these molecules and proper crystal class is crucial, especially for χ (2) NLO effects (e.g., SHG), where the absence of an inversion centre is the symmetry requirement.
Crystal engineering of structures containing polar organic molecules mostly utilises hydrogen bonding as the dominant inter-molecular interaction for adequate crystal packing, and the role of this bonding has been widely discussed in previous studies [12][13][14][15][16][17]. Furthermore, implementing a suitable hydrogen-bonded counter molecule/ion can counteract undesirable tendencies of the selected organic molecules (cations) to form centrosymmetric pairs. The concept [18], based on the use of chains/networks formed by inorganic hydrogenbonded anions (i.e., anions acting as hydrogen bonds donors), was successfully applied to several salts and cocrystals prepared in our laboratory, e.g., references [19][20][21][22]. However, the presence of deprotonated anions acting as strict hydrogen bond acceptors can also lead to the formation of promising non-centrosymmetric phases, e.g., references [19,[23][24][25].
Biguanides and their derivatives are an interesting class of compounds not only for their extensive medicinal applications but also for their ability to form stable metal complexes and intermediates useful for heterocycles synthesis and organocatalysis. Their applications and structural characteristics have been recently reviewed [26]. Furthermore, biguanides and their derivatives, such as N-phenylbiguanides, have been considered promising building blocks in crystal engineering. Crystals containing their cations have been studied as potential proton transfer materials [27,28] and prospective materials for NLO [24,29,30]. DFT computations on the molecular level indicate that biguanidium and N-phenylbiguanidium cations exhibit promising values of hyperpolarisability [29][30][31], as confirmed by the experimental determination of SHG efficiency, especially for inorganic salts of biguanide [24].

Crystal Structures
The crystal structures were solved using the single-crystal X-ray diffraction methodbasic crystallographic data of all studied salts are summarised in Tables 1 and 2. Measurement and refinement details are presented in Tables S1 and S2, Supplementary Materials. Selected bond lengths and angles, including those of hydrogen bonds, are presented in Tables S3-S9, Supplementary Materials. All figures were created using the PLATON software [35].
N-phenylbiguanidium(1+) nitrate, phbiguaNO 3 , crystallises in the monoclinic system (C2/c space group). Atom numbering is shown in Figure S1, Supplementary Materials. The crystal structure is based on chains containing alternating cations and anions interconnected via N-H...O hydrogen bonds, which form double layers, as shown in Figure 1.   The chains in these layers can be described [36,37] by the graph set descriptor C(7). Phenyl rings of the cations are oriented outwards from these layers. All nitrate anions are hydrogen-bonded to four N-phenylbiguanidium(1+) cations. Two of these cations are involved in eight-membered hydrogen-bonded rings-graph set descriptor R N-phenylbiguanidium(1+) perchlorate, phbiguaClO 4 , crystallises in the triclinic system (space group P-1). Atom numbering is shown in Figure S2  N-phenylbiguanidium(1+) hydrogen carbonate, phbiguaHCO 3 , crystallises in the monoclinic system (space group P2 1 /c). The atom numbering is shown in Figure S3  Bis(N-phenylbiguanidium(1+)) sulfate, phbigua 2 SO 4 , crystallises in the monoclinic system (non-centrosymmetric space group C2). Atom numbering is shown in Figure S4, Supplementary Materials. The crystal structure is based on layers (parallel to the ab plane) which consist of alternating cations and anions (see  The resulting tubular crystal packing with phenyl rings oriented towards cavity axes is shown in Figure 5. Each hydrogen phosphate anion interacts with three water molecules and five cations, in addition to interacting with its neighbouring anion. The resulting hydrogen-bonded rings can be described by graph set descriptors R 2 1 (4), R 1 2 (6), R 2 2 (8), R 2 3 (8), and R 3 4 (10). Bis(N-phenylbiguanidium(1+)) phosphite, phbigua 2 HPO 3 , crystallises in the monoclinic system (non-centrosymmetric space group P2 1 ). Atom numbering is shown in Figure  S5  Bis(N-phenylbiguanidium(1+)) phosphite dihydrate, phbigua 2 HPO 3 2H 2 O, crystallises in the monoclinic system (space group P2 1 /n). Atom numbering is shown in Figure S6, Supplementary Materials. The cations form chains along the b axis through intermolecular N-H...N hydrogen bonds (i.e., N12-H12B . . . N13 a and N22-H22B . . . N23 d with N...N distances of 3.169(2) and 3.172(2) Å, respectively), which can be described by graph set descriptor C(4). Pairs of anions form a cluster with four water molecules (see Figure S7  From a crystal engineering standpoint, the anions that are present in the structures can be divided into several sub-groups, according to different criteria. Considering their hydrogen bonding potential, we can distinguish strict H-bond acceptors (nitrates, perchlorates, sulfates, and phosphites) and mixed H-bond donors/acceptors (hydrogen carbonates and hydrogen phosphates). By considering only general symmetry (i.e., disregarding the position of their hydrogen atoms), we can differentiate planar (nitrates and hydrogen carbonates), trigonal pyramidal (hydrogen phosphates and phosphites), and tetrahedral (perchlorates and sulfates) anions.
Surprisingly, these criteria have only a limited effect on the final crystal packing. Planar anions present in phbiguaNO 3 and phbiguaHCO 3 crystals, as expected, allow the formation of layered structures with the phenyl ring oriented outwards from the double layers. Unexpectedly, the tetrahedral sulfate anions are involved in a very similar crystal packing in the non-centrosymmetric phbigua 2 SO 4 .
The trigonal pyramidal hydrogen phosphates and phosphites are involved in tubular structures of phbigua 2 HPO 4 1.5H 2 O and phbigua 2 HPO 3 2H 2 O, with their phenyl rings oriented towards the axes of the cavities. In addition, the presence of water molecules likely contributes to the formation of this type of crystal packing.
The crystal structure of phbiguaClO 4 , containing tetrahedral perchlorate anions, can be considered, to some extent, a crossover between a layered structure (with layers interconnected through intermolecular H-bonds forming centrosymmetric pairs of cations) and the 3D crystal packing that is observed in non-centrosymmetric phbigua 2 HPO 3 , containing pyramidal phosphite anions.
The conformation of N-phenylbiguanidium(1+) cations in the crystal structures is further stabilised by one (in phbigua 2 SO 4 and in phbigua 2 HPO 3 -one of the cations) or two N-H...N intramolecular hydrogen bonds (see Tables S3-S9, Supplementary Materials) described by graph set descriptor S(6). A very similar situation was also observed in the previously studied biguanide salts [24,29,32,33].
Most of these structures also contain at least one weak intramolecular C-H . . . N interaction (except for phbiguaClO 4 ), which involves a phenyl ring C-H group and the nonprotonated nitrogen N3 (or NY3, where Y = 1-8)-see Tables S3-S9 All aforementioned crystal structures of N-phenylbiguanidium(1+) salts contain an extensive network of hydrogen bonds involving nitrogen atoms. These interactions significantly affect the geometry of the cations. Cation conformation can be discussed based on two main parameters: (a) the interplanar angle of two plains of "guanidine" sub-skeletons forming the biguanide part of the cation and (b) the interplanar angle between the phenyl ring plane and the neighbouring "guanidine" plane.
The values of interplanar angles in the biguanide part of the cations range from 36 to 50 • . To compare the values determined in previously studied biguanide and Nphenylbiguanide salts, see references [24,29,30,32,33]. The values of the interplanar angles of the phenyl ring plane and neighbouring "guanidine" plane range from 15 to 71 • . Monitoring these conformation parameters has proved crucial for explaining the isostructural phase transformation mechanism of bis(N-phenylbiguanidium(1+)) oxalate [30].

Vibrational Spectra
Vibrational spectra of the inorganic salts of N-phenylbiguanidium(1+) are presented in Figures 8 and 9, and the list of the recorded maxima is provided in Supplementary Materials (see Listing S1, Supplementary Materials). The IR spectra are dominated by structured, strong-to-medium intensity bands (3500-2500 cm   Raman spectra of all seven salts contain bands of C-H stretching modes of an aromatic ring (3080-3060 cm −1 region). The intensity and shape of these bands can be correlated with the involvement of phenyl rings in C-H...X (X = N or O) hydrogen bonds (see Figure 8 and Tables S3-S9, Supplementary Materials). Strong symmetric bands were observed in phbiguaClO 4 and phbiguaHCO 3 with no or only single, weak intramolecular C-H...N interactions, respectively. The formation of slightly stronger intramolecular and/or intermolecular C-H...N or intermolecular C-H...O interactions is reflected by medium-intensity asymmetric doublets (phbigua 2 SO 4 , phbiguaNO 3 , phbigua 2 HPO 3 , and phbigua 2 HPO 4 1.5H 2 O) or even by the quadruplet (phbigua 2 HPO 3 2H 2 O). Other intense bands typical for phenyl groups were recorded at~1600 cm −1 and~1000 cm −1 , corresponding to stretching C-C vibrations and "breathing" modes of phenyl rings, respectively. The sharp and very strong bands at 1056, 935, and 980 cm −1 in the Raman spectra of phbiguaNO 3 , phbiguaClO 4 , and phbigua 2 SO 4 can be assigned to symmetric stretching modes of NO 3 − , ClO 4 − , and SO 4 2− anions, respectively. The solid-state vibrational modes of phbigua 2 SO 4 and phbigua 2 HPO 3 crystals were devised (see Figures S8-S11, Supplementary Materials) within the computational characterisation of these non-centrosymmetric materials. In addition to the standard B3LYP functional, the advanced CAM-B3LYP functional (incorporating additional long-range corrections) was used to improve the agreement between computed vibrational modes and the recorded spectra of phbigua 2 SO 4 . The latter approach was particularly successful in improving the frequencies of the calculated modes (mainly of those involving contributions from the antisymmetric stretching vibrations of SO 4 2− anions) in the 1150-1000 cm −1 region of the IR spectrum (see Figure S8, Supplementary Materials) and the intensities calculated below 950 cm −1 in the Raman spectrum (see Figure S9, Supplementary Materials).
These preliminary conclusions about recorded spectra are based on the literature on vibrational manifestations of inorganic anions [38] and on previous spectroscopic studies focused on bis(N-phenylbiguanidium(1+)) oxalate polymorphs, including N-phenylbiguanidium(1+) chloride [30]. The spectra of N-phenylbiguanide inorganic salts will be assigned in detail in a following paper.

Thermal Behaviour
To evaluate the thermal stability of the materials prepared in this study, DSC measurements were performed for all compounds in the temperature range from 93 K to a temperature close to their melting points. No thermal anomalies were recorded in the low-temperature region, and crystals of phbiguaNO 3

Linear and Nonlinear Optical Properties of phbigua 2 SO 4 and phbigua 2 HPO 3
Based on the results from crystal structure determination and solid-state DFT computations, the non-centrosymmetric crystals of phbigua 2 SO 4 and phbigua 2 HPO 3 are optically anisotropic biaxial negative and biaxial positive materials, respectively. The values of the calculated refractive indices and static (λ = ∞) nonlinear susceptibility tensor χ (2) components are presented in Table 3. The results indicate that phbigua 2 SO 4 crystals have larger birefringence than phbigua 2 HPO 3 crystals. The refractive indices of phbigua 2 SO 4 assessed using the CAM-B3LYP functional were systematically lower than those determined using the B3LYP functional. The results of the main χ (2) component values of phbigua 2 SO 4 were, nevertheless, quite similar when comparing the functionals that were used in this study (except for the opposite sign of the of χ (2) yyy component). The SHG measurements of powdered samples were performed for phbigua 2 SO 4 and phbigua 2 HPO 3 crystals, using powdered KDP as the reference material. These measurements yielded a relative SHG efficiency d phbigua2SO4 = 5.57 d KDP and d phbigua2HPO3 = 1.60 d KDP at 800 nm laser irradiation. The value of SHG efficiency determined for phbigua 2 HPO 3 is slightly underestimated due to the presence of approximately 5% (as determined by powder XRD) of the centrosymmetric phbigua 2 HPO 3 2H 2 O phase. The phase-matchability of the more promising phbigua 2 SO 4 was studied by particle-size-dependent measurements with an 800 nm laser line (see Figure 10). The results clearly confirm the possibility of SHG phase matching for this salt. The small discrepancy between the two smallest fractions can be explained by problems in preparing powders with isometric particles and effective control of the particle size in the family of soft molecular crystals, which are quite common (e.g., [22,39,40]). Figure 10. Phase-matching curve (i.e., particle size vs. SHG intensity) of phbigua 2 SO 4 (800 nm fundamental wavelength); the red curve drawn is to guide the eye and is not fitted to the data. Black horizontal segments represent particle size intervals.
For further optical characterisation of phbigua 2 SO 4 , UV-Vis spectra of its powdered sample were recorded (see Figure S12, Supplementary Materials), confirming the wide transparency range of this material, at least from the Vis region to 320 nm.
We have also preliminarily determined the second-order nonlinear optical tensor [d SHG ijk ] of the phbigua 2 SO 4 crystal using the Maker-fringe method [41] at 1064 nm laser irradiation. The tensor of monoclinic phbigua 2 SO 4 (space group C2) has eight independent components and in the contracted form (Voigt's notation using d mn components), it reads as [42]: Unfortunately, accurate experimentally determined values of refractive indices remain unknown, and the size of the bulk single crystals grown was insufficient for their precise determination. Therefore, we were only able to estimate the magnitudes of d mn s. The results indicate that the component d 21 is larger than the component d 36 in reference KDP [43].

Methods
The collection of X-ray diffraction data of suitable single crystals was performed on the following experimental devices: KappaCCD (Bruker Nonius, Billerica, MA, USA) diffractometer (MoK α radiation, graphite monochromator, COLLECT [44], and DENZO [45] software) for samples phbiguaNO 3 , phbiguaClO 4 , phbiguaHCO 3 , phbigua 2 SO 4 , phbigua 2 HPO 4 1.5H 2 O, and phbigua 2 HPO 3 2H 2 O; KM4 CCD (Oxford Diffraction Ltd., Abingdon, UK) diffractometer (CrysAlis CCD [46] software) for phbigua 2 HPO 3 2H 2 O and D8 VENTURE Kappa Duo (Bruker, Billerica, MA, USA) diffractometer with PHOTON III CMOS detector and IµS micro-source with MoK α radiation and SAINT V8.40B [47] software for phbigua 2 HPO 3 . The temperature of the crystal samples was controlled by a liquid nitrogen Cryostream Cooler (Oxford Cryosystems, Oxford, UK). The phase problem was solved by direct methods (with corresponding diffractometer software [48][49][50]) and the non-hydrogen atoms were refined anisotropically, using the full-matrix least-squares procedure (appropriate version of SHELX software [50,51]). Hydrogen atoms attached to carbon atoms were calculated in a geometrically idealised position with Csp 2 -H = 0.93 Å, and constrained to ride on their parent atoms, with U iso (H) = 1.5U eq (C). The positions of the hydrogen atoms attached to the oxygen and nitrogen atoms were localised on difference Fourier maps and were fixed during refinement using rigid body approximation with assigned displacement parameters equal to 1.2U iso (pivot atom). A hydrogen atom attached to a phosphorus atom was localised on a difference Fourier map and refined isotropically.
The determination of absolute structure of the non-centrosymmetric crystals was based on anomalous-dispersion effects of either the phosphorus or the sulphur atoms; the resulting chirality parameters [52] were 0.01(3) and −0.01(2) for phbigua 2 HPO 3 and phbigua 2 SO 4 , respectively.
Powder X-ray diffraction patterns of prepared polycrystalline products (see Tables S10-S16, Supplementary Materials) were collected at room temperature using a Philips/PANalytical (Royston, UK) X'pert PRO MPD diffractometer (Bragg-Brentano geometry, ultrafast X'Celerator detector, and Cu Kα radiation, λ = 1.5418 Å). Theoretical diffraction patterns for confirmation of phase purity of prepared products were calculated using the FullProf software [53].
Infrared spectra of powdered samples were recorded by DRIFTS technique (samples mixed with KBr at~1:50 ratio and grinded in agate mortar) on a Thermo Fisher Scientific (Madison, WI, USA) Nicolet 6700 FTIR spectrometer in the 400-4000 cm −1 region (2 cm −1 resolution, KBr beamsplitter, Happ-Genzel apodization). Raman spectra of powdered samples were recorded on a Nicolet 6700 FTIR spectrometer equipped with a Thermo Fisher Scientific (Madison, WI, USA) Nicolet Nexus FT Raman module (2 cm −1 resolution, Happ-Genzel apodization, 1064 nm Nd:YVO 4 laser excitation,~400 mW power at the sample) in the 100-3700 cm −1 region. Raman spectra of microcrystalline samples were also collected on a dispersive confocal Raman microscope MonoVista CRS+ (Spectroscopy & Imaging GmbH, Warstein, Germany) interfaced to an Olympus microscope (objectives 20x and 50x) using a 785 nm diode excitation laser (10 mW laser power, 40-3800 cm −1 spectral range, 300 lines/mm grating). The spectrometer was wavelength-and intensity-calibrated using a software-controlled auto-alignment and calibration procedure with mercury and Ne-Ar lamps.
The UV-Vis spectra of phbigua 2 SO 4 powdered samples were recorded using an Agilent Technologies (Santa Clara, CA, USA) Cary Series 4000 UV-Vis spectrometer equipped with internal DRA-900 accessory (diffuse reflectance with integrating sphere) in the 190-850 nm spectral region with a 0.5 nm spectral resolution.
DSC measurements were carried out on Perkin Elmer (Waltham, MA, USA) Pyris Diamond DSC and DSC 7 instruments in the 93-473 K temperature region (20 mL/min nitrogen or helium flow above and below 298 K, respectively). A heating rate of 10 K/min was selected to measure approximately 20 mg of the finely ground sample placed in hermetically sealed aluminium pans.
SHG measurements of powdered phbigua 2 HPO 3 and phbigua 2 SO 4 were performed using the modified Kurtz-Perry method [54]. The samples were irradiated with 90 fs laser pulses generated at an 82 MHz repetition rate by a Spectra Physics (Santa Clara, CA, USA) Tsunami Ti:sapphire laser at 800 nm. For quantitative determination of SHG efficiency, the intensity of the back-scattered laser light generated in the sample at 400 nm was measured on a grating spectrograph with a diode array (InstaSpect II, Oriel Newport Corporation, Irvine, CA, USA), and the signal was compared with that produced by a potassium dihydrogen phosphate (KDP) standard. The initial experiments were performed on a powdered sample (75-150 µm particle size) loaded into a 5 mm glass cell using a mechanical vibrator. The measurements were repeated on different areas of the same sample, and the results were averaged. This experimental procedure was employed to minimise signal fluctuations induced by sample packing. Lastly, the measurements of phbigua 2 SO 4 were performed also with size-fractioned samples (particle size: 25-45, 45-63, 63-75, 75-100, 100-125, and 125-150 µm).
The SHG study focused on the determination of particular d mn coefficients of the second-order NLO tensor that were performed on cuts of phbigua 2 SO 4 single crystals using the standard Maker-fringe method [41]; with a BMI (Evry, France) Q-switched Nd-YAG laser (1064 nm,~6 ns pulses of less than 0.1 mJ at 20 Hz) as a light source. For details of the applied measurement strategy (monoclinic crystals of point group 2), see e.g., [31]. Single crystals were grown by isothermal crystallisation of aqueous solution at 296 K in a dark room for several months.
Solid-state DFT studies of phbigua 2 SO 4 and phbigua 2 HPO 3 were carried out using the CRYSTAL17 program (version 1.0.2.) [55]. Primary computations employed the B3LYP functional; the 6-31+G(d,p) basis set was used for all oxygen and the 6-31G(d) basis set for all other atoms. DFT integrations utilised the large integration grid implemented in the program. The Pack-Monkhorst and Gilat nets used for sampling the Brillouin zone consisted each of 8 points. Derivatives required to compute IR and Raman spectra were obtained by the coupled-perturbed Kohn-Sham analytical approach [56,57]. Unit cell parameters were not optimised. The starting atomic coordinates were those obtained from diffraction experiments, which were subjected to optimisation that used tightened convergence criteria. Spectral and optical properties were subsequently computed on the stationary-point geometries obtained. Theoretical Raman spectra were corrected empirically for experimental wavelength 1064 nm and temperature 293 K. Vibrational frequencies obtained from computations were scaled uniformly by an empirical factor of 0.97. To improve the fit of computed vibrational frequencies and recorded spectra of phbigua 2 SO 4 , the CAM-B3LYP functional and the triple-zeta POB-TZVP basis set [58] for all atoms was also used. The Pack-Monkhorst and Gilat nets used for sampling the Brillouin zone in this case consisted each of 12 points. Tolerance criteria for coulomb and exchange sums were tightened by using the respective values 7, 7, 7, 9, and 30 for the TOLINTEG input keyword. Vibrational frequencies obtained from computations were scaled uniformly by an empirical factor of 0.96.

Conclusions
Seven novel salts were prepared and characterised in studied systems containing N-phenylbiguanide and selected inorganic oxyacids. The involved inorganic anions differ in their symmetry and in the role that they play in hydrogen bond networks. The crystal packing of these salts ranges from layered (phbiguaNO 3 , phbiguaHCO 3 , and phbigua 2 SO 4 ) to tubular (phbigua 2 HPO 4 1.5H 2 O and phbigua 2 HPO 3 2H 2 O) structures, which are enabled not only by the suitable symmetry of their anions but also by the involvement of water molecules in the H-bonding patterns. Two salts-phbigua 2 SO 4 and phbigua 2 HPO 3 -met our expectations of a non-centrosymmetric crystal arrangement, necessary for observing χ (2) NLO effects.
The results from the recorded vibrational spectra match the results of crystal structure determination, i.e., crystals based on hydrogen-bonded networks containing Nphenylbiguanidium(1+) cations and inorganic anions (eventually water molecules). The formation of various hydrogen bonds is reflected in vibrational spectra by the presence of strong wide bands (IR spectra) of stretching vibrations of N-H and O-H groups involved in these bonds, by the different positions of P-H stretching modes in phbigua 2 HPO 3 2H 2 O and phbigua 2 HPO 3 (both spectra) and by the shape of the bands of C-H stretching vibrations in phbigua 2 SO 4 , phbiguaNO 3 , phbigua 2 HPO 3 , phbigua 2 HPO 4 1.5H 2 O, and phbigua 2 HPO 3 2H 2 O (Raman spectra).
Powder SHG measurements of phbigua 2 SO 4 coupled with preliminary single-crystal optical measurements confirmed that this compound especially is a promising (SHG efficient, phase-matchable, transparent down to 320 nm, and thermally stable to 483 K) novel material for NLO. In addition to favourable χ (2) NLO properties, both non-centrosymmetric salts could also be considered as potential SRS-active materials for Raman laser converters. Especially intense sharp Raman bands of sulfate symmetric stretching vibration of phbigua 2 SO 4 and phosphite stretching P-H vibration of phbigua 2 HPO 3 could, similarly as in relative guanidine-based compounds [10,11], act as SRS-promoting vibrational modes. Unfortunately, problems related to bulk single-crystals growth complicate broader technical applications of phbigua 2 SO 4 crystals.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.