ORIGINAL RESEARCH article

Front. Chem., 13 November 2018

Sec. Inorganic Chemistry

Volume 6 - 2018 | https://doi.org/10.3389/fchem.2018.00552

Field-Induced Dysprosium Single-Molecule Magnet Based on a Redox-Active Fused 1,10-Phenanthroline-Tetrathiafulvalene-1,10-Phenanthroline Bridging Triad

  • Univ Rennes, CNRS, ISCR (Institut des Sciences Chimiques de Rennes)-UMR 6226, Rennes, France

Article metrics

View details

11

Citations

4k

Views

1,4k

Downloads

Abstract

Tetrathiafulvalene and 1,10-phenanthroline moieties present, respectively remarkable redox-active and complexation activities. In this work, we investigated the coordination reaction between the bis(1,10-phenanthro[5,6-b])tetrathiafulvalene triad (L) and the Dy(hfac)3·2H2O metallo precursor. The resulting {[Dy2(hfac)6(L)]·CH2Cl2·C6H14}3 (1) dinuclear complex showed a crystal structure in which the triad L bridged two terminal Dy(hfac)3 units and the supramolecular co-planar arrangement of the triads is driven by donor-acceptor interactions. The frequency dependence of the out-of-phase component of the magnetic susceptibility highlights three distinct maxima under a 2000 Oe static applied magnetic field, a sign that 1 displays a Single-Molecule Magnet (SMM) behavior with multiple magnetic relaxations. Ab initio calculations rationalized the Ising character of the magnetic anisotropy of the DyIII ions and showed that the main anisotropy axes are perpendicular to the co-planar arrangement of the triads. Single-crystal rotating magnetometry confirms the orientation of the main magnetic axis. Finally combining structural analysis and probability of magnetic relaxation pathways through Quantum Tunneling of the Magnetization (QTM) vs. excited states (Orbach), each DyIII center has been attributed to one of the three observed magnetic relaxation times. Such coordination compound can be considered as an ideal candidate to perform redox-magnetic switching.

Introduction

Single-Molecule Magnets (SMMs) are intensively studied for more than 25 years due to their potential applications in high density data storage, quantum computing, and spintronics (Leuenberger and Loss, 2001; Gatteschi et al., 2006; Lehmann et al., 2007; Bogani and Wernsdorfer, 2008; Mannini et al., 2009; Ganzhorn et al., 2013). Especially the interest for lanthanide ions has quickly grown during the last decade due to their high magnetic moment and strong magnetic anisotropy making them potential candidates for the elaboration of SMMs (Benelli and Gatteschi, 2002; Sessoli and Powell, 2009; Rinehart and Long, 2011; Woodruff et al., 2013; Liddle and van Slageren, 2015; Goodwin et al., 2017; Guo et al., 2017; Pointillart et al., 2017; Gupta and Murugavel, 2018). Since most lanthanide centers are luminescent (Piguet and Bünzli, 1999; Comby and Bünzli, 2007), the use of such element pushes up the interest of both chemist and physicist communities for luminescent SMMs with proper correlation between optical and magnetic properties (Cucinotta et al., 2012; Long et al., 2012; Ehama et al., 2013; Yamashita et al., 2013; Ren et al., 2014; Yi et al., 2014; Pointillart et al., 2015c). In a general manner, the interest is focused on multi-properties SMM (Long et al., 2015; Ou-Yang et al., 2016), and thus redox active SMMs could be designed using tetrathiafulvalene(TTF)-based ligands (da Cunha et al., 2013; Gao et al., 2013, 2014; Pointillart et al., 2013b; Soussi et al., 2015). Obviously the design of such molecular objects implies modifications of the TTF core with various acceptor groups (Sessoli and Powell, 2009; Pointillart et al., 2013a) in order to guaranty the coordination reaction with the lanthanide ions. Among the plethora of possible decorations of the TTF fragment, the 1,10-phenanthroline (phen) is an excellent choice to construct donor-acceptor (D-A) systems (Jia et al., 2007; Keniley et al., 2010, 2013; Qin et al., 2011). Their association with metallic precursors as led already to the observation of auspicious optical (Dupont et al., 2013; Lapadula et al., 2015) and magnetic (Pointillart et al., 2015b, 2016) properties. Even more promising are the acceptor-donor-acceptor (A-D-A) triads, such as bis-(quinoxaline)-TTF (Nishida et al., 2006), bis(tetracyanoquinodimethane)-TTF (Otón et al., 2011a), bis(naphthoquinone)-TTF (Otón et al., 2011b), bis(pyrazine)-TTF (Schuler et al., 2014), bis-(dipyrido[3,2-a:2′,3′-c]phenazine)-TTF (Jia et al., 2014), and bis(bistert-butyl-o-quinone)-TTF (Kuropatov et al., 2010) because the majority of them offers the possibility to bridge two metallic units which could be optically and magnetically active. In other words, A-D-A triads where A is a phen fragment and D a TTF moiety could play the role of a bridging unique redox-switch between two multi-properties units. Few triads involving the phen moiety are already available in the literature such as the bis-(dipyrido[3,2-a:2′,3′-c]phenazine)-TTF (Jia et al., 2014) and Bis(1,10-phenanthro[5,6-b])tetrathiafulvalene (Chen et al., 2016).

In these lines, the latter triad was selected to bridge two SMM units as it could be readily obtained in a few synthetic steps. Thus the coordination reaction between the Dy(hfac)3·2H2O units and the bis(1,10-phenanthro[5,6-b])tetrathiafulvalene triad (L) leads to the formation of a dinuclear complex of formula {[Dy2(hfac)6(L)]·CH2Cl2·C6H14}3 (1). The crystallographic structure highlights the role of bridge of the triad between the two Dy(hfac)3 units while the dynamic magnetic measurements show a field-induced multi-SMM behavior due to the six crystallographically independent DyIII centers. The magnetic properties are fully rationalized using SA-CASSCF/SI-SO calculations.

Materials and methods

General

The precursor Dy(hfac)3·2H2O (hfac = 1,1,1,5,5,5-hexafluoroacetylacetonate anion) (Richardson et al., 1968) and the 1,3-dithiole-2-thione[4,5-f][1,10-phenanthroline] (Qin et al., 2011) compound were synthesized following previously reported methods. All other reagents were purchased from Aldrich Co., Ltd. and used without further purification.

Single crystal of 1 was mounted on a APEXIII D8 VENTURE Bruker-AXS diffractometer for data collection (MoKα radiation source, λ = 0.71073 Å), from the Centre de Diffractométrie (CDIFX), Université de Rennes 1, France (Table 1). Structure was solved with a direct method using the SHELXT program (Sheldrick, 2015a) and refined with a full matrix least-squares method on F using the SHELXL-14/7 program (Sheldrick, 2015b). Complete crystal structure results as a CIF file including bond lengths, angles, and atomic coordinates are deposited as Supporting Information.

Table 1

Compounds{[Dy2(hfac)6(L)]·CH2Cl2·C6H14}3 (1)
FormulaC189H102Cl6Dy6F108N12O36S12
Fw/g.mol−16741.24
Crystal systemtriclinic
Space groupP-1 (N°2)
Cell parametersa = 20.507(8) Å b = 20.492(7) Å c = 33.957(13) Å α = 93.123(12)°β = 96.248(13)°γ = 119.224(10)°
Volume/Å312284.0(80)
Z2
T/K150 (2)
radiationMo Kα
2θ range/°4.10 ≤ 2θ ≤ 56.37
ρcalc/g.cm−31.823
μ/mm−12.117
Number of reflections375,353
Independent reflections53,657
Rint0.1276
Fo2 > 2σ(Fo)239,102
Number of variables2,896
GOOF1.129
R1, ωR20.1396, 0.3716

X-ray crystallographic data for 1 (CCDC 1826960).

The elemental analyses of the compounds were performed at the Centre Régional de Mesures Physiques de l'Ouest, Rennes. Cyclic voltammetry was carried out in dried and degassed CH2Cl2 solution, containing 0.1 M N(C4H9)4PF6 as supporting electrolyte. Voltammograms were recorded at 100 mVs−1 at a platinum disk electrode. The potentials were measured vs. a saturated calomel electrode (SCE). The dc magnetic susceptibility measurements were performed on solid polycrystalline samples with a Quantum Design MPMS-XL SQUID magnetometer between 2 and 300 K under an applied magnetic field of 0.2 kOe for temperatures of 2–20 K, 2 kOe between 20 and 80 K and 10 kOe above. These measurements were all corrected for the diamagnetic contribution as calculated with Pascal's constants. The ac magnetic susceptibility measurements were performed on both Quantum Design MPMS-XL SQUID and Quantum Design PPMS magnetometers.

Wavefunction-based calculations were carried out on the partially-optimized structure of a mononuclear model complex (vide infra) by using the SA-CASSCF/RASSI-SO approach, as implemented in the MOLCAS quantum chemistry package (versions 8.0; Aquilante et al., 2016). In this approach, the relativistic effects are treated in two steps on the basis of the Douglas–Kroll Hamiltonian. First, the scalar terms were included in the basis-set generation and were used to determine the spin-free wave functions and energies in the complete active space self-consistent field (CASSCF) method (Roos et al., 1980). Next, spin-orbit coupling was added within the restricted-active-space state-interaction (RASSI-SO) method, which uses the spin-free wave functions as basis states (Malmqvist and Roos, 1989; Malmqvist et al., 2002). The resulting wave functions and energies are used to compute the magnetic properties and g-tensors of the lowest states from the energy spectrum by using the pseudospin S = 1/2 formalism implemented in the SINGLE-ANISO routine (Chibotaru et al., 2008; Chibotaru and Ungur, 2012). Cholesky decomposition of the bielectronic integrals was employed to save disk space and speed-up the calculations (Aquilante et al., 2008). The active space in the CASSCF calculation consisted of the nine 4f electrons of the DyIII ion spanning the seven 4f orbitals, i.e., CAS(9,7)SCF. State-averaged CASSCF calculations were performed for all of the sextets (21 roots), all of the quadruplets (224 roots), and 300 out of the 490 doublets (due to software limitations) of the Dy(III) ion. Twenty one sextets, 128 quadruplets, and 107 doublets were mixed through spin–orbit coupling in RASSI-SO. All atoms were described by ANO-RCC basis sets (Roos et al., 2004, 2005, 2008). The following contractions were used: [8s7p4d3f2g1h] for Dy, [4s3p2d1f] for the O and N atoms [3s2p1d] for the C and F atoms, [4s3p1d] for the S atoms and [2s] for the H atoms. DFT geometry optimization has been performed on model complexes for the 6 asymmetric Dy centers. The atomic positions of the molecule were extracted from the X-ray crystal structure. For each Dy center, the corresponding dimer unit has been cut in half at the bridging C-C double bond, the corresponding truncated C atom being replaced by a H atom (see Figure S10). In the optimization, the Dy(III) center was replaced by Y(III). Only the H and F positions were optimized while all the other positions were kept frozen. The calculations were carried using the PBE0 (Perdew et al., 1996; Adamo and Barone, 1999) hybrid functional implemented in the Gaussian 09 (revision D.01) package (Frisch et al., 2013). The “Stuttgart/Dresden” basis sets and effective core potentials were used to describe the yttrium atom (Dolg et al., 1993) while the other atoms were described with SVP basis sets (Weigend and Ahlrichs, 2005).

Synthesis

Bis(1,10-phenanthro[5,6-b])tetrathiafulvalene triad (L)

L was synthesized by modifying the published method (Chen et al., 2016). The 1,3-dithiole-2-thione[4,5-f][1,10-phenanthroline] was used instead of the 1,3-dithiole-2-one[4,5-f][1,10-phenanthroline] compound in order to cancel the oxidation step of the thione in one derivative. Yield: 45%.

{[Dy2(hfac)6(L)]·CH2Cl2·C6H14}3 (1)

Eleven milligram of L (0.022 mmol) were added to 20 ml of 1,2-dichloroethane. The suspension was heated to reflux and then a solution of 10 mL of 1,2-dichloroethane containing 35 mg of Dy(hfac)3·2H2O (0.043 mmol) was added. After 6 h of reflux, the 1,2-dichloroethane was eliminated under vacuum and the residue was dissolved in CH2Cl2. n-hexane was layered on the CH2Cl2 solution of 1 in the dark to give red single crystals suitable for X-ray diffraction study. I.R. bands (KBr): 2965, 1651, 1564, 1538, 1499, 1465, 1260, 1216, 1147, 1100, 806, 662, and 588 cm−1. Anal. Calcd (%) for C189H102Cl6Dy6F108N12O36S12 (1): C 33.64, H 1.51, N 2.49; found: C 33.71, H 1.66 N, 2.37.

Results and discussion

Structural analysis

The X-ray crystallographic data for compound 1 are given in Table 1. It crystallizes in the P-1 (N°2) triclinic space group. The asymmetric unit is composed by two dinuclear complexes of formula [Dy2(hfac)6(L)], two half dinuclear complexes, three dichloromethane and three n-hexane molecules of crystallization. In other words, six crystallographically independent DyIII centers are present in the molecular structure (Figure S1). All the DyIII centers are linked to three hfac anions and one L ligand giving an N2O6 surrounding with a D4d symmetry (square antiprism) for Dy1, D2d symmetry for Dy2-Dy5 and an intermediate symmetry between D2d and C2v for Dy6. The distortion is visualized by continuous shape measures performed with SHAPE 2.1 (Table S1; Llunell et al., 2013). The two 1,10-phenanthroline coordination sites are occupied by Dy(hfac)3 units and L plays the role of bridge between two Dy(hfac)3 units (Figure 1).

Figure 1

As we expect for such large and extended aromatic ligands, L is planar offering optimal stacks in the crystals along the c unit cell direction. The molecular planes in the stack are coplanar (top part of Figure 2) but they are slightly shifted (or/and rotated of an angle ranging from 83 to 102°) relative to each other to guaranty efficient π-π interactions despite of the steric hindrance of the Dy(hfac)3 units (bottom part of Figure 2). No bent structures are obtained and the charge sensitive central C = C bond is found to be in average 1.335 Å, corresponding to neutral form of the TTF ligand part. Complementary interactions such as S···N or S···S might also play a role in the crystal packing. The dithiol sulfur atoms from the upper molecule are localized above the nitrogen atoms of the phenanthroline moiety of the lower one. The S···N distances range from 3.681 to 3.812 Å which is close to the sum of the van der Waals radii (3.550 Å). Few S···S contacts are also identified with distances ranging from 3.871 to 4.141 Å. The minimal interplanar distance in the stack is 3.301 Å, suggesting the driving forces for the formation of stacks could be these acceptor-donor intermolecular interactions (bottom part of Figure 2).

Figure 2

To investigate sample homogeneity, PXRD was performed on crystalline powder of 1 (see Figure S2 for comparison between simulated and experimental PXRD pattern). Unfortunately, the powder showed a drastic loss of crystallinity character, probably due to solvents evaporation during the measurement (vide infra).

Electrochemical properties

In contrast with the poor solubility of L in most common organic solvents, complex 1 is well soluble in CH2Cl2 and its redox properties could be investigated by cyclic voltammetry (Figure S3). The cyclic voltammogram shows two mono-electronic oxidations at 0.87 V and 1.22 V for the first and second oxidations, corresponding to the formation of a radical cation and a dication TTF fragment, respectively (Table 2).

Table 2

E/VE/V
OxERedEOxERedE
10.9440.7901.3101.120

Oxidation potentials (V vs SCE, nBu4NPF6, 0.1 M in CH2Cl2 at 100 mV.s−1) of the ligand L and complex [Dy(hfac)3(L)]·CH2Cl2.

These values are found very close to those reported by Zuo et al (Chen et al., 2016). The reversibility of the electrochemical properties is a clear evidence of the complex stability, meaning that almost no de-coordination seems to occur once L is converted to its radical cationic form, under these conditions. In addition, such results attest the reversibility of the oxidation potentials and the redox-activity of ligand (L) after complexation.

Magnetic properties

The static magnetic properties were determined by measuring the thermal dependence of the magnetic susceptibility (χM) between 2 and 300 K (Figure 3A). The room temperature of the χMT product is equal to 27.86 cm3 K mol−1 i.e., an average value of 13.93 cm3 K mol−1 per Dy(III) center. Such value is in agreement with the expected value 14.17 cm3 K mol−1 for a Dy(III) free ion (6H15/2, g = 4/3; Kahn, 1993). The χMT(T) curve monotonically decreases due to the depopulation of the MJ states until 2 K where the χMT product reaches the value of 22.50 cm3 K mol−1. At 2 K, the field dependence of the magnetization shows a classic behavior with a value of 10.54 Nβ (31.61 Nβ for the six DyIII, Figure S3) at 50 kOe which is close to the expected value of 10 Nβ for an Ising ground state of the dinuclear complex.

Figure 3

The dynamic magnetic properties were probed by measuring the frequency dependence of the magnetic susceptibility. Such measurements only show an out-of-phase component of the magnetic susceptibility at high frequency and thus no maxima were observed in the experimental window at zero applied field. Consequently a scan field of the magnetic susceptibility at 2 K between 0 and 2,000 Oe was performed (Figure S4) in order to determine the optimal field to cancel the zero-field Quantum Tunneling of the Magnetization (QTM; Gatteschi and Sessoli, 2003; Ishikawa et al., 2005a) and both dipolar and hyperfine interactions (Ishikawa et al., 2005b; Guo et al., 2011; Pointillart et al., 2015a). Under an applied magnetic field, few out-of phase magnetic susceptibility maxima shift in the experimental window and could be attributed to the crystallographically different Dy(III) centers. At the optimal DC field of 2,000 Oe, three maxima were observed at 2 K (3 Hz, 315 Hz and 5,000 Hz, Figure 3B) whereas six different DyIII ions have been identified in the crystal structure. It is worth noticing that the three sets of data have almost the same intensity meaning that the same number of Dy(III) are involved in each of them. The frequency dependence of the magnetic susceptibility under 2,000 Oe field and at various temperatures (Figure 3B and Figures S5–S7) confirmed the multiple relaxation behavior of 1. Such a multi-peak magnetic response renders the fitting of the Arrhenius plots with an extended Debye model quite problematic (Figure S8). The dependence of the relaxation time (τ) with the temperature was represented considering the maxima of the out-of-phase susceptibility curves at each temperature. We manually selected the frequencies maxima to plot log(τ) vs. T (Figure 3C). The Arrhenius data was fully extracted for the intermediate relaxation time (blue plots, called τ1) while the Arrhenius plot was extracted at low temperature for the metal center(s) which display(s) the slowest magnetic relaxation time (green plots, called τ2) and finally no Arrhenius was extracted for the metal center(s) which display(s) the fastest magnetic relaxation time because the maxima shift at too high frequency for T >2 K. The Arrhenius plot deviates from the linearity at low temperature suggesting the coexistence of more than one relaxation process (Guo et al., 2010; Zhang et al., 2015, 2016; Lu et al., 2017; Zhu et al., 2019). The Arrhenius plot was fitted using a combination of Orbach relaxation (exp(-Δ/T) + ) (Cole and Cole, 1941), which provided an energy barrier of Δ = 48(2) cm−1 and τ(1)0 = 2.37(10) × 10−7 s and remaining QTM. A tentative fit of the Arrhenius plot for τ2 gave Δ = 32(3) cm−1 and τ(2)0 = 3.12(46) × 10−7 s and remaining QTM. Nevertheless, the existence of Raman (Orbach, 1961a,b; Scott and Jeffries, 1962; Abragam and Bleaney, 1970; Shrivastava, 1983), and/or Direct (Orbach, 1961a,b; Scott and Jeffries, 1962; Abragam and Bleaney, 1970; Shrivastava, 1983) processes cannot be excluded.

Rotating single-crystal magnetometry allows to determine experimentally the average susceptibility tensor. The angular dependence of the magnetization was measured at 2 K in three orthogonal planes of an oriented single crystal. The molar magnetic susceptibility was then fitted with

where α and β are the directions X, Y, and Z (Figure S9) in a cyclic permutation and θ is the angle between H and α. In the effective spin 1/2 formalism, the largest principal value of the g- tensor is equal to 16.53, close to the tabulated value (20.00) for a purely axial magnetic moment. In order to rationalize the experimental observations, ab-initio calculations (SA-CASSCF/SI-SO) were performed (see Computational details, Figure S11). The crystal field splitting, g-tensor components and wavefunction composition for each Kramers's Doublet of the ground-state multiplet (6H15/2) for the six Dy(III) ions were determined (Table 3 and Tables S2–S7).

Table 3

Energy (cm−1)gXgYgZWavefunction composition*
Dy1GD00.080.1518.990.86|±15/2> +0.11|±11/2>
1st ES96.41.192.5213.820.56|±13/2> +0.24|±9/2> +0.09|±5/2>
Dy2GD00.010.0119.460.93|±15/2> +0.07|±11/2>
1st ES146.00.170.2515.640.75|±13/2> +0.19|±9/2>
Dy3GD00.010.0219.330.91|±15/2> +0.08|±11/2>
1st ES139.00.160.2614.920.60|±13/2> +0.27|±9/2> +0.08|±5/2>
Dy4GD00.020.0519.230.89|±15/2> +0.10|±11/2>
1st ES111.10.280.5414.600.60|±13/2> +0.26|±9/2> +0.07|±5/2>
Dy5GD00.010.0219.270.90|±15/2> +0.09|±11/2>
1st ES133.40.410.8314.160.58|±13/2> +0.27|±9/2> +0.13|±5/2>
Dy6GD00.010.0119.400.92|±15/2> +0.07|±11/2>
1st ES128.40.190.3815.860.77|±13/2> +0.15|±9/2> +0.05|±11/2>

Computed energies, g-tensor components and wavefunction composition for the ground doublet state (GD) and 1st Excited State (ES) of the ground-state multiplet for the six Dy(III) centers in 1.

*Only the contributions ≥ 5% are given.

The calculated crystal field splittings of the ground multiplet state lead to a calculated χMT product, averaged over the six positions, which is perfectly in agreement with the experimental data (Figure S10). The composition of the ground doublet reproduces finely the experimental magnetization at 2 K. As expected, an Ising type anisotropy with the largest g value ranging from 18.99 to 19.46 (Table 3) was found with the orientation of the main magnetic axis perpendicular to the plane formed by the nitrogen atoms (Figure 4) and the metallic center. In other words, the anisotropy axis direction is positioned along the most charged direction of the coordination polyhedral as already observed in several N2O6 environments (da Cunha et al., 2013; Pointillart et al., 2015b; Ou-Yang et al., 2016; Speed et al., 2017). These are in agreement with the experimental data for the orientation of the anisotropy axes, in which the g- value is slightly smaller due to the contribution of all the Dy atoms forming the unit cell. The ground doublet state is mainly composed of the MJ = ±15/2 doublet with small MJ = ±11/2 contributions suggesting the possibility to have magnetic relaxation through QTM. The calculated energy gaps between the ground doublet state and the first excited doublet state range from 96.4 to 146 cm−1. Thus the calculations allow to rationalize the SMM behavior (at least under an applied magnetic field) but the activated energy barrier is overestimated probably because the intermolecular and hyperfine interactions are not included in the calculations (Vignesh et al., 2017). Of utmost importance, (i) the coplanar arrangement of the TTF-based ligands is driven by the supramolecular stack through both π-π and donor-acceptor interactions, (ii) previous works on oblate Dy(III) ion in N2O6 environment demonstrated the perpendicular orientation of the main anisotropy axis compared to the plan formed by the two nitrogen atoms and thus permit to predict to some extent the parallel direction of all the anisotropy axes.

Figure 4

In order to rationalize the experimental ac magnetic measurements and to identify the three families of Dy centers, the transition matrix elements of the magnetic moment were computed. The strength of these matrix elements allows to sketch the magnetic relaxation pathway for each Dy ion (Figure 5). It appears that the two Dy1 and Dy4 centers displayed the most efficient QTM and might correspond to the maxima at 5,000 Hz at 2 K. Dy2 and Dy6 are the ions with the most calculated Ising type magnetic anisotropy and the magnetic relaxation may involve the first excited state. They should display the slowest relaxation time with maxima of the out-of-phase component of the magnetic susceptibility centered at 3 Hz. Finally, the intermediate maxima at 315 Hz could be attributed to the Dy3 and Dy5 ions.

Figure 5

Conclusions

In summary, the redox-active Bis(1,10-phenanthro[5,6-b])tetrathiafulvalene triad (L) was used to bridge two Dy(hfac)3 units giving the dinuclear complex of formula [Dy2(hfac)6(L)]·CH2Cl2·C6H14. The supramolecular organization of the dinuclear species in the crystal driven by the donor-acceptor interactions leads to co-planar stacks of the triads. The coordinated ligands lead to a N2O6 environment with a crystal field splitting suitable for the observation of slow magnetic relaxation under an applied magnetic field. In a crystallographic point of view six independent Dy(III) ions have been identified resulting in three families of metal centers displaying different relaxation time of their magnetization in the frequency windows of 1–10 kHz at 2 K. Experimental and the ab-initio study confirmed the Ising nature of the magnetic anisotropy for the Dy(III) ions with the main anisotropy axes perpendicular to the co-planar stacks of the triads. Finally, based on structural analysis, composition of the ground doublet state and probability of QTM, a tentative attribution of the Dy(III) ions to the three different families of field-induced mononuclear SMM was achieved, as out-of-phase signal of the magnetic susceptibility for Dy2/Dy6, Dy3/Dy5, Dy1/Dy4, respectively centered at 3, 315, and 5,000 Hz. The redox-activity of the triad is promising in the perspective to use it as a switch between the mononuclear SMM.

Statements

Author contributions

BL and FP made the synthesis; OG and OC performed the magnetic measurements and interpreted them; JF performed single crystal rotating magnetometry; BLG and VM performed the ab initio calculations. All authors participated in the writing process of the manuscript. All authors read and approved the final version of the manuscript.

Funding

This work was supported by Région Bretagne, Rennes Métropole, CNRS, Université de Rennes 1, FEDER and the European Research Council (ERC, MULTIPROSMM project N°725184).

Acknowledgments

BLG and VM thank the French GENCI/IDRIS-CINES center for high-performance computing resources.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2018.00552/full#supplementary-material

References

  • 1

    AbragamA.BleaneyB. (1970). Electron Paramagnetic Resonance of Transition Ions.Oxford: Clarendon Press.

  • 2

    AdamoC.BaroneV. (1999). Toward reliable density functional methods without adjustable parameters: the PBE0 model. J. Chem. Phys.110, 61586170. 10.1063/1.478522

  • 3

    AquilanteF.AutschbachJ.CarlsonR. K.ChibotaruL. F.DelceyM. G.De VicoL.et al. (2016). MOLCAS 8: new capabilities for multiconfigurational quantum chemical calculations across the periodic table. J. Comput. Chem.37, 506541. 10.1002/jcc.24221

  • 4

    AquilanteF.MalmqvistP.-A.PedersenT.-B.GhoshA.RoosB. O. (2008). Cholesky decomposition-based multiconfiguration second-order perturbation theory (CD-CASPT2): application to the spin-state energetics of CoIII(diiminato)(NPh). J. Chem. Theory Comput.4, 694702. 10.1021/ct700263h

  • 5

    BenelliC.GatteschiD. (2002). Magnetism of lanthanides in molecular materials with transition-metal ions and organic radicals. Chem. Rev.102, 23692388. 10.1021/cr010303r

  • 6

    BoganiL.WernsdorferW. (2008). Molecular spintronics using single-molecule magnets. Nat. Mater.7, 179186. 10.1038/nmat2133

  • 7

    ChenB.LvZ.-P.HuaC. F.LeongC.TunaF. M.D'AlessandroD.et al. (2016). Dinuclear ruthenium complex based on a π-extended bridging ligand with redox-active tetrathiafulvalene and 1,10-phenanthroline units. Inorg. Chem.55, 46064615. 10.1021/acs.inorgchem.6b00437

  • 8

    ChibotaruL. F.UngurL. (2012). Ab initio calculation of anisotropic magnetic properties of complexes. I. unique definition of pseudospin Hamiltonians and their derivation. J. Chem. Phys.137, 064112064122. 10.1016/1.4739763

  • 9

    ChibotaruL. F.UngurL.SonciniA. (2008). The origin of nonmagnetic kramers doublets in the ground state of dysprosium triangles: evidence for a toroidal magnetic moment. Angew. Chem. Int. Ed.47, 41264129. 10.1002/anie.200800283

  • 10

    ColeK. S.ColeR. H. (1941). Dispersion and absorption in dielectrics I. alternating current characteristics. J. Chem. Phys.9, 341351. 10.1063/1.1750906

  • 11

    CombyS.BünzliJ.-C. G. (2007). Handbook on the Physics and Chemistry of Rare Earths. Amsterdam: Elsevier BV.

  • 12

    CucinottaG.PerfettiM.LuzonJ.EtienneM.CarP. E.CaneschiA.et al. (2012). Magnetic anisotropy in a dysprosium/DOTA single-molecule magnet: beyond simple magneto-structural correlations. Angew. Chem. Int. Ed.51, 16061610. 10.1002/anie.201107453

  • 13

    da CunhaT. T.JungJ.BoulonM. E.CampoG.PointillartF.PereiraC. L. M.et al. (2013). Magnetic poles determinations and robustness of memory effect upon solubilization in a DyIII-based single ion magnet. J. Am. Chem. Soc.135, 1633216335. 10.1021/ja4089956

  • 14

    DolgM.StollH.PreussH. (1993). A combination of quasirelativistic pseudopotential and ligand field calculations for lanthanoid compounds. Theor. Chim. Acta85, 441450.

  • 15

    DupontN.RanY.-F.LiuS.-X.GriljJ.VautheyE.DecurtinsS.et al. (2013). A donor-acceptor tetrathiafulvalene ligand complexed to iron(II): synthesis, electrochemistry, and spectroscopy of [Fe(phen)3-TTF-dppz)](PF6)2. Inorg. Chem.52, 306312. 10.1021/ic3019277

  • 16

    EhamaK.OhmichiY.SakamotoS.FujinamiT.MatsumotoN.MochidaN.et al. (2013). Synthesis, structure, luminescent, and magnetic properties of carbonato-bridged ZnLn complexes [(μ4-CO3)2{ZnIILnLnIII(NO3)}2] (LnIII = GdIII, TbIII, DyIII; L1 = N,N′-Bis(3-methoxy-2-oxybenzylidene)-1,3-propanediaminato, L2 = N,N′-Bis(3-ethoxy-2-oxybenzylidene)-1,3-propanediaminato). Inorg. Chem.52, 1282812841. 10.1021/ic4022273

  • 17

    FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.et al. (2013). Gaussian 09, Revision D.01.Wallingford, CT: Gaussian Inc.

  • 18

    GanzhornM.KlyatskayaS.RubenM.WernsdorferW. (2013). Strong spin-phonon coupling between a single-molecule magnet and a carbon nanotube nanoelectromechanical system. Nat. Nanotechnol.8, 165169. 10.1038/nnano.2012.258

  • 19

    GaoF.CuiL.LiuW.HuL.ZhongY. W.LiY. Z.et al. (2013). Seven-coordinate lanthanide sandwich-type complexes with a tetrathiafulvalene-fused schiff base ligand. Inorg. Chem.52, 1116411172. 10.1021/ic401421h

  • 20

    GaoF.ZhangX. M.CuiL.DengK.ZengQ. D.ZuoJ. L. (2014). Tetrathiafulvalene-supported triple-decker phthalocyaninato dysprosium(III). complex: synthesis, properties and surface assembly. Sci. Rep.4, 59285935. 10.1038/srep05928

  • 21

    GatteschiD.SessoliR. (2003)Quantum tunneling of magnetization related phenomena in molecular materials. Angew. Chemie Int. Ed.42, 268297. 10.1002/anie.200390099

  • 22

    GatteschiD.SessoliR.VillainJ. (2006). Molecular Nanomagnets.New York, NY: Oxford University Press.

  • 23

    GoodwinC. A. P.OrtuF.RetaD.ChiltonN. F.MillsD. P. (2017). Molecular magnetic hysteresis at 60 kevlin in dysprosocenium. Nature548, 439442. 10.1038/nature23447

  • 24

    GuoF. S.Day benjaminM.ChenY. C.TongM. L.MansikkamäkiA.LayfieldR. A. (2017). A dysprosium metallocene single-molecule magnet functioning at the axial limit. Angew. Chem. Int. Ed.56, 1144511449. 10.1002/anie.201705426

  • 25

    GuoY. N.XuG. F.GamezP.ZhaoL.LinS. Y.DengR.et al. (2010). Two-step relaxation in a linear tetranuclear dysprosium(III). aggregate showing single-molecule magnet behavior. J. Am. Chem. Soc.132, 85388539. 10.1021/ja103018m

  • 26

    GuoY. N.XuG. F.WernsdorferW.UngurL.GuoY.TangJ.et al. (2011). Strong axiality and ising exchange interaction suppress zero-field tunneling of magnetization of an asymmetric Dy2 single-molecule magnet. J. Am. Chem. Soc.133, 1194811951. 10.1021/ja205035g

  • 27

    GuptaS. K.MurugavelR. (2018). Enricing Lanthanide single-ion magnetism through symmetry and axiality. Chem. Commun.54, 36853696. 10.1039/C7CC09956H

  • 28

    IshikawaN.SugitaM.WernsdorferW. (2005a). Quantum tunneling of magnetization in lanthanide single-molecule magnets: Bis(phthalocyaninato)terbium and Bis(phthalocyaninato)dysprosium Anions. Angew. Chem. Int. Ed.44, 29312935. 10.1002/anie.200462638

  • 29

    IshikawaN.SugitaM.WernsdorferW. (2005b). Nuclear spin driven quantum tunneling of magnetization in a new lanthanide single-molecule magnet: Bis(Phthalocyaninato)holmium anion. J. Am. Chem. Soc.127, 36503651. 10.1021/ja0428661

  • 30

    JiaC.LiuS. X.TannerC.LeiggenerC.NeelsA.SanguinetL.et al. (2007). An experimental and computational study on intramolecular charge transfer: a tetrathiafulvalene-fused dipyridophenazine molecule. Chem. Eur. J.13, 38043812. 10.1002/chem.200601561

  • 31

    JiaH.DingJ.HauserA.DecurtinsS.LiuS.-X. (2014). Large π-conjugated chromophores derived from tetrathiafulvalene. Asian J. Org. Chem.3, 198202. 10.1002/ajoc.201300144

  • 32

    KahnO. (1993). Molecular Magnetism.Weinhem: VCH.

  • 33

    KenileyL.JrDupontN.RayL.DingJ.KovnirK.HoytM. J.et al. (2013). Complexes with redox-active ligands: synthesis, structure, and electrochemical and photophysical behavior of the Ru(II). compex with TTF-annulated phenanthroline. Inorg. Chem.52, 80408052. 10.1021/ic4006949

  • 34

    KenileyL.Jr.RayL.KovnirK. A.ShatrukM. (2010). TTF-annulated phenanthroline and unexpected oxidation cleavage of the C = C bond in its ruthenium(II). complex. Inorg. Chem.49, 13071309. 10.1021/ic902230f

  • 35

    KuropatovV.KlementievaS.FukinG.MitinA.KetlovS.BudnikovaY.et al. (2010). Novel method for the synthesis of functionalized tetrathiafulvalenes, an acceptor-donor-acceptor molecule comprising of two o-quinone moieties linked by a TTF bridge. Tetrahedron66, 76057611. 10.1016/j.tet.2010.07.038

  • 36

    LapadulaG.TrummerD.ConleyM. P.SteinmannM.RanY.-F.BrasseletS.et al. (2015). One-photon near-infrared sensitization of well-defined Yb(III). surface complexes for NIR-to-NIR single nanoparticle imaging. Chem. Mater.27, 20332039. 10.1021/acs.chemmater.5b00306

  • 37

    LehmannJ.Gaita-ArinoA.CoronadoE.LossD. (2007). Spin qubits with electrically gated polyoxometalate molecules. Nat. Nanotechnol.2, 312317. 10.1038/nnano.2007.110

  • 38

    LeuenbergerM. N.LossD. (2001). Quantum computing in molecular magnets. Nature410, 789793. 10.1038/35071024

  • 39

    LiddleS. T.van SlagerenJ. (2015). Improving f-element single molecule magnets. Chem. Soc. Rev.44, 66556669. 10.1039/C5CS00222B

  • 40

    LlunellM.CasanovaD.CireraJ.BofillJ. M.AlemanyP.AlvarezS. (2013). SHAPE (v. 2.1).Barcelona.

  • 41

    LongJ.RouquetteJ.ThibaudJ. M.FerreiraR. A. S.CarlosL. D.DonnadieuB.et al. (2015). A high-temperature molecular ferroelectric Zn/Dy complex exhibiting single-ion-magnet behavior and lanthanide luminescence. Angew. Chem. Int. Ed.54, 22362240. 10.1002/anie.201410523

  • 42

    LongJ.VallatR.FerreiraR. A.CarlosL. D.PazF. A. A.GuariY.et al. (2012). A bifunctional luminescent single-ion magnet: towards correlation between luminescence studies and magnetic slow relaxation processes. Chem. Commun.48, 99749976. 10.1039/C2CC35321K

  • 43

    LuJ.GuoM.TangJ. (2017). Recent developments in lanthanide single-molecule magnets. Chem. Asian. J.12, 27722779. 10.1002/asia.201701032

  • 44

    MalmqvistP. A.RoosB. O. (1989). The CASSCF state interaction method. Chem. Phys. Lett.155, 189194. 10.1016/0009-2614(89)85347-3

  • 45

    MalmqvistP. A.RoosB. O.SchimmelpfennigB. (2002). The restricted active space (RAS). state interaction approach with spin-orbit coupling. Chem. Phys. Lett.357, 230240. 10.1016/S0009-2614(02)00498-0

  • 46

    ManniniM.PineiderF.SainctavitP.DanieliC.OteroE.SciancaleporeC.et al. (2009). Magnetic memory of a single-molecule quantum magnet wired to a gold surface. Nat. Mater.8, 194197. 10.1038/nmat2374

  • 47

    NishidaJ. I.KumakiD.TokitoS.YamashitaY. (2006). High performance n- and p-type field-effect transitors based on tetrathiafulvalene derivatives. J. Am. Chem. Soc.128, 95989599. 10.1021/ja0630083

  • 48

    OrbachR. (1961a). Spin-lattice relaxation in rare-earth salts. Proc. R. Soc. Lond. A Math. Phys. Eng. Sci.264, 458484. 10.1098/rspa.1961.0211

  • 49

    OrbachR. (1961b). On the theory of spin-lattice relaxation in paramagnetic salts. Proc. Phys. Soc.77, 821826.

  • 50

    OtónF.LloverasV.Mas-TorrentM.Vidal-GancedoJ.VecianaJ.RoviraC. (2011a). Coupling tetracyanoquinodimethane to tetrathiafulvalene: a fused TCNQ-TTF-TCNQ triad. Angew. Chem. Int. Ed.50, 1090210906. 10.1002/anie.201104841

  • 51

    OtónF.PfattnerR.PavlicaE.OlivierY.BratinaG.CornilJ.et al. (2011b). Electronic and structural characterization of a tetrathiafulvalene compound as a potential candidate for ambipolar transport properties. CrystEngComm13, 65976600. 10.1039/C1CE05559C

  • 52

    Ou-YangJ. K.SalehN.Fernandez GarciaG.NorelL.PointillartF.GuizouarnT.et al. (2016). Improved slow magnetic relaxation in optically pure helicene-based DyIII single molecule magnets. Chem. Commun.52, 1447414477. 10.1039/C6CC08638A

  • 53

    PerdewJ. P.BurkeK.ErnzerhofM. (1996). Generalized gradient approximation made simple. Phys. Rev. Lett.77, 38653868. 10.1103/physrevlett.77.3865

  • 54

    PiguetC.BünzliJ.-C. G. (1999). Mono- and polymetallic lanthanide-containing functional assemblies: a field between tradition and novelty. Chem. Soc. Rev.28, 347358. 10.1039/A804240C

  • 55

    PointillartF.BernotK.GolhenS.Le GuennicB.GuizouarnT.OuahabL.et al. (2015a). Magnetic memory in an isotopically enriched and magnetically isolated manonuclear dysprosium complex. Angew. Chem. Int. Ed.54, 15041507. 10.1002/anie.201409887

  • 56

    PointillartF.CadorO.Le GuennicB.OuahabL. (2017). Uncommon Lanthanide ions in purely 4f single molecule magnets. Coord. Chem. Rev.346, 150175. 10.1016/j.ccr.2016.12.017

  • 57

    PointillartF.GolhenS.CadorO.OuahabL. (2013a). Paramagnetic 3d coordination complexes involving redox-active tetrathiafulvalene derivatives: an efficient approach to elaborate multi-properties materials. Dalton Trans.42, 19491960. 10.1039/c2dt32150e

  • 58

    PointillartF.JungJ.Berraud-PacheR.Le GuennicB.DorcetV.GolhenS.et al. (2015b). L. luminescence and single-molecule magnet behavior in lanthanide complexes involving a tetrathiafulvalene-fused dipyridophenazine ligand. Inorg. Chem.54, 53845397. 10.1021/acs.inorgchem.5b00441

  • 59

    PointillartF.Le GuennicB.CadorO.MauryO.OuahabL. (2015c). Lanthanide ion and tetrathiafulvalene-based ligand as a “magic” Couple toward luminescence, single molecule magnets, and magnetostructural correlations. Acc. Chem. Res.48, 28342842. 10.1021/acs.accounts.5b00296

  • 60

    PointillartF.Le GuennicB.GolhenS.CadorO.MauryO.OuahabL. (2013b). A redox-active luminescent ytterbium based single molecule magnet. Chem. Commun.49, 615617. 10.1039/C2CC37635K

  • 61

    PointillartF.LiuX.KepenekianM.Le GuennicB.GolhenS.DorcetV.et al. (2016). Thermal and near-infrared light induced spin crossover in a mononuclear iron(II). complex with a tetrathiafulvalene-fused dipyridophenazine ligand. Dalton Trans.45, 1126711271. 10.1039/C6DT00920D

  • 62

    QinJ.HuL.LiG. N.WangX. S.WangX. S.XuY.et al. (2011). Synthesis, characterization, and properties of rhenium(I). tricarbonyl complexes with tetrathiafulvalene-fused phenanthroline ligands. Organometallics30, 21732179. 10.1021/om101141h

  • 63

    RenM.BaoS. S.FerreiraR. A.ZhengL. M.CarlosL. D. (2014). A layered erbium phosphonate in pseudo-D5h symmetry exhibiting field-tunable magnetic relaxation and optical correlation. Chem. Commun.50, 76217624. 10.1039/C4CC02085E

  • 64

    RichardsonM. F.WagnerW. F.SandsD. E. (1968). Rare-earth trishexafluoroacetylacetonates and related compounds. J. Inorg. Nucl. Chem.30, 12751289. 10.1016/0022-1902(68)80557-3

  • 65

    RinehartJ. D.LongJ. R. (2011). Exploiting single-ion anisotropy in the design of f-element single-molecule magnets. Chem. Sci.2, 20782085. 10.1039/C1SC00513H

  • 66

    RoosB. O.LindhR.MalmqvistP.-A.VeryazovV.WidmarkP.-O. (2004). Main group atoms and dimers studied with a new relativistic ANO basis set. J. Phys. Chem. A108, 28512858. 10.1021/jp031064+

  • 67

    RoosB. O.LindhR.MalmqvistP.-A.VeryazovV.WidmarkP.-O. (2005). New relativistic ANO basis sets for transition metal atoms. J. Phys. Chem. A109, 65756579. 10.1021/jp0581126

  • 68

    RoosB. O.LindhR.MalmqvistP. A.VeryazovV.WidmarkP. O.BorinA. C. (2008). New relavistic atomic natural orbital basis sets for lanthanide atoms with applications to the Ce diatom and LuF3. J. Phys. Chem. A112, 1143111435. 10.1021/jp803213j

  • 69

    RoosB. O.TaylorP. R.SiegbahnP. E. M. (1980). A complete active space SCF method (CASSCF). using a density matrix formulated super-CI approach. Chem. Phys.48, 157288. 10.1016/0301-0104(80)80045-0

  • 70

    SchulerB.LiuS. X.GengY.DecurtinsS.MeyerG.GrossL. (2014). Contrast formation in kelvin probe force microscopy of single π-conjugated molecules. Nano Lett.14, 33423346. 10.1021/nl500805x

  • 71

    ScottP. L.JeffriesC. D. (1962). Spin-lattice relaxation in some rare-earth salts at helium temperatures; observation of the phonon bottleneck. Phys. Rev.127, 3251. 10.1103/PhysRev.127.32

  • 72

    SessoliR.PowellA. K. (2009). Strategies towards single molecule magnets based on lanthanide ions. Coord. Chem. Rev.253, 23282341. 10.1016/j.ccr.2008.12.014

  • 73

    SheldrickG. M. (2015a). SHELXT – Integrated space-group and crystal-structure determination. Acta Crystallogr. Sect. A Found Adv.71, 38. 10.1107/S2053273314026370

  • 74

    SheldrickG. M. (2015b). Crystal structure refinement with SHELXL. Acta Crystallogr. Sect. C71, 38. 10.1107/S2053229614024218

  • 75

    ShrivastavaK. N. (1983). Theory of spin–lattice relaxation. Phys. Status Solidi117, 437458. 10.1002/pssb.2221170202

  • 76

    SoussiK.JungJ.PointillartF.Le GuennicB.LefeuvreB.GolhenS.et al. (2015). Magnetic and photo-physical investigations into DyIII and YbIII complexes involving tetrathiafulvalene ligand. Inorg. Chem. Front.2, 11051117. 10.1039/C5QI00087D

  • 77

    SpeedS.FengM.Fernandez GarciaG.PointillartF.LefeuvreB.RiobéF.et al. (2017). Lanthanide complexes involving multichelating TTF-based ligands. Inorg. Chem. Front.4, 604617. 10.1039/c6qi00546b

  • 78

    VigneshK. R.LangleyS. K.MurrayK. S.RajaramanG. (2017). Exploring the influence of diamagnetic ions on the mechanism of magnetization relaxation in {CoLn} (Ln = Dy, Tb, Ho). “butterfly” complexes. Inorg. Chem.56, 25182532. 10.1021/acs.inorgchem.6b02720

  • 79

    WeigendF.AhlrichsR. (2005). 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.7, 32973305. 10.1039/b508541a

  • 80

    WoodruffD. N.WinpennyR. E. P.LayfieldR. A. (2013). Lanthanide single-molecule magnets. Chem. Rev.113, 51105148. 10.1021/cr400018d

  • 81

    YamashitaK.MiyazakiR.KataokaY.NakanishiT.HasegawaY.NakanoM.et al. (2013). A luminescent single-molecule magnet: observation of magnetic anisotropy using emission as a probe. Dalton Trans.42, 19871990. 10.1039/C2DT32785F

  • 82

    YiX.BernotK.Le CorreV.CalvezG.PointillartF.CadorO.Le GuennicB.et al. (2014). Unraveling the crystal structure of lanthanide-murexide complexes: use of an ancient complexometry indicator as a near-infrared-emitting single-ion magnet. Chem. Eur. J.20, 15691576. 10.1002/chem.201303833

  • 83

    ZhangL.JungJ.ZhangP.GuoM.ZhaoL.TangJ.et al. (2016). Site-resolved two-step relaxation process in an asymmetric Dy2 single-molecule magnet. Chem. Eur. J.22, 13921398. 10.1002/chem.201503422

  • 84

    ZhangP.ZhangL.TangJ. (2015). Lanthanide single molecule magnets: progress and perspective. Dalton Trans.44, 39233929. 10.1039/C4DT03329A

  • 85

    ZhuZ.GuoM.LiX.-L.TangJ. (2019). Molecular magnetism of lanthanide: advances and perspectives. Coord. Chem. Rev.378, 350364. 10.1016/j.ccr.2017.10.030

Summary

Keywords

tetrathiafulvalene, triads, dysprosium, single-molecule magnet, ab initio calculations

Citation

Lefeuvre B, Galangau O, Gonzalez JF, Montigaud V, Dorcet V, Ouahab L, Le Guennic B, Cador O and Pointillart F (2018) Field-Induced Dysprosium Single-Molecule Magnet Based on a Redox-Active Fused 1,10-Phenanthroline-Tetrathiafulvalene-1,10-Phenanthroline Bridging Triad. Front. Chem. 6:552. doi: 10.3389/fchem.2018.00552

Received

25 May 2018

Accepted

25 October 2018

Published

13 November 2018

Volume

6 - 2018

Edited by

Hitoshi Ishida, Kitasato University, Japan

Reviewed by

Sukhendu Mandal, Indian Institute of Science Education and Research, Thiruvananthapuram, India; Juan Manuel Herrera, Universidad de Granada, Spain

Updates

Copyright

*Correspondence: Olivier Cador Fabrice Pointillart

This article was submitted to Inorganic Chemistry, a section of the journal Frontiers in Chemistry

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics