Gholam Hossien Rounaghi*a,
Mostafa Gholizadeha,
Fatemeh Moosavia,
Iman Razavipanaha,
Hossein Azizi-Toupkanloob and
Mohammad Reza Salavatia
aDepartment of Chemistry, Faculty of Sciences, Ferdowsi University of Mashhad, Mashhad, Iran. E-mail: ghrounaghi@yahoo.com; ronaghi@um.ac.ir; Fax: +98-51-38796416; Tel: +98-51-37626388
bDepartment of Chemistry, Faculty of Sciences, University of Neyshabur, Neyshabur, Iran
First published on 15th December 2015
The complexation reaction between the La3+ cation and a macrocyclic ligand (kryptofix 22DD) in ordinary methanol (MeOH) and extraordinary methanol solvent (magnetized methanol) exposed to a magnetic field at different temperatures has been investigated using the conductometric method, in order to determine the effect of magnetic solvents on the thermodynamics of the complexation process between the metal cation and the ligand. The stoichiometry of the complex formed between the La3+ cation and the ligand is 1:
1 [M
:
L] in both solvent systems. Comparison of the corresponding molar conductivity versus [L]t/[M]t for (kryptofix 22DD·La)3+ complex shows that the complex formed between the metal cation and kryptofix 22DD in the magnetized methanol is weaker compared to the ordinary methanol. The results of the density functional theory shows that addition of the metal ion to the ligand completely deforms the ligand from its planar shape and the orientation of two carbon chains attached to the ring changes in comparison with the free ligand. In addition, the molecular dynamics simulation can shed light on the influence of the magnetic field, from the molecular point of view, on the properties of solvents and is able to explain the variation of hydrogen bond and transport properties. These changes may be the reason for the interaction of the solvent and the solute as well as the stability constant of the complex formed between the La3+ cation and kryptofix 22DD in solution. As the strength of the magnetic field increases, the number of the hydrogen bonds formed between the methanol molecules increases and the methanol self-diffusion coefficient decreases. These changes restrict the movement of the molecules which results in increasing the viscosity of the methanol and therefore the stability constant of the (kryptofix 22DD·La)3+ complex decreases in magnetized methanol with respect to the ordinary methanol.
Lanthanides are widely distributed in low concentrations throughout the Earth's crust.10 La3+ cation with potentially high and variable coordination numbers and flexible coordination environments behaves as a hard acid with a strong affinity for a hard base containing neutral and negatively charged oxygen atoms.11–15
In solutions, the metal ion, the macrocyclic ligand, and also the metal ion–macrocyclic ligand complex would be surrounded by the solvent molecules. Because of different physicochemical properties, each solvent can effect on the stability of metal ion complexes in different ways. Therefore, it is of great importance to study the effect of solvent on the complexation of metal ions with the macrocyclic ligands. One approach to achieve a detailed understanding of these interactions is to study the influence of magnetic field on the physicochemical properties of solvent and as a consequence on the formation of metal ion–macrocyclic complexes in solutions.
As the macrocyclic compounds are good models for the antibiotics and some of the other drugs treatment in biological systems,16 it is of great interest to investigate the effect of magnetic field upon complexation of such ligands with the metal ions. To our knowledge, the data about the effect of magnetic field upon formation of metal ion-macrocyclic complexes in non-aqueous solvents are very sparse. The first goal of the present study is to investigate how the magnetic field can influence on the physicochemical properties of a non-aqueous solvent passing through a magnetic field and as a result, the stability of the metal–ion complex at different temperatures.
The changes occurring in the structure of solvent under the effect of an external magnetic field are important in various applications such as chemical processes. It has been found that various aspects of the most common solvent (water) properties, including the size of the cluster,17–19 hydrogen bonding,20,21 electrical conductivity,22 dissolution rate into water,23,24 van der Waals bonding between molecules,25 and its melting point21 change when it is exposed to a magnetic field. Special attention has been devoted towards the structure of alloy clusters26 as well as their growth.27 Attempts to generate perpendicular magnetic anisotropy have been performed to check the shape anisotropy due to dipolar interactions, experimentally and theoretically. In addition, the use of a magnetic field is of considerable interest in a number of practical applications due to the self-diffusion coefficient alteration under a magnetic field. From molecular point of view, temperature influences on the properties of organic solvents in the presence of a magnetic field, the point that has not yet been under investigation.
Molecular dynamics (MD) simulations provide a powerful means of investigating the enhanced hydrogen bonding mechanism from an atomic viewpoint. Moreover, it is able to explore how surrounded ions by the solvent molecules may be retarded providing complex. In other words, solvent can effect on metal ion complex stability; as a result, achieving detailed insight on the influence of magnetic field on the physicochemical properties of solvent is of great importance. Temperature effect of external magnetic field on the properties of organic solvents has not yet been reported. Therefore, the temperature effect on external magnetic field change with the properties of methanol solvent is also under consideration to answer the above points in line with experimental results.
The static magnetic field in a compact form, a unit called “AQUA CORRECT”, was used. The equipment has a coaxial static magnetic system of 6000 G field strength, and was imported from Germany (H.P.S Co.) [DN = 20, 3/4 in, flow 2 m3 h−1] for the experiments. The equipment was connected from one end to the liquid pump and the other end to the pipelines of the solvent reservoir. The solvent had to flow through a coaxial static magnetic gap and came back to the solvent reservoir. Therefore, the solvent could pass through the magnetic field for many times in a closed cycle.
![]() | (1) |
Neighborhood lists were updated every 10 time steps with a distance of 15 Å. Full electrostatic evaluations were calculated. The coupling methods for pressure and temperature were applied by Nose–Hoover thermostat–barostat every 1.0 and 0.1 ps, respectively.36,37 Using Verlet leapfrog scheme,38 the positions and velocities of particles were calculated.
The external magnetic field aligned along the Z direction at different temperatures may shed light on the thermodynamic, transport, and structural properties of the solvent using MD computation. All MD simulations were conducted using DL_POLY version 2.17 simulation package.39,40 Notice that all mentioned details of simulation were applied in the case that the complex is present; mole fraction of the complex is 0.002 and the effect of temperature and magnetic field was studied. In this case, the cut off distance was set to 16 Å.
In this study, a static and constant magnetic field parallel to the Z-direction of the simulation box was applied. According to the Lorentz force, the output of the magnetic field depends on the velocity of solvent flow, magnetic force, number of particles with positive or negative charges in solvent, and the angle between particle velocity and applied magnetic field, which is always upright.
An accurate practical means for the analysis of the structure is provided by pair correlation function, g(r), defined as:
![]() | (2) |
![]() | (3) |
Besides, the molecular mobility was investigated by mean-square displacement.
MSD = 〈|rj(t) − rj(0)|〉2 | (4) |
![]() | ||
Fig. 1 Molar conductance–mole ratio plots for (kryptofix 22DD·La)3+ complex in methanol solution at different temperatures: 25 °C (open circle), 35 °C (asterisk), 45 °C (open triangle). |
As is shown in Fig. 1, the slope of the corresponding molar conductivity versus [L]t/[M]t changes at the point where the ligand to cation mole ratio is about one which is an evidence for formation of a relatively stable 1:
1 [M
:
L] complex between kryptofix 22DD and La3+ cation in ordinary methanol solvent. On the other hand, the plot of corresponding molar conductivity versus [L]t/[M]t does not show any considerable change in slope emphasizing the formation of relatively weak 1
:
1 [M
:
L] complex in the magnetized methanol solvent (see Fig. 2).
In order to make the 1:
1 [M
:
L] complex model more clear, the fitting and the experimental curves for the (kryptofix 22DD·La)3+ complex in magnetized methanol solvent at 298.15 K is shown in Fig. 3. As it is clearly evident in this Figure, there is an accurate acceptable agreement between experimental and least square regression.
![]() | ||
Fig. 3 The experimental data and fitting curve for (kryptofix 22DD·La)3+ complex in magnetic methanolic solution at 298.15 K. |
To gather accurate information about the conformational change of kryptofix 22DD upon complexation with lanthanum(III) cation, the molecular structures of the uncomplexed ligand and its 1:
1 La3+ complex were computed with quantum calculations. The structure of the free ligand was optimized at B3LYP/Lanl2dz level of theory. The optimized structure of the ligand was then applied to find out the initial structure of its 1
:
1 lanthanum(III) complex (at the same level of theory). The optimized structures of the uncomplexed ligand (i.e. kryptofix 22DD) and its 1
:
1 complex with La3+ cation (i.e. kryptofix 22DD·La)3+ are presented in Fig. 4. In order to study the effect of solvent upon complexation, we compared the order of stability of the complexes formed between kryptofix 22DD and La3+ cation in the gas phase which were obtained based on binding energies with the order of stability of the complexes formed between kryptofix 22DD and La3+ cation in the presence of methanol. Fig. 4A and B shows the optimized structures of free ligand and its complex with La3+ cation in the gas phase. It is obvious from Fig. 4A that the ligand forms a more or less planar but interestingly, in the optimized structure of the 1
:
1 (kryptofix 22DD·La)3+ complex (Fig. 4B), lanthanum(III) ion is well incorporated inside the twisted macrocyclic ligand and coordinated to six donating atoms of the ligand, including four oxygen atoms and two nitrogen. It is fascinating to note that when the metal ion is added to the ligand, the ligand was twisted and completely deformed from its planar shape. Moreover, the orientation of two carbon chains attached to the ring changes in comparison with the free ligand. The optimized structures of free ligand and its complex with La3+ in the presence of methanol were also studied. For the calculations in the presence of methanol as a solvent, the Polarizable Continuum Model (PCM) was chosen. In this model, the cavity is created via a series of overlapping spheres. The PCM was applied within the B3LYP/LanL2DZ level to predict the solvent effects on the structures and interaction energy of kryptofix 22DD ligand and [kryptofix 22DD·La]3+. The results show there is no major changes occurred in the structures of free ligand and its complex in compression to the gas phase calculation. For this reason, the structures of free ligand and its complex are not shown in the presence of methanol. Moreover, the effect of solvent on the binding energy of complex was investigated. The pair wise binding energy ΔE between kryptofix 22DD ligand and La3+ is estimated as the difference between the energy of the complex [kryptofix 22DD·La]3+ and the energies of isolated partners:
ΔE = E[kryptofix 22DD·La]3+ − (ELa3+ + Ekryptofix 22DD) |
The binding energy for the [kryptofix 22DD·La]3+ complex were obtained −53.1 and −481.2 kcal mol−1 in the gas phase and in the presence of methanol, respectively. These results show that the stable complex is formed in the presence of methanol as a solvent.
Now, the question is how the magnetic field impacts on the solvent and how this magnetized solvent can change the thermodynamic of the complexation process between the kryptofix 22DD and La3+ cation. Chang and Weng17 have investigated the effects of the magnetic field on the hydrogen-bonded structure of water and they found that the number of the hydrogen bonds increases when it is exposed to a magnetic field. Since water molecules have two donors and two acceptor sites, liquid water forms a continuous network of hydrogen bonds. The presence of spatial network of hydrogen bonds is the main property of the liquid water.41 Furthermore, some researchers have carried out a series of studies to know if there is any relationship between the hydrogen bonding and temperature of water.42,43 It has been recognized that the hydrogen bonding becomes weaker with increasing the thermal motion of the atoms involved.44 Consequently, the water structure breaks down as the temperature increases. It has been suggested that the effects of a magnetic field on water hydrogen bonds are similar with the effects of decreasing the temperature.45 The theory of cluster models points out that the average size of water cluster decreases with increasing temperature.46,47 Normally, the viscosity and also the dielectric constant of liquid water decreases as the temperature increases.48,49 From the other side of view, methanol is a polar solvent and its molecules are able to constitute hydrogen bonds through their hydrogen and oxygen atoms; therefore, it is expected a similar behavior is observed for this organic solvent as water. It has been shown that the hydrogen bond network in methanol molecules increases as the temperature decreases.50 In addition, the viscosity and the dielectric constant of the liquid methanol increase as the temperature decreases.49,51
The values of the stability constant (logKf) of (kryptofix 22DD·La)3+ complex are given in Table 1. It can be seen that by applying magnetic field, the stability of the (kryptofix 22DD·La)3+ complex decreases at all studied temperatures. Similar results have been observed by the authors in the previous study which further corroborate the results reported in Table 1.52 This result may be attributed to the higher dielectric constant and possibly changing the donor number of the magnetized methanol solvent than that of the ordinary methanol. The higher dielectric constant and change in the donicity of the magnetized solvent with respect to the ordinary methanol may lessen the electrostatic interactions between the ligand and La3+ cation. As a result, the complex formation between kryptofix and La3+ cation weakens by imposing magnetic field on the solvent. Furthermore, it seems that the higher viscosity of the magnetized methanol solvent, leads to a decrease in the rate of the complexation reaction in this extraordinary solvent compared to the ordinary methanol solvent.
In order to shed light on this suggestion, MD was handled and the results are given in the following paragraphs.
![]() | ||
Fig. 5 Comparison of the target solvent H–O distance in the presence and absence of magnetic field at different temperatures. |
Commonly, radial distribution function (RDF) examines the structure of a liquid; in the case of methanol, gO–H and gO–O are considered to investigate the structural dependence on the external magnetic field. The other correlations as C–C, C–H, and H–H show insignificant revolution with magnetic field that are not shown. In other words, the correlations due to hydrogen bond enjoy considerable variations that are under study here. The structural differences of pure liquid with and without the application of an external magnetic field are depicted for gO–H in Fig. 6, at T = 298.15, 308.15, and 318.15 K. The first peak occurs at 2.575 Å, at the absence of the magnetic field and at T = 298.15 K, which corresponds to the oxygen–hydrogen distance between the two hydrogen-bonding methanol molecules, which reduces to 2.525 Å, at T = 298.15 K and the presence of the magnetic field. It should be noted that this distance is greater than an O–H bond length obtained by DFT computation (0.961 Å), i.e., the aforementioned distance does not correspond to the O–H bond length. The tetrahedral structure of the near neighbors is related to the second peak of O–H RDF belonging to two molecules both hydrogen bonded to a third methanol.
![]() | ||
Fig. 6 Comparison of g(r) O–H in the presence and absence of the magnetic field at different temperatures. |
In order to shed light on the compactness of methanol molecules because of magnetic field, O–O RDF, the distribution of the oxygen atoms of the methanol molecules, at T = 298.15 K, 308.15 K, and 318.15 K was also considered, Fig. 7. Surprisingly, the distance between two methanol molecules in the absence and presence of the magnetic field is unchanged (3.475 Å) and the intensity of the peak just only increases sharply. However, the height of the second peak shows a slight change indicating the distance between the two oxygen atoms of two next-nearest neighbors is not influenced by the magnetic field. From the other side of view, the average tetrahedral hydrogen bonding length may not be affected by external magnetic field. The increased height of the first and second valleys under external magnetic field demonstrates that more molecules exist between methanol shells. The slight increase in the height of the first valley indicates that more molecules are located between the shells. However, magnetic field reforms the structure of the liquid methanol, forces more methanol molecules between the shells, and enhances the connectivity between molecules to progress the stability of the solvent–solvent network. Consequently, a small increase in the number of hydrogen bonds is apparent under an external magnetic field. According to the results presented in Fig. 7, the methanol molecules tend to form more stable connections with other molecules in all directions.
![]() | ||
Fig. 7 Comparison of g(r) O–O in the presence and absence of the magnetic field at different temperatures. |
In summary, the results presented in Fig. 6 and 7 demonstrate that the magnetic field enhances the binding between the methanol molecules and stabilize the structure of the liquid methanol.
In this study, we also investigated the self-diffusion coefficients of the solvent molecules under a magnetic field. Study of the transport properties of a liquid is an important topic, both in fundamental sciences and also in its applications. The mobility of methanol molecules is indicated by the value of mean-square displacement (MSD), which depends on the temperature, pressure, structure, and so on. The value of the self-diffusion coefficients can be obtained from the Einstein relation based on MSD function. The present results show that a static magnetic field enhances the stability of methanol molecules and hence influences their mobility. Calculating the self-diffusion coefficient of the molecules provides a clearer understanding of their change in mobility. Fig. 8 presents the profiles of MSD under various temperatures and magnetic strengths. It is noticeable that the self-diffusion coefficient reduces as the magnetic field is applied. The decreasing self-diffusion coefficient by applying magnetic field points to the mobility decrease of methanol molecules. If the molecule mobility changes, the physical properties such as viscosity, thermal conductivity, and their melting point also change. Although the decrease of the self-diffusion coefficient caused by the magnetic field is not large, it indicates a change in the properties of the liquid. In general, based on the above finding, it seems that the changes of structure of the methanol solvent under the influence of the magnetic field are probably effective in changing the (kryptofix 22DD·Y)3+ complex stabilization in solution.
![]() | ||
Fig. 8 Comparison of the mean-square displacement of the solvent in the presence and absence of the magnetic field at different temperatures. |
Now it is the point to emphasize on the second target of this computation. The goal is to investigate the effect of magnetic field on the complex stabilization since a molecular knowledge about the solvent and magnetic field is present. For this aim, a system containing the complex and solvent was simulated. The results of atomic RDF between the solvent, hydrogen atom connected to oxygen atom of complex represented by H(O), and the complex, N atom, at the presence of magnetic field and its absence are presented in Fig. 9. As the figure demonstrates, the complex interaction with the solvent increases at a distance about 3 Å with magnetic field without any substantial variation of distance between solvent and complex. Interestingly, the complex structure does not witness change showing that the complex is stable at the presence of magnetic field. However, the observation was quantified by computing coordination number where at different temperatures is compared between the system at the presence of magnetic field and its absence in Table 2. As the table tabulates, the previous statement is confirmed by coordination number values.
![]() | ||
Fig. 9 Comparison of g(r) N⋯H(O) in the presence and absence of the magnetic field at different temperatures. The inset is g(r) in the absence of magnetic field at three studied temperatures. |
CN value | 298.15 K | 308.15 K | 318.15 K | |
---|---|---|---|---|
Pure methanol | Magnetic field off | 0.89 | 0.86 | 0.84 |
Magnetic field on | 1.26 | 1.25 | 1.25 | |
Solution | Magnetic field off | 0.88 | 0.86 | 0.83 |
Magnetic field on | 1.23 | 1.22 | 1.27 |
Based on Table 2, it can be understood that coordination number (CN) of methanol at the presence of magnetic field enhances that is in good agreement with the decreasing temperature. However, the complex presence leads to a very slight decrease of coordination number. Additionally, the intermolecular interaction between the complex and the solvent may lessen the magnetic field effect on the solvent; as a result, the CN is less than its corresponding value for the pure methanol and the presence of magnetic field. In other words, the strong interaction between complex and solvent at magnetic field causes retard of complex provide.
In addition, molecular dynamics simulation was applied to examine the effect of a static magnetic field on the liquid methanol at a various temperatures of 298.15, 308.15, and 318.15 K. It has been shown that an external magnetic field influences the number of hydrogen bonds, its strength, the structure of the liquid methanol, as well as its molecular mobility. The increasing number of the hydrogen bonds, due to the external magnetic field, indicates the formation of larger methanol molecule clusters. The magnetic field induces a tighter bonding between the methanol molecules and improves the stability of this liquid. Under the effect of the magnetic field, the structure of the liquid changes and more molecules are forced between the methanol shells decreasing the stability of the complex.
These molecules connect the shells and hence create a more stable methanol–methanol network. The transport properties of the molecules, as indicated by MSD, are of considerable interest in many applications. The current simulation results have shown that the molecular mobility reduces when a magnetic field is applied to the methanol solvent and, hence, resulting in a reduction in the self-diffusion coefficient of this organic solvent. In other words, the magnetic field constrains the movement of the methanol molecules, and hence changes both the thermal conduction and the viscosity in the liquid state. This change is in a good agreement with decreasing the stability constant of the (kryptofix 22DD·La)3+ complex in magnetized methanol compared to the ordinary methanol solvent.
This journal is © The Royal Society of Chemistry 2016 |