Open Access Article
This Open Access Article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence

Direct observation of Mg2+ complexes in ionic liquid solutions by 31Mg β-NMR spectroscopy

Daniel Szunyogh a, Ryan M. L. McFadden bc, Victoria L. Karner bc, Aris Chatzichristos cd, Thomas Day Goodacre e, Martin H. Dehn cd, Lia Formenti f, Derek Fujimoto cd, Alexander Gottberg e, Evan Kallenberg d, Ildikó Kálomista a, Robert F. Kiefl cd, Flemming H. Larsen g, Jens Lassen e, C. D. Philip Levy e, Ruohong Li e, W. Andrew MacFarlane bc, Iain McKenzie hi, Gerald D. Morris h, Stavroula Pallada a, Matthew R. Pearson h, Stephan P. A. Sauer a, Paul Schaffer j, Peter W. Thulstrup a, Lars Hemmingsen *a and Monika Stachura *j
aDepartment of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 København Ø, Denmark. E-mail: lhe@chem.ku.dk
bDepartment of Chemistry, The University of British Columbia, 2036 Main Mall, Vancouver, BC V6T 1Z1, Canada
cStewart Blusson Quantum Matter Institute, The University of British Columbia, 2355 East Mall, Vancouver, BC V6T 1Z4, Canada
dDepartment of Physics and Astronomy, The University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada
eAccelerator Division, TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada
fDepartment of Physics, University of Guelph, 50 Stone Road E., Guelph, ON N1G 2 W1, Canada
gDepartment of Food Science, University of Copenhagen, Rolighedsvej 26, 1958 Frederiksberg C, Denmark
hPhysical Sciences Division, TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada
iDepartment of Chemistry, Simon Fraser University, 8888 University Drive, Burnaby, BC V5A 1S6, Canada
jLife Sciences Division, TRIUMF, 4004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada. E-mail: mstachura@triumf.ca

Received 7th June 2018 , Accepted 12th September 2018

First published on 3rd October 2018


Abstract

NMR spectra of Mg2+ ions in ionic liquids were recorded using a highly sensitive variant of NMR spectroscopy known as β-NMR. The β-NMR spectra of MgCl2 in EMIM-Ac and EMIM-DCA compare favourably with conventional NMR, and exhibit linewidths of ∼3 ppm, allowing for discrimination of species with oxygen and nitrogen coordination.


Over the last few decades, β-detected nuclear magnetic resonance (β-NMR) has enjoyed use as a technique to study nuclear physics1,2 and condensed matter,3,4 the latter mainly focusing on solids, but also with applications to liquids.5–11 More recently, several groups have pursued advancing modern incarnations of the technique to study the solution chemistry of β-NMR probes,12–15 motivated by its potential to solve chemical and biochemical problems.16–19 With the advent of high-intensity radioactive ion beam (RIB) facilities, many new nuclei are available in quantities sufficient for the study of condensed matter. The main challenge to the use of RIBs as probes in liquids is the incompatibility of typical solutions with high vacuum usually applied to achieve optimal RIB transport. Though specialized setups12,13 aim to circumvent this limitation, room temperature ionic liquids are readily compatible with high vacuum, as they exhibit virtually zero vapour pressure.20 Here we demonstrate that Mg2+ chemical shift measurements in ionic liquids can be made with the precision necessary to discriminate between oxygen and nitrogen coordination environments using 31Mg β-NMR. This resolution, required for chemical and biochemical studies, indicates that the prospect for future use in these fields is promising, and may be extended to other metal ions.

Mg2+ was selected as a test case because 31Mg exhibits appropriate nuclear properties (I = ½, T½ = 236 ms, and asymmetric β-emission21,22) for β-NMR spectroscopy. Moreover, the closed shell Mg2+ is difficult to observe by most spectroscopic techniques, and it is an important metal ion in biological systems, involved in, e.g. ATP chemistry, nucleic acid folding, and as an integral component of chlorophyll.23 Finally, there is considerable and rising interest in characterization of Mg2+ chemistry in non-aqueous solutions, such as ionic liquids, because magnesium-based batteries may outperform lithium ion batteries.24 In this work, two ionic liquids were used as solvents, 1-ethyl-3-methylimidazolium acetate (EMIM-Ac), and 1-ethyl-3-methylimidazolium dicyanamide (EMIM-DCA), aiming to probe differential chemical shifts of Mg2+ complexes with typical oxygen and nitrogen containing ligands, see Fig. 1. NMR spectra recorded by 31Mg β-NMR and conventional 25Mg NMR are presented in Fig. 2. In these experiments, anhydrous MgCl2 was dissolved in the ionic liquids to a final concentration of 25 mM, to achieve comparable conditions in conventional NMR and β-NMR. The most important observation, and the key result of this work, is that the chemical shift difference for Mg2+ in the two ionic liquids is very similar using the two techniques. The similarity is expected if the experiments reflect the same Mg2+ species, and if these species exhibit rapid rotational diffusion. This establishes that 31Mg β-NMR can be applied to record chemical shifts of Mg2+ species in solution, and that the technique provides sufficient spectral resolution to discriminate between typical coordination environments of this metal ion. A priori, it was not obvious that the 40 keV beam of 31Mg+ implanted into the ionic liquids would rapidly (within the 236 ms half-life) be oxidized to Mg2+, and reach the same chemical equilibrium of complexation as observed in conventional 25Mg NMR, but the spectra in Fig. 2 indicate that this is indeed the case.


image file: c8dt02350f-f1.tif
Fig. 1 Chemical structure of the two ionic liquids used as solvents in this work. Left: 1-ethyl-3-methylimidazolium acetate (EMIM-Ac); right: 1-ethyl-3-methylimidazolium dicyanamide (EMIM-DCA).

image file: c8dt02350f-f2.tif
Fig. 2 31Mg β-NMR spectra (top) and conventional 25Mg NMR spectra (bottom) of 25 mM MgCl2 in EMIM-Ac (red) and EMIM-DCA (blue). The β-NMR experiments were conducted at 22 °C, ultra-high vacuum (<10−9 mbar), and an external magnetic field of 3.41 T (45.9 MHz) using crystalline MgO as chemical shift reference. The conventional NMR experiments were conducted at 70 °C to improve signal-to-noise, ambient pressure, and an external magnetic field of 11.7 T (500 MHz) using 11 M MgCl2 as chemical shift reference. See the text and ESI for details.

Another central observation is the high resolution of the 31Mg β-NMR spectrum exhibiting linewidths of ∼3 ppm (∼140 Hz at a resonance frequency of 45.9 MHz). This is analogous to the narrow lines observed in conventional solution (as opposed to solid state) NMR spectroscopy.25 Given the small span of chemical shifts typically observed in 25Mg NMR,26 high resolution is essential, to discriminate different Mg2+ complexes. For additional discussion of the origin of the linewidth of the β-NMR resonances, see the ESI.

The 31Mg β-NMR spectrum of Mg2+ in EMIM-DCA displays a major resonance at −60.2 ± 2.4 ppm with respect to MgO and a weaker signal at −52.0 ± 2.4 ppm, see Fig. 2. For EMIM-Ac the main resonance is observed at −31.9 ± 2.4 ppm, with a shoulder consistent with another unresolved resonance at −38.1 ± 2.4 ppm. Thus, co-existence of two Mg2+ complexes is observed in EMIM-DCA and possibly in EMIM-Ac. The fact that both resonances are observed indicates no or slow exchange between them, i.e. the lifetime of each species is longer than milliseconds. In the conventional 25Mg NMR spectra speciation is not resolved. This may either be due to the poorer signal-to-noise, faster dynamics at the higher temperature (70 °C), or nuclear quadrupole interactions (NQIs) giving rise to line broadening for 25Mg (I = 5/2), while NQIs do not contribute to the line width for 31Mg (I = 1/2), vide infra.

To demonstrate the extraordinary sensitivity of 31Mg β-NMR, additional experiments were conducted at ultra-trace concentrations, with only (∼108 cumulative) magnesium ions originating from the 31Mg+ RIB implanted into the sample, see Fig. 3. The spectra presented in Fig. 2 and 3 differ significantly, because 25 mM non-radioactive MgCl2 was added for the data displayed in Fig. 2, while none was added for the data displayed in Fig. 3. Thus, the speciation observed in Fig. 3 may originate from trace amounts of high affinity ligands (such as water) for Mg2+, while these impurity sites are essentially saturated in Fig. 2, shifting the spectrum towards pure DCA coordination. Given the current and energy of the 31Mg+ beam, the half-life (236 ms), the beam-sample geometry, and sample volume, the 31Mg concentration may be estimated to be 0.1 pM–1000 pM, depending on the diffusion coefficient of Mg2+ in the ionic liquid. The actual concentration of Mg2+ is presumably higher and dictated by the trace amounts that may be present in the ionic liquids (<3.5 μM, see Table S1). Interestingly, the two resonances observed with 25 mM Mg2+ in EMIM-DCA, Fig. 2, still appear to be present in the spectrum recorded with ∼pM 31Mg, Fig. 3, presumably reflecting the same two species. The minor species observed at −51.7 ppm however, becomes more prominent, and an additional (third) resonance is present above the noise level. The fact that the relative intensity of the two latter resonances increases when no MgCl2 is dissolved in the ionic liquid prior to implantation of 31Mg+ support the notion that they may originate from coordination by water or other dilute impurities with higher affinity for Mg2+ than DCA in the ionic liquid. Moreover, they are observed at chemical shifts closer to that of oxygen coordination (in EMIM-Ac), corroborating that they originate from water replacing DCA in the first coordination sphere. The resonance of 31Mg in EMIM-Ac with ultra-trace amounts of Mg2+ is slightly shifted and somewhat broadened as compared to the signal observed with 25 mM MgCl2, Fig. 2 and 3. This may also reflect coordination by water or other impurities, and the observation that only a minor change occurs, agrees with the expectation that replacing one oxygen coordinating ligand (acetate) with another (water) leads only to small changes in chemical shift.


image file: c8dt02350f-f3.tif
Fig. 3 31Mg β-NMR spectra recorded with ultra-trace amounts (∼pM) of 31Mg in EMIM-Ac (red) and EMIM-DCA (blue). The chemical shift reference is crystalline MgO. See the text for additional details.

Complementary experiments using 1H-, 13C-, 15N-, 25Mg-, and 35Cl-NMR spectroscopy, and FTIR spectroscopy, see Fig. S1–S9, and Tables S1–S6, indicate that the ionic liquid anions (acetate and DCA respectively) coordinate to Mg2+ in rapid exchange with free anions, that the DCA anions coordinate with the terminal nitrogen, and that trace amounts of impurities such as water in the ionic liquids may give rise to the observed speciation, see Fig. 2 and 3. These observations agree with the μs residence time of water molecules reported at 25 °C for the aqua ion of Mg2+,27 with crystal structures of Mg(DCA)2(s) and Mg(DCA)2(H2O)4(s),28,29 and with hard–soft-acid–base theory.30,31 To further substantiate the structural interpretation of the recorded 31Mg β-NMR data, we present calculated NMR properties for selected Mg2+ complexes in Table 1.

Table 1 Calculated spectroscopic properties of selected complexes. Calculations were done at the B3LYP/pcSseg-2//B3LYP/pc-2 level with an iefpcm solvent model and a static dielectric constant of 11.7 corresponding to EMIM-DCA.32–36 Isotropic shielding, σ, chemical shift δ = (σrefσ)/(1 − σref), arbitrarily with respect to [Mg(H2O)6]2+, and electric field gradient, Vzz and η. See the ESI for details
Complex (in solvent with ε = 11.7) σ/ppm δ/ppm 25Mga 31Mga V zz η
δ expt/ppm δ expt/ppm a.u.
a Experimental chemical shifts, assuming the indicated assignment (the absolute values of the chemical shifts differ because different chemical shift references were used).
[Mg(H2O)6]2+ 577 0 0.00
[Mg(DCA)2(H2O)4] 578 −1 0.18 0.04
[Mg(DCA)3(H2O)3] 581 −4 −0.12 0.87
[Mg(DCA)4(H2O)2]2− 585 −8 −43.2 −0.12 0.00
[Mg(DCA)5(H2O)]3− 589 −12 −52.0 −0.04 0.22
[Mg(DCA)6]4− 597 −20 −13.3 −60.2 0.00
[MgCl2(DCA)2(H2O)2]2− 573 4 0.16 0.66
[Mg(Ac)2(H2O)4] 570 7 −31.9 0.18 0.40
[Mg(Ac)4(H2O)2]2− 573 4 10.2 −38.1 −0.01 0.41


In agreement with the 31Mg β-NMR experimental data, nitrogen coordination (DCA) gives rise to shielding as compared to oxygen coordination. Accordingly, the consecutive replacement of DCA anions with water molecules gives rise to deshielding, and a tentative interpretation that accounts for the main observed β-NMR resonances for Mg2+ in EMIM-DCA is [Mg(DCA)6]4− (−60.2 ppm), [Mg(DCA)5(H2O)]3− (−52.0 ppm), and [Mg(DCA)4(H2O)2]2− (−43.2 ppm), see Fig. 2, 3 and Table 1. Assuming that this assignment is correct, it is possible to estimate the equilibrium constants in EMIM-DCA for:

image file: c8dt02350f-t1.tif

image file: c8dt02350f-t2.tif
consistently reproducing the data in Fig. 2 and 3. For Mg2+ in EMIM-Ac, only minor differences in chemical shift are observed experimentally, Fig. 2 and 3, in agreement with the small changes calculated when replacing acetate with water in the first coordination sphere, see Table 1. Thus, the resonances observed for EMIM-Ac may originate from [Mg(Ac)4(H2O)2]2− (−38.1 ppm), and [Mg(Ac)2(H2O)4] (−31.9 ppm) or other complexes with mixed acetate and water coordination.

Interestingly, the calculated electric field gradients, Vzz, are very small for the [Mg(H2O)6]2+, [Mg(Ac)4(H2O)2]2− and [Mg(DCA)6]4− complexes. For [Mg(H2O)6]2+ and [Mg(DCA)6]4−, this is expected due to the high symmetry of the complex, i.e. these complexes may essentially escape quadrupolar line broadening. The remaining complexes are likely to experience quadrupolar line broadening of the conventional 25Mg NMR resonances. Although there are other contributions to the line width, this suggests that only the [Mg(Ac)4(H2O)2]2− and [Mg(DCA)6]4− complexes are readily observable in 25Mg NMR, in agreement with the fact that only two resonances are observed. Moreover the observed chemical shift difference is ∼23 ppm, see Fig. 2, Table 1, and Table S5.

In summary, we have successfully recorded 31Mg β-NMR spectra for Mg2+ in solution. Contrary to 25Mg (I = 5/2), 31Mg (I = 1/2) NMR spectra do not suffer from line broadening due to nuclear quadrupole interactions, and the sensitivity of 31Mg β-NMR is far superior to that of 25Mg NMR illustrating that β-NMR allows for interrogation of elements that are otherwise difficult or impossible to observe by NMR spectroscopy. The excellent resolution of the 31Mg β-NMR spectra, allows for the discrimination of various complexes with the typical oxygen and nitrogen coordination formed by Mg2+, and consequently opens for further characterization of the coordination chemistry of this spectroscopically elusive metal ion in solution. Finally, the results indicate that similar success may be achieved with application of β-NMR to other metal ions. Even formation of transient species and surface properties of liquids may be explored, as the time from implantation to measurement is of the order of a second, and the implantation depth can be controlled by the energy of the ion beam if diffusion is slow. Therefore, broader applications of β-NMR in inorganic and bioinorganic chemistry may be anticipated.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

TRIUMF receives federal funding via a contribution agreement through the National Research Council of Canada and the Natural Sciences and Engineering Research Council of Canada (RGPIN-2018-04030). LH thanks The Danish Council for Independent Research | Natural Sciences, the Agency for Science, Technology and Innovation under the Ministry of Higher Education and Science, Denmark, for financial support. RMLM and AC acknowledge additional support from their NSERC CREATE IsoSiM fellowship. MD and DF acknowledge additional support from a QuEST fellowship. SP acknowledges the Greek State Scholarships Foundation (I. K. Y.), and Project scholarships I. K. Y. from resources of “Education and Lifelong Learning” of the European Social Fund (ESF) of NSFR 2007–2013.

Notes and references

  1. K. Asahi and K. Matsuta, Nucl. Phys. A, 2001, 693, 1 CrossRef.
  2. G. Neyens, Rep. Prog. Phys., 2003, 66(4), 633 CrossRef CAS.
  3. J. Christiansen and M. Ackermann, Hyperfine interactions of radioactive nuclei, Springer-Verlag, Berlin, New York, 1983 Search PubMed.
  4. W. A. MacFarlane, Solid State Nucl. Magn. Reson., 2015, 68–69, 1 CrossRef CAS PubMed.
  5. P. Heitjans, G. Kiese, H. Ackermann, B. Bader, W. Buttler, K. Dörr, F. Fujara, H. Grupp, A. Körblein and H.-J. Stöckmann, J. Phys. Colloq., 1980, 41(C8), C8–409 Search PubMed.
  6. Ionic liquids, molten salts, and polyelectrolytes: proceedings of the international conference held in Berlin (West), ed. K.-H. Bennemann, F. Brouers and D. Quitmann, 1982 Search PubMed.
  7. P. Heitjans, G. Kiese, C. Vandermarel, H. Ackermann, B. Bader, P. Freilander and H. J. Stockmann, Hyperfine Interact., 1983, 16(1–4), 569 CrossRef CAS.
  8. C. Vandermarel, P. Heitjans, H. Ackermann, B. Bader, P. Freilander, G. Kiese and H. J. Stockmann, J. Non-Cryst. Solids, 1984, 61–2(Jan), 213 CrossRef.
  9. C. Vandermarel, P. Heitjans, H. Ackermann, B. Bader, P. Freilander, A. Schirmer and H. J. Stockmann, Z. Phys. Chem., Neue Folge, 1988, 156, 629 CrossRef CAS.
  10. F. Hilmer, A. Schirmer, H. Ackermann, P. Heitjans and H. J. Stockmann, Z. Naturforsch., A: Phys. Sci., 1994, 49(1–2), 27 CAS.
  11. W. Geertsma and C. Vandermarel, J. Phys.: Condens. Matter, 1995, 7(47), 8867 CrossRef CAS.
  12. A. Gottberg, M. Stachura, M. Kowalska, M. L. Bissell, V. Arcisauskaite, K. Blaum, A. Helmke, K. Johnston, K. Kreim, F. H. Larsen, R. Neugart, G. Neyens, R. F. Garcia Ruiz, D. Szunyogh, P. W. Thulstrup, D. T. Yordanov and L. Hemmingsen, Chemphyschem, 2014, 15(18), 3929 CrossRef CAS PubMed.
  13. T. Sugihara, M. Mihara, J. Shimaya, K. Matsuta, M. Fukuda, J. Ohno, M. Tanaka, S. Yamaoka, K. Watanabe, S. Iwakiri, R. Yanagihara, Y. Tanaka, H. Du, K. Onishi, S. Kambayashi, T. Minamisono, D. Nishimura, T. Izumikawa, A. Ozawa, Y. Ishibashi, A. Kitagawa, S. Sato, M. Torikoshi and S. Momota, Hyperfine Interact., 2017, 238, 20 CrossRef.
  14. M. Stachura, R. M. L. McFadden, A. Chatzichristos, M. H. Dehn, A. Gottberg, L. Hemmingsen, A. Jancso, V. L. Karner, R. F. Kiefl, F. H. Larsen, J. Lassen, C. D. P. Levy, R. Li, W. A. MacFarlane, G. D. Morris, S. Pallada, M. R. Pearson, D. Szunyogh, P. W. Thulstrup and A. Voss, Hyperfine Interact., 2017, 238, 8 CrossRef.
  15. M. Kowalska, P. Aschenbrenner, M. Baranowski, M. L. Bissell, W. Gins, R. D. Harding, H. Heylen, G. Neyens, S. Pallada, N. Severijns, V. Ph, M. Walczak, F. Wienholtz, Z. Y. Xu, X. F. Yang and D. Zakoucky, J. Phys. G: Nucl. Part. Phys., 2017, 44(8), 084005 CrossRef.
  16. M. Stachura, A. Gottberg, M. Kowalska, K. Johnston and L. Hemmingsen, Nucl. Phys. News, 2015, 25(2), 25 CrossRef.
  17. R. F. G. Ruiz, M. L. Bissell, A. Gottberg, M. Stachura, L. Hemmingsen, G. Neyens and N. Severijns, EPJ Web Conf., 2015, 93, 07004 CrossRef.
  18. M. Stachura, A. Gottberg, K. Johnston, M. L. Bissell, R. F. G. Ruiz, J. M. Correia, A. R. G. Costa, M. Dehn, M. Deicher, A. Fenta, L. Hemmingsen, T. E. Molholt, M. Munch, G. Neyens, S. Pallada, M. R. Silva and D. Zakoucky, Nucl. Instrum. Methods B, 2016, 376, 369 CrossRef CAS.
  19. A. Jancso, J. G. Correia, A. Gottberg, J. Schell, M. Stachura, D. Szunyogh, S. Pallada, D. C. Lupascu, M. Kowalska and L. Hemmingsen, J. Phys. G: Nucl. Part. Phys., 2017, 44(6), 064003 CrossRef.
  20. E. F. Smith, I. J. Garcia, D. Briggs and P. Licence, Chem. Commun., 2005, 45, 5633 RSC.
  21. M. Kowalska, D. T. Yordanov, K. Blaum, P. Himpe, P. Lievens, S. Mallion, R. Neugart, G. Neyens and N. Vermeulen, Phys. Rev. C: Nucl. Phys., 2008, 77(3), 034307 CrossRef.
  22. G. Neyens, M. Kowalska, D. Yordanov, K. Blaum, P. Himpe, P. Lievens, S. Mallion, R. Neugart, N. Vermeulen, Y. Utsuno and T. Otsuka, Phys. Rev. Lett., 2005, 94(2), 022501 CrossRef CAS PubMed.
  23. I. G. Bertini, B. Harry, E. I. Stiefel and J. S. Valentine, Biological inorganic chemistry: structure and reactivity, University Science Books, Sausalito, Calif., 2007, pp. 739 Search PubMed.
  24. J. Muldoon, C. B. Bucur and T. Gregory, Chem. Rev., 2014, 114(23), 11683 CrossRef CAS PubMed.
  25. L. Simeral and G. E. Maciel, J. Phys. Chem., 1976, 80(5), 552 CrossRef CAS.
  26. J. C. C. Freitas and M. E. Smith, in Annual Reports on NMR Spectroscopy, ed. G. A. Webb, Academic Press, 2012, vol. 75, pp. 25–114 Search PubMed.
  27. L. Helm and A. E. Merbach, Coord. Chem. Rev., 1999, 187(1), 151 CrossRef CAS.
  28. B. Jurgens, E. Irran and W. Schnick, J. Solid State Chem., 2001, 157(2), 241 CrossRef CAS.
  29. E. Irran, B. Jurgens, S. Schmid and W. Schnick, Z. Anorg. Allg. Chem., 2005, 631(8), 1512 CrossRef CAS.
  30. R. G. Parr and R. G. Pearson, Absolute, J. Am. Chem. Soc., 1983, 105(26), 7512 CrossRef CAS.
  31. R. G. Pearson, J. Am. Chem. Soc., 1963, 85(22), 3533 CrossRef CAS.
  32. F. Jensen, J. Chem. Phys., 2001, 115(20), 9113 CrossRef CAS.
  33. F. Jensen, J. Phys. Chem. A, 2007, 111(44), 11198 CrossRef CAS PubMed.
  34. F. Jensen, J. Chem. Theory Comput., 2015, 11(1), 132 CrossRef CAS PubMed.
  35. J. Tomasi, B. Mennucci and E. Cances, J. Mol. Struct. (THEOCHEM), 1999, 464(1–3), 211 CrossRef CAS.
  36. M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone, T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega, G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, T. Vreven, K. Throssell, J. A. Montgomery Jr., J. E. Peralta, F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin, V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand, K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar, J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi, J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas, J. B. Foresman and D. J. Fox, Gaussian 16, Wallingford, CT, 2016 Search PubMed.

Footnotes

Electronic supplementary information (ESI) available. See DOI: 10.1039/c8dt02350f
These authors contributed equally to this work.

This journal is © The Royal Society of Chemistry 2018
Click here to see how this site uses Cookies. View our privacy policy here.