Extending NMR Quantum Computation Systems by Employing Compounds with Several Heavy Metals as Qubits

: Nuclear magnetic resonance (NMR) is a spectroscopic method that can be applied to several areas. Currently, this technique is also being used as an experimental quantum simulator, where nuclear spins are employed as quantum bits or qubits. The present work is devoted to studying heavy metal complexes as possible candidates to act as qubit molecules. Nuclei such 113 Cd, 199 Hg, 125 Te, and 77 Se assembled with the most common employed nuclei in NMR-QIP implementations ( 1 H, 13 C, 19 F, 29 Si, and 31 P) could potentially be used in heteronuclear systems for NMR-QIP implementations. Hence, aiming to contribute to the development of future scalable heteronuclear spin systems, we specially designed four complexes, based on the auspicious qubit systems proposed in our work ( J. Phys. Chem. A 2020 , 124 , 4946–4955), which will be explored by quantum chemical calculations of their NMR parameters and proposed as suitable qubit molecules. Chemical shifts and spin–spin coupling constants in four complexes were examined using the spin–orbit zeroth-order regular approximation (ZORA) at the density functional theory (DFT) level, as well as the relaxation parameters ( T 1 and T 2 ). Examining the required spectral properties of NMR-QIP, all the designed complexes were found to be promising candidates for qubit molecules.


Introduction
Nuclear resonance magnetic (NMR) is a celebrated technique with plenty of applications in chemistry, physics, medicine, and structural biology. Besides its well-known uses, such as identifying molecules and determining their structure and dynamics, NMR has also been used as a experimental quantum simulator [1], where nuclear spins are employed as quantum bits or qubits-the fundamental information unit used in quantum information processing (QIP) [2].
Meanwhile, the advancement of the aimed at large-scale NMR quantum processors depends mainly on two challenges: first, the accurate and efficient control of the quantum system states; second, the construction of systems containing a large number of qubits [3]. The former have been augmented by the development of algorithms that are able to manipulate, with high fidelity, the quantum states of relatively large systems [4]. Recently, Peterson et al. [5] designed a fast and scalable algorithm for controlling the quantum states of systems containing 100 qubits.
These recent algorithms combined with modern NMR equipment and sophisticated theoretical techniques make possible the manipulation of the nuclear spin system states through magnetic radiofrequency (RF) pulses, which are used to implement quantum gates, the basic computational step [5].
Regarding the second challenge, chemistry can leverage the development of NMR large-scale quantum systems, as it depends especially on the design of a qubit molecule containing a large number of suitable spin-1/2 nuclei, which act as the qubits. The properties that determine whether a molecule is suitable as a QIP molecule are-apart from of course the stability of the molecule-its NMR parameters such as the chemical shifts and the indirect nuclear spin-spin coupling constants [6]. These two molecular properties can currently accurately and efficiently be calculated by various quantum chemical methods [7][8][9][10][11][12]. The most widely employed of them is certainly density functional theory (DFT). The natural sequence of designing molecules is today therefore to screen potential candidate molecules by quantum chemical calculations of their relevant properties prior to an actual synthesis. In this way, theoretical chemistry can speed up the sharpening of the potential QIP molecule concerning its physical requirements [13], towards a specially synthesized system, giving resolvable spectra to the upcoming large-scale NMR quantum processor.
Nuclear spins-1/2 are natural qubits and are therefore widely used in NMR-QIP, where systems containing 1 H, 13 C, 19 F, 29 Si, and 31 P nuclei are the most employed ones [3,14,15]. Furthermore, other spin-1/2 nuclei, such as 113 Cd,199 Hg, 125 Te, and 77 Se nuclei, have been investigated as qubits, and assembled with the most commonly employed nuclei, they could potentially be used in heteronuclear systems for NMR-QIP implementations.
NMR-QIP is usually implemented in an ensemble of identical and non-interacting molecules at room temperature [3]. An NMR quantum computer with n qubits requires a molecule with n coupled spin-1/2 [16]. This qubit molecule should present appropriate NMR parameters for QIP's efficient implementation, these being [14,16,17]: (i) large spinspin coupling constants (J) between directly coupled spins, as the typical time required for two-qubit quantum gates is inversely proportional to the size of the coupling between the two coupled nuclei; (ii) a large range of chemical shifts; the frequency differences of the nuclei should be as large as possible to allow the selective manipulation of the individual spins, which is a basic computational step; (iii) the relaxation time of the nuclei should also be large enough to perform a huge number of logic quantum gates in a given algorithm.
In our previous studies, we investigated suitable qubit molecules, as well as suggested structural modifications of them, with the aim to enhance their NMR parameters for the QIP efficient implementation. Lino et al. [16] theoretically explored benzyldene-2,3-dihydro-1H- [1,3]diphosphole (BDF) derivatives in order to develop new promising phosphorus heterocycle compounds for QIP. The BDF-NO 2 derivative was found to be the best candidate for NMR quantum computations among the studied phosphorus heterocycles containing π-conjugated molecular skeletons.
In our following work [18], exceptionally large through-space (TS) P-P spin-spin coupling constants observed in 1,8-diphosphanaphthalenes (PPN) and in naphtho[1,8-cd]-1,2-dithiole phenylphosphines (NTP) were proposed and investigated to provide more accurate control within large-scale NMR-QIP. We explored the large TS J 31 P-31 P couplings to provide tighter, more accurate control through large-scale NMR-QIP, proposing the application of TS J( 31 P, 31 P) coupling as a resource for universal logic in NMR quantum computers. From our results, PPN o -F, PPN o -ethyl, and PPN o -NH 2 were the best candidates for NMR-QIP, in which the large TS SSCCs could face the need for long-time quantum gates' implementations. The following work [19] was devoted to evaluating and calculating directly the relaxation parameters of the suitable qubit molecules proposed in the last mentioned work. In NMR-QIP, decoherence is a key problem, and considering the unique qubits, the time of the coherence phase can be well measured by the time of transversal relaxation (T 2 ), wherein the longer the relaxation time, the better the information processing is. The relaxation parameters (T 1 and T 2 ) of the molecules PPN o -F, PPN o -ethyl, and PPN o -NH 2 were calculated, and the results supported and confirmed that they are suitable qubit molecules, which could improve the control accuracy through large-scale NMR-QIP, highlighted by the promising longer coherence time, avoiding a high decoherence rate.
Along the same lines, our latest work [20] focused on heteronuclear systems with heavy metals such as 113 Cd,199 Hg, 77 Se, and 125 Te as qubits for NMR-QIP. We examined the NMR parameters of metal complexes with phosphine chalcogenide ligands (called MRE) using spin-orbit ZORA and four-component relativistic methods. We developed a computational design strategy for prescreening molecules that could enable many and heteronuclear qubits for NMR-QIP implementations. Particularly, the influence of different conformers, basis sets, functionals, and methods to treat the relativistic, as well as solvent effects was studied. The MRE complexes were found to be multiple-spin systems with Larmor frequencies appropriately dispersed, so well-defined qubits, allowing qubit addressability, together with an exceptionally large spin-spin coupling between the pair of spins, which enables the two-qubit operations.
In this work, we aimed to go a step further in relation to the use of heavy metals ( 113 Cd, 199 Hg, 77 Se, and 125 Te) as qubits. Supported by the findings of our last mentioned work [20], the use of heavy metals combined with the most frequently used qubits ( 1 H, 13 C, 19 F, 29 Si, and 31 P) can boost the emergent scalable heteronuclear spin system in NMR-QIP. An NMR computer can be programmed electronically analogously to a quantum computer, but also, it can be implemented at room pressure and temperature using macroscopic liquid samples [21]. Thereby, we specially designed four complexes, still based on the auspicious qubit systems proposed in our work [20], which will be explored by quantum chemical calculations of their NMR parameters in order to investigate their suitability as qubit molecules.

Computational Details
Geometry optimizations and vibrational frequency calculations of all compounds were performed with the ORCA 4.2 program [22] at the DFT level in the framework of the two-component relativistic zeroth-order regular approximation (ZORA) [23] with the PBE0 exchange-correlation functional and ZORA-def2-TZVP basis set [24]. The SARC-ZORA-TZVP [24] basis set was used for Hg and the old-ZORA-TZVP was used for Cd, Te, and Se. The segmented all-electron relativistically contracted (SARC/J) basis sets were used in all ORCA calculations as the auxiliary basis sets. The SARC/J auxiliary basis set is a decontracted def2/J auxiliary. Corresponding frequency calculations confirmed that the geometries were indeed minima. The Cartesian coordinates of all optimized structures of the isolated molecule and vibrational frequencies are given in the Supplementary Material.
The NMR parameters of the studied complexes were investigated for the isolated molecule at the DFT level also in the ZORA framework, including spin-orbit effects and the Gaussian finite nucleus, as implemented in the Amsterdam density functional (ADF) program [25]. The calculations of the absolute nuclear magnetic shielding constants σ were carried out with the PBE0 exchange-correlation functional together with the triple-ζ (TZ2P) Slater-type basis set [26]. Gauge including atomic orbitals were employed in these calculations. The indirect spin-spin coupling constants (SSCC or J) were evaluated also with the PBE0 functional, but with the TZ2P-J Slater-type basis set [27], which is specially designed for NMR SSCC calculations. Those two theoretical levels were studied in our previous work and exhibited the best performance within the performed benchmarking [20,[28][29][30][31]. One should note, however, that in a recent study, it was shown that Dyall's energy-optimized Gaussian-type basis sets had to be extended with additional f-type functions in order to obtain convergence for spin-spin coupling constants in correlated wave function calculations at the SOPPA(CCSD) level [32]. Whether this also would be the case for the Slater-type basis sets and DFT calculations needs to be seen.
For the present work, the two-component spin-orbit ZORA method was used. This approach has been successfully applied in several systems involving transition metals, including all atoms of the systems used in this work [29,31,[33][34][35]. From a previous study, the concomitant changes in the NMR spectroscopy parameters were, in general, small going from the ZORA calculation to the four-component DKS approach [20].
NMR chemical shifts (δ) were calculated relative to: tetramethylsilane (Si(CH 3 ) 4 , TMS) for 1 H, 13 C, and 29 Si; CFCl 3 for 19 F; PH 3 for 31 P; SeMe 2 for 77 Se; TeMe 2 for 125 Te; CdMe 2 for 113 Cd; HgMe 2 for 199 Hg. All these references compounds were optimized and calculated at the same level as mentioned before. Table 1 shows the individual absolute shieldings. The chemical shifts (δ in ppm) of the nucleus (N) were then calculated as: The 31 P and 113 Cd chemical shifts were subsequently corrected by the gas phase experimental value of CdMe 2 , −706.15 ppm [36], and PH 3 , −266.1 ppm [37], as they are not the reference compounds of the standard chemical shift scale for Cd and P. Comparing our calculated absolute shielding constant for the carbon in TMS, 191.3 ppm, with a value of 188.1 ppm from a semi-experimentally derived absolute shielding scale [38] showed good agreement, taking into account that our value is for the equilibrium geometry, while the experimentally derived value is still for 300 K. In addition, one should recall that calculated chemical shifts or differences in Larmor frequencies, as we study here, are normally more accurate due to error cancellations. Whilst the herein employed theoretical model provides insightful information on the relationship between structure and spectroscopic properties, the clear practical utility lies in its ability to predict the performance of unique molecule classes, thereby directing future efforts. In this line, quantum chemical calculations performed on isolated molecules are useful, as a first approach, to rationalize qubit molecules. Later on, more sophisticated calculations can be employed to introduce thermal and solvent or environment effects, thus making the model more realistic.
Molecular dynamics (MD) simulations were carried out in the ORCA 4.2 program [22], using the functional PBE with the basis set ZORA-def2-TZVP; the temperature of 300 K was included in all MD simulations. The simulation consisted of a thermalization stage of 1 ps, obtaining a total of 1000 conformations. After the MD simulations, we obtained the relaxation times (T 1 and T 2 ) and relaxation rates (r 1 and r 2 ) of the atoms by means of the spectral density. This methodology has already been used successfully in other works [19,39,40]; it uses the distance between two atoms for each conformation selected by the procedure of statistical inefficiency. That way, the parameters T 1 and T 2 were obtained by the spectral density method in the MATLAB 7.6 software [41]. For more datails, please see [19,39].

Results and Discussion
The complexes proposed and investigated in this work are HgF 2 [((CH 3 ) 3 SiCH 2 )(CH 3 3 )HPSe], which are denoted as complex 1, complex 2, complex 3, and complex 4 throughout the article (Figure 1). These complexes were theoretically designed based on the complexes and heavy atoms studied in our previous work [20] and pointed out as auspicious qubit molecules. In this work, structure modifications were performed in a way that all the basic requirements for NMR-QIP were fulfilled. These modifications were based on fragments of chemical structures previously synthesized [42][43][44]. Although the proposed molecules studied in this work have not been synthesized yet, they are based on structures already constructed and could, in principal, be synthesized as well.

Conformational Flexibility
In our previous work [20], we investigated the influence of different conformers in the calculations of the NMR parameters. We carried out a conformer search using Open Babel 2.4.1 [45], and up to 50 conformers were generated in order to identify other conformations. From the Boltzmann averaging, we observed only small variations in both δ and J and concluded, for this reason, that there is no need to perform calculations of the NMR parameters for all the main conformers. Therefore, we followed the trends here in this work and kept working with only one optimized geometry. Therefore, starting from an optimized geometry of each complex, the NMR parameters were calculated at the ZORA/PBE0/TZ2P level for chemical shifts and ZORA/PBE0/TZ2P-J for spin-spin couplings.

Spectroscopic Parameters: Chemical Shift Values
The designed structure of the complexes proposed in this work could overcome the restriction regarding the number of suitable qubits for the forthcoming NMR quantum computer. The chemical shifts of the studied molecules are shown in Table 2. As expected, the complexes presented a huge range of chemical shifts, when all the nuclei were taken into account, as we dealt with heterogeneous systems. Complex 1 exhibited the biggest chemical shift range, 1567 ppm, varying from −1513.02 ppm for Hg1 to 53.53 ppm for C15. Complex 2 presented the second biggest range, which was 1447 ppm, varying from −1387.93 for Hg1 to 58.59 for Se3. The following range is given by complex 3: the range of chemical shifts was 1027 ppm, varying from −700.37 for Te3 to 326.71 for Hg1. The smallest range was seen in complex 4, which was 554 ppm, varying from −250.82 for Te3 to 303.00 for Hg1. In summary, Hg complexes showed the widest chemical shift range; these results are as expected, as Hg spans a wider chemical shift range compared to Cd.
The different chemical shifts afforded to each nucleus in the complexes are a reflection of the different electronic environments in the molecules. When nuclei have distinct gyromagnetic ratios γ; they provide qubits that can easily be individually addressed through RF pulses [46]. Hence, using different spin-1/2 nuclei in differing chemical environments is a kick off for scaling up NMR quantum computers.

Spectroscopic Parameters: Spin-Spin Coupling Constant Values
For multiple-spin systems, while homonuclear spins are often individually addressed by the distinct δ due to different local environments, heteronuclear spins are easily distinguished due to the distinct nuclear gyromagnetic ratio γ and, thus, very different Larmor frequencies ω 0 [4]. Taking this into account, further in the paper, we discuss the NMR-QIP prerequisite of the large chemical shift range in terms of the frequency differences of the coupled nuclei. The nuclei Larmor frequencies of the studied complexes are shown in Tables S3-S6.
To evaluate the relative magnitudes of the spin-spin coupling constant in consonance with the spread in the Larmor frequencies, one needs to come back to the foundations of NMR quantum computing [17]. Quantum logic gates were implemented through RF pulses based on a scalar coupling Hamiltonian [47]. The Hamiltonian for scalar J-coupling of a molecule containing N spin-1/2 coupled nuclei is given by: where J KL is the coupling strength between the spins K and L. In Equation (2), the system is said to be strongly coupled. When the frequency separation between the spins ω K 0 − ω L 0 is large compared to their coupling strength, i.e., when |J KL | |ω K 0 − ω L 0 |, one can simplify Equation (2) to When the condition |J KL | |ω K 0 − ω L 0 | applies, the spectra are also said to be first order and the system is said to be weakly coupled; this is also known as the minimal coupling approximation. This condition is easily satisfied for heteronuclear spins and for small homonuclear molecules when the chemical shift between the nuclei is sufficiently large. Mawhinney and Schreckenbach [17] demonstrated that the differences in chemical shifts and the spin-spin coupling between nuclei are not required to be extremely large. Therefore, they evaluated the required spectral properties and concluded that what is actually required for a first-order spectrum is a relationship of ∆ω 0 J ≥ 10 [48]. Hence, in compliance with minimal coupling approximation and in line with the criteria proposed by Mawhinney and Schreckenbach [17], we examined the required spectral properties for NMR-QIP, evaluating the ratio ∆ω 0 J ≥ 10 for all the designed complexes. The ratio analysis of the studied complexes is shown in Tables S7-S10. We first investigated the directed coupled nuclei, because even though the nuclei do not need to be mutually coupled, they must form a contiguous network of couplings [17]. The spin-spin coupling constant for directed coupled nuclei and the difference in the Larmor frequencies (for a 16.5 T magnet) between coupled nuclei are displayed in Table 3. Taking into account all four complexes, all the pairs of nuclei exhibited a ratio grater than 10. The smallest ratio was 5658 for the pair of nuclei Hg1-Se3 in complex 2, and the largest ratio was 5,172,004 for the pair P4-C14 in complex 1. Analyzing the individual complexes and starting with the Hg complexes, for complex 1 with Te, the smallest ratio, 17, was observed for the pair of nuclei C16-C17, which exhibited a J-coupling constant of 23.53 Hz and a Larmor frequency difference of 410 Hz. At the other end of the range, the largest ratio, 5,172,004, was found for the pair of nuclei P4-C14, which had a J-coupling constant of 20.86 Hz and a Larmor frequency difference of 107,887,999 Hz. For complex 2 with Hg and Se, the smallest ratio, 13, was predicted for the pair of nuclei C16-C17, which exhibited a J-coupling constant of 23.53 Hz and the smallest Larmor frequency difference of 300 Hz. The largest ratio, 4,714,635, was found for the pair of nuclei C12-H27 for which the J-coupling constant was 111.54 Hz and the Larmor frequency difference was 525,851,552 Hz. Concerning the Cd complexes, complex 3 with Te exhibited also for the pair of nuclei C16-C17 the smallest ratio, 11, with a J-coupling constant of 23.64 Hz and a Larmor frequency difference of 268 Hz. The biggest ratio, 4,801,382, was shown by P4-C14, with a J-coupling constant of 22.47 Hz and a Larmor frequency difference of 107,887,060 Hz. The last complex, complex 4 with Cd and Se, exhibited the smallest ratio of 6 for the pair of nuclei C16-C17, for which the J-coupling constant was 23.76 Hz and the Larmor frequency difference was 138 Hz. The largest ratio was presented by C11-H25, with a J-coupling constant of 110.23 Hz and a Larmor frequency of 525,851,481 Hz. As mentioned above, the pair of nuclei C16-C17 presented the smallest ∆ω 0 J ratio for all four complexes. However, we still could propose the nuclei C16 and C17 as suitable qubits, because even though the pair of nuclei C16-C17 exhibited a ratio smaller than 10 for complex 4, for example, the C16 and C17 nuclei are involved in other couplings, forming a contiguous network of couplings. In complex 4, the pair of nuclei C16-H39 presented a ratio of 4,337,957, and the pair of nuclei C17-H42 exhibited a ratio of 4,487,478, which fulfills the requirement outlined by Mawhinney and Schreckenbach [17].

Correlation Time and Spectral Density
The nuclear spin-1/2 systems can be precisely manipulated due to their coherence (lifetime of arbitrary superposition states). Despite this, they decohere slowly on the time scale of the spin interactions, the so-called relaxation time. Furthermore, NMR coherences vanish at the time of the spin-spin relaxation (T 2 ), so the NMR-QIP implementations are typically designed as pulse sequences shorter in time than T 2 [49].
Recently, the combination of NMR and MD investigations has been shown to provide a comprehensive description of the fast conformational dynamics of small molecules, as well as proteins and organic compounds [50]. However, MD simulations generate thousands of structures, and thus, to perform the subsequent quantum calculations, taking into account all of them is computationally expensive and infeasible, which makes the selection of the main conformations important. Thereby, in the present work, the selection of the main MD conformations was carried out by the statistical inefficiency (SI) method, which was previously developed and validated by the group of Coutinho and Canuto [51,52]. This method uses the statistical interval obtained from the energy autocorrelation; thus, it is possible to calculate the interval of the uncorrelated configurations, τ, by integrating C(n) (Equation (5)) from zero to infinity. In light of this, the configurations that are separated by an interval of 2τ or more are pointed out as uncorrelated configurations.
The number of structures selected from the SI method are shown in Table 4. The Hg complexes presented the biggest number of conformations selected, 28 for complex 1 and 24 for complex 2. The molecule used in the validation of the theoretical methodology, TFE, exhibited 19 conformations selected. The smallest number of conformations selected was provided by Cd complexes, eight for both complex 3 and 4. In fact, the higher the correlation time, the smaller the number of selected structures was. After selecting the structures, the next step is to obtain the relaxation times for the atoms of each structure. To obtain these parameters, we used another methodology called spectral density, which has already been successfully applied in other works [19,39]; therefore, the calculation of T 1 and T 2 (through the spectral density) is given by the following equations.
In Equation (6), the rotational correlation time is given by τ C , which is proportional to the inverse of the rotation diffusion constant, S 2 is the order parameter, and τ e is the effective correlation time (τ −1 = τ −1 C + τ −1 e )). From the spectral density, the relaxation rate R 1 and R 2 can be calculated by Equations (7) and (8).
In the above equations, I is the spin quantum number, ω 0 is the Larmor angular frequency of the atoms, µ 0 is the vacuum permeability,h is the Planck constant, and φ is the average distance between the atoms of interest. The Larmor angular frequency refers to the precession rate of the proton magnetic moment around the external magnetic field (B 0 ) and the gyromagnetic radius of the atoms (γ), given by ω 0 = γB 0 [19]. The following topic shows the validation of the spectral density methodology, which is important, as it shows the agreement between the data obtained theoretically and the experimental results.

Validation of the Theoretical Methodology
The compound iodotrifluoroethylene (TFE) (Figure 2) was employed to validate the methodology used. In this step, the distances between the F1-F3 and F1-F2 atoms were taken into account to calculate the relaxation times values. Table 5 shows the theoretical and experimental results of the relaxation times for this molecule. The experimental values were taken from [53]. For F1-F3 atoms, the differences between theoretical and experimental values were: for T 1 , 0.08 s, and for T 2 , 0.05 s. For the F1-F2 atoms, the differences between theoretical and experimental values were 0.31 s for T 1 and 0.03 s for T 2 . Consequently, the differences between the theoretical and experimental values were found to be very small, showing that the methodology used is effective in obtaining the relaxation parameters. Furthermore, the use of DFT offered an efficient way to determine the geometries; this conclusion is in agreement with [54][55][56], which used DFT to successfully calculate relaxation times and hyperfine coupling constants with a relatively low computational cost.

Spectroscopic Parameters: Relaxation Times
Qubit molecules composed of spin-1/2 nuclei are broadly accepted as one of the most powerful systems against environmental noise, which affords them with relatively longer coherence time and allows the accurate control of quantum gate operations [57]. The relaxation process occurs when an arbitrary ensemble evolves toward its equilibrium state. This process may lead the system to decoherence, destroying qubit superposition states [47]. In larger qubit systems, the relaxation rate for each added spin contributes to multiple quantum coherences of the system. However, on the other hand, these individual contributions could destroy qubit superposition states, leading the system to decoherence [58]. For single qubits such as nuclear spin states, a good measure of the lifetime of arbitrary superposition states is the 'longitudinal' (T 1 ) and 'transverse' (T 2 ) relaxation times [58].
The relaxation times (T 1 and T 2 ) values and, consequently, the R 1 and R 2 values (R 1 = 1 T 1 and R 2 = 1 T 2 ) for all nuclei of each complex were calculated (Table 6). Taking into account all the studied complexes, the Si nucleus presented the highest T 1 and T 2 values, 44.80 s for T 1 and 9.97 s for T 2 , while the H the smallest, 0.95 s for T 1 and 0.35 s for T 2 . Consequently, Si had the smallest R 1 and R 2 values, 0.02 s −1 for R 1 and 0.10 s −1 for R 2 , while H the highest, 1.05 s −1 for R 1 and 2.85 s −1 for R 2 . Analyzing the complexes individually, complex 1 provided the highest values of T 1 and T 2 for the Si nucleus, 44.8 and 3.52 s, respectively. For complex 2, the highest T 1 was 32.77 s for Hg and the highest T 2 was 9.34 s for Si. Concerning the Cd complexes, for complex 3, the nucleus with the highest T 1 was Si, 13.99 s and the highest T 2 was Te, 2.05 s. For complex 4, the highest T 1 was for Se, 21.56 s and the highest T 2 was for Si, 9.97 s.
As already mentioned, the longer the relaxation time of the atoms, the better the coherence phase of the system. In this sense, it is important to understand the relaxation time increasing in the systems. From Table 6, the Hg T 1 value is much higher for complex 2 when compared to complex 1, which can be explained based on the electronegativity of atoms directly bonded to Hg. In complex 1, the Hg atom is directly bonded to the Te atom, while in complex 2, the Hg atom is directly bonded to the Se atom. The Se atom is more electronegative than the Te atom, and this causes in complex 2 a higher electron density, which can contribute to the increase in the relaxation time values. The same reasoning can be applied to the other systems. For the complexes with Te, the Te relaxation times were longer for complex 3 than for complex 1, because the Te atom is directly bonded to Hg in complex 1, while Te is directly bonded to Cd in complex 3, and as the Cd atom is more electronegative than Hg, the electron density in complex 3 becomes greater than in complex 1.
As expected, the T 1 is longer than T 2 . Therefore, for single qubits, the coherence phase time could be well measured by the transverse (or spin-spin state) relaxation time, which is the minimum lifetime of arbitrary superposition states [3]. In this way, the NMR relaxation timescale T 2 must be as long as possible in each nucleus of the qubit molecule, thus avoiding a high decoherence rate [59].
In light of this and considering that the time needed for basic quantum operations is of the order of 100 µs [60], the T 2 values presented by all the nuclei from the studied complexes are thus suitable to be used as qubits for NMR-QIP, which implies that the operational accuracy is in principle close to what is needed for a scalable fault-tolerant NMR quantum processor [60].
Thus, the effects of the modifications on the relaxation times are significantly important for NMR-QIP, as the longer relaxation time in each nucleus of the qubit molecule provides the system a relatively longer coherence time, allowing accurate control on the quantum gate operations in a large-scale NMR-QIP.

Conclusions
All the specially designed complexes were found to be promising candidates as qubit molecules. Examining the required spectral properties for NMR-QIP, the suitable qubit molecules studied in this work presented appropriate NMR parameters for efficient NMR-QIP implementation. Following the requirements for NMR quantum computing and a first-order spectrum, the four studied complexes presented: (i) spin-spin coupling constants large enough, enabling two-qubit operations, which could preserve the coherence and reduce the time of quantum gate operations; (ii) an appropriate range of chemical shifts, where the frequency differences of the nuclei are ample, allowing the selective manipulation of the individual spins, the qubit addressability; (iii) nuclei relaxation times large enough to perform the logic quantum gates in a certain algorithm; as in NMR-QIP, the decoherence is a key problem the longer the relaxation time in each nucleus of the qubit molecule, providing the system a relatively longer coherence time and better information processing.
Hence, the structural modifications proposed for the complexes gathered appropriate NMR parameters for the NMR-QIP, which could overcome the restriction regarding the number of suitable qubits for the forthcoming NMR quantum computers and what is needed for a scalable fault-tolerant NMR quantum processor.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data that support the findings of this study are available within the article and the Supporting Information. Further data are available from the corresponding author upon reasonable request.