Molecular nanostructures may constitute the fabric of future quantum technologies, if their degrees of freedom can be fully harnessed. Ideally one might use nuclear spins as low-decoherence qubits and optical excitations for fast controllable interactions. Here, we present a method for entangling two nuclear spins through their mutual coupling to a transient optically excited electron spin, and investigate its feasibility through density-functional theory and experiments on a test molecule. From our calculations we identify the specific molecular properties that permit high entangling power gates under simple optical and microwave pulses; synthesis of such molecules is possible with established techniques. DOI: 10.1103/PhysRevLett.104.200501 PACS numbers: 03.67.Bg, 03.67.Lx, 31.15.EÀ, 33.40.+f Molecules are promising building blocks for quantum technologies, due to their reproducible nature and ability to self-assemble into complex structures. However, the need to control interactions between adjacent qubits represents a key challenge [1][2][3][4]. We here describe a method for optical control of a nuclear spin-spin interaction that presents several advantages over conventional NMR quantum processors: the spin-spin interaction can be switched, the gates are faster, and the larger energy transitions facilitate polarization transfer onto the nuclear spins. After explaining the theory behind our method, we present an experimental study of a test molecule, and show with density-functional theory that an entangling gate could be achieved.We consider two nuclear spin qubits labeled n and n 0 and one mediator e. The mediator possesses a paramagnetic excited state jei with spin one character and a diamagnetic, spinless ground state j0i. The two nuclear qubits do not interact with each other directly, but both couple to the excitation via an isotropic hyperfine (HF) coupling with generally unequal strengths A and A 0 , see Fig. 1(a). The Hamiltonian in a magnetic field is given by (@ ¼ 1):where S z;i and S i are the Pauli spin operators and ! i denotes the respective Zeeman splittings (i ¼ n, n 0 , e). D is the zero-field-splitting (ZFS). ! 0 denotes the optical frequency corresponding to the creation energy of the excitation. Let us first analyze the case of a symmetric (homonuclear) system with ! n ¼ ! n 0 and A ¼ A 0 . Using degenerate perturbation theory and assuming that the electronic Zeeman splitting is much larger than that of the nuclei, ! n ( ! e and the ZFS, D ( ! e , we obtain an effective Hamiltonian by approximating H sym : ¼ hejHjei; thus,where V ¼ ðjc 1 ijc 2 i Á Á Á jc 12 iÞ is the matrix of the approximate eigenvectors up to first order and the E i (i ¼ 1; . . . ; 12) are the eigenenergies up to second order. This reveals that the time evolution of the entire system can be