Lorena E.
Rosaleny
and
Alejandro
Gaita-Ariño
*
Instituto de Ciencia Molecular, Universitat de València, c/ Catedrático José Beltrán, 2, 4980, Paterna, Spain. E-mail: alejandro.gaita@uv.es; Fax: +34 96354 3273; Tel: +34 96354 4421
First published on 20th October 2015
Lanthanoid complexes are amongst the most promising compounds both in single ion magnetism and as molecular spin qubits, but their organization remains an open problem. We propose to combine Lanthanide Binding Tags (LBTs) with recombinant proteins as a path for an extremely specific and spatially-resolved organisation of lanthanoid ions as spin qubits. We develop a new computational subroutine for the freely available code SIMPRE that allows an inexpensive estimate of quantum decoherence times and qubit–qubit interaction strengths. We use this subroutine to evaluate our proposal theoretically for 63 different systems. We evaluate their behavior as single ion magnets and estimate both decoherence caused by the nuclear spin bath and the interqubit interaction strength by dipolar coupling. We conclude that Dy^{3+} LBT complexes are expected to behave as SIMs, but Yb^{3+} derivatives should be better spin qubits.
This bio-nanotechnological strategy will eventually be applied for quantum computing purposes, where the challenging goal of scalability requires the ability to organize different kinds of quantum building blocks. Obviously, the use of biopolymers to control quantum effects in complex organized systems is still a long-term goal. Nevertheless, the nascent field of quantum biology, devoted to the study of coherent quantum effects in processes as diverse as photosynthesis in plants,^{5} geolocation in birds^{6} and possibly smell in insects,^{7} shows that this complex organization of quantum coherent processes already takes place in nature. The challenge is then to achieve this artificially.
Spin-carrying metalloproteins, which are already being studied by manipulating their quantum states via pulsed EPR,^{8} are promising systems for this synthetic quantum biology. We will focus on magnetic lanthanoid complexes because of the interest they awaken both as Single Ion Magnets (SIMs) and as spin qubits, that is, because of their favourable magnetic and quantum properties.^{9}
A key experimental advance in this context is Lanthanide Binding Tags (LBTs). These are oligopeptides based on calcium-binding motifs of EF-hand proteins^{10} that have been designed to interact very specifically with lanthanoids.^{11}
These new building blocks constitute a key advance towards an interdisciplinary region, since an LBT can be seen not only as a small part of a protein but also as a standard coordination complex. LBTs are easily incorporated at the DNA level into any recombinant protein, a potential route to highly complex organisation, e.g. via histones (Fig. 1).^{14}
Fig. 1 (a) Coordination environment created by a 17 aminoacid-long LBT. Only the coordinating lateral residues are shown, α-carbons are labeled. (b) Polypeptidic chain of the LBT without lateral residues, the α-carbons of aminoacids involved in the coordination sphere are labeled. (c) Each nucleosome (8 histones) could organise up to 8 LBTs (4 shown). (d) Nucleosome Positioning Sequences^{12} could organise a scalable sequence of nucleosomes. (c) & (d) Reproduced with permission from ref. 13. |
The first goal of this work is the development of a new computational subroutine to allow an inexpensive estimate of decoherence times and interaction strengths. A tool of this kind is lacking in the field: while wavefunctions are routinely calculated, only general arguments are given concerning the tunneling gap and its relation to decoherence, but detailed numbers are rarely offered. A second purpose will be to theoretically explore the possibility of using LBT for organising the lanthanoid ions, either for their use as SIMs or as spin qubits. This necessarily includes obtaining realistic estimates for decoherence times and interaction strengths that pave the way for the first experimental studies.
First we analyzed the aptitude of these systems as single ion magnets. Next, we started from this set of energy levels and wavefunctions and used a specially crafted version of SIMPRE to calculate the expected quantum decoherence due to interaction of the electron spin qubit with the nuclear spins in the biomolecule, both in the native form and after deuteration (see Methods: Calculations). Finally we calculated interqubit coupling, which needs to be strong enough for two-qubit quantum gates to happen within the coherent time window.
Fig. 2 Box plots with the full distribution of expectation values 〈J_{z}〉 for diverse LBT structures substituted by Dy (left) and Er (right). |
Fig. 3 Box plots with the full distribution of Ω/Δ ratios for diverse LBT structures substituted by Yb (left) and Dy (right). |
In all the structures studied, the expectation values 〈J_{z}〉 stay close to the maximum theoretical values for both Dy and Er, that is, an almost pure Ising behavior is obtained for those two metals. This contrasts with the rest of the series, where a dispersion of behaviors is obtained. The second relevant parameter is the energy level scheme, here summarized in the Ω/Δ ratio. From this point of view, Yb and Nd are unrivalled, with values that routinely go over a 100 and never below 10, and Dy is an acceptable compromise, with a rather narrow distribution around Ω/Δ = 10. This means that, within the expected variation of coordination environments in LBT complexes, only Dy is consistently expected to show single molecule magnet behavior.
Fig. 4 Box plots with the full distribution of estimates decoherence times τ for diverse LBT structures substituted by Tb (left) and Yb (right). |
Because of the nature of dipolar coupling, the relative orientation between the magnetic axes and the field establishes a vanishing lower bound for the coupling between two ions. Therefore, we estimate here the upper bounds for interqubit coupling in double-lanthanide-binding tags, assuming an optimal alignment between the magnetic axes of two neighbouring magnetic complexes. We do this for two extreme examples: Tm and Dy. Both ions present an adequate energetic isolation of the ground doublet: (Ω/Δ)_{Tm} = 7.97, (Ω/Δ)_{Dy} = 26.69, but, as discussed below, are practically opposite in the nature of their ground doublets.
In the case of Tm[2OJR], the ground doublet has an easy plane character, and as a result the expectation values are dominated by the x orientation (see Table 1). In turn, these result in the following differences (Δ_{01}H_{α}) between the fields (H_{α}) created by the two qubit states |Ψ_{0}〉, |Ψ_{1}〉 of a Tm site on a neighbouring Tm site: Δ_{01}H_{x} = 1.10 mT, Δ_{01}H_{y} = 0.99 mT, Δ_{01}H_{z} = 0.29 mT. This means an upper limit of 0.198 μeV for the interqubit coupling.
Main M_{J} components (weight) | 〈J〉 | |
---|---|---|
|0, Dy〉 | |−15/2〉 (95%) | 〈J_{xy}〉 = 0.59 |
〈J_{z}〉 = 7.29 | ||
|1, Dy〉 | |+15/2〉 (95%) | 〈J_{xy}〉 = 0.61 |
〈J_{z}〉 = −7.28 | ||
|0, Tm〉 | |−3〉 (10%), |0〉 (58%), |+3〉 (10%) | 〈J_{xy}〉 = 3.29 |
〈J_{z}〉 = 0.03 | ||
|1, Tm〉 | |−4〉 (8%), |−2〉 (32%), |+2〉 (32%), |+4〉 (8%) | 〈J_{xy}〉 = 6.40 |
〈J_{z}〉 = −0.01 |
For Dy[2OJR], the ground doublet has a marked easy axis character, with almost maximal expectation values 〈J_{z}〉 (see Table 1). These relatively high magnetic moments result in correspondingly larger differences Δ_{01}H_{α} between the field created by the two qubit states |Ψ_{0}〉, |Ψ_{1}〉 of a Dy site on a neighbouring Dy site: Δ_{01}H_{x} = 4.20 mT, Δ_{01}H_{y} = 2.17 mT, Δ_{01}H_{z} = 1.49 mT. This means an upper limit of 3.54 μeV for the interqubit coupling.
A dipolar interaction that cannot be switched off (“always-on”) in the order of the μeV means times for swap operations in the order of the nanosecond, which is also the order of magnitude for pulsed EPR operations. It is also within technologically accessible limits (1–100 GHz, i.e. 4–400 μeV or 0.05–5 K). According to this estimate and considering decoherence times as calculated above, this approach is then theoretically feasible. The ability of proteins to produce an on-demand spatial distribution of qubits means any conceivable scheme of dipolar couplings is available. Of course, in order to actually exploit polypeptide-organized qubits for a scalable quantum information processor, a new operating scheme would need to be developed. Lloyd's proposal, based on a periodic organisation of three different qubit types, would probably be a good start for this, as it is based on energetic, rather than spatial, addressing of the qubits.^{15}
1TJB is a 2.0 Å resolution X-ray crystal structure of a 17-residue lanthanide-binding peptide, complexed with Tb^{3+}, which excludes water molecules from the primary coordination sphere.^{17} There are two crystallographically independent metal sites, corresponding to two separate copies of the same LBT.
2OJR is a construct of a double-lanthanide-binding tag as an N-terminal fusion of ubiquitin complexed with Tb^{3+}, with a 2.60 Å resolution.^{18}
3LTQ is a 2.1 Å X-ray crystal structure of a construct containing an LBT insert between the middle S-loop residues of interleukin-1β complexed with Tb^{3+}.^{19} There is only one crystallographically independent metal site.
3VDZ is a 2.4 Å X-ray crystal structure of a modified dLBT-ubiquitin chimera complexed with Gd^{3+}. The LBT sequence was modified to (a) enhance the exchange of water molecules in the vicinity of the magnetic site, (b) keep a high affinity for the lanthanoid and (c) favour crystallization. There are four crystallographically independent metal sites, corresponding to two separate copies of a “dinuclear” peptide.^{20}
For all metal sites in 1TJB, 2OJR and 3LTQ, we considered eight oxygen atoms in the coordination sphere: two bidentate carboxylates, three monodentate carboxylates, and a carbonyl group belonging to the LBT backbone.
In contrast, the coordination spheres in 3VDZ are more diverse and less well-defined, up to the point where the considered coordination number is somewhat arbitrary. A r < 3.5 Å criterion results in eight oxygen atoms in the coordination sphere. One of these always belongs to a backbone carbonyl group, while the rest, depending on the cases, are from three to five carboxilate groups and either zero or one water molecule.
Hydrogen atoms were not resolved crystallographically, and instead their positions were estimated with Mercury software. As our purpose is to estimate the order of magnitude effect of a hydrogen cloud in peptide-coordinated lanthanoid ions, this lack of precision in the coordinates of each particular H atom is not considered critical.
Lanthanoid complexes are commonly isostructural to each other, with the metal–ligand bond distance being the main structural parameter that varies with the nature of the metal. Thus, we adapted the coordination environment from the original Ln = Tb/Gd structures to the complete series Ln = Nd, Tb, Dy, Ho, Er, Tm, Yb by changing the radial coordinates in the coordination sphere according to the variation in the ionic radii (see Table S1†).
Note also that the current version of SIMPRE automatically chooses the orientation of the coordinate axes that correspond to the most simple expression of the wavefunction, and applies the magnetic field in this z direction. We maintain the standard definition of Δ as an extradiagonal term in the qubit basis, meaning that we redefine the qubit states and that they do not necessarily correspond to the spin being aligned with the easy axis of magnetization. In turn, this results in sub-optimal decoherence times. As a test case, we chose the case of Yb[3VDZ4], as Yb^{3+} is the ion for which we calculated the longest decoherence times and Yb[3VDZ4] is the LBT complex where Yb^{3+} is expected to present the most marked Ising behavior: M_{J} = +7/2 accounts for 92% of the wavefunction, resulting in an expectation value 〈J_{z}〉 = 3.2 and a decoherence time due to the nuclear spin bath τ = 5.2 × 10^{−5} s, which is among the lowest calculated for Yb[LBT] complexes. After a 90° rotation of the molecule, the field is along a hard axis and therefore there is a quenching of the expectation value 〈J_{z}〉 = 0.33. As a consequence, the calculated decoherence time rises to τ = 5.7 × 10^{−4} s which is among the highest calculated for Yb[LBT] complexes. This is comparable with Yb[3VDZ5], the LBT complex described as an easy-plane behavior, with M_{J} = +1/2 accounts for 92% of the wavefunction, resulting in an expectation value 〈J_{z}〉 = 0.3 and decoherence time due to the nuclear spin bath τ = 4.7 × 10^{−4} s.
SIMPRE was further adapted to extract the expectation values of 〈J_{α}〉 (with α = x, y, z) from the wavefunctions, using the Pauli matrices σ_{α}:
〈J_{α}〉 = 〈ψ|σ_{α}|ψ〉 | (1) |
As we mainly intend to distinguish between Ising and non-Ising character here, in Table 1 we present 〈J_{xy}〉 = (〈J_{x}〉^{2} + 〈J_{y}〉^{2})^{1/2}.
Moreover, this specially crafted version of SIMPRE also takes the coordinates of the hydrogen atoms as an input. Of course, there are an effectively infinite number of hydrogen nuclei in a crystal structure. A cutoff radius for the hydrogen nuclei to be included in our calculation is needed. We neglect every hydrogen nucleus which, on average, is expected to produce 1/100^{th} of the effect produced by the hydrogen nucleus closest to the metal. As the hyperfine interaction falls with the third power of the distance, this means the cutoff radius is a factor of farther away than the nearest hydrogen atom.
From the expectation values of 〈J_{x}〉, 〈J_{y}〉, 〈J_{z}〉 and these coordinates, the dipolar magnetic field (H) felt by each nucleus and the hyperfine interaction energy (E) can be trivially calculated, for each of the two states of the qubit.
(2) |
Of course, by including the nearest magnetic ion, this procedure can be immediately used to estimate the (dipolar) interqubit interaction strength. From the set of hyperfine interactions, we also estimate the nuclear spin bath decoherence time using the standard equation (eqn (3)).^{23} This estimate of decoherence depends on the sum of the energy differences, for each proton i between the two qubit states |0〉, |1〉. The decoherence time is then estimated as a function of the tunneling gap Δ and this energy sum (ω_{i} = E_{0} − E_{1}):
(3) |
Note that this model is only valid at low temperatures and in the cases where the tunneling splitting is much larger than the hyperfine couplings, something that is generally verified for lanthanoid ions, and fields of the order of hundreds of mT.^{25}
For the representation of the results, box plots^{24} were generated with Wavemetrics’ IgorPro and include the full range of values.
Footnote |
† Electronic supplementary information (ESI) available: For all Ln^{3+} not included in Fig. 1, 2 and 3, results are in Fig. S1, S2 & S3. Parameters relevant to the calculations are in Tables S1 & S2. See DOI: 10.1039/c5qi00127g |
This journal is © the Partner Organisations 2016 |