PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of ncommsLink to Publisher's site
 
Nat Commun. 2017; 8: 16038.
Published online 2017 July 24. doi:  10.1038/ncomms16038
PMCID: PMC5527286

Coherent superpositions of three states for phosphorous donors in silicon prepared using THz radiation

Abstract

Superposition of orbital eigenstates is crucial to quantum technology utilizing atoms, such as atomic clocks and quantum computers, and control over the interaction between atoms and their neighbours is an essential ingredient for both gating and readout. The simplest coherent wavefunction control uses a two-eigenstate admixture, but more control over the spatial distribution of the wavefunction can be obtained by increasing the number of states in the wavepacket. Here we demonstrate THz laser pulse control of Si:P orbitals using multiple orbital state admixtures, observing beat patterns produced by Zeeman splitting. The beats are an observable signature of the ability to control the path of the electron, which implies we can now control the strength and duration of the interaction of the atom with different neighbours. This could simplify surface code networks which require spatially controlled interaction between atoms, and we propose an architecture that might take advantage of this.

In the silicon donor impurity phosphorus, the outermost electron is bound to the ion at low temperature with hydrogen-like states and a level spectrum according to the Kohn–Luttinger (KL) model1. The energy and length scales are determined by the dielectric constant and the effective mass, so the transitions occur in the region of 10 THz and the effective Bohr radius is An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m1.jpg=3 nm. There is a long history of incoherent spectroscopy of shallow impurities (P, Bi and so on) in silicon, normally utilizing the Fourier Transform Infrared technique (FTIR) (Si:Bi2; Si:P2; and more recently our own high field (30 T) work4 in Si:P). Si:P can now be positioned and detected with atomic precision5 to make devices that demonstrate control over coherent quantum interactions between neighbouring donors6,7 and semiclassical interaction with contacts8. Optical control of Rydberg orbits of trapped atoms in vacuum has provided high-fidelity gate interactions9 and fine control over the electron trajectory10, with applications in atomic clocks and quantum computing11,12. Similar experiments have been carried out on optical transitions in semiconductors13, but it is known that in free space atoms one is able to gain better control by including many orbital eigenstates in the wavepacket10,13. The importance of control over neighbouring qubit interactions in silicon has been shown separately6,7,9,13, and it is also appreciated that control over interactions is important for the implementation of surface code error quantum correction14,15,16,17. A convenient radiation source for coherent studies in the THz spectral region is the free-electron laser. Particularly relevant in this regard are Rabi oscillations reported on the 1s→2p+ system of hydrogenic donors in GaAs in a magnetic field18 and the observation of the excitonic Autler–Townes effect, utilizing InGaAs/AlGaAs quantum wells19.

The present work bridges these previous studies by investigating three-eigenstate control of phosphorous orbitals in silicon. An interferometric experiment is implemented using a free-electron laser as the source, targeting the 1s(A1)−2p± orbital transitions under a small magnetic field. A three-eigenstate superposition is generated, and its time evolution demonstrated through direct observation of beats in the interferogram. Measurements in the Si:P system are compared with the theoretical results for both the simple hydrogenic sytem and the multi-valley effective mass approximation appropriate for the crystal. We find that while the multi-valley effective mass theory is complicated, given simplifications appropriate for our experimental conditions the results are directly comparable to the simple and intuitive hydrogenic picture.

Results

Hydrogen three-level superpositions

Orbital quantum control of the hydrogenic states of isolated Si:P donors with THz pulses has so far been restricted to simple two-level superpositions20,21 as in the example of 1s and 2p0 states. We illustrate this with the simpler hydrogen 1s and 2p0, and their superpositions shown in Fig. 1a,b; see next section for donor wavefunctions in Si. In a superposed state, the relative phase of the two contributing hydrogen eigenstates evolves as a function of time resulting in oscillation of the centre of mass. Oscillations of different symmetry may be obtained, for example, with superposition of 1s and 2p states by a circularly polarized pulse, shown in Fig. 1c. Building wavepackets out of multiple states provides a greater variety of possible overlaps and spatio-temporal evolution. For example the 1s, 2p+ and 2p states form a three-eigenstate system that allows freedom of control over the probability density in the x–y plane. To produce an arbitrary superposition state with a single pulse requires complete control of the polarisation and pulse shape. The control can be simplified significantly if the degeneracy of the 2p± states is lifted with a small applied magnetic field, in which case a pair of linearly polarised pulses with variable delay is sufficient for most purposes. The superposition amplitudes, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m2.jpg, may be found analytically in the case of long dephasing times. Consider a Zeeman splitting of the excited state 2ħΔ about a central frequency, ω, and a coherent, resonant, x-polarized pulse with Rabi frequency Ω. If the splitting is small (An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m3.jpg,1/tp, where tp is the pulse length) then at time td in the dark after the pulse

Figure 1
Control of hydrogen orbital wavefunctions using superposition states.
An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m4.jpg

where, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m5.jpg, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m6.jpg. We applied the rotating wave approximation and ignored a global phase and all dephasing/relaxation. Figure 1d shows the evolution of |ψright angle bracket in the dark for the case where s=1/[check]2. The fast oscillation at frequency ω is shown across each row on Fig. 1d, and the slower evolution at frequency Δ, from row to row. As is clear from Fig. 1d, the transverse component of the dipole moment of the superposition is An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m7.jpg, where μ0 is a constant, so the intensity of the transverse luminescence at frequency ω oscillates slowly in the manner known as quantum beating22, often used for extraction of Δ through its Fourier transform (FT)23. Greater flexibility in state preparation is provided by a second pulse. If it is identical to the first and arrives at time td, it produces

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m8.jpg

where An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m9.jpg, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m10.jpg, and t is the time in the dark after the second pulse. In this case the quantum beat may be observed in the absorption, which is proportional to

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m11.jpg

As an aside, we note that the usual transverse luminescence observation of quantum beats has the same form as the second term. It allows both extraction of Δ, and may also be used for discrimination between quantum beats and polarisation interference24 between two independent oscillators in situations, where the system is not well understood (not the case for the donors here). Quantum beating where the difference frequency is of the order 1 THz has been observed in semiconductor wells and dots using interband (near-IR) pumps13,25,26,27. Terahertz pulses can also induce other quantum coherent effects such as the Autler–Townes splitting28, and ac Stark effect29. Here, we investigate a GHz beat induced by THz pump pulses, in this case using neutral donor orbital Rydberg states in silicon, thus extending the analogy with quantum beat experiments carried out for free atomic vapours and atoms in vacuum traps.

Our aim is to demonstrate the flexibility in state preparation provided a pair of pulses. The state (equation (3)) is described by the phases An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m12.jpg, Δtd, ωtd, Δt, ωt. Even if Δ and ω are fixed, the fact that An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m13.jpg means we still have near complete freedom to choose An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m14.jpg independent of ωtd and An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m15.jpg independent of ωt. This is sufficient control to take an atom from the 1s ground state to most of the Hilbert space spanned by 1s and 2p± at well-determined times (the value of |ψ0|2 can be varied from 0 to 1, but the maximum |ψ±|2 is 27/32). Complete control may be gained by allowing variable polarisation and amplitude of each pulse in the pair (see Methods).

Theory of donor wavefunctions in silicon

Silicon is an indirect gap semiconductor with 6 anisotropic conduction band minima lying near the X-points at [±kx, 0, 0], [0, ±ky, 0], [0, 0, ±kz], where k≈0.85 2π/a and a is the lattice constant. The KL approach to the silicon donor1,30 is to begin with the solutions of the Schrödinger equation for a single minimum at k=0. The single valley 1s, 2p0 and 2p± state wavefunctions are, to a good approximation, hydrogen-like with extent approximately equal to the effective Bohr radius An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m16.jpg, but display a contraction along the valley axis due to the anisotropy in the effective mass (see Methods for more details). This anisotropy lifts the degeneracy between the 2p0 and 2p±.

The quantum defect of the donor atom and local strains produce a central cell correction (CCC) potential. The CCC makes the 1s state more tightly bound, that is, shrinks its extent. Although its precise form is unknown, the CCC has tetrahedral symmetry, which produces inter-valley mixing of the 1sx,y,z, and in Si:P (and most other shallow donors in silicon) one resulting multi-valley linear combination (called 1s(A1) by KL and referred to here as |0right angle bracket for brevity) has significantly lower energy than the other five:

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m17.jpg

where

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m18.jpg

is the ground state function translated from k=0 to the X-points along the j axis, where j runs over x, y, z. There is a similar set of functions where the crystal part of the envelope is sin(kj), but as they only appear in the higher energy multi-valley components of 1s (called 1s(T2) by KL), we ignore them here.

The excited states that have non-zero dipole moment overlap with |0right angle bracket, that is, those with odd parity, are almost unaffected by the CCC due to its short range. We let

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m19.jpg

for the excited state where i and j both run over x, y, z, so that for example, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m20.jpg is the 2px state in the z-valley and so on. Again, there is a similar set of states with crystal part sin(kj) that are degenerate with our set An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m21.jpg, but they are dipole forbidden for transitions from An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m22.jpg and again we ignore them here.

Since k is large, the crystal part of the wavefunction oscillates rapidly compared with the slowly varying part of the envelope, f j(r), even in the case of the 1s state. This simplifies calculation of matrix elements of slowly varying operators; inter-valley matrix elements are negligible, and intra-valley matrix elements such as the dipole moment matrix elements may be approximated thus:

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m23.jpg

in the example of a transition between 1s and the 2px state in the z-valley and so on. We find d±=0.65a*0. For light polarised with electric field along e=[ex, ey, ez], the multi-valley 2p± state that couples to 1s(A1) is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m24.jpg

as shown in Fig. 2 for e=[001]. The appropriate Rabi frequency for any e is An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m25.jpg, where F is the electric field amplitude of the light wave.

Figure 2
The silicon multi-valley 2p± state.

In the presence of a magnetic field the static Hamiltonian is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m26.jpg

where, H0 is the zero field Hamiltonian that includes the anisotropic mass effects. An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m27.jpg=/2m*=0.30 meV T−1 is the effective Bohr magneton and An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m28.jpg=19.9 meV is the effective Rydberg energy. In our experiment, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m29.jpg/An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m30.jpg (~66 T) so the quadratic terms are negligible. The 2p± excited states have their angular momentum, L, quantized along the valley axis (along left angle bracket001right angle bracket directions) and so for B parallel to [111] the linear B term is the same for all valleys: the 2p+ states in all valleys obtain the same Zeeman energy ħΔ, while the 2p states all obtain energy –ħΔ. Figure 3a shows the calculated level spectrum in the low field limit with B along [111] using the Lanczos method described elsewhere4. These calculations produce Δ/2πB=42 GHz T−1. Other field directions give rise to a more complex spectrum (see for example, ref. 4).

Figure 3
Time evolution of multi-valley donor states in silicon.

For the three levels produced by B along [111], we find (see Methods) that the time dependence of the amplitudes of the states may be taken directly from the preceding hydrogen theory equations (1, 2, 3), and Fig. 3b, which shows the time evolution according the equation (3), is appropriate for either situation. The 1s ground state of hydrogen is replaced by the silicon 1s(A1) state, and the excited states by the silicon equivalents of the 2p+ and 2p. The wavepacket components are

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m31.jpg

where, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m32.jpg. Figure 3c shows the evolution of the density after a pulse polarised along [11–2] appropriate for our experiment. The density has been averaged over a volume (2π/k)3 corresponding to one oscillation period of the crystal part of the wavefunction,

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m33.jpg

Figure 1 shows a hydrogen atom, and the Si:P donor Fig. 3c has a larger scale as already mentioned but is otherwise equivalent—differing essentially by some cubic distortions (the [111] plane shown has hexagonal symmetry) and a relatively more compact ground state wavefunction due to the CCC. The change in ground state size due to the CCC is behind the relatively low contrast in changes along rows in Fig. 3c compared to Fig. 1d.

Evolution of superposition wavefunctions in silicon

The first term in equation (3) is normally very difficult to access, but for Si:P ω is in the THz region of the spectrum, and thus we are easily able to produce control over ωtd with high precision (0.05 radians was our instrument limit) using ordinary stepper motors. We are therefore not restricted to investigation solely of dynamics on the scale of Δtd, and so are able to probe the beat directly via the first term in equation (3).

The parameter-space is large so we restrict our investigation to the td dependence, which exhibits variation in two of the five phases in equation (2). We chose to use weak pulses (Ωtp~π/20) to avoid multi-photon effects, since we are anyway interested in the s2 term in equation (3) rather than the s4 terms. The latter are more susceptible to intensity jitter and were unobserved. Reducing s clearly changes only the amplitudes of the td-dependent terms, equation (3) (Fig. 3b).

Experiment

The transmitted light through the sample was measured as a function of the pulse delay for different magnetic field strengths. The FT of equation (3) (with respect to td) contains a doublet at ω±Δ (as well as a zero frequency term and a term at 2Δ) and Fig. 4 shows the FTs of the experimental interferograms in which the doublet is clearly observed, with Δ/2πB=42 GHz T−1 as expected from the level splitting calculated in Fig. 3b.

Figure 4
Frequency domain quantum beats in transmission.

These experiments are similar to a variety of pulsed-laser spectroscopies where a short optical pulse is used to excite an oscillation, and correlations between copies of these oscillations are measured by taking advantage of interference effects. Our technique is distinct primarily in that the free-electron laser (FEL) source used allows THz pulses with a selective, narrower band corresponding to a few hundred coherent field oscillations. As for other experiments of this type, and in contrast to incoherent frequency domain spectroscopies such as FTIR, the interference is most intuitively understood in terms of the time-domain behaviour: the second pulse interferes with the long-lived excitation produced by the first pulse. The time-domain interference pattern is the autocorrelation of the sample impulse response, so its Fourier transform is its frequency spectrum.

The characteristic resolving power of our time-domain interferometer was tested by measuring the transmission fringes induced by an optically polished slab of high purity float-zone silicon at room temperature, shown in Fig. 5. The results were compared with the transmission measured by high-resolution FTIR with resolution 0.01 cm−1 (0.0003 THz), and excellent agreement was produced.

Figure 5
Calibration of the experimental setup using a Si optical flat.

The features in the transmission are superimposed on a background due to the system response (red data on Fig. 4), which includes the bandwidth-limited laser spectrum31 and an absorption line from some residual water vapour in the beam path, containing a few cm of air between the evacuated optical setup and the cryostat. This means that the raw time-domain data include oscillations with several periodicities and are therefore not interpretable by direct inspection. It is preferable to expose the quantum level splittings by examination of the proportion of excited atoms within the sample rather than the transmitted field pattern at a point far from the sample. This is for the simple reason that the transmitted light spectrum is the product of the transmission spectrum with the system response spectrum (as in Fig. 4), which contains a complicated background, while the product of the absorption spectrum and the system response spectrum is background-free, that is, it contains just the atomic transition features. One method, contactless photothermal ionization spectroscopy (PTIS), relies on a small degree of phonon-induced ionisation from the excited state, which generates a small number of free charge carriers. These conduction electrons may be detected via the AC dielectric properties of the sample32,33. The number of donors ionized is proportional to the average population in the excited states21.

We show in Fig. 6 a set of contactless PTIS interferograms which display a beat in the excited state probability. Again, we kept the intensity low (Ωtp~π/20) to avoid multi-photon effects. When there is no applied field and the excited states are degenerate (Δ=0), a single lobe is observed in the interferogram at zero time delay, which decays on a time-scale determined by the inverse of the inhomogeneous line width. The single lobe becomes a characteristic beat pattern upon the application of a magnetic field. Figure 6d–f also shows the FTs of these interferograms, the peak positions of which correlate directly to the absorptions observed in Fig. 4 and predicted by equation (3). We performed fits in the frequency domain using the FT of equation (3) convolved with a gaussian broadening, and the results are shown in grey. The inverse FTs of the fits are also shown in the time domain. The beat period is clearly controllable with the magnetic field, Fig. 7, hence the beating behaviour is induced by the energy level splitting in the atoms according to equations (2) and (3).

Figure 6
Quantum beats in absorption.
Figure 7
Splitting of the 1s to 2p± transition measured in absorption.

Figure 6a shows the laser autocorrelation interferogram aquired simultaneously with the zero field PTIS. It dies off much more rapidly than the PTIS (in spite of a slight lengthening in the time-domain due to the sharp water vapour line); in other words we have shown that the coherent excitation in the sample is significantly longer lived than the light pulse duration. Our experiment of Fig. 6c demonstrates control over the relative phase and amplitude of the 2p+ and 2p states using a pair of identical linearly x-polarized pulses, which we illustrate by comparison to Fig. 3b. The antinodes of the beat pattern in Fig. 6c are points at which the second pulse arrives when both excited states are in phase or 180° out of phase with each other (first and third rows of Figs 1d and and2a);2a); they are equally doubled (maxima in the rapid oscillations at the average frequency) or destroyed (minima). The nodes in Fig. 6c, on the other hand, are points where the second pulse arrives when both excited states are 90° out of phase (second and fourth rows of Figs 1d and and3a);3a); the second pulse transfers amplitude from one component to the other with no change in the total excited state population.

Discussion

We have used a coherent time-domain analogue of FTIR to examine the magnetic field split 1s–2p± transition in Si:P. We have made a calculation for donor wavefunctions in Si in small B along the [111] crystal axis; simultaneous pumping of two Zeeman split transitions induces beating, thus demonstrating control over a three-state admixture. This control provides an additional tuning parameter, that is, the relative admixture of the three states, for production of Si:donor-based quantum logic gates. Utilization of the orbital motion in the xy plane could have great advantages for surface code networks14,15,16,17. Suppose there are physical (spin) qubits positioned at the red cross positions in Fig. 1, and that we position a fiducial donor qubit of a different species (P in our case) at the centre; the radius of the wavefunction for the spin qubits is much less than the separation between them so as to avoid direct overlap in the ground state. The central donor can be excited to the 2p+ state by a THz pulse, while the surrounding qubits remain unexcited because the associated Rydberg levels are at different energies. The excited donor electron now spends time in the vicinity of the logical qubits during which it experiences exchange interaction, and acquires a phase depending on their state. Projecting the fiducial qubit back onto z at an appropriate time and then observing it is equivalent to evaluating a weighted sum of surrounding qubits, which is the required operation for error detection. The strength of the exchange interaction available in the excited state has been considered by ref. 34. Controlling the orbital trajectory by controlling superpositions of 2p+ and 2p gives additional control, which could allow a bigger variety of fiducial measurement strategies. Finally, we note that our THz quantum beat spectroscopy technique opens up the possibility of high-resolution spectroscopies for GHz splittings of excited states.

Methods

Experimental methods

The magnetic field was provided by a Newport Type C water cooled electromagnet, which was calibrated using a Hall probe. The maximum field measured was 0.285 T at an applied current of 2 A. Our Si:P sample was ordinary commercial float-zone grown silicon (natural isotope composition), with growth direction <110>. The dopant density was 2 × 104 cm−3. It was mounted in a an Oxford Instruments Microstat, in vacuum, cooled by liquid helium to ~10 K. The cryostat containing the sample was placed between the poles of the magnet, such that the sample surface normal was approximately perpendicular to the magnetic field, which was parallel to a <111> direction.

The sample was suspended strain-free, between a pair of capacitor plates formed using pockets of copper foil which did not wrap tightly around the sample, and a PCB. No DC electrical conduction was measurable through the system either at room temperature or under cryogenic conditions. The electronic measurement was obtained by supplying a 1 V peak–peak sine wave (15 MHz) generated using a HAMES 50 MHz HMF2550 arbitrary waveform generator, and measuring the supplied current to the circuit. The current was amplified using a FEMTO DHCPA-100 transimpedance amplifier with 102 V A−1 gain and 200 MHz bandwidth, and measured through a National Instruments PXIe-5162 oscilloscope. Figure 4 shows the r.m.s. current for |B|=0 and 0.285 T, a fuller field dependence (at taken with lower resolution due to laser beam-time restraints) is shown in Fig. 7.

Laser pulses from the Dutch terahertz Free-Electron Laser FELIX were passed through a variable delay Mach–Zender interferometer arrangement. Figure 4e shows an optical layout diagram of the experiment, with the location of each detector and the sample. A pellicle beam-splitter was used to split the FELIX beam into a pair of equal Ramsey pulses, as in ref. 21, and recombined using a pellicle beam-splitter in a Mach–Zender arrangement. The laser beams were linearly polarized perpendicular to the plane of the bench, and normally incident on the sample with electric field perpendicular to the static magnetic field (Voigt geometry). To measure an interferogram, the optical path difference between the pulses was varied using a computer controlled delay stage. At each value of path difference, the electronic waveform from each detector was recorded over one FELIX macropulse (~110 micropulses). Datapoints were obtained from the mean or RMS of the waveform over the stable part of the macropulse, to minimize noise.

The laser was tuned to resonance with the 1s(A1)–2p± orbital transition, that is, optical frequency 9.48 THz. The FELIX pulses had a spectral linewidth FWHM 0.47% and are bandwidth limited, corresponding to a temporal power FWHM of 9.8 ps. The measured micropulse energy at the sample per beam was 30–50 μJ before the cryostat window. The light intensity transmitted through the sample was measured using a liquid helium cooled Ga:Ge detector. The laser autocorrelation spectrum was taken simultaneously using wasted light at the other output of the downstream beam-splitter (Fig. 4e) using a pyroelectric detector. The beam-splitters are not 50:50 and the fringe contrast of the autocorrelation signal is poorer than for the output directed towards the sample, but it has the same temporal profile.

Optical signals were averaged over the FELIX macropulse, and the r.m.s. was measured for the electrically detected signal. The interferograms of Fig. 6 are unprocessed. Further data processing was applied to the interferograms before their Fourier transform was taken to reduce machine precision errors, including: running median background subtraction, normalization and correction for systematic errors along the time axis (on the order of 0.02%). To minimize data collection time, under sampling was employed such that the Nyquist frequency was 1 THz, and high frequency information was extracted by unfolding the spectrum in a standard manner.

The characteristic resolving power of our interferometer was tested by measuring the Fabry–Perot transmission fringes induced by an optically polished slab of high purity float-zone silicon at room temperature. The slab was of thickness 3.05 mm and a calibration standard was set by transmission measurement in a high-resolution FTIR IFS125HR with resolution 0.01 cm−1 (0.0003 THz). We compare the FT of the transmission interferogram of the FELIX pulse through the slab, with the product of the FELIX autocorrelation spectrum and the reference spectrum in Fig. 5. The results compare favourably, so we set an upper bound for the resolution at 0.0146 THz, which is given by the peak spacing of the fringes. The lower bound is 0.0003 THz, limited by the length of the delay line34.

The fidelity of the fringes in Fig. 6 was calculated as (median−minimum)/median in the centre 10% of the interferogram data set. It was 95% for 0.000 T and 82% for 0.285 T. The readout fidelity in these measurements is limited most strongly by the under sampling of the interferogram, and a more densely sampled measurement would produce a stronger fidelity.

Theory of three-level quantum optics

The state vector ψ with components |2pright angle bracket, |1sright angle bracket, |2p+right angle bracket, and with amplitudes ψ, ψ0, ψ+,

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m34.jpg

evolves according to the time-dependent Schrödinger equation Hψ=i[variant Planck's over 2pi][partial differential]ψ/[partial differential]t. The off-diagonal matrix elements of the perturbation for a wave of frequncy ω and phase θ, polarized along x is V=[variant Planck's over 2pi]Ωcos(ωtθ), where the Rabi frequency is Ω=μxFx/[variant Planck's over 2pi], the electric field amplitude is Fx, and μx is the dipole moment, and the Hamiltonian for a basis order of 2p, 1s, 2p+ is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m35.jpg

We make a transformation ψ=Uχ to a rotating frame in which

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m36.jpg

and

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m37.jpg

The time-dependent Schrödinger equation in the rotating frame is [partial differential]χ/[partial differential]t=−iWχ, where

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m38.jpg

after use of the rotating wave approximation. Here, 2ħΔ=E+E is the Zeeman splitting of the excited state. The drive frequency is assumed to be half way between the 1s to 2p± transitions, but for short pulses this restriction is unimportant. The result of integrating [partial differential]χ/[partial differential]t is χfinal=T(t)χinitial where

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m39.jpg

and An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m40.jpg. The result for a delta-function pulse (pulse duration An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m41.jpg) on an atom starting in the ground state is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m42.jpg

where, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m43.jpg, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m44.jpg. In the dark, a time td after the pulse, the result is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m45.jpg

and returning to the laboratory frame, ψ=Uχ produces equation (1) (ignoring the global phase factor in front of U). The effect of a second, identical pulse (arriving at time td) and a further wait in the dark of time t is equation (2). The state is described by the phases Ωtp/An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m46.jpg, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m47.jpg, ωtd, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m48.jpg, ωt. Even if Δ and ω are fixed, the fact that An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m49.jpg means we still have near complete freedom to choose An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m50.jpg independent of ωtd and An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m51.jpg independent of ωt.

Theory of donor wavefunctions in silicon

The KL1 multi-valley wavefunction (excluding the atomic part of the Bloch function) is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m52.jpg

where for valley v, kv is the wavevector at the valley minimum, thus kv.r runs over ±kx, ±ky, ±kz, where k=0.85 2π/a and a is the lattice constant. The single-valley envelope function fv(r) is the same for valleys in opposite directions.

As mentioned in the main text, the single valley odd parity excited state wavefunctions are, to a very good approximation, hydrogen-like with a length scale, L, that is approximately equal to the effective Bohr radius An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m53.jpg, and with a contraction along the valley axis by a factor An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m54.jpg, where Γγ=mt/ml, the ratio of effective masses transverse and longitudinally along the valley axis (Table 1). This anisotropy lifts the degeneracy between the 2p0 and 2p±. The wavefunctions may be obtained from a Lanczos procedure4 or variational calculations1,30. The z-valley variational envelope functions are given in Table 1.

Table 1
KL wavefunctions and parameter.

For a superposition the density is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m55.jpg

where An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m56.jpg is the spatial profile of state An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m57.jpg. Neglecting the inter-valley terms

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m58.jpg

Averaging over a volume (2π/k)3 corresponding to one oscillation period of the crystal part of the wavefunction:

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m59.jpg

The An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m60.jpg ground state contains contributions from all valleys, and in general has non-zero dipole moment matrix element with the excited states in more than one valley. The dipole moment matrix elements are of the formAn external file that holds a picture, illustration, etc.
Object name is ncomms16038-m61.jpg with i, j, α, β running over x, y, z. From symmetry arguments

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m62.jpg
An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m63.jpg

For an electric field F=F [ex, ey, ez], where [ex, ey, ez] is the normalized polarization direction and [ex, ey, ezF=|F|=F, we can construct a normalized excited state An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m64.jpg whose dipole moment matrix element An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m65.jpg is parallel to F. In the case of the 2p0:

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m66.jpg
An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m67.jpg

so the Rabi frequency is Ω=d0F/An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m68.jpg[variant Planck's over 2pi], while for the 2p± it is

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m69.jpg
An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m70.jpg

so the Rabi frequency is An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m71.jpg. It is easy to show that An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m72.jpg for all other excited states An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m73.jpg such that An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m74.jpg, that is, they are not involved in dipole excitation and can be ignored. Figure 2 shows the 2p± excited state superposition appropriate for z-polarized light, that is,:

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m75.jpg

For B parallel to [111], all the 2p+ states in all the valleys obtain the same Zeeman energy ħΔ, while the 2p states all obtain energy –ħΔ, so to find the time evolution we resolve the superposition into its 2p+ and 2p components. We therefore change to a basis of 2p+ and 2p states in each valley:

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m76.jpg

and cyclic permutations of x, y, z, and resolve the excited state into two components |+right angle bracket and |–right angle bracket:

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m77.jpg

At this point, all the three-level quantum optics theory from the preceding section now follows as before, so after a delta-function pulse and a wait in the dark, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m78.jpg, where the amplitudes are given by ψ0=c, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m79.jpg. To find the density, we require the amplitudes An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m80.jpg. It is convenient to change basis again

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m81.jpg

so that An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m82.jpg

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m83.jpg

and now

An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m84.jpg

where, An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m85.jpg. Figure 3a shows the evolution of An external file that holds a picture, illustration, etc.
Object name is ncomms16038-m86.jpg after a pulse polarized along [11–2] appropriate for our experiment.

Data availability

The data for Figs 3a, 4a–d, ,5,5, 6a–f and 7a–d are available online35. Data in more verbose formats are available upon request to the corresponding author.

Additional information

How to cite this article: Chick, S. et al. Coherent superpositions of three states for phosphorous donors in silicon prepared using THz radiation. Nat. Commun. 8, 16038 doi: 10.1038/ncomms16038 (2017).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Peer Review File:

Acknowledgments

We acknowledge financial support from the UK Engineering and Physical Sciences Research Council [COMPASSS/ADDRFSS, Grant No. EP/M009564/1], and from the research programme of the Foundation for ‘Fundamental Research on Matter' (FOM), which is part of the ‘Netherlands Organisation for Scientific Research' (NWO). B.N.M. is grateful for a Royal Society Wolfson Research Merit Award.

Footnotes

The authors declare no competing financial interests.

Author contributions B.N.M. initiated the work; S.C., N.S., K.S., B.R., G.M. and M.N. performed the experiments; S.C. and N.S. analysed the data; B.N.M. and P.T.G. provided theoretical tools; S.C., C.R.P., G.A. and B.N.M. wrote the paper.

References

  • Kohn W. & Luttinger J. M. Theory of donor states in silicon. Phys. Rev. 98, 915–922 (1955).
  • Zwerdling S., Button K. J. & Lax B. Zeeman effect of impurities in silicon. Phys. Rev. 118, 975–986 (1960).
  • Mu Y.-M., Peng J.-P., Liu P.-L. & Shen S.-C. Effects of an external magnetic field on shallow donor levels in semiconductors. Phys. Rev. B 48, 10864–10869 (1993). [PubMed]
  • Murdin B. N. et al. . Si:P as a laboratory analogue for hydrogen on high magnetic field white dwarf stars. Nat. Commun. 4, 1469 (2012). [PubMed]
  • Oberbeck L. et al. . Imaging of buried phosphorus nanostructures in silicon using scanning tunnelling microscopy. Appl. Phys. Lett. 104, 253102 (2014).
  • Veldhorst M. et al. . A two-qubit logic gate in silicon. Nature 526, 410–414 (2015). [PubMed]
  • Dehollain J. P. et al. . Bell’s inequality violation with spins in silicon. Nat. Nanotechnol. 11, 242–246 (2016). [PubMed]
  • Fuechsle M. et al. . A single atom transistor. Nat. Nanotechnol. 7, 242–246 (2012). [PubMed]
  • Jau Y.-Y. et al. . Entangling atomic spins with a Rydberg-dressed spin-flip blockade. Nat. Phys. 12, 71–74 (2016).
  • ten Wolde A., Noordam L. D., Lagendijk A. & van den Heuvell H. B. Observation of radially localized atomic electron wave packets. Phys. Rev. Lett. 61, 2099–2101 (1988). [PubMed]
  • Ludlow A. D., Boyd M. M. & Ye J. Optical atomic clocks. Rev. Mod. Phys. 87, 637–701 (2015).
  • Muralidharan S. et al. . Ultrafast and fault-tolerant quantum communication across long distances. Phys. Rev. Lett. 112, 250501 (2014). [PubMed]
  • Lee J. P. et al. . Ramsey interference in a multilevel quantum system. Phys. Rev. B 93, 085407 (2016).
  • Dennis E., Kitaev A., Landahl A. & Preskill J. Topological quantum memory. J. Math. Phys. 43, 4452–4505 (2002).
  • Fowler A. G., Mariantoni M., Martinis J. M. & Cleland A. N. Surface codes: Towards practical large-scale quantum computation. Phys. Rev. A 86, 032324 (2012).
  • Hill C. D. et al. . A surface code quantum computer in silicon. Sci Adv 1, e1500707 (2015). [PMC free article] [PubMed]
  • Pica G. et al. . Surface code architecture for donors and dots in silicon with imprecise and nonuniform qubit couplings. Phys. Rev. B 93, 035306 (2016).
  • Cole B. E. et al. . Coherent manipulation of semiconductor quantum bits with terahertz radiation. Nature 410, 60–63 (2001). [PubMed]
  • Carter S. G. et al. . Quantum coherence in an optical modulator. Science 310, 651–653 (2005). [PubMed]
  • Greenland P. T. et al. . Coherent control of Rydberg states in silicon. Nature 465, 1057–1062 (2010). [PubMed]
  • Litvinenko K. et al. . Coherent creation and destruction of orbital wavepackets in Si:P with electrical and optical read-out. Nat. Commun. 6, 6549 (2015). [PMC free article] [PubMed]
  • Haroche S., Shimoda K. (eds) Quantum Beats and Time-Resolved Fluorescence Spectroscopy (Springer-Verlag, 1976).
  • Yei W., Sieradzan A. & Havey M. D. Delayed-detection measurement of atomic Na 3p2P3/2 hyperfine structure using polarization quantum-beat spectroscopy. Phys. Rev. A 48, 1909–1915 (1993). [PubMed]
  • Chini M., Wang X., Cheng Y. & Chang Z. Resonance effects and quantum beats in attosecond transient absorption of helium. J. Phys. B At. Mol. Opt. Phys. 47, 124009 (2014).
  • Stolz H., Langer V., Permogorov S. & von der Osten W. Picosecond quantum-beat spectroscopy of bound excitons in CdS. Phys. Rev. Lett. 67, 679–682 (1991). [PubMed]
  • Razbirin B. S. et al. . Bound biexcitons in II-VI semiconductors,’. Solid State Commun. 93, 65–70 (1995).
  • Sekine N., Hirakawa K. & Arakawa Y. Terahertz emission from quantum beats in coupled quantum wells. Jpn J. Appl. Phys. 37, 1643–1645 (1998).
  • Serapiglia G. B. et al. . Laser-induced quantum coherence in a semiconductor quantum well. Phys. Rev. Lett. 84, 1019–1022 (2000). [PubMed]
  • Wagner M., Teich M., Helm M. & Stehr D. Temperature dependence of the intraexcitonic AC Stark effect in semiconductor quantum wells. Appl. Phys. Lett. 100, 051109 (2012).
  • Thomas G. A. et al. . Optical study of interacting donors in semiconductors. Phys. Rev. B 23, 5472–5494 (1981).
  • Knippels G. M. H. et al. . Generation and complete electric-field characterization of intense ultrashort tunable far-infrared laser pulses. Phys. Rev. Lett. 83, 1578–1581 (1999).
  • Andreev B. A. et al. . Contactless photo thermal ionization spectroscopy of shallow defects in semiconductors. Mater. Sci. Forums 143–147, 1365–1370 (1994).
  • Saeedi K. et al. . Room-temperature quantum bit storage exceeding 39 minutes using ionized donors in silicon-28. Science 342, 830–833 (2014). [PubMed]
  • Stoneham A. M., Fisher A. J. & Greenland P. T. Optically driven silicon-based quantum gates with potential for high-temperature operation. J. Phys.: Condens. Matter 15, L447–L451 (2003).
  • Chick S. Supporting Data for ‘Coherent Superpositions Of Three States for Phosphorous Donors in Silicon Prepared using THz Radiation’ [Dataset]. Zenodo. doi:10.5281/zenodo.545834 (2017). [PMC free article] [PubMed]

Articles from Nature Communications are provided here courtesy of Nature Publishing Group