|Home | About | Journals | Submit | Contact Us | Français|
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 10THz and the effective Bohr radius is =3nm. 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 (30T) 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.
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, , 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 (,1/tp, where tp is the pulse length) then at time td in the dark after the pulse
where, , . We applied the rotating wave approximation and ignored a global phase and all dephasing/relaxation. Figure 1d shows the evolution of |ψ in the dark for the case where s=1/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 , 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
where , , 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
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 1THz 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 , Δtd, ωtd, Δt, ωt. Even if Δ and ω are fixed, the fact that means we still have near complete freedom to choose independent of ωtd and 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).
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 , 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 |0 for brevity) has significantly lower energy than the other five:
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 |0, that is, those with odd parity, are almost unaffected by the CCC due to its short range. We let
for the excited state where i and j both run over x, y, z, so that for example, 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 , but they are dipole forbidden for transitions from 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, fj(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:
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
as shown in Fig. 2 for e=. The appropriate Rabi frequency for any e is , where F is the electric field amplitude of the light wave.
In the presence of a magnetic field the static Hamiltonian is
where, H0 is the zero field Hamiltonian that includes the anisotropic mass effects. =eħ/2m*=0.30meVT−1 is the effective Bohr magneton and =19.9meV is the effective Rydberg energy. In our experiment, / (~66T) so the quadratic terms are negligible. The 2p± excited states have their angular momentum, L, quantized along the valley axis (along 001 directions) and so for B parallel to  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  using the Lanczos method described elsewhere4. These calculations produce Δ/2πB=42GHzT−1. Other field directions give rise to a more complex spectrum (see for example, ref. 4).
For the three levels produced by B along , 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
where, . 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,
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  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.
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).
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=42GHzT−1 as expected from the level splitting calculated in Fig. 3b.
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.01cm−1 (0.0003THz), and excellent agreement was produced.
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 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.
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  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 x–y 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.
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.285T at an applied current of 2A. Our Si:P sample was ordinary commercial float-zone grown silicon (natural isotope composition), with growth direction <110>. The dopant density was 2 × 104cm−3. It was mounted in a an Oxford Instruments Microstat, in vacuum, cooled by liquid helium to ~10K. 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 1V peak–peak sine wave (15MHz) generated using a HAMES 50MHz HMF2550 arbitrary waveform generator, and measuring the supplied current to the circuit. The current was amplified using a FEMTO DHCPA-100 transimpedance amplifier with 102VA−1 gain and 200MHz bandwidth, and measured through a National Instruments PXIe-5162 oscilloscope. Figure 4 shows the r.m.s. current for |B|=0 and 0.285T, 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.48THz. The FELIX pulses had a spectral linewidth FWHM 0.47% and are bandwidth limited, corresponding to a temporal power FWHM of 9.8ps. 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 1THz, 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.05mm and a calibration standard was set by transmission measurement in a high-resolution FTIR IFS125HR with resolution 0.01cm−1 (0.0003THz). 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.0146THz, which is given by the peak spacing of the fringes. The lower bound is 0.0003THz, 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.000T and 82% for 0.285T. 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.
The state vector ψ with components |2p−, |1s, |2p+, and with amplitudes ψ−, ψ0, ψ+,
evolves according to the time-dependent Schrödinger equation Hψ=iψ/t. The off-diagonal matrix elements of the perturbation for a wave of frequncy ω and phase θ, polarized along x is V=Ωcos(ωt−θ), where the Rabi frequency is Ω=μxFx/, the electric field amplitude is Fx, and μx is the dipole moment, and the Hamiltonian for a basis order of 2p−, 1s, 2p+ is
We make a transformation ψ=Uχ to a rotating frame in which
The time-dependent Schrödinger equation in the rotating frame is χ/t=−iWχ, where
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 χ/t is χfinal=T(t)χinitial where
and . The result for a delta-function pulse (pulse duration ) on an atom starting in the ground state is
where, , . In the dark, a time td after the pulse, the result is
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/, , ωtd, , ωt. Even if Δ and ω are fixed, the fact that means we still have near complete freedom to choose independent of ωtd and independent of ωt.
The KL1 multi-valley wavefunction (excluding the atomic part of the Bloch function) is
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 , and with a contraction along the valley axis by a factor , 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.
For a superposition the density is
where is the spatial profile of state . Neglecting the inter-valley terms
Averaging over a volume (2π/k)3 corresponding to one oscillation period of the crystal part of the wavefunction:
The 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 form with i, j, α, β running over x, y, z. From symmetry arguments
For an electric field F=F [ex, ey, ez], where [ex, ey, ez] is the normalized polarization direction and [ex, ey, ez]·F=|F|=F, we can construct a normalized excited state whose dipole moment matrix element is parallel to F. In the case of the 2p0:
so the Rabi frequency is Ω=d0F/, while for the 2p± it is
so the Rabi frequency is . It is easy to show that for all other excited states such that , 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,:
For B parallel to , 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:
and cyclic permutations of x, y, z, and resolve the excited state into two components |+ and |–:
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, , where the amplitudes are given by ψ0=c, . To find the density, we require the amplitudes . It is convenient to change basis again
where, . Figure 3a shows the evolution of after a pulse polarized along [11–2] appropriate for our experiment.
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.
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.
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.