Open Access Article
Steen H. Hansen*a,
Christian D. Buch
a,
Bela E. Bode
*b and
Stergios Piligkos
*a
aDepartment of Chemistry, University of Copenhagen, Universitetsparken 5, 2100, Denmark. E-mail: piligkos@chem.ku.dk
bEaStCHEM School of Chemistry, Biomedical Sciences Research Complex, and Centre of Magnetic Resonance, University of St Andrews, North Haugh, St Andrews KY16 9ST, UK. E-mail: beb2@st-andrews.ac.uk
First published on 1st September 2025
We herein demonstrate the synthesis of a pair of enantiomerically pure YbIII complexes by post-functionalisation of the parent YbIII complex via condensation with an enantiomerically pure chiral amine. The enantiomeric pair is structurally characterised by single crystal and powder X-ray diffraction, showing that it crystalises in the P212121 Sohncke space group with Flack parameters close to zero, which confirms their enantiopurity. Circular Dichroism (CD) and absorption spectroscopies in the NIR reveal sharp 2F7/2 → 2F5/2 f–f transitions, with gabs values up to 0.07, indicating a chiral environment for the ytterbium centre. Furthermore, a dynamic mechanism with mixing of ligand states is shown to contribute to the CD intensity. X-band pulse Electron Paramagnetic Resonance spectroscopy, on a magnetically dilute single crystal containing 1% of YbIII complexes within the isostructural YIII diamagnetic host, reveals a phase memory time, Tm, of the electronic spin of 600 ns and that it can be coherently manipulated by microwave pulses, as evidenced by Rabi nutations.
Coherently addressing and manipulating the state of a spin qubit, such as a molecular spin, is typically achieved by use of magnetic dipole transitions induced by the oscillating magnetic field component of microwave pulses in an externally applied magnetic field.17 However, a very interesting alternative avenue is the potential substitution of microwave pulses with the application of electric fields coupled to the spin qubit,18 the magnetic dipole transition in this case being mediated by the magneto-electric coupling.19 Furthermore, electric fields can be used to tune the resonance frequency of molecular qubits and could be used to selectively bring specific qubits on and off resonance within multiqubit processor setups.20 The main advantages of using electric rather than magnetic fields include the ability to control electronics on a nanosecond timescale, the highly precise electronics instrumentation already developed and the very power-efficient operation of electronic circuits.21–23 Efficient coupling of molecular spins to electric fields requires absence of an inversion centre at the molecular level.20,24 Hence, designing new chiral molecules with relevant properties for quantum information processing or other QTs is of great interest.25 Since Yb(trensal) has been shown to display numerous properties of interest for QTs, an enantiomeric pure analogue of Yb(trensal) would provide an interesting avenue for studying chiral effects related to quantum information processing. From a preparative perspective, there is only a very limited number of options available for the synthesis of enantiomerically pure lanthanide complexes where chirality is due to the helicity of the coordination sphere environment. This is due to the very small stabilisation effect of the ligand field, resulting in near instantaneous racemization in these complexes in solution.26 Therefore, the synthesis of chiral lanthanide complexes often involves reacting a simple lanthanide salt or nonchiral coordination complex with a previously synthesized enantiomerically pure ligand,27–35 or by precipitation of an ionic complex with a chiral counterion.27 We have previously shown that post-functionalization of the complex LLn (Fig. 1), with Ln = Gd, Tb, Dy, Ho, Er, Tm, Yb or Lu, with primary amines provides a direct route towards adding desired properties to Ln-based coordination complexes.36 In these previous studies we have only used non-chiral amines.37,38 Herein, we demonstrate that condensation of LYb with an enantiomerically pure primary amine (S- or R-methyl-benzylamine) provides the means to form chiral ytterbium complexes exhibiting coordination sphere helicity (ΔYb or ΛYb), an avenue that has not been previously explored for trensal-like ligands. These studies affirm the versatility of our post-functionalization approach to synthesize functional coordination compounds. The chiroptical properties of the synthesized ytterbium complexes are examined, and their origin investigated. Furthermore, these chiral materials are examined by a combination of continuous wave (c.w.) and pulse Electron Paramagnetic Resonance (EPR) spectroscopy demonstrating that they retain a suitable phase memory time and the ability of their electronic spin to be coherently manipulated by microwave pulses.
![]() | ||
| Fig. 1 Reaction scheme for post-functionalisation of LLn with methyl-benzylamine. Red asterisks indicate the stereogenic centre. | ||
:
1 chloroform/methanol solution (20 ml) and S- or R-methyl-benzylamine (0.4 ml, 3 mmol) was added resulting in ΔYb or ΛYb, respectively (vide infra). The reaction was stirred for 1 hour after which a clear solution was obtained. 200 ml diethylether were added to the solution which was covered with a glass cover and left to crystalize for 2–4 days.
Yield: around 160 mg, 58%. Compositional and phase purity was confirmed by elemental analysis (Table S1), IR (Fig. S11 and S12), and X-ray powder diffraction (Fig. S15).
1H-NMR (Fig. S16–S18) was obtained using a Bruker 500 MHz instrument equipped with a cryoprobe. For 1H-NMR calibration was done against solvent signals from the deuterated solvent. Positive-ion mode MALDI mass spectrometry (Fig. S19–S22) was performed on a Bruker Solarix XR 7T ESI/MALDI FT-ICR MS instrument at the Department of Chemistry, University of Copenhagen. Elemental (C, H, and N) analyses were performed on a FlashEA 1112 instrument at The Microanalytical Laboratory at the Department of Chemistry, University of Copenhagen.
Single crystals were measured on a BRUKER D8 VENTURE diffractometer equipped with a Mo Kα High-brilliance IμS 53 radiation source (λ = 0.71073 Å). A PHOTON 100 CMOS detector was employed. The instrumentation was controlled using APEX2. The sample was cooled to 120 K using an Oxford cryosystem. The resulting data were modelled using SHELXT with intrinsic phasing and refined using SHELXL (Least squares).39 Visualisation of the refinement process was done using OLEX2.40 Hydrogens were added using the “Add H” function in Olex2 and refined isotopically, while all non-hydrogen atoms were refined anisotropically.40
Absorption spectroscopy was performed on a PerkinElmer Lambda 2 UV/Vis spectrometer with a 2 nm slit width. Circular dichroism (CD) was measured on a Jasco J1700 equipped with an InGaAs detector for NIR with a 5 nm slit. Both Absorption and CD (Fig. 3, S13 and S14) were measured in approximately 15 mg ml−1 solutions in chloroform in 10.0 mm path quartz cuvettes. Absorption and CD were measured back to back on the same cuvette. Infrared spectra were measured on polycrystalline samples on an Agilent Technologies Cary 630 FTIR.
EPR spectra were measured on a magnetically dilute sample of ΛYb in the isostructural diamagnetic ΛY host lattice at a concentration of 1% (ΛY0.99Yb0.01), as determined by Inductively Coupled Plasma Mass Spectrometry (ICP-MS). X-band c.w.-EPR data were measured on a polycrystalline powder sample using a BRUKER E500 EPR spectrometer fitted with a Bruker X-band ER4122 SHQE cavity resonator. The measurements were performed at 20 K using an Oxford Instruments helium flow cooling system.
Single crystal pulse EPR was measured on a BRUKER X-band E580 EPR spectrometer fitted with an MS3 resonator appropriate for smaller samples and offering higher B1 excitation fields, on single crystals of approximate size 0.2 × 0.2 × 3 mm.3 A cryogen free variable temperature cryostat from Cryogenics Ltd was used. Echo-Detected Field-Swept (EDFS) EPR spectra were measured using a standard Hahn echo sequence (π/2–τ–π–τ-echo) with π/2 = 16 ns. Tm was measured by recording the time evolution of the Hahn echo, with π/2 = 128 ns to decrease ESEEM (electron spin echo envelope modulation). Spin lattice relaxation, T1, was measured using a standard inversion recovery sequence (π–T–π/2–τ–π–τ-echo) with π/2 = 16 ns. Both T1 and Tm were modelled as mono-exponential decays. Transient nutation was measured using a θ–T–π/2–τ–π–τ-echo sequence with θ the nutation pulse angle and T = 5000 ns, at 6 different power levels between 6–18 dB.
The c.w.- and EDFS EPR spectra were fitted using Pepper in Easyspin (version 6).41 EDFS spectra can be assimilated to zero-th derivative c.w. spectra in cases where Tm is the same for all observed transitions and ESEEM or other dynamic modulation effects affecting individual echo intensities are negligible. The single crystal orientation was determined by fitting the sample orientation and molecular frame using the spin Hamiltonian parameters obtained for the polycrystalline powder.
![]() | ||
| Fig. 2 Crystal structure of ΔYb and ΛYb seen from the side showing the different handedness of the structures. Colour code: Yb (green), N (Blue), O (red). Hydrogen atoms were omitted for clarity. | ||
| Napi–Ln/Å | Nimi–Ln/Å | O–Ln/Å | (Napi–Ln–O)/° | |
|---|---|---|---|---|
| LYb | 2.609(3) | 2.420(2) | 2.160(2) | 127 |
| ΔYb/ΛYb | 2.597(3) | 2.412(2) | 2.156(2) | 125 |
Chiral molecules interact differently with left and right circularly polarised light because they break space parity symmetry. ΔYb or ΛYb display several absorption bands between 900–1050 nm belonging to the 2F7/2 → 2F5/2 transition, with the three major bands (marked a–c in Fig. 3) displaying fine structure which is assigned as vibrational fine structure (Fig. 3). The intensity of these transitions is very weak (ε < 20 L mol−1 cm−1) which is characteristic for lanthanides.42 Two minor absorption bands (marked dh) are assigned as hot-bands due to their emergence at lower energy and being of very weak intensity, in agreement with previous observations for Yb(tensal).12 The nearly identical absorption spectra in the solid state and in solution indicate dissolution of complexes without major structural changes or decomposition. This is supported by 1H-NMR (Fig. S16–S18) and MALDI-MS (Fig. S19–S22). The CD spectra of the 2F7/2 → 2F5/2 transitions are mirror images of each other with altering phase for various transients. The largest CD signal is observed at transition a which shows a Δε = 0.7 L mol−1 cm−1 and the largest dissymmetry factor, gabs = 0.07. Although the theoretical maximum absolute value of gabs is 2, typical values for chiral organic molecules or transition metal complexes lie in the range 10−4 to 10−2, with lanthanide complexes displaying in general values above this range, sometimes even close to the theoretical maximum.43–45 Thus, the dissymmetry displayed by ΔYb or ΛYb is sizeable but not exceptional for f–f transitions. This sizeable dissymmetry indicates that the relevant Yb centre feels a chiral electric potential. Thus, a strong magneto-electric response to electric fields is possible. The origin of the CD signal can either be via a static (mixing with ytterbium-based 5d orbitals) or dynamic (mixing with exited states of the ligand) mechanism. For a purely Ln-centred transition to have a transition moment, the transition must contain in part some 5d orbital component to be parity allowed. This is mediated by the dissymmetric ligand allowing this mixing. Due to sum rules the integral of the whole band is hence zero.46 On the other hand, a dynamic mechanism involves polarisation from the ligand which allows for transitions to be slightly allowed. However, due to mixing with the ligand, the integral of the transition can be non-zero.46 The integral of the CD spectra of ΔYb/ΛYb show clear symmetric divergence from zero (Fig. S13), suggesting that some of the intensity of the CD spectra originates from polarisation from the ligand. Hence, designing ligands which facilitate mixing of the lanthanide orbitals with exited states of the ligand could provide a path towards very sensitive electric field coupling.
![]() | ||
| Fig. 3 Absorption and CD spectra of ΔYb and ΛYb in chloroform. A solid-state absorption spectrum of ΔYb is also shown. Labels are defined in the main text. | ||
Having shown that the YbIII centre experiences a chiral environment due to the dissymmetry induced by the chiral ligand, we now examine whether condensation with a chiral amine perturbs its susceptibility to coherent manipulation by microwave pulses, as compared to LYb. Detailed information about the g- and hyperfine coupling tensor, A, of ΔYb/ΛYb was obtained from c.w.-EPR spectra on polycrystalline ΛY0.99Yb0.01 (Fig. 4). The c.w.-EPR spectrum was modelled as originating from an effective S = 1/2 system, due to the large zero field splitting of the eigenstates of the 2F7/2 term, resulting in a thermally isolated ground Kramers doublet. The following Hamiltonian was used to model the obtained c.w.-EPR spectra:
The fit (χ2 = 0.0163) matches the experimental data very well (Fig. 4) and results to the following best-fit parameters: gx = 2.574(2), gy = 3.274(0), gz = 3.697(2), Ax = 523(4) MHz, Ay = 664(2) MHz, Az = 766(4) MHz.
The resulting parameters are within the range observed for other similar systems based on the Yb(trensal) motive.12,36,47 In contrast to these other systems, ΔYb/ΛYb display rhombic g- and hyperfine exchange tensors, even though to a good approximation the molecule possesses trigonal symmetry. However, unlike similar molecules, ΔYb/ΛYb do not possess strict crystallographic trigonal symmetry, hence allowing the g- and hyperfine exchange tensors to be rhombic. These observations illustrate the importance of, and the necessity for, strict crystallographic symmetry for obtaining systems characterised by Hamiltonians of accordingly high symmetry.
Pulse-EPR measurements were conducted on a single crystal of ΛY0.99Yb0.01 for which orientation selectivity between the four different magnetically inequivalent sites of the unit cell was possible without overlapping contributions, as would be the case for a polycrystalline sample. This also results to long Rabi nutation traces possible due to the larger B1-field homogeneity over the single crystal, as compared to the one for polycrystalline powder samples which are intrinsically less homogeneous. The specific orientation was chosen such that a large magnetic field splitting between the observed lines was obtained (Fig. 5). The EDFS spectrum reveals 4 major lines, labelled A–D in Fig. 5, each corresponding to one of the 4 different magnetically inequivalent orientations of crystallographically identical YbIII sites within the crystal possessing no nuclear spin (I = 0). The smaller peaks observed originate from hyperfine coupling to the less abundant 173Yb and 171Yb isotopes possessing a nuclear spin, I = 5/2 and I = 1/2, respectively. The dynamics of the A–D lines (Fig. S23–S30 and Tables S7 and S8) originating from I = 0 isotopes were measured for all 4 inequivalent orientations within the sample (Fig. 5 and 6). Slightly longer Tm's than for Yb(trensal) were observed at similar doping level (1%), likely due to the methyl-benzyl arms increasing the average distance between paramagnetic centres. T1 shows a steep temperature dependence, T1 ∝ T−4.4, which is also observed for similar ytterbium compounds due to the large orbital angular momentum of the ground state resulting in strong coupling to the lattice. Consequently, Tm slowly decreases with increasing temperatures, with 14 K being the highest temperature at which an echo is detectable, due to T1 limiting Tm at higher temperatures. This limit is however slightly lower than for Yb(trensal) which is T1-limited at around 20 K. This difference could originate from the increase of the number of phonons resulting from the introduction of the extra benzyl amine moieties.
![]() | ||
| Fig. 6 Temperature dependence of T1 and Tm for the four main transitions A–D of ΛY0.99Yb0.01. A model of T1 ∝ T−4.4 is shown as a black line. | ||
Transient nutation experiments display clear oscillations of the echo intensity (Rabi oscillations) with nutation pulse duration (Fig. 7 and S31–S34). As expected, the Rabi frequencies, extracted by Fourier Transform of transient nutation traces (Fig. S35–S38), show a linear dependence to microwave power (B1). The Rabi frequencies at a given B1 come in pairs with the A and B resonances in general showing slower nutation frequencies than the C and D ones, due to the anisotropic g-tensor. This indicates that the orientation of the sample plays an important role for the Rabi frequency, for systems with large g-tensor anisotropy. Transitions A–D display more than 30 coherent oscillations indicative of a homogenous B1 field across the sample. Observation of a large number of Rabi oscillations is indicative of potential for high fidelity of operations to be implemented on the system.
![]() | ||
| Fig. 7 Left: Rabi frequencies for the four transitions A–D of ΛY0.99Yb0.01 versus B1. Right: Rabi nutation at various power levels for transition A. | ||
Supplementary information is available. See DOI: https://doi.org/10.1039/d5qi00977d.
CCDC 2428735 (ΔYb), 2428867 (ΛYb), 2428550 (ΔY) and 2428518 (ΛY) contains the supplementary crystallographic data for this paper.48a–d
All structures were measured at 120 K.
| This journal is © the Partner Organisations 2025 |