Atomtronics is an emerging field of ultra-cold atomic physics that focuses on quantum circuits that operate as atom analogues to electrical circuits, where an atom current takes the place of an electron current and the chemical potential takes the place of a voltage. In their design, atomtronic circuits are drastically different from their electronic circuit counterparts, as they involve manipulating ensembles of ultra-cold atoms to flow through a variety of potential energy landscapes, such as optical potentials produced from laser beams and/or harmonic magnetic potentials, e.g., created with atom chips.
The electronic transistor is the fundamental building block of modern electronic devices, and is ubiquitous in modern electronic systems. An electronic transistor is a semiconductor device used to amplify or electronic signals and electrical power. In a field effect transistor (FET), an electric field (or lack thereof) at a transistor gate can be used to control a current from a source terminal to a drain terminal.
An analogous atomtronic transistor can include a source (chemical potential) well, a gate well, and a drain well. The source well and the gate well are separated by a source-gate barrier, while the gate is separated from the drain by a gate-drain barrier. Atoms with thermal energy in excess of the chemical potentials associated with the source-gate barrier and the gate-drain barrier can flow from source to drain.
Most of the drawings herein have counterparts in the form of illustrations in the 2018 Ph.D. Thesis of Brad A. Dinardo, “Single Atom Delivery into a Bottle Beam Trap Using an Optical Conveyor Belt and an Ultracold Gain Mechanism in a Matterwave Transistor Oscillator”, incorporated by reference herein. This document is available online at https://scholar.colorado.edu/phys_gradetds/. Accordingly, in the listing below, illustration numbers in parentheses following figure numbers identify the illustration number in the thesis that corresponds to the referenced figure in this application. Similarly, formula numbers in parenthesis following formula numbers for the application are indicative of the formulas in the thesis.
The present invention provides a quantum tunneling matterwave transistor in which the flow of particles (e.g., atoms and molecules) through the transistor is a result of resonant tunneling from a source well, through a gate well and into a drain well (as opposed to being a result of collisions, as in a classical atomtronic transistor). The transistor current of matterwave particles can be controlled as a function of the breadth of resonant tunneling conditions of the gate well. For example, the resonant tunneling conditions of a gate well that does not include a dipole-oscillating Bose-Einstein condensate (DOBEC) can be broadened by including a DOBEC in the gate well. Similarly, the breadth of resonant tunneling conditions of the gate well can be changed by changing the particle population of a DOBEC in the gate well.
As shown in
As a comparison of
Herein, ultra-cold means temperatures at or below 120 nanoKelvin (nk). In one embodiment, the mean temperature of the source well particles is initially 78 nK and rises to 109 nK once particles begin to flow into the gate well. The source well temperature then starts cooling as particles flow into the drain well. In this embodiment, the flux is not a Bose-Einstein condensate (BEC). In another embodiment, the initial temperature of the source well is 54 nK and the resulting flux is or includes a BEC.
The wells and barriers of transistor 100 are formed using: 1) optical potentials produced from laser beams; and 2) harmonic magnetic potentials created using an atom chip. As shown in
Particle source 412 introduces atoms or other particles into cooling chamber. Laser system 404 and magnetic system 406 cooperate to produce two-dimensional and three-dimensional magneto-optical traps used to cool the atoms in cooling chamber 414. Electronics 410 provides radio-frequency signals used to modulate laser beams to implement evaporative cooling in the 3D MOT to achieve ultra-cold temperatures. Laser system 404 cooperates with harmonic magnetic potentials generated by atom chip 418 to produce the chemical potentials that define transistor 100. Electronics 408 can modulate laser beams so that the positions of the source-gate and gate-drain barriers can be adjusted and varied.
A matter-wave transistor process 500 is flow charted in
At 502, in effect, the source well is populated with ultra-cold particles. In practice, the ultra-cold particles can be loaded before the source-drain barrier is instantiated. The source-drain barrier can then be formed over the eventual position of the gate well, the gate-drain barrier, or the drain well. The source-gate drain barrier is then swept into its final position, moving particles on its way so that they end up densely packed against the source-gate barrier by the time it reaches its eventual position.
At 503, while the transistor is in a first state in which the gate well is empty, that is, does not contain a DOBEC, the source-to-drain matter-wave flux is at a relatively low level, e.g., zero or relatively small.
At 504, the transistor is switched to a second state in which the gate well is populated with a DOBEC. Due to the degeneracy between eigenstate transitions, this DOBEC has the effect of broadening the resonance boundary conditions of the gate, increasing the likelihood of source-particle energies that can tunnel into the gate well via the source-drain barrier and tunnel out of the gate well and into the drain well via the gate-drain barrier. This increase in likelihood is represented in the graph of
At 506, changes in DOBEC numerosity (number of particles constituting the DOBEC) positively correlate with changes in flux. In other words, flux positively correlates with DOBEC numerosity. Thus, the addition of particles to the DOBEC results in increased flux.
This description details: 1) transporting and loading single neutral cesium atoms into optical bottle beam traps in order to facilitate a continuously operating quantum computer; and 2) developing a purely quantum thermodynamic model and experiment for a matter-wave transistor oscillator that demonstrates the existence of a matter-wave gain mechanism.
In the first section below, a single atom can be loaded from a reservoir into a blue-detuned crossed vortex bottle beam trap using a dynamic one-dimensional optical lattice. The lattice beams are frequency chirped using acousto-optic modulators, which causes the lattice to move along its axial direction and behave like an optical conveyor belt. A stationary lattice is initially loaded with approximately 6000 atoms from a reservoir, and the conveyor belt transports them 1.1 millimeters (mm) from the reservoir to a bottle beam trap, where a single atom is loaded via light-assisted collisions. Photon counting data confirm that an atom can be delivered and loaded into the bottle beam trap 13.1% of the time.
Further below, the second section addresses the field of atomtronics and aims to demonstrate the existence of a matterwave gain mechanism in an atomtronic transistor oscillator for the eventual use in ultracold atomic experiments. An ultracold gain mechanism for a matter-wave transistor oscillator is described. The description starts with a well-established semi-classical model of an atomtronic transistor and model the system with a purely quantum mechanical formalism. The quantum model predicts dynamics differently from the semi-classical results only when the atoms owing through the transistor oscillator have sufficiently low enough temperatures such that the motional state of a DOBEC, placed in the transistor itself, couples atom transitions between high lying transistor energy eigenstates. In this ultra-cold temperature regime, interesting physics not observed in the semi-classical theory or experiments includes the presence of a gain mechanism that broadens the resonant tunneling conditions of a matter-wave current flowing from the source-well, through the gate-well, and into the drain-well. Experiments indicate that the transistor gain mechanism is coherent and increases the range of matter-wave energies (and by extension, momenta) that flow into the transistor drain-well.
First proposed by Richard Feynman at Caltech in 1959, quantum computing is a general, all-encompassing term that describes storing information in a two-state quantum system in order to perform algorithms analogous to how classical computers perform algorithms using binary transistor states, or bits. The two-state quantum system, called a qubit, could be the spin of an atomic ground state, a ground state and Rydberg state atom pair, electronic states of an ion, a Josephson junction, or a SQUID to name a few examples. Irrespective of the type of qubit, the ambition of quantum computing is to perform algorithms exponentially faster than classical computers. A major milestone in the field was achieved in 1995, when quantum computing was first demonstrated by using a single trapped ion to serve as a qubit to perform quantum logic. Since then, multi-qubit quantum logic gate operations have been performed experimentally with multiple trapped ions, numerous superconductors, collections of photonic qubits, and neutral Rydberg atom pairs. Just like a classical computer, the number of qubits in a quantum computer directly influences its productivity. Most recently, using superconducting qubits, IBM has steadily increased their commercially available quantum computer from two qubits in 2016, to sixteen qubits in 2017, and finally to using a 20-qubit register in 2018.
Further progress calls for increasing numbers of qubits. Large scale quantum computers that would compute with tens to hundreds of qubits would exceed the capabilities of classical computers by being able to solve integer factorization, such as using Shor's algorithm, or inverting functions by using Grover's algorithm, faster than any classical computer. While experiments with trapped ions and superconducting qubits have already succeeded in realizing high fidelity many-qubit quantum gates, using an array of neutral atom qubits also remains an attractive and viable solution to scaling a qubit register to an arbitrarily large amount.
As shown in
Herein, qubit gates are based on exciting neutral atoms to Rydberg states. Using neutral atoms, one can theoretically scale the amount of qubits N to as high as one desires with very minimal crosstalk between the remaining N−1 qubits. This is because neutral atoms differ from ions with respect to their state dependent interaction properties. This is evident in the graph 800 of
First, ions separated by a distance r interact predominantly with the 1/r Coulomb potential. Because the Coulomb potential only decreases with a single power of r, it remains strongly interacting, even at long distances, which are great attributes for qubits. The drawback is the Coulomb interaction is always on, making it difficult to create an arbitrarily large register of ion qubits. Next, the interaction between two ground state neutral atoms is dominated by the 1/r6 van der Waals potential at short distances and by the 1/r3 magnetic dipole-dipole interaction at large distances. With a separation of 10 microns (m), this 1/r6 ground state interaction is very weak at approximately U/ℏ≈10−5 Hz in frequency units. However, at this same separation, two Rydberg state atoms have a very strong r3 electric dipole-dipole interaction of approximately 107 Hz.
This result gives the core principle of operation of the neutral atom qubit array: selectively turning on the Rydberg interaction (by exciting an atom to such a state) allows one to selectively increase an atom-atom interaction by 12 orders of magnitude. The ability to control the interaction (2 Here, long distances refers to any atom separation r 30 nm.) strength by a factor of 1012 over a wide spatial range (e.g., >30 nanometers (nm)) appears to be unique to the Rydberg system and makes an array of Rydberg-mediated quantum gates a viable solution to a many-qubit quantum computer. In fact, recently, a controlled-NOT (CNOT) gate has been demonstrated between two Rydberg state qubits. The measured truth table 900 for this CNOT gate is shown in
A drawback of using neutral atom qubits, as opposed to ion qubits and superconducting qubits, is that neutral atoms must be held in optical traps enclosed inside ultra-high vacuum chambers. This results in unavoidable atom losses due to finite lifetimes on all of the neutral atom qubits trapped in the array shown in
A solution to atom loss includes transporting Atoms to Reload Distant Unoccupied Array Optical Traps. Site selective atom reloading of a qubit array can be accomplished by using a dynamic, one-dimensional optical lattice (also called an optical conveyor belt) to deliver a single atom to an unoccupied array site. The concept of such an atom reloading scheme is shown in
Moreover, the atom transport system can be scalable such that it is capable of addressing and delivering to an array of sites and not just to a single site. Accordingly, the atom transport system must be “steerable” in that it can deliver an atom to one site and then be repositioned to deliver an atom to another array site. Finally, the qubit array in the AQuA-49 project is located in a miniaturized and portable vacuum chamber which makes optical isolation an intricate problem. Any stray light from the atom transportation lasers can destroy the qubit array coherence and will disrupt any ongoing quantum computations. Thus, to transport atoms from a reservoir to an empty qubit array site without disrupting any occupied array sites, the issue of optical isolation presents a compelling problem.
A system 1000 providing for such an atom reloading is shown in
System 1000 includes a vacuum cell 1002 with a hexagonal cross section. A qubit array 1004 is centrally located within the cell. A face 1006 is used for transmitting laser beams used for readout, cooling, optical pumping, AC Stark compensation, Raman lasers, and excitation lasers used for Rydberg atoms. An opposing face 1008 is used for counter-propagating beams used for readout, cooling, and Rydberg excitation. Faces 1010 and 1012 are used to admit counter-propagating beams that interfere to produce an optical conveyor belt for transporting ultra-cold atoms from a reservoir 1014 to vacated trap sites of array 1004. Face 1016 is used to admit 532 nm laser light for forming bottle beams to help guide transported atoms into their target trap sites. Face 1018 is used as an exit for readout of array 1004. A dichroic mirror 1020, external to cell 1002, reflects the light used to form the bottle beams to remove that light from the readout and direct it to a 532 nm beam dump 1022. As shown in the detail to the right in
System 1000,
Starting from the Hamiltonian of a two-level atom interacting with a laser field leads to forces that give rise to optical dipole trapping, magneto-optical trapping, and polarization gradient cooling; all of which are absolutely essential to the experiments in this dissertation and can cool an initially room temperature ensemble of 10 atoms to about 10 to 20 μK in about one second. The resulting pre-cooled atoms are subsequently trapped in magnetic potentials that are produced on an atom chip, which allows further cooling with forced radio-frequency evaporation to produce a Bose-Einstein condensate with a temperature as low as 50 nK.
Consider a two-level atom with a ground state |g and an excited state |e, with energies Eg and Ee, respectively. By placing the atom in an oscillating electric field E produced by a laser beam with frequency) the atoms will gain an induced electric dipole moment −{circumflex over (d)}. In such a configuration, the atom experiences an energy shift given by the interaction Hamiltonian
H′=−d·E (1)(2.1)
Because H′ is a perturbation to the energy of the atom, the associated change in the atom's ground state energy Eg is given by standard perturbation theory as
where the first term is the energy shift due to the atom absorbing a photon of energy ℏω from the laser and the second term represents the energy change due to the atom emitting a photon of energy ℏω. In both the single atom loading and transistor oscillator experiments described in this dissertation, the cooling laser frequency ω is very close to resonance. That is Ee−E≈ℏω. Thus, the absorption contribution to the ground state energy shift (equation 2) is negligible compared to the emission contribution. Using this approximation, the ground state energy shift for an atom in the laser field can be approximated as
where {circumflex over (∈)} is a unit vector indicating the direction of the electric field with magnitude E produced by the laser.
The ground state energy shift (equation 3) can also be examined using classical electrodynamics. Using that the induced dipole moment is proportional to the electric field by means of the atom polarizability, the shift to the ground state energy of the atom in the laser beam can also be expressed as
where the relation {circumflex over (d)}=α(ω)E was used with the atomic polarizability α(ω)) defined as a function of the frequency ω of the laser electric field E. By equating equations 3 and 4, one finds that the atomic polarizability for a two-level atom in the laser electric field is
In its current state, the result of equation 5 intrinsically assumes an infinite lifetime for an atom in the excited state |e. The infinite lifetime arises because the energy Ec of the excited state |e is purely real. Thus, the probability of finding the atom in the excited state is a stationary state and will never decay, since
This, of course, is false since an excited state decays exponentially with time constant τ=1/Γe, where Γe is the linewidth of the transition |e →|g. The polarizability in equation 5 can be corrected to reflect this finite lifetime by including the complex energy shift of iℏΓe/2 due to transitions to the ground state. That is Ee→Ee−iℏΓe/2=ΔEe. By including the imaginary energy shift iℏΓe/2 to the excited state energy, the probability amplitude correctly reflects a decay with time constant τ=1/Γe due to the finite lifetime 1/Γe of the excited state.
Applying this energy shift in equation 5, the atomic polarizability becomes:
where the laser detuning is δ=ωe−ωg−ω with ωc=Ee/ℏ and ωg=Eg/ℏbeing the energy eigenvalues of the |e and |g states, respectively, as expressed in frequency units.
This argument is also extended to the ground state |g where in a very general case, the ground state energy shift ΔEg in a laser beam is also written as a complex quantity
ΔEg=Vg−iℏΓg/2. (9)(2.7)
The real component Vg represents the energy shift of the atom ground state due to the laser electric field E and the imaginary component iℏΓR accounts for the finite lifetime τ=1/Γg of |g due to transitions |g →|e that are induced by the atom absorbing the laser radiation field. Equation 8 can be inserted into equation 4 and equate the result to equation 9 to formulate an expression representing the ground state energy shift of an atom due to a laser field:
where we have introduced the Rabi frequency Ω=|g|{circumflex over (d)} ·E|e|/ℏ, which is merely the magnitude of the original perturbation H′=−d·E (equation 1) expressed in frequency units. Equation 10 is a very important result for the laser cooling and trapping used herein. By equating the purely real components, equation 10 gives the atom energy shift due to scattering photons from the laser electric field and directly leads to the optical dipole force that is used to trap and repel atoms. Additionally, by equating the purely imaginary components of equation 10, one arrives at the atom energy shift due to absorbing photons from the laser field and leads immediately to the radiation force, which permits the cooling of atoms in an optical molasses and leads to more effective cooling by means of a magneto-optical trap.
An expression for the effective potential Vg can be attained by equating the purely real components of both sides of equation 10:
This effective potential Vg is called the optical dipole potential. Equation 11 can be rewritten to make it clear that Vg permits atom trapping at the maxima of an intense laser beam:
where c is the velocity of light, ω0=ωe−ωg, and ∈0 is the permittivity of free space. Some identities can be used to further simplify this result. The intensity I(x) of an electromagnetic wave with amplitude E can be expressed as I=c∈0R2. Additionally, we use that linewidth of the excited state |e can be alternatively expressed as
Γe=ω03|g|{circumflex over (d)}·{circumflex over (∈)}|e|2/3ℏc3π∈0 (13)
Finally, the lasers used for dipole trapping in the experiments described herein are detuned far from resonance; hence δ2>>Γe2. Examples described further below use 1064 nm and 760 nm lasers for dipole trapping and repelling of 133Cs and 87Rb atoms, respectively. That corresponds to detunings of order 1012 Hz for 133Cs and 109 Hz for 87Rb. For both types of atoms, the linewidths are of order 106 Hz, thus validating δ2>>Γe2. Using these identities and the approximation δ2>>Γe2, the optical dipole potential imparted onto an atom by a far-detuned laser beam can be written as:
Equation 14 shows that the optical dipole potential can be attractive (and trap atoms) or repulsive (and repel atoms) depending on the sign of the laser detuning. As shown in
The result Fdip 1120 of equation 15 is called the optical dipole force and is a very useful tool to trap (or repel) neutral atoms that have been previously cooled to milliKelvin temperatures or lower. The experimental consequences of equation 15 are used for facilitating the optical trapping and transport of single neutral 133 Cs atoms and for constructing and preparing atoms in a matter-wave transistor oscillator potential.
Equating the purely imaginary components of equation 10 yields an expression for Γg, or the rate that atoms absorb a photon from the laser beam:
where δ is the detuning of the laser beam from resonance. Now, for simplicity, consider a two-level atom 1202 placed in a one-dimensional light field produced by a single pair of counter propagating laser beams 1204 and 1206 along the {circumflex over (x)} direction as shown in
which is called the radiation force. Equation 17 (2.14) is only valid when the atom in the laser field is at rest. Atoms move with some non-zero velocity. An atom moving along the +{circumflex over (x)} direction with velocity ν observes the frequency of each laser beam (with wavenumber k) to be Doppler shifted by an amount k·v. (This Doppler shift of the laser beams is observed by the atom in its rest frame.) Thus, the radiation forces on an atom moving at velocity ν due to a pair of counter-propagating laser beams is
where the subscripts ±x denote the force from each laser beam propagating in the ±x direction. For a red-detuned laser beam (δ>0), equation 18 (2.15) imparts a velocity dependent, frictional force on the atoms with the dependence on δ shown in
Physically, the radiation force occurs due to the Doppler effect: the atom absorbs (and scatters) more photons from the laser that propagates against the atom's motion because this light, if red-detuned by some amount δ, is subsequently blue-shifted in the atom frame back onto (or near) resonance with the |g →|e transition. Consequently, the atoms absorb far less light from the other laser beam (which propagates in the same direction of the atom) since in the atom frame this laser is Doppler shifted out of resonance. By red-detuning the lasers to the optimal δ (
While the radiation force given in equation 18 (2.15) has the capabilities of cooling atoms to a few hundred microKelvin (μK), atom trapping with it is unstable as the radiation force Frad is divergenceless. Because the electric-field intensity (encapsulated in the Rabi frequency) from the laser beam is Gaussian, the radiation force can be described more generally in cylindrical coordinates as:
and taking the divergence yields
Because ∇Frad=0, the radiation force cannot possibly be center-seeking at all points in space and consequently fails as a stable means to trap atoms. The force given in equation 18 (2.15) can be made stable by introducing a positionally dependent magnetic force in addition to the radiation force by placing the entire atom-radiation field system into an inhomogeneous magnetic field.
In the one-dimensional MOT 1300 of
Consider a one-dimensional case where an atom is placed into a magnetic quadrupole field generated by anti-Helmholtz coils
ΔE=gμbm B′x (22)(2.16)
where g is the Lande g factor, μb is the Bohr magneton, and m is the magnetic quantum number of the atom. The important result of equation 22 (2.16) is that magnetic field establishes a positionally dependent energy shift for the atom. Given an atom moving with velocity ν in an inhomogeneous magnetic field, the atomic transition |g →|e experiences an effective Doppler shift (in frequency units) of Δνdop=±kv as well as an additional Zeeman shift (in frequency units) of Δνzee=g μb×B′/ℏ. Therefore, when inserting values for the total detuning δ for the radiation force Frad given in equation 18 (2.15), the effects of both the Doppler shift of the laser light and the Zeeman shift of the atomic transition. In this case, the total radiation force on an atom becomes
With the result of equation 23 (2.17), the radiation force Frad is now proportional to both the velocity and position of the atom. In this instance, the divergence of Frad is non-zero and in cases where k2>>g μb B′/ℏ (which is common in cold atom experiments), the motion of an atom subjected this force is that of an overdamped harmonic oscillator (
To ensure that atoms remain trapped, the radiation force must be directed towards the B-field minimum (located at x=0). This can occur by giving the counter-propagating lasers orthogonal, circular polarizations σ+ and σ− (
Using the laser cooling methods, a MOT of approximately 109 133Cs or 87Rb atoms can be created but only cooled to their respective Doppler temperature limits of 125 μK and 145.57 μK. To further cool the atoms to between 10-20 μK, kinetic energy can be continuously removed from the atoms by utilizing polarization gradient cooling (PGC).
Real atoms have multiple ground states. As a simple example, consider an atom with a triplet ground state (Jg=1) coupled to a quintuplet excited state (Je=2) in a region with no external magnetic field. Since the cooling lasers are red-detuned from resonance, atoms absorb more light from the laser beam they are propagating against. Thus, in the one-dimensional setup of
The opposite effect occurs when the atoms move with negative velocity ν<0, illustrated in
Laser cooling and trapping methods are used to cool ensembles of both 133Cs and 87Rb to about 10-20 μK. This temperature is sufficiently cold enough for the IARPA single atom transport and loading, but is still too warm by orders of magnitude for an ultra-cold matter-wave transistor. Ensembles of 87Rb atoms are cooled, forced radio-frequency (RF) evaporative cooling, to approximately 50 nK, which is beyond the capabilities of the laser cooling methods. Furthermore, the phase space density of the atoms must be high enough such that the 87Rb atoms transition into a Bose-Einstein condensate (BEC). To reach BEC, pre-cooled atoms (by means of PGC described above) are confined to a harmonic magnetic potential produced on an atom chip and then cooled to degeneracy by performing forced radio-frequency (RF) evaporative cooling.
Consider a neutral atom placed in a magnetic field. The potential 9 energy between the atom's magnetic moment y and the external magnetic field B is given by the Zeeman effect as
U=−{circumflex over (μ)}·B=mƒgFμBB (24)(2.18)
where mƒ is the z-component of the atomic angular momentum (or magnetic quantum number), gF is the Landé-g factor, μB is the Bohr magneton, the Bohr magneton is μB=eℏ/2 me where e and me are the charge and mass of the electron. While this interaction is very weak, magnetically trapping atoms provides a very efficient method to compress and manipulate previously cooled atoms into a Bose-Einstein condensate. If the magnetic field is spatially varying, then the atom experiences a force given by the gradient of equation 24 (2.18)
F=−∇U=∇({circumflex over (μ)}·B)=−mƒgƒμB∇B (25)(2.19)
The magnetic force can be attractive (thus strong-field seeking) or repulsive (weak-field seeking) depending on the sign of mf gF. This is very important for carrying out magnetic trapping experiments as all magnetic fields are divergenceless (∇·B=0), which means they do not feature a local maximum. This renders strong-field seeking magnetic potentials useless as a means to magnetically trap neutral atoms. Thus, the potential in equation 24 2.18 must be made to be a repulsive, weak-field seeking potential. This can be achieved by optically pumping the atoms into a magnetic sublevel where mf gF<0. With U=−{circumflex over (μ)}·B<0 fulfilled, atoms can be magnetically trapped but only if the external magnetic field B features a local minimum.
An ensemble of atoms at a temperature of approximately 10−20 μK can be prepared in a harmonic magnetic trap. To further cool the atoms such that their phase space density is high enough to form a BEC, a process called “radio frequency (RF) evaporation” can be used to continuously remove the warmest atoms from the magnetic trap, leaving only the “coldest of the cold” to remain in the magnetic potential. An ensemble of 87Rb atoms can be optically pumped to a magnetic trapping sublevel, such as |F=3, mƒ=2. To lower the temperature of this trapped atom ensemble, RF radiation with frequency ωRF irradiates the atoms. Since the magnetic field strengths are small, the energy difference between adjacent Zeeman sublevels in the magnetic potential is approximately ΔE=gF μBB/ℏ, RF photons that are tuned to the frequency
ωRF=gFμBB/ℏ (26)(2.20)
induce RF spin ips where atoms undergo multi-photon transitions from the |m=+2 trapped state to the |mf=−2 anti-trapped state and are subsequently scattered from the trap.
Consequently, by sweeping the RF frequency ωRF from infinity down to some arbitrary ωflip, atoms with energy above ℏωflip are pumped from the |mF=+2 →|mF=−2 magnetic sublevel and subsequently expelled from the trap. This process permits selectively removes warmer atoms above a certain energy in the magnetic potential without changing the tightness of the trap, allowing for efficient rethermalization. The expelled atoms carry away energy, leaving the atom ensemble at a reduced temperature. By continuously ramping ωRF lower and lower, forced RF evaporation can maintain a high cooling rate by regularly removing the hottest atoms in the trap as long as the RF sweeping speed is slow enough to allow the remaining cooler atoms time to rethermalize. RF evaporation and rethermalization is repeated until the atom phase space density is high enough for a BEC. RF evaporation to cool 20 μK atoms to energies 10-20 kHz above the trap bottom, which then condense into a 50 nK BEC. The phenomenon of runaway evaporation, which is the maximum efficiency of RF evaporative cooling, occurs when the time per collision, τcol decreases as the atom ensemble cools. Runaway evaporation emerges when the trap tightness and RF sweeping rate are such that τtrap=τt is about equal to or greater than 100 where τtrap is the 1/e lifetime of the magnetic trap.
An RF evaporation process is illustrated in
RF evaporation is the final stage to condensing atoms into a BEC. For the transition from a gas to a BEC to occur during evaporation, the de Broglie wavelength
of each atom with momentum p and temperature T must become larger than the inter-atomic spacing of the gas. That is, the phase space density, PSD must be greater than unity
PSD≡λ3n≥1 (28)(2.22)
where n is the number of atoms per unit volume. Equation 2.21 reveals that as the gas is cooled, the de Broglie wavelength of each atom increases, eventually reaching a critical temperature Tc whereby the de Broglie waves of each atom begin to overlap. By cooling the gas even further to the ultra-cold regime (T<<Tc), the de Broglie wavelength of each atom becomes so large that all of the atom de Broglie waves completely overlap, rendering each atom of the gas indistinguishable. At this temperature, nearly the entire gas has condensed into a BEC that must be treated quantum mechanically and is described en masse according to the Bose-Einstein distribution. (For temperatures where the de Broglie wavelengths of each atom are much smaller than the inter-atomic spacing d between each atom, the particles are distinguishable from each other and hence the gas behaves classically and is described with Maxwell-Boltzmann statistics).
Consider an ensemble of atoms in a state of thermal equilibrium at temperature T and chemical potential p in a magnetic potential with the simplification that the atoms are not interacting with each other. The average occupancy
Moreover, in the matter-wave transistor detailed further below, the cold atom gas is confined to an anisotropic harmonic oscillator potential given as
V(x,y,z)=½m(wx2x2+ωy2y2+ωz2z2) (30)(2.24)
where m is the mass of each atom and the ωi represents the atom oscillation frequencies along the i={circumflex over (x)}, ŷ, {circumflex over (z)}, axis of the potential. The corresponding energy levels of the anisotropic potential are
where nx, ny and nz are integers. Each atom thus has an energy ∈x, ∈y and ∈z due to oscillation along the respective axis. The energy of this system (equation 31 (2.25)) is plotted in
The integral η(∈x), and hence the number of available energy states in the anisotropic potential, is merely the volume of first octant of
It immediately follows that the density of states, or the number of available energy states between ∈ and ∈+d∈ is:
Knowing the density of states (equation 33 (2.27)), as well as the Bose-Einstein distribution (equation 29 (2.23)), one can readily determine the experimental conditions on atom temperature T and trapped atom number N where all atoms of the gas are in an excited state (meaning the ensemble temperature is exactly at criticality Tc). For this simplified case of a non-interacting gas, the number of particles in the excited state Nex reaches a maximum when the chemical potential μ=0. That is
where Γ(n) and ι(n) are the Gamma and Riemann zeta functions, respectively. Rearranging equation 34 (2.28) and solving for Tc gives the critical temperature
where for T>Tc, all atoms in the gas are in an excited state. It immediately follows that if the atom ensemble is cooled lower than Tc, atoms begin to abruptly occupy the anisotropic harmonic potential ground state, which is the onset of Bose-Einstein condensation. By cooling the atom temperature T<<Tc, one achieves a macroscopic occupation of the harmonic oscillator ground state and can produce a pure BEC.
In the matter-wave transistor oscillator, interactions between low energy, ultra-cold atoms are not only omnipresent, but they drive the matter-wave gain mechanism. Specifically, the interactions that we must address are binary elastic scattering events between pairs of atoms with very small total energy. By “small” is meant that the scattering process is dominated by the s-wave contribution to the wavefunction: that is components of each atom's wavefunction with orbital angular momentum l=0 are considered. A s a result, the total scattering cross section a drastically simplifies to
where δl is the phase shift experienced by the l=0 partial wave due to the scattering event.
For s-wave scattering of identical bosons, the scattering cross section is σ=8πa2 where a is the scattering length given in the Born approximation as
where mr is the reduced mass of the two atoms, V(r) is the effective two-body interaction and where k and k′ are the wave vectors of the atom before scattering and after scattering, respectively. Because the interactions here are due to s-wave scattering, the momentum transfer ℏk′−ℏk due to the scattering process is approximately zero. Thus, equation 37 (2.30) can be approximated as
The effective two-body interaction V (r; r′) between identical atoms at points r′ and r in equation 38 (2.31) is treated as a discrete contact interaction
V(r′,r)=gδ(r′−r) (39)(2.32)
where g denotes the energy of the interaction (defined shortly) and δ(r′−r) is the Dirac delta function.
Inserting equation 39 (2.32) into equation 38 (2.31), one arrives at the scattering length between two atoms in the BEC as
Using equation 40 (2.33), the result, important for transistor oscillator theory, that the effective interaction between two colliding atoms at low energy in the BEC is constant and is given by
Given that there are N interacting atoms in the BEC, the BEC is described with a many bodied Hamiltonian
Using this many-bodied Hamiltonian in the Schrödinger equation yields the Gross-Pitaevskii equation
where the eigenvalue 1@=@E/@N is the BEC chemical potential and Ã(r) is the total wavefunction for a system of N indistinguishable bosons making up the BEC, which is given as
ψ(r)=ψ1(r1)ψ2(r2) . . . ψN(rN) (44)(2.37)
where ri refers to the position of the ith boson in the BEC. The Gross-Pitaevskii equation is a form of the Schrödinger equation that consists of a linear contribution from the external potential (r) confining the atoms and a non-linear contribution g|Ã(r)|2 modeling the averaged, or mean field, interaction with all of the other atoms in the BEC.
At equilibrium, the mean-field pseudopotential and the external potential are both proportional to R2, while the kinetic energy is proportional to 1/R2. Additionally, the mean-field increases linearly with atom number N. Thus, the kinetic energy of atoms in the BEC is negligible compared to V and g.
An important approximation can be made to equation 43 (2.36) when the kinetic energy of the atoms is much less than the mean-field. In such cases, the kinetic energy operators −ℏ2∇i2/2 m are ignored, which reduces the Gross-Pitaevskii equation to the much simpler
[V(r)+g|ψ(r)|2]ψ(r)=μψ(r). (4s)(2.38)
This result is called the Thomas-Fermi approximation and equation 45 (2.38) has the solution
Equation 46(2.39) gives the atom density profile of a BEC. It states that in the Thomas-Fermi approximation, atoms condensed into a BEC fill the external trapping potential V(r) uniformly up to the condensate chemical potential.
The external trapping potential is harmonic along the {circumflex over (x)}, ŷ, and {circumflex over (z)} axes. From equation 46 (2.39), the atomic density n of the BEC is an inverted parabola 1800 (
The number of atoms in the condensate can be found by inserting the harmonic potential into the solution to the Gross-Pitaevskii equation 46 (equation 2.39) giving
Integration over the radial coordinate from the origin to the Thomas-Fermi radius R gives the total atom number in the BEC as a function of the chemical potential 1.
where in the final line of equation 49 (2.42), the definitions of the Thomas-Fermi radius (equation 47 (2.40)) and the scattering length (equation 40 (2.33)) are used. Rearranging this result and solving for yields the Thomas-Fermi approximation for the BEC chemical potential:
Vacuum chambers described herein can be compact (small enough to carry on a commercial airplane) and provide UHV (pressures below 10−10 Torr) environments suitable for MOTs and BECs. A vacuum chamber for atom transport minimizes unwanted magnetic fields that could introduce noise by containing no metal1 nuts, bolts, or screws and being completely held together with optical and anodic bonds.
Transport and delivery of single atoms into optical bottle beam traps is conducted in a hex-cell vacuum system 1900, shown in
To achieve the above listed performances, the hex-cell system design is based on of a double MOT vacuum configuration featuring a 133Cs source region, a 2D MOT region, a 3D MOT region, and a hexagonal chamber featuring six windows with anti-reflection (AR) coating on both ambient and vacuum faces. The hex cell was fabricated with no metal nuts, bolts, or screws and is entirely held together by over 50 optical contact bonds and anodic bonds. Even the ion pumps themselves are anodically bonded to the cell. Using a much larger vacuum pump in conjunction with a bake-out process, the hex cell can be pumped down to UHV pressures.
Cell manifold 1908 serves as a mechanical foundation for hex-cell vacuum system 1900. Manifold 1908 is a 1 cm thick, machined Pyrex™ block that is polished on all faces. Channels are machined throughout the interior of manifold 1908 to connect the 3D MOT chamber 1906 to Hex cell 1910 and also to connect ion pumps 1917 to their respective pumping loads. Two Pyrex™ support rods 1920 are fixed to the bottom of the 2D MOT chamber 1902 and are connected to manifold 1908. A vacuum seal 1922 is completed by anodically bonding a Pyrex window in-vacuo (while still connected to a bake-out station) onto manifold 1908.
As shown in
Shown in
Anodically bonded directly above the cesium source tube 1902 is the 2D MOT chamber 1904, which is a 1 cm×1 cm×4.5 cm Pyrex vacuum chamber designed to support a 2D+ MOT from the background 133Cs gas produced in the source. Because this chamber is filled with warm, cesium vapor, the pressure in this chamber is higher than UHV at approximately 10-100 nanoTorr. 2D MOT chamber 1904 is capped with a 500 1m thick silicon wafer that features a 750 1m diameter pinhole drilled in the center. A 2D+ MOT 2032 (
Anodically bonded to the top of cylinder 2014 is the 3D MOT chamber 1906, which has the same dimensions as 2D MOT chamber 1904 and is also made of Pyrex. In order to have MOT lifetimes of 10-20 seconds, the 3D MOT chamber is kept at a lower pressure of approximately 0.1 nanoTorr, which is orders of magnitude lower than the adjacent connecting 2D MOT chamber pressure. Sufficient isolation between the high and low pressure MOT chambers is provided by the 2D-3D pinhole, as the 0.024 mL/sec conductance of pinhole 1916 is negligible compared to the 0.4 L/sec pumping rate of the ion pumps described further below.
Because of the very low background pressure of 0.1 nanoTorr in the 3D chamber, a 3D MOT cannot be loaded from a background gas. Instead, a stream of pre-cooled 2D+ MOT atoms are pushed through the 2D-3D pinhole 1916 and into 3D MOT chamber 1906, where the atoms are used to source the 3D MOT. This design of using a 2D MOT to directly source a 3D MOT in an entirely separate vacuum chamber is called a double-MOT cell and it improves the 3D MOT lifetime by decoupling the atom lifetime (which is set by the 3D chamber vapor pressure) from the MOT load time (which is set by the 2D chamber vapor pressure). Similar to 2D MOT chamber 1904, 3D MOT chamber 1906 is capped with a silicon disc but features a larger pinhole (aperture 1918) with a 3.18 mm diameter. This pinhole is larger so that the 3D MOT itself can be transported vertically into Hex cell 1910, located above 31) MOT chamber. 1906. This entire source-2D-3D MOT vacuum ensemble is then anodically bonded to manifold 1908. The magnetic fields that are required to make the 2D+ and 3D MOTs are produced with permanent magnets 2042 (
Located above the 3D MOT chamber is the six sided Hex Cell 1910 (3.3), which provides an UHV environment of less than 0.1 nanoTorr with six faces of optical access. The Hex cell walls are AR coated on all inside and outside faces for the 532, 780, 852, 1038, and 1064 nm wavelength lasers used herein. Hex cell 1910 is optically contact bonded to the top of the Pyrex manifold 1908 (
System 1900 features two 0.4 L/sec ion pumps 1917, one of which shown in
Body 2308 of each ion pump 1917 is a machined Pyrex block that features ta vertically oriented bore hole 2316 and a horizontally oriented bore hole 2318. A titanium cylinder 2020, which serves as an anode 2022, extends most of the way through vertically oriented bore hole 2216. A conductive rod 2324 is placed inside horizontally oriented hole 2318 and is welded to titanium cylinder 2320. This rod 2324 extends all the way to the very left-most side of ion pump 1917. Horizontal bore hole 2318 is sealed at the distal end of rod 2324 with a metalized conductive silicon cap 2326. Electrical contact between cap 2326 and rod 2324 is made with a metal spring to form an anode feedthrough.
Metalized silicon cathodes 2302 and 2314 respectively cap the top and the bottom of ion pump 1917. Small Pyrex rings anodically respectively bonded to each cathode 2302, 2314. Placed in between the Pyrex rings and the cathodes are titanium discs 2304, 2312 that serve as cathode feedthroughs. Titanium disks 2304 and 2312 are slightly bowed, giving the discs spring-like properties. Connectivity between the anode and cathode is provided by the spring force present in the titanium discs. Due to its slight radius of curvature, the disc is constantly pushing on both the cathode and the titanium cylinder, which maintains constant electrical connectivity. Neodymium magnets 2328 and 2330 are mounted, respectively, on top and bottom cathodes 2302 and 2314 to produce the magnetic field required for ion pump 1917 to function. A custom Mu Metal box placed over the ion pumps provides magnetic shielding. The bottoms of ion pumps 1917 are anodically bonded directly to Pyrex manifold 1908 (
An atomtronics vacuum system 2400 is shown in
Bases of the 2D and 3D MOT chambers 2404 and 2408 are attached to ConFlat flanges via glass-to-metal anodic bonds 2412 and 2414. Atomtronics system 2400 is then held together by bolting each MOT chamber 2404, 2408 onto a spherical cube 2416. Isolation between the 2D and 3D MOT chambers 2404 and 2408 is provided by a 750 1m pinhole 2418, which is drilled into a silicon wafer placed at the top of the 2D MOT chamber 2404. The pinhole permits differential pumping between the 2D and 3D chambers, resulting in a background pressure 0.1 nanoTorr in 3D chamber 2408.
An Agilent 2 L/sec ion pump is also bolted to the spherical cube 2416 to maintain UHV. When bolting each chamber or ion pump onto the spherical cube 2416, a copper gasket is placed between the base of the ConFlat and the spherical cube. A vacuum seal is established when knife edges, found at the base of the ConFlat, as well as on the spherical cube, are tightly pressed into the copper gaskets, creating a knife-edge seal.
Atomtronics system 2400 is initially pumped down to UHV by bolting the entire vacuum system to a much larger vacuum station via a copper pinch off tube 2418 and following a standard bake-out procedure at approximately 300° C. Following bake-out, the atomtronics system was separated from the bake-out station by “pinching” pinch off tube 2418 with a set of hydraulic jaws. After removal, the system is soft baked for 5-7 days at 120° C. to remove any water moisture that erroneously entered the cell 2400 during pinch-off.
As best seen in
Formation of harmonic potentials for the transistor uses guide wires 2406 and H-wires 2408 positioned at the centrally located chip window. A single pair of parallel guide-wires 2606 is shown in the detail of
Consider running a current I, through one of the chip guide-wires 2406. Compared to the size of the atoms as well as the dimensions of a BEC, the guide-wire is considered to be infinitely long. (Herein, typical BEC dimensions are 50-100 1m in the longitudinal direction.) From Ampere's law, the current I, running through the guide-wire produces a magnetic “guide-wire” field
with corresponding field gradient B′ and curvature B″ of
where 10 is the permeability of free space and r is the distance from the guide wire. The corresponding field lines for the magnetic field 2710 generated using the guide wire are shown in
Since the magnetic field 2730 is symmetric along the axis of the guide-wire (in and out of the page with
By choosing to orient the guide-wire along the {circumflex over (x)} axis, the local minima of the quadrupole field is found by simply taking the limit of equation 54 as y→0 and z→d. In this limit, the guide-wire fields can be approximated:
With equations 55a and 55b, one of the requirements for magnetic trapping is satisfied: a localized minimum for a magnetic field has been produced. However, the quadrupole fields from equations 55a and 55b produced by the atom chip can be problematic for atom trapping and reaching BEC. Plotting the fields By and Bz (
Using the magnetic field from equation 56 3.6, we can form a magnetic potential that can confine atoms along the radial direction but cannot trap atoms along the longitudinal direction, since ∇xB=0. To confine atoms along the {circumflex over (x)}, or longitudinal direction, a current IH can be run through a pair of H-wire 2602 (
By combining the magnetic fields from the guide-wire (equation 56 (3.6)) and the H-wire (equation 3.7), the total field forms a three-dimensional harmonic magnetic potential capable of trapping atoms in a local field minimum. This configuration is called an Ioffe-Pritchard trap. Accordingly,
The magnetic fields used to radially confine ultra-cold atoms (equation 56 (3.6)) was formed by running current Ix through a single guide-wire. This places the magnetic trap a distance d (equation 53 (3.3)) below the wire, which blocks optical access. In order to gain optical access to ultra-cold atoms trapped in the magnetic potential, the trap must be located out from under the guide wire. This can be accomplished by adiabatically ramping up an oppositely flowing current through the second guide-wire (running parallel to the first, separated by distance l;
and is plotted in
Similar to the single guide-wire trap, a local minimum in equation 58 (3.8) can be formed by applying a bias field, this time oriented in the {circumflex over (z)} direction with magnitude Bz (
Physically, the quantity D represents the z value where the z-bias field, βz, cancels the {circumflex over (z)} component of the split-guide wire field (equation 58 (3.8)) along the entire line y=0 (
The local minima of the total magnetic field that forms the split guide-wire Ioffe-Pritchard trap is found by taking the limit of equation 60 (3.10) as y→0 and z→d, giving
Recall from equation 59 (3.9) and the sentence that follows it, that B, vanishes when z=D. Therefore, equation 61a (3.11a) is not used in the magnitude of the total field at the location z=D. From equations 61b (3.11b) and 61c (3.11c), the magnitude of the total field at the location z=D is
While this expression gives the total B-field located a distance D under the chip window, it can be more useful to also only consider the field at small y values, since, after all, atoms are localized to y≈0 when trapped in the field minimum. Therefore, the field is expanded in a Maclaurin series
The trapping forces on an atom in this spatially varying magnetic field are readily calculated by evaluating the first and second spatial derivatives of equation 63 (3.13), giving corresponding field gradient B′ and curvature B″ of
where the definition of D given in equation 59 3.9 is used in solving |B″|. Using equations 64a (3.14a) and 64b (3.14b), the radial trapping force imparted onto atoms in the split guide wire trap is:
which features a trapping frequency in the radial direction of
The trapping potential along the longitudinal direction is still provided by the same H-wire field given in equation 57 (3.7). Thus the trapping force and corresponding trapping frequencies in the longitudinal direction are found by evaluating the first and second spatial derivatives of equation 57 3.7 and are
Advances in neutral atom quantum computing have led to the development and demonstration of a 49-element atomic qubit array, where quantum information is stored in the ground states of neutral cesium atoms and two-qubit gates implemented with Rydberg states. As shown in
Due to collisions with thermal background atoms 3202,
Initially, a few thousand atoms are loaded into a red-detuned 1D optical lattice from a background cold atom reservoir. Precise frequency control over the lattice beams causes the lattice to move back-and-forth along its axial direction, making it function as an “optical conveyor belt.” This optical conveyor belt can be used to transport the small atom ensemble approximately 1.1 mm to an empty optical BoB trap followed by transferring exactly one atom at most into the BoB trap.
The single atom transport and reloading experiment requires the use of three sets of lasers with various detunings from the 133Cs D2 line (λ=852.347 nm, ν=351.725 THz): 1) three 852 nm lasers that are on resonance with the 133Cs D2 line; 2) one 780 nm laser for BoB traps that is blue-detuned from resonance; and 3) one 1064 nm laser for atom transport that is red-detuned from resonance.
852 nm resonant laser light can be produced using Vescent Photonics D2-100-DBR-852 distributed Bragg reactor (DBR) diode lasers. Atom transport and reloading uses three individual DBRs: one to provide cooling light, one to provide repump light, and one to act as a probe laser to perform absorption imaging. In addition, a New Focus VAMP TA-7616 tapered amplifier (TA) can be used to amplify the light emitted from the cooling DBR from 70 mW to 500 mW of total usable laser power. Blue detuned light at 780 nm, used to form optical bottle beam traps and bottle beam arrays, is produced with a Sanyo DL 7140-201 S laser diode and is amplified by a Toptica BoosTA Pro to 3 W of laser power. Finally, the red detuned 1064 nm laser that is used for optical dipole trapping is a 30 W, IPG Photonics YLR-30-1064-LP-SF laser. A block diagram of the 852 nm, 780 nm, and 1064 nm laser systems are shown in
The 852 nm cooling, repump, and probe laser setup 3300 (
To keep atoms contained within the cycling transition, repump laser 3304 is locked resonant to the |F=3 →|F′=4 transition
Finally, the probe laser 3402 (
The setup 3600 for the 1064 nm atom transport laser 3602 is shown in
Δν1=2×78=156 MHz and the other by a variable amount of
Δν2=2×(78+s)=156+26 MHz, where the amount δ is tunable. To minimize phase jitter between the two beams, both AWGs are frequency locked and phase synchronized via their internal 10 MHz clocks. The frequency shifted light emerging from the double pass AOMs is then combined and spatially overlapped at a final PBS. This setup produces two co-propagating 1064 nm beams with orthogonal linear polarizations with a mutual frequency detuning of Δν=2δ. The pair of beams is then sent to the Hex-cell vacuum 1900 (
Shown in
The objective of the crossed vortex BoB setup is to generate two co-propagating Laguerre-Gaussian (LG01) beams with equal intensity and orthogonal polarizations. Shown in
As shown in
For every complete 2π rotation of the azimuthal angle ϕ of the TEM00 mode, the total phase delay Δϕ imparted on the light with wavelength λ passing through SPP 3714 is
where d and n are the thickness and refractive index of SPP 3714, respectively. As a result of the phase delay, SPP 3714 transforms the laser electric field from a Gaussian mode into a Laguerre-Gaussian mode LG, described by
where p is the radial index, l is the topological charge, Lpl(x) are the associated Laguerre polynomials, and R(z) is the radius of curvature of the beam, given as R(z)=z+zR2/z. The phase profile of SPP 3714 is shown in 38A (4.4a) and it was carefully manufactured to impart a topological charge of exactly l=+1 on the incoming TEM00 mode of 780 nm light, which transforms it into an LG01 mode beam in the far-field (
where P is the total laser power in the LG01 mode. The LG01 beam is then split into two co-propagating beams separated by 2.5 mm with equal intensity and orthogonal polarizations by passing the beam through a calcite beam displacer 3716. Upon exit from the calcite, the two LG01 beams pass through a 2× telescope designed to increase both the beam separation to 5 mm and the beam waist to 1.61 mm. This configuration of two co-propagating “vortex” LG01 beams with orthogonal polarization is then sent to the hex cell 1910 (
Returning to the 780 nm laser system shown in
A set of telescopes 3724, 3726, and 3728 and additional calcite beam displacers 3722 (
After passing through the calcite beam displacer 3722, the array of eight Gaussian beams is further de-magnified by a 200:60 telescope 3724 and then demagnified again with a 500:23.125 telescope 3726
The 852 nm cooling and MOT lasers 3302 and 3304 (
BoB traps 4012 are formed at the focus of the 780 nm BoB lasers 4022 and are offset from the 1064 nm atom transport lasers (4024 functioning as an optical conveyor belt) by an angle of θ=60°. Additionally, the cooling lasers 4020 that form a 3D MOT intersect at an angle of θ=60° (as opposed to the traditional 90° angle from most MOT setups).
A 3D magneto-optical trap of cesium atoms located in the hex cell provides a continuous reservoir of cold atoms for replenishing the qubit array. Atoms that source the 3D MOT are initially trapped and cooled in a 2D+ MOT located in a chamber below the Hex cell and then sent upwards into the 3D MOT chamber via a “push beam”.
An optical setup 4100, shown in
A second branch from PBS 4106 is split further by a PBS 4108 after passing through a half-wave plate λ/2. A first branch from PBS 4108 passes through a half-wave plate λ/2 and a quarter-wave plate λ/4 before entering hex-cell vacuum chamber 4000. The beam exits chamber 4000 only to be reflected back into chamber 4000. passing through a quarter-wave plate λ/4 coming and going. A second branch from PBS 4108 passes through a half-wave plate λ/2 and a quarter-wave plate λ/4 before entering hex-cell vacuum chamber 4000. This branch then exits chamber 4000 and is then reflected back into chamber 4000.
The probe beam exits fiber 1104, is shaped by a lens with a 50 mm focal length, and rotated by a quarter-wave plate λ/4 before entering chamber 4000 to acquire an image of bottle-beam array 4012 (
The inset for
The final orientation of the 3D MOT beams as they appear in the actual experiment in the lab is shown in the inset of
With the laser setup completed, a pair of coils 4130, where each is made with 100 turns of 20 gauge wire, is oriented along the {circumflex over (z)} direction of the hex cell 1910. Running 1 A of current through the coils when the coil pair is oriented in the anti-Helmholtz configuration generates a positionally dependent magnetic field B(x) with gradient magnitude of |B′|=16 G/cm. A flux of atoms is sent into Hex cell 1910 from the 2D+ MOT formed in a chamber 1906 below Hex cell 1910 and sources the 3D MOT 4302. Using this described 2D+/3D double MOT setup, a 3D MOT of approximately 10′ atoms is created within 1 second of loading time.
An optical setup 4200 to form a 2D+ MOT in the 2D chamber is shown in
A beam exiting PBS 4208 is expanded using a 30 mm lens and a 75 mm lens. The result is polarization rotated using a half-wave plate λ/2 and then split using a PBS 4210. Both outputs are processed by respective quarter-wave plates λ/4 and directed through adjacent walls of 2D MOT chamber 1904. Each of the beams passes through chamber 1904, processed by a respective quarter-wave plate λ/4, reflected back through the respective quarter-wave plate λ/4 on its way back into chamber 1904. The two direct σ+ beams and their respective reflection beams σ− are shown with the 2D MOT chamber 1904 in the upper inset of
A second 2D cooling beam output from PBS 4206 passes through a 100 mm lens and then a 30 mm lens to be used as a push beam to move atoms into and through 2D MOT chamber 1904, as indicated in the lower inset of
Shown in
A telescope made from a pair of cylindrical lenses then enlarges the beam and gives it an elliptical profile with a beam waist of ω=0.5 mm along the semi-minor axis and a waist of ω=1.0 mm along the semi-major axis. Using a PBS 4210, the light is then split into two beams with equal intensity. A pair of λ/4 waveplates for each beam controls the circular polarization such that all beams overlap at an angle of 90° at the center of the 2D chamber (upper inset of
Additionally, four permanent magnets located outside the 2D chamber create a 2D quadrupole field with gradient 15 G/cm. A 2D+ MOT is then created by red-detuning the 2D cooling light by δ=2:9Γ. The resulting MOT force provides strong confinement of atoms in the {circumflex over (x)} and ŷ directions, but no confinement in the {circumflex over (z)} direction. A vertically oriented push beam overlaps the 2D MOT along the loose {circumflex over (z)} axis (inset 2 of
After loading 3D MOT 4302, the cooling beams can be further detuned to −20Γ which serves to act as an effective compressed MOT (CMOT) stage, since only atoms at the most intense region of the MOT beams will consistently absorb the laser light. After the CMOT stage, time-of-flight absorption imaging can be performed to measure the temperature of the atoms to be Tx=80.1 μK and Ty=67.7 μK along the {circumflex over (x)} and ŷ directions, respectively, giving an average 3D MOT temperature of T=73.9 μK.
Polarization gradient cooling (PGC) is used to obtain maximum single atom loading probability and transport efficiency in the BoB trap and optical conveyor belt (after forming the 3D MOT to decrease the MOT temperature and increase the MOT atom density). To perform PGC, the MOT coils 4128 (
The coil system for generating the larger magnetic quadrupole field for the 3D MOT and the considerably smaller shim fields needed for PGC can provide shim coils for all six hex cell windows in order to bias fields along the optical axis of all cooling beams. Moreover, the system must also provide shim coils for the {circumflex over (z)} direction. The coil system includes a set of anti-Helmholtz coils for producing the 3D MOT magnetic fields. The mounting scheme includes a “pagoda” coil mount 4702 as shown in
Once each MOT coil 4704 and each shim coil 4706 is wrapped on its respective face of the pagoda mount, each coil can be permanently fastened to the mount using Araldite epoxy. The finished coil system can then be mounted over the hex cell as shown in
The merged output from PBS 4812 is split by PBS 4814 into first and second branches, each of which is input to a respective one of opposing faces of hex-cell vacuum chamber 4000, resulting in counter-propagating beams that interfere to form one-dimensional lattice 4840. The two branches differ in frequency by a variable amount δ. Varying δ causes interference fringes from the two branches to move. Atoms trapped by the fringes move with the fringes, so the interference pattern serves as an optical conveyor belt for the atoms.
The first branch from PBS 4814 passes through a 300 mm lens on its way into chamber 4000, while the second branch passes through a half-wave plate λ/2 and a 300 mm lens on its way into chamber 4000. Light exiting chamber 4000 returns to PBS 4814, which diverts some of the returned light to a photodiode 4850 via another 300 mm lens. A resulting electric signal from photodiode 4850 is conducted over a conductor 4852 to a radio-frequency spectrum analyzer 4854, which measures the velocity of the optical conveyor belt 4840 by measuring Δω of the lattice beams.
Optical conveyor belt 4840 is a 1D optical lattice potential that can move atoms back and forth along the axial direction of the lattice. The optical conveyor belt is created inside the hex cell by overlapping two P=5 W, λ=1064 nm, counter-propagating Gaussian laser beams 4860 and 4862 with frequencies ω1 and ω2 and equal intensities of I=3.75×101 mW/cm2. The resulting standing wave potential is
where ω0 is the laser beam waist size, ω(z) is the Gaussian beam waist, k=2π/λ is the wavenumber and Δω=ωi−ω2 is the detuning between the two 1064 nm beams.
The quantity U0 in equation 72 is the maximum dipole potential and is defined as
where ˜ℏ=1.05×10−34, J is the reduced Planck's constant, P and I0 are the total laser power and saturation intensity at the atoms, Δ=70 THz is the detuning between the 1064 nm beam and the cesium D2 transition frequency, and δ=2π×5.22 MHz is the natural linewidth of the cesium D2 transition. Both λ=1064 nm lasers have an equal power of P=5 and intensity of I=3.75×108 mW/cm2. From equation 734 4.5, focusing the beams to a waist of ω0=30 μm at the center of the hex cell produces a 1D lattice with a maximum trapping depth of U0=kB=−1.88 mK. where kB, =1.38×1023, J/K is the Boltzmann constant. The axial and radial harmonic oscillation frequencies ∩i and Ωrad of atoms trapped at the focus of the 1D lattice are
where the numerical values are computed using the values for U0 and ω0 listed above.
The velocity ν of the lattice potential U(r,z,t) along the axial direction of the lattice is controlled by changing the lattice beam detuning Δω. To understand why detuning the lattice beams induces motion of the lattice, realize that there exists a reference frame moving along the axial direction of the lattice at a velocity ν where both lattice beams are Doppler shifted to the same frequency, that is Δω=0 in this frame. From the lab frame of reference, this is consistent with the lattice propagating along the axial direction with the same velocity ν. Consequently, atoms that are trapped in this 1D lattice potential U(r,z,t) are also translated with the same velocity ν. The velocity of the 1D lattice as a function of laser detuning is
The value of Δω, and hence ν is set by passing each lattice beam through an acousto-optic modulator (AOM) 4806, 4808 initially driven at 78 MHz in a double-pass configuration as shown in
It is important to note that the lattice detuning Δω, and hence conveyor belt velocity ν, cannot just be “snapped on” instantaneously, as the resulting acceleration can heat the atoms out of the lattice. Rather, motion of the lattice is achieved by adiabatically frequency chirping one lattice beam by Δω. Frequency chirping the beams at an adiabatic rate is important since the total potential experienced by atoms trapped in the optical conveyor belt becomes a tilted lattice potential. In the atom reference frame, the complete optical conveyor belt potential along the direction of the lattice acceleration becomes
where a is the acceleration of the atoms.
The total potential is tilted upward or downward when the respective conveyor belt acceleration is positive or negative. During any atom transport sequence, as long as the 1D lattice potential is greater than the maz contributions from atom accelerations, then atoms will remain trapped in the optical conveyor belt. To solve for the maximum permitted acceleration in the tilted, moving potential, the local minimum of equation 77 (4.9) as a function of acceleration a is computed by equating the spatial derivative of Utot to zero:
and solving for a
Since the maximum value of cos(kz) sin(kz)=0.5, it immediately follows that the maximum of equation 79 4.11, and hence maximum optical conveyor belt acceleration amax, is
where the numerical value is computed using the values for U0, k, and m. With the result in equation 80 (4.12), there is an absolute speed limit for how fast the atoms can be transported using the optical conveyor belt. However, this does not tell us how to accelerate the atoms. During transport, a minimal amount of heating can occur when frequency chirping one of the lattice beams in a fashion that minimizes the rate of change of acceleration of the atoms.
Using equation 76 4.8, the velocity of the optical conveyor belt as a function of the mutual lattice beam detuning is
where ω=2πc/λ is the initial, un-shifted frequency of the λ=1064 nm transport laser. It immediately follows that the time rate of change of acceleration of the optical conveyor belt is expressed as
The lattice chirping Δω can be a linear ramp. Thus, the d2Δω/dt2 term vanishes. Additionally, the lattice axial trap frequency Ω2>>Δω, and thus we can approximate that Δω/Ω3≈0. Using these approximations, the magnitude of the time rate of change of the conveyor belt acceleration is
In order to minimize this function, the lattice chirping must satisfy
Equation 84 4.16 can be used to define a dimensionless adiabatic chirping parameter
whereby ξ<<1 must be satisfied to transport atoms with minimal heating imparted from changes in acceleration of the optical conveyor belt. The desired adiabaticity parameter is ξ≈10−3.
A composite time-lapse image of an entire transport process 5100 is shown in
As the optical conveyor belt motion abruptly starts and stops, large forces are imparted to atoms trapped in the 1D lattice. Consequently, cold atoms in the conveyor belt can be heated and even ejected from the lattice. This heating effect is achieved by chirping the moving lattice at a much lower rate than the square of the lattice axial oscillation frequency. For motion of the optical conveyor belt the inequality holds as Ωz2˜1012 s−2 and dΔω/dt˜1018 s−2. Additionally, this transport process preserves approximately 87% of the atoms trapped in the optical conveyor belt. This atom survival probability is approximately equal to the atom survival in the stationary lattice over the same time period.
Because the process of transferring a single atom from the conveyor belt and loading it into a BoB trap relies heavily on the atom density in the lattice, 400 kHz frequency chirping is used as it permitted the farthest and fastest transport with the highest remaining atom density after transport. Faster transport over greater distances can be achieved by using higher frequency chirps and longer transportation times. While greater transport distances of up to 5.32 mm can be achieved in this way, the more intense accelerations resulted in the atom survival rate rapidly declining to below 50 percent. More importantly, for these more intense accelerations, the remaining atom density after transport drops far below the remaining atom density of the stationary lattice over the same time durations, showing that heating due to these greater accelerations is outweighing the natural heating rate of the lattice.
Successful atom transport with the optical conveyor belt depends critically on the degree of overlap between the two 1064 nm lattice beams and their relative polarization. The amplitude of the beat note Δω of the 1D lattice beams is a maximum when the beams are entirely overlapped with parallel polarization. Moreover, the line width of the beat note is used to quantify the phase jitter, and hence the axial position fluctuations of the conveyor belt. The beat note is measured by focusing the overlapped 1D lattice beams onto a photodetector that is connected to RF spectrum analyzer 4804 (
As shown in
Atoms in chamber 4000 are subjected to 852 nm molasses light, that causes 852 nm atom fluorescence, some which is directed on a return path through lens 5410, dichroic 5408, a 60 mm lens, and an 852 nm filter 5420 to be collected by a 10 mm core single-mode fiber 5424 for transfer to a photon counter 5426.
Additionally, to avoid beam distortions on the 780 nm BoB when passing through the 7.3 mm thick cell windows, the front face of 0.4 NA lens 5410 must be aligned normal to the cell walls. Otherwise, the BoB does not completely close in the axial direction, giving atoms an opportunistic path to escape the trap. To have precise control over the position of 0.4 NA lens 5410, the lens is mounted with a compact lens positioner that provides adjustments along the x; y; and z axes in addition to control over θx and θy.
To deliver and load an atom into the BoB trap, the foci of the 780 nm and 1064 nm beams must overlap. As shown in
Performing absorption imaging provides visual proof of atoms ejecting from the trap and approximates the degree of overlap between the conveyor belt and BoB trap. Course adjustments of the conveyor belt position can be made until only atoms at the focus of the conveyor belt are observed to be ejected. Fine alignment of the conveyor belt position was made by turning on the 3D MOT lasers after transporting atoms to the BoB site and counting the photons emitted from the atoms while still held in the lattice. Micron-scale adjustments of the conveyor belt position can be made to optimize the fluorescence signal. With the photon counting signal at a maximum, there exists adequate overlap between the foci of the BoB trap, photon detection optics, and conveyor belt lasers. However, this degree of overlap contains the largest source of error for this technique, due to the lack of a distinct, physical reference point for the microscopic BoB trap.
An optical bottle beam (BoB) is an electric field distribution that features a region of low optical intensity completely surrounded by high optical intensity. Herein, the laser light used to form the BoB intensity profiles has wavelength λ=780 nm and is blue-detuned from resonance. From equations 14 (2.11) and 15 (2.12), by using blue-detuned light, the regions of the BoB with high intensity produce a repulsive potential, forcing the atoms to seek regions of low optical intensity. Thus, atoms can become trapped in the central, dark region of the BoB. While a red-detuned dipole trap is considerably easier to implement for trapping atoms, the upside of the blue-detuned BoB trapping scheme is that it is very useful for confining single cesium atoms for use in quantum computing. Atoms in the dark, blue-detuned traps subsequently scatter less photons than if they were trapped in a red-detuned dipole trap of comparable depth. Since the qubit decoherence time is proportional to the photon scattering rate of the trapped qubit, trapping in blue-detuned BoBs can result in coherence times 300 times longer than those observed in red-detuned dipole traps with similar trapping depths
The intensity profiles of the BoBs can generate radial confinement (along {circumflex over (x)} and ŷ) but not axial confinement along the direction of beam propagation {circumflex over (z)}. To provide confinement along the axial direction, additional blue-detuned light must close the ends of the BoB (i.e., “cap” the bottle beam). A BoB trapping potential can be formed using a high numerical aperture (0.4 NA) lens to tightly cross two co-propagating Laguerre-Gaussian LG01 beams with orthogonal linear polarizations at a tight focus. Because the LG01 beams are sometimes called “vortex beams” this optical configuration gives the trap the appropriate name of “Crossed Vortex Bottle Beam Trap.”
The complete, detailed schematic for generating the co-propagating pair of orthogonally polarized LG01 beams with identical beam waists of w=1:61 mm and separated by distance d=5 mm is shown back in
As the beams are focused with lens 5410, they are rotated by
Because the rotation is small (θ<10°) the changes in polarization to each LG01 beam can be neglected. The resulting crossed vortex intensity profile ICV is also shown in
ICV(r,z)=Il,p(r+,z+)+Il,p(r−,z−)=I1,0(r+,z+)+I1,0(r−,z−) (87)(4.19)
where the intensity of each LG beam, Il,p(r,z) is defined in equation 71 (4.3). For the intensities in equation 4.19, the radial coordinates rx are defined as
r±=√{square root over (y2+(x cos ↓±z sin θ)2)} (88)(4.20)
and the axial coordinates z± are
z±=z cos θ∓x sin θ. (89)(4.21)
where the ± subscripts refer to the π, and π linearly polarized LG01 beams, respectively.
The total optical trapping potential of the crossed vortex BoB trap is then numerically calculated by inserting equation 87 (4.1) into equation 14 (2.11). The calculated trapping profiles along the orthogonal xz 5602, yz 5604, and xy 5606 planes at the focus of the crossed vortex BoB are plotted in
Finally, in order to get an intuitive understanding of how the final size of the crossed vortex BoB affects the trapping depth, the equations for the trapping potential can be vastly simplified by choosing to perform a Taylor series expansion of the intensity (equation 87 4.19) about the origin. Such expansion of the trap depth about (x=0; y=0; z=0) gives
Equations 90A (4.22a), 90B (4.22b), and 90C (4.22c) convey a very important result that the crossed vortex BoB trap potential scales as ω−4, so meticulous care must be made when focusing the LG01 beams though the thick hex cell walls to ensure that as little distortions as possible occur such that the final beam waist of ω0=5.5 μm is correct.
Fluorescence imaging is used to detect the presence of a single atom loaded into the crossed vortex BoB trap. When a single atom is loaded, a 3D optical molasses is applied to the trapped atom by turning on the 3D MOT lasers and detuning them −30 MHz. Repump light is also turned on. Additionally, the MOT coils are turned off and shim coils are used to eliminate any residual magnetic fields. This induces atom fluorescence while continuously cooling the trapped single atom.
As shown in
The SPCM 5426 works alongside an avalanche photodiode to count individual photons emitted from the trapped atom over an arbitrary time interval. When exposed to optical molasses 5412, the trapped atom scatters on the order of 10′ photons over a 50 ms interrogation time. Since the fluorescence collecting lens has an effective focal length of 23.125 mm and an aperture of 18 mm, the collection angle is Ω/4π=4.67%. Total transmission of the photon collecting optics is approximately 50% and the quantum efficiency of SPCM 5406 at 852 nm is 54%. Thus, the maximum photon collection efficiency is approximately 1:1%. Additionally, noise of approximately 100 counts/sec of 780 nm light and approximately 25 dark counts/sec is recorded using the photon counter.
To verify that the BoB is capable of trapping a single atom, the BoB trap can be loaded directly from a 3D MOT. The MOT was loaded for 1 second followed by PGC for 10 ms. During the final 5 ms of PGC, the BoB is snapped on in 300 μs and overlaps the central and most dense region of the 3D MOT. If multiple atoms are initially loaded into the BoB trap, a useful process between pairs of atoms called light assisted collisions reduces the trap occupation to zero or one atom.
As atoms confined in the BoB approach each other in the presence of the MOT laser light, there exists a critical inter-atomic separation where the laser light excites a long-range, attractive dipole-dipole interaction between pairs of atoms, which results in the atom pair gaining considerable kinetic energy leading to collisions that expel both atoms from the BoB trap. Consequently, if an even number of atoms occupy the BoB trap, they will all pair-collide until no atoms remain in the trap, whereas an odd number of atoms will pair-collide until only one atom remains in the trap.
After loading the BoB trap and allowing time for light assisted collisions, all MOT lasers are turned off for 20 ms to let any untrapped atoms fall away due to gravity. After this fall time, the MOT beams are turned back on and fluorescence imaging on the trapped single atom is performed. A histogram 5700 of the photon counting data recorded by the SPCM 5426 during fluorescence imaging of the BoB trap is shown in
The photon counting statistics for 1002 independent atom loading attempts show two Poisson distributions corresponding to either 0 or 1 atom confined in the BoB trap after any given loading. When 0 atoms are present in the BoB trap, a mean of nbg=125 background photons are counted over 50 ms. This signal arises from stray 852 nm molasses light and 780 nm trapping light being collected by the detection optics. During runs when one atom is loaded in the BoB trap, the number of photons detected over 50 ms increases to a mean of ni=215 counts due to the added fluorescence from the atom trapped in the BoB. Photon counting statistics show that a single atom is loaded into the BoB trap in 45.2% of all loading trials. The near complete absence of multi-atom loading events is due to light assisted collisions and the collisional blockade effect.
The photon counting data is fitted to a compound Poisson fit
For random loading following a Poisson distribution
the probability of loading a single atom (assuming that on average we load one atom) has an upper bound of P(1)=e−1=0.368. Thus, following random loading, P(1)≤0.368. However, from
Each run begins with using the 3D MOT system to create a reservoir of 106 atoms at T≈75 K within 1 second. During the last 480 ms of the 3D MOT loading time, the 1064 nm optical conveyor belt lasers with no mutual detuning (Δω=0) are un-shuttered and overlap the MOT. Approximately 5000 to 6000 of the 3D MOT atoms are loaded into the stationary 1D lattice. Polarization gradient cooling is applied throughout the final 10 ms of the conveyor belt loading time, cooling the loaded atoms to T≈25 1K and increasing the atom density to p≈1010 atoms/cm3.
The probability of loading a single atom into the BoB trap is directly proportional to the atom density trapped in the optical conveyor belt. With the BoB approximated as an ellipsoid, an upper bound for the BoB trap volume is
V=4/3πrrad2rz (93)(4.25)
where rrad is the radius along the radial direction and rz is the axial radius.
After loading the optimal atom density into the optical conveyor belt, atoms are transported 1.1 mm to the BoB trap site by following the frequency chirping sequence previously outlined where one lattice beam is linearly frequency chirped by 400 kHz over a 5 ms time interval and accelerated to a maximum velocity of 212 mm/s. The lattice beam is then linearly chirped down −400 kHz over another 5 ms period to instantaneously stop the conveyor belt before repeating the chirping sequence in reverse to bring the atoms back to the starting point. This series of lattice beam frequency chirps first transports atoms a distance of 0.530 mm along the axial direction followed by an additional 0.530 mm in the reverse direction for a total transport distance of 1.1 mm in 20 ms. Time-of-flight temperature measurements performed on the transported atoms show that atom temperatures rise to 70 1 K after transport in the conveyor belt. Ballistic expansion in the radial direction of the conveyor belt was chosen for time-of-flight temperature measurements since it expands faster and therefore shows greater temperature dependence than in the axial direction. To load a single atom in the BoB trap from the conveyor belt and also to successfully image the loaded atom, the transported atom temperature should optimally be a factor of 50 colder than the BoB trap depth. To cool the transported atoms to such an optimal temperature, an additional 10 ms of PGC is applied to the atoms immediately after transport while still held in the 1D optical lattice. The intense light shift imparted by the lattice (equation 14 (2.11)) detunes the atoms an additional −8δ, giving a total detuning of −14.6Γ. When applying PGC after transport the atoms have an average temperature of 15 1K as compared to a temperature of 70′K when no extra PGC is applied (
During the final 5 ms of the additional PGC, the BoB trap is snapped on and overlaps the focus of the optical conveyor belt. This permits an atom to load from the transported 1D lattice into the BoB trap. After loading, the optical conveyor belt is turned off while still holding the BoB trap to allow all untrapped atoms to fall away. At this point, the BoB is loaded with either zero or one atom. A complete synopsis of the timing scheme for single atom loading via an optical conveyor belt is presented in
Fluorescence imaging is performed on the BoB trap for 100 ms to observe if zero or a single atom is successfully delivered and loaded into the trap.
This technique provides a method for correcting atom loss errors for single atom BoB traps used in neutral atom quantum computing systems. To provide an optimal vacuum environment for trapped qubits and qubit gate operations, the 3D MOT that supplies the source of cold atoms for atom reloading should be located in a separate vacuum chamber positioned adjacent to the chamber that houses the BoB traps. This also permits the 3D MOT to be on continuously, as the MOT footprint does not disturb any previously trapped atoms. Using an uninterrupted 3D MOT, the presented technique for single atom delivery and loading takes a total of 680 ms, as the procedure involves loading the conveyor belt (530 ms), transporting and loading the atom into a BoB trap (50 ms), and then finally imaging to verify the presence or absence of a single atom (100 ms). The total time for delivering and loading an atom can be reduced to a little as 55 ms total, as far-off resonance dipole traps, and hence the optical conveyor belt, can be loaded in as little as 10 ms, depending on the trap detuning, intensity, and the vacuum background pressure, while single atom signals can be observed using fluorescence imaging with exposure times as little as 5-10 ms. Additionally, the BoB hold time can be decreased to 10 ms. Even shorter hold times are possible but at the expense of introducing fluorescence noise from untrapped atoms that may be collected during imaging.
With a 13.1% probability of successful single atom delivery and loading, 33 loading trials will need to run in order to have single atom loading with 99% confidence. Working in the lower limit of a 55 ms loading time, the procedure calls for 1.82 seconds to reload an atom if all 33 attempts are executed, which is within the necessary interval of time of 2 seconds to keep a 49 element qubit array continuously occupied. Additionally, with a 13.1% probability of successful atom reloading, 50 loading trials are required to reload an atom with 99.9% confidence, which consumes 2.75 seconds if all attempts are executed, which is slightly greater than the necessary time to reload a site.
Finally, the single atom loading efficiency of this technique may be improved by incorporating methods that suppress the inherit stochastic characteristics present with using red-detuned, light-assisted collisions to prepare a single atom in the BoB trap. One such method uses a feed-back mechanism to rapidly prepare a 3D MOT with a single atom, while another method uses additional blue-detuned lasers, rather than MOT lasers, to induce repulsive light-assisted collisions to prepare single atoms in dipole traps. These methods have the capabilities to prepare single atoms directly from a MOT within 20-200 ms with average efficiencies ranging from 80% to 90%. If implemented after atoms are delivered to the BoB trap via transport, they may increase the single atom loading efficiency.
The optical conveyor belt is initially loaded with a few thousand atoms. Using the optical conveyor belt, atoms are transported a distance of 1.1 mm to a crossed vortex BoB trap. After transport, a single atom is then loaded into the BoB trap via light assisted collisions. Photon counting data confirms that an atom can be delivered and loaded into the bottle beam trap 13.1% of the time.
Implementing this work to reload a functional qubit array calls for deterministically loading any arbitrary site of an array of BoB traps. An additional pair of acousto-optic modulators for each of the 1D lattice beams can reposition the optical conveyor belt to permit loading of any BoB array site at an arbitrary location. Modifications can be made to the atom transport and loading apparatus to “steer” the optical conveyor belt with the eventual goal of loading any site of an arbitrarily large array of bottle beam traps. Steering the optical conveyor belt involves giving the atom transport system the capabilities of delivering an atom to one BoB trap and the repositioning the conveyor belt and delivering an atom to another array BoB trap.
The optical conveyor belt must be “steerable” to transport and load atoms into specific sites in an array of multiple bottle beam (BoB) traps. The conveyor belt can be transported, not only along its axial direction, but also along the {circumflex over (x)} and ŷ directions by utilizing an additional pair of acousto-optic modulators for each lattice beam. This allows one to reposition the optical conveyor belt to permit loading of a BoB trap (e.g., within an array of BoB traps) at an arbitrary (x; y; z) location.
The process of making an array of crossed vortex BoB traps involves taking a pair of LG01 beams and using a diffractive optic element (DOE) to replicate or “copy” the beams into an array of LG01 beams on a conjugate plane of the image plane of the final 0.4 NA lens used to focus the BoBs into the hex cell. While this setup concept is relatively simple to construct, it creates extensive problems with the resulting array of crossed vortex BoBs. The calcite beam displacer and DOE result in the tightly crossing BoBs at the image plane having the same polarization, rather than orthogonal polarizations. As a result, considerable phase dependent interference effects commonly arise, making it very difficult to implement as a qubit array capable of performing sensitive quantum computations. Accordingly, the illustrated embodiment employs a “Gaussian-beam array” (GBA) rather than a crossed vortex array. The GBA is used to for an array of BoB traps.
The optical system 3700 for creating the Gaussian beam array (GBA) with two BoB trapping sites is shown formally in
A 66:35 telescope 3728 then images the 2×2 array pattern through a 353 μm thick calcite crystal 3722 to replicate the 2×2 array into an arrangement of eight identical Gaussian beams separated from their nearest neighbor by d=353 μm. Equal spacing of the eight array beams is achieved by orienting the fast axis of the calcite 45 with respect to the diffractive beam splitter axis. After passing through the calcite 3722, the array of eight Gaussian beams is demagnified by a 200:60 telescope 3724 and then further demagnified and imaged through the walls of the hex cell 4000 (
After passing through all of the optics required to produce the array, the total optical power remaining in the array is 750 mW, which gives a uniform distribution of 47 mW per beam in the array setup. An image 6400 showing the optical intensity of the GBA as it is swept through the focus of the 200:60 telescope is shown in
I(r)=I0e**(−2r2/ω202). (94)(5.1)
where ω0 is the focused beam waist and I0 is the maximum intensity that occurs at r=0. As shown in
where d is the separation between neighboring Gaussian beams. The lowest optical intensity in the central, dark region ID, occurs at a distance r=d/√{square root over (2)} from the middle of each focused beam, and is equal to the sum of intensities from all four Gaussian beams at this central point r=d/√{square root over (2)}:
Comparing equations 95 (5.2) and 96 (5.3), there exists a non-zero, optical intensity gradient between the bright saddle points Isp and the dark central region 1D. Since the laser light used to create these intensity profiles is blue-detuned from resonance, equations 2.11 and 2.12 state that a potential energy well is created that permits 2D atom confinement at the dark regions; that is, the four beams together form a single 2D BoB trap. Furthermore, from equation 2.11, the total trapping depth is proportional to the difference of bright optical intensity (equation 95 (5.2)) and dark optical intensity (equation 96 (5.3)):
where s≡d/ω0 is the aspect ratio of the array. Equation 97 (5.4) can be expressed in terms of the average optical intensity of the entire 2D trap (shown in the inset of
where P is the total power of all four Gaussian beams combined.
Thus, we can use equation 98 (5.5) to express I0 in terms of the average intensity with
Inserting equation 5.6 into equation 5.4 gives the total 2D trapping intensity as
The total optical intensity in equation 100 (5.7) results in 2D atom confinement in the radial direction but provides no confinement in the loose, axial direction of the array. For the array BoBs to trap atoms in three dimensions, a repulsive optical potential in the axial direction must also be formed. Similar to how tightly overlapping two LG01 beams “caps” off the crossed vortex BoB trap with high light intensity, the overlap of the Gaussian beams as they rapidly diverge when propagating out of the GBA focus (
When finally focusing the array into the hex cell to its desired size (using with the f=23.125 mm lens), the GBA has focused beam waists of ω0=2.21 μm and beam separations of d=4.95 μm (
which is very close to optimal aspect ratio of s=2.20. Note that the focal plane of the array is carefully aligned to overlap with the focus of the optical conveyor belt (
To get an explicit calculation for the trap depth as a function of the atom position in each BoB, the intensities from each diagonally opposite pairs of Gaussian beams are added:
where cis the velocity of light, ∈0 is permittivity of free space, and E(x; y;z) is the electric field of a Gaussian laser beam given as
Equation 5.9 can be plugged into equation 2.11 to numerically calculate a dipole trapping depth of 547 μK along the radial (x,y) directions and a trap depth of 578 μK along the axial z direction when the total optical power in the entire array is 750 mW. The radial 6810 and axial 6820 trapping potentials for each BoB trap in the array are plotted in
The array can be overlapped with the center of a 3D MOT for 500 ms. After loading the array from a background MOT, polarization gradient cooling can be performed on the atoms for 5 ms and then all MOT beams can be extinguished for an additional 25 ms to allow un-trapped atoms to fall away. Fluorescence imaging can be performed using the setup 6900 of
Results of the fluorescence imaging are shown in
A setup 7100 for steering the optical conveyor belt is shown in
One branch beam from PBS 7114 is input to a pair of AOMs, that is an “X1” AOM 7116 and a “Y1” AOM 7118. Another branch beam from PBS 7114 is input to another pair of AOMs, namely, “X2” AOM 7120 and “Y2” AOM 7122. The four AOMs 7116, 7118, 7120, and 7120 are used to steer the optical convey belt so that atoms can be delivered to specific traps in the bottle-beam array. The modulated branches from PBS 7114 are focused into respective opposing faces of hex-cell vacuum chamber 4000 by respective multi-element lenses 7124 and 7126. The branch beams are thus counter-propagating and form interference fringes in the form of a 1D optical lattice 7130 within chamber 4000. The fringes can be moved by changing S in AOM 7108; atoms attracted to bright fringes thus can be “conveyed” by the moving fringes. Return beams can be detected by a photodiode 7140, the output of which (via 7142) can be analyzed using radio-frequency (rf) Spectrum Analyzer 7142.
Setup 7100 results from modifying the dynamic 1D lattice setup shown in
An optical relay system 7200 shown in
The 1064 nm beams are split into two beam paths. Each beam passes through an x-AOM 7206, 7230, which alters the angle of the beam in the x dimension. A 1:1 telescope relay 7210, 7234 images the center of the x-AOM 7206, 7230 on the center of the y-AOM 7208, 7232. This y-AOM 7208, 7208 changes the beam angle in the y dimension. A second 1:1 telescope 7214, 7236 relays the center of the y-AOM onto the center of a 5× telescope 7216, 7236, which magnifies the beam to 1 cm in diameter. The beam is then focused with a 0.4 NA custom lens 7228, 7240 onto the qubit array 7250 in hex cell 4000. At the qubit array plane 7250, the x and y positions of the beams can be laterally translated by δx and 6y by varying the x and y angular deflections of the AOMs.
For each lattice beam, the relay system, shown in
As a result of this relay system, by driving the beam deflecting AOMs, each 1D lattice beam can be displaced by δx and δy at a resolution 1.06 μm/MHz. Since the drivers can be controlled to a resolution of δν=10.1 MHz, the beam steering setup can laterally displace the optical conveyor belt by some distance at a resolution of 0:1 μm. With this setup, for both lattice beams to remain overlapped after repositioning, the x-AOM drivers are set to the same driving frequencies but opposite sign (that is x-AOM1=+δνx and x-AOM2=−δνx) while the y-AOM drivers are set to the same frequency (y-AOM1=y-AOM2=δνy) in the steerable conveyor belt setup. This can cause a non-zero mutual detuning δ≠0 between the transport beams, thus giving it an unwanted velocity. To correct this error, one must simply adjust the detuning δ in the original frequency chirping AOMs. While each AOM is driven with an independent RF driver, they are all phase locked to the same 10 MHz clock reference signal in order to minimize phase jitter of the lattice beams.
To demonstrate the single BoB site addressing capabilities with the steerable optical conveyor belt, the focused BoB array was imaged onto an Andor iXon EMCCD using the same imaging system used in
On this image, light from the focused optical conveyor belt can also be seen. By only changing the driving frequencies of the lattice deflecting x-AOMs and y-AOMs, the focus of the optical conveyor belt is repositioned such that it overlaps the center of either of the two BoB traps in the GBA. Images 7310, 7320, 7330, and 7340 demonstrate that the beam steering system is capable of addressing both BoB trapping sites in the GBA. These images were taken by only changing the 1064 nm beam deflecting x-AOM and y-AOM driving frequencies with no manual realignment of any beams. In the image 7310 at the upper left of
The sequence of delivering atoms to either of the array BoB traps begins by first choosing a BoB trap to deliver atoms to. With the spatial coordinates of both BoB traps known, the driving frequencies of the beam deflecting AOMs can be set to the proper value such that the focus of the optical conveyor belt is repositioned to overlap the center of the target BoB trap. Specifically, in this setup shown in
Once a BoB destination is chosen, a 3D MOT is loaded into the hex cell for 1 second. During the final 530 ms of MOT loading, the 1064 nm transport beams are un-shuttered allowing the optical conveyor belt to overlap the MOT where five atoms on average are loaded into the optical conveyor belt. This relatively low average atom payload may be explained by the appearance of the collisional blockade effect due to reducing the lattice beam waists combined with any light assisted collisions from the cooling light, which prevents more than a single atom occupying in each potential well of the 1D lattice potential. After loading the conveyor belt, the MOT was extinguished by shuttering the 852 nm cooling and repump beams as well as turning off the 3D MOT magnetic fields. The conveyor belt can be held in place for 20 ms to allow the remaining untrapped atoms to fall away due to gravity.
A setup 7400 for fluorescence imaging of a 1D optical lattice 7401 is shown in
An optical molasses 7420 is applied to atoms trapped in the 1D lattice 7401 and a fraction of the induced fluorescence is collected with the 0.4 NA custom lens 7406 and focused onto an Andor iXon EMCCD 7412. The occupancy of a conveyor belt can be measured by illuminating it with a 3D optical molasses by turning on 3D MOT lasers and detuning them −30 MHz. Repump light can also be turned on. Additionally, shim coils are used to eliminate any residual magnetic fields. As a result, the 3D molasses light induces atom fluorescence from any atoms loaded into the optical conveyor belt. The 852 nm scattered fluorescence light is collimated with one of the 0.4 NA, f=23.125 mm lenses 7406 used to focus the conveyor belt and separated from the 1064 nm trapping light 7412 with a 900 nm short pass dichroic mirror 7408. Using an f=500 mm achromatic lens 7409, the fluorescence light is then focused onto the Andor iXon EMCCD 7412 with single photon sensitivity where photon counting is performed over a 10 ms time interval. The overall imaging system 7400 features a magnification of 21.6 and results in a 0.740 μm/pixel ratio for EMCCD 7412. A set of three 852 nm laser line filters 7410 are used to eliminate any stray 1064 nm lattice light.
The statistics show that the most probable loading number is five atoms. The appearance of this compound probability distribution is due to the combination of: 1) the random, discrete number of atoms loading into the optical conveyor belt; and 2) the Poissonian spread of observed photon counts from each of the loaded atoms. The fitting result is μ=0:5, N=5, n=10, nbg=120, and n1=120. The appearance of multiple peaks in the photon counting histogram (as compared to 2 peaks in
Unlike the crossed vortex BoB trap loading, the BoB traps formed by the GBA are never turned off during the atom transport sequence. Turning off the GBA in order to reload a new atom to an empty BoB trap would ruin any previously loaded traps and ongoing quantum computations would also be ruined, Because the GBA and this optical conveyor belt are focused through the same 0.4 NA, f=23.125 mm lens, the conveyor belt will pierce atoms through the weaker, axial confining potential of the array BoBs. During this process, atoms confined to the conveyor belt are not expelled, as the axial trapping frequency of the conveyor belt, Ωz=21*450 kHz is considerably greater than the axial potential of the array BoB trap Ωz=2T*2.58 kHz.
To transport the small number of loaded atoms to the desired array BoB trap, one lattice beam is linearly frequency chirped by S=20 kHz over 5 ms. Because the frequency chirping AOMs 7102 and 7104 in
Demands for maintaining a satisfactory vacuum environment for trapped qubits and qubit gate operations call for locating the 3D MOT which supplies the source of cold atoms for atom reloading in a separate vacuum chamber positioned adjacent to the cell that houses the BoB array. Being in a separate chamber permits the 3D MOT to be on continuously, as the MOT footprint does not disturb any previously trapped atoms. One practical atom reloading sequence transports atoms vertically from the 3D MOT, using what is called a Moving Molasses MOT, to a position in a hex cell that is off-axis from the BoB array.
In the setup 7700 of
By making a 6-Beam MOT with three pairs of counter-propagating cooling beams in a 2D magnetic field, the entire MOT can be propelled vertically, in what a Moving Molasses MOT, by mutually detuning the vertically oriented MOT beams by δω. Similar to the physics behind the optical conveyor belt, this detuning creates a reference frame where both pairs of cooling beams are Doppler shifted to the same frequency. This frame moves vertically (as an atom fountain) with velocity ν=√{square root over (2)}2λδω. A 12 W, 1064 nm dipole trapping beam helps to guide the atoms upwards. Once the atoms reach a maximum height in the Hex cell, they are trapped in a crossed dipole trap which serves as secondary cold atom reservoir 7704. The atoms from this secondary reservoir are then horizontally transported to a BoB site using an optical conveyor belt. Using this transport scheme, while considerably more difficult, maintains a better vacuum in the Hex cell, and solves isolation issues by keeping the cold atom reservoir from disrupting any concurrently operating qubit gates taking place in the array, which allows for a continuous atom reloading.
A optical system 7800 for steering an optical conveyor belt to address six BoB traps is shown in
Optical system 7800 provides fully capable atom reloading system for a qubit array without disturbing neighboring, previously loaded BoB traps. A 2×2 array of identical Gaussian beams is created by Holo/Or MS-248-X-Y-A diffractive beam splitter 7806. The 100:35 telescope 7810 images the arrays in between two calcites 7812 and 7814. The first calcite 7812 replicates the 2×2 array into an array of eight identical Gaussian beams separated by 353 μm, while the second calcite 7808 duplicates the 8 beams into a full 16-beam array with each beam separated by 250 μm.
In order to have equal spacing of each array site, the fast axis of each of the calcites is rotated by 45° with respect to each other. After passing through both calcites, the 2×3 array of Gaussian beams is demagnified by a 200:60 telescope to give each beam a waist of 34 μm and beam separation of 74 μm. The array is then demagnified one final time with a 500:23.125 telescope to image the array through the hex cell. At the focus of the 0.4 NA, f=23.125 mm lens, each beam has a separation of d=4.67 μm and each beam waist is ω=2.22 μm, giving an aspect ratio of d=w=2.10.
An Ultra-cold Gain Mechanism for a Matterwave Transistor Oscillator introduces a new gain mechanism for a matterwave transistor oscillator. Atomtronics is an emerging field of ultra-cold atomic physics that focuses on quantum circuits that operate as atom analogues to electrical circuits, where an atom current takes the place of an electron current and the chemical potential takes the place of a voltage. In their design, atomtronic circuits are drastically different from their electronic circuit counterparts, as they involve manipulating ensembles of ultra-cold atoms to flow through a variety of potential energy landscapes, such as: 1) optical potentials produced from laser beams; 2) harmonic magnetic potentials created with atom chips; and 3) a hybrid combination of optical and magnetic potentials, such as the transistor potential used throughout.
At the heart of all atomtronic circuits is a Bose-Einstein condensate (BEC), which supplies a source of coherent, ultra-cold atoms that are controlled to flow throughout a network of quantum potentials. Some recently demonstrated atomtronic circuits of note include BECs flowing throughout double-well and triple-well potentials, as well as BECs owing through 1D optical lattices and optical ring lattices. Irrespective of the circuit potential geometry, the objective of atomtronics remains the same: to manipulate a current of ultra-cold atoms in some tailored potential to study fundamental physics or to create purely quantum mechanical technologies.
Accompanied with such high aspirations of atomtronics comes a difficult, inherent challenge ingrained in the field: such ultra-cold atom circuits are inherently a many-bodied, open quantum system, since atoms and BECs interact and dissipate as they flow throughout the circuit from some source to a sink or drain. While this immensely complicates the theoretical description of atomtronic circuits, complete descriptions are usually simplified by invoking various approximations, such as: 1) Neglecting dissipation and treating the atomtronic circuit as a closed system; 2) reducing the many-bodied system from a very large number to a three or four-body problem; 3) ignoring the BEC mean-field interaction effects.
The transistor oscillator model described herein uses all three of the above bulleted approximations and justifies the validity of each when applied to this transistor oscillator model. Despite the necessity of using approximations just to solve for something, atomtronics has made unwavering progress in ultra-cold atom technologies such as realizing quantum simulators, improving precision measurements, creating an atomtronic analogue to a SQUID (called an AQUID, or atomtronic quantum interference device), and realizing an atom. Moreover, if used in conjunction with well-established matter-wave interferometry, atomtronics can help enhance inertial sensing and gravimetry sensitivities by up to 10 orders of magnitude as compared to their light wave interferometry counterparts. Finally, atomtronics can conceive new types of quantum technologies by combining various atomtronic circuit elements into an integrated ultra-cold atom circuit. One such integrated ultra-cold circuities a matter-wave transistor oscillator.
Fundamental to many technological developments, the transistor is a three-terminal semiconductor electronic device that uses an applied voltage or current on one terminal to amplify or switch electrical signals on another terminal. An atomtronic transistor behaves analogously to its semiconductor counterpart by switching or amplifying atom, or matter-wave, currents rather than electrical currents. Moreover, in place of three electrical terminals (commonly known as the source, gate and drain terminals) found on semiconductor transistors, an atomtronic transistor features three quantum potential energy wells that are capable of confining ultra-cold atoms, as shown in
Most importantly, for such source-gate-drain system to be called an atomtronic transistor, it must display a gain mechanism for a matter-wave or atom current. The overall concept of an atom current gain mechanism in an atomtronic transistor is shown in
In a semi-classical formalism, transistor-like behavior of ultracold atoms flowing through networks of quantum potential energy wells involve creating matter-wave analogues of Kirchoff's current and voltage laws for atoms flowing across a double Gaussian barrier potential. Compared to the quantum model, the semi-classical picture is defined as the set of conditions where the thermal energy of the atoms flowing through the transistor potential is much larger than the energy level spacing of the transistor gate-well. In the semi-classical model atom currents are driven by collisions and their flow is described using classical Boltzmann statistics. Atoms that have enough thermal energy to traverse the Gaussian barriers flow into the drain and establish an atom current Is. A non-intuitive result of the semi-classical model is that given an initially empty gate-well, an atom population quickly grows in the gate and rapidly condenses into a Bose-Einstein condensate. The BEC steady state population influences the atom current ISD flowing from source to drain, which is characteristic of a gain mechanism.
As shown in
There are major differences between the semi-classical and the quantum transistor models. Namely, in the quantum model: 1) the thermal energy of the source-well atom ensemble is comparable (ideally equal) to the gate-well energy level spacing; 2) atoms flowing through the transistor potential are described as matter waves; 3) the gate-well condensate is treated as a coherent state; the ow of atoms through the transistor is a result of resonant tunneling from the source well, through the gate-well, and into the drain-well (as opposed to the collisions primarily driving the semi-classical atom current); and 4) when the atom current reaches low enough energies, a new interaction, not observed in the semi-classical model or experiment, arises in the gate-well between the gate-well BEC and the atom current.
This new interaction, which arises due to the motional state of an oscillating gate-well BEC coupling high lying gate-well energy eigenstates, is similar to the Rabi model interaction between two-level atoms and a photon field. The interaction serves to increase the atom current by ultimately amplifying the matter-wave probability amplitude in the transistor drain-well. This introduces an additional gain mechanism that only appears when the atom current is at sufficiently low enough temperatures for the BEC coupling (described above) to occur.)
The matter-wave transistor oscillator described herein is constructed by overlapping a harmonic, magnetic potential with two repulsive Gaussian barriers. This setup forms three potential energy wells which are labeled the source well, gate well, and drain well (
Given a thermal atom ensemble originating in the source-well, a flux Φ of atoms per unit time flows through the transistor potential and into the drain-well by resonantly tunneling through the transistor gate-well. Our model operates with the assumption that only the phenomenon of tunneling facilitates the flow of atoms throughout the transistor potential and does not allow atoms to classically traverse the SG or GD barriers. Thus, the atoms are hereafter described as matter waves. This definition leads to a first principle of our matter-wave transistor model: “Only the source matter-waves satisfying the resonant tunneling boundary conditions of the gate-well successfully flow with flux Φ into the drain-well (
This leads to the second principle of our matter-wave transistor model: “by placing a dipole oscillating BEC in the gate-well, the matter-wave current flowing into the drain-well with flux Φ′ is amplified due to the BEC broadening the resonant tunneling boundary conditions (
The mechanism responsible for this increase of matter-wave current is the emergence of an atom-BEC interaction that occurs only when the motional state of a gate-well BEC couples high lying gate-well energy eigenstates. An interesting consequence is the atom-BEC interaction allows matter-waves that do not initially match a resonance condition of the gate-well to absorb energy from or deposit energy into the interaction, allowing off-resonant matter-waves to shift their energy and satisfy the resonant boundary conditions of the gate-well. The now-resonant matter-waves are then permitted to resonate and build up inside the gate-well followed by subsequent flow out into the drain. Combining the first and second principles yields: “when a dipole oscillating BEC occupies the transistor gate well, the matter-wave current flowing into the drain with flux Φ′ is greater than the flux D when the BEC is absent. This characteristic shows the presence of a gain mechanism, where gain g is defined as
=Φ′/Φ>1 (104)(7.1)
VS(x)=VD(x)=0 (105)(7.2)
where matter-waves are described for free particles and the source and drain wave functions, ΨS(x; t) and ΨD(x,t) are given as plane wave solutions
ΨS(x,t)=ASe±i(k
ΨD(x,t)=ADe±i(k
where wavenumber ki=√{square root over (2 mEi)}/ℏ for kinetic energy Ei(where i=S, D for source and drain) and mass m. The centrally located gate-well is separated from the source-well and drain-well by the SG and GD repulsive Gaussian barriers featuring equal 1/e waists of ω0=2 m and asymmetric heights 1 of VSG/ℏ=30 kHz and VGD/ℏ=33 kHz. The asymmetry in the barrier heights is chosen such that the ratio ν≡(VSG−VGD)/Ts, called the “feedback parameter” of the transistor oscillator, is equal to ν=1. When barrier heights VSG, VGD and source-well ensemble temperature Ts give a feedback parameter ν=1, there is an optimal balance of particle and energy currents flowing into and out of the gate-well to yield a stable atom ensemble in the gate well with high enough phase-space density to condense in steady state. This remarkable behavior is equivalent to a resistor with negative resistance since the source atoms at temperature Ts that enter the gate and condense into a BEC are cooled to a lower temperature. If one were to flip the barrier heights so that VSG=33 kHz and VGD=30 kHz, the atom current would actually heat up upon entering the gate and no BEC would form. This latter case is equivalent to a resistor with positive resistance. Thus, the harrier height asymmetry in the model allows for a steady-state BEC to collect in the gate-well as source atoms flow through the gate and into the drain-well.
By separating the SG and GD barriers by a distance d=4:85 μm, the gate-well becomes a slightly anharmonic, asymmetric oscillator potential capable of confining ultra-cold atoms. The anharmonicity establishes a degeneracy amongst transitions between the gate-well ground and first excited states and the highest lying pair of bound energy eigenstates.
Finally, despite the slight anharmonicity of the gate-well, the gate-well wavefuctions n(x; t) can be approximated as solutions to the quantum harmonic oscillator with
where n corresponds to the nth energy eigenstate of the gate well and Hn(x) is the nth Hermite polynomial. With the gate-well eigenfunctions assumed to be harmonic oscillator solutions, as well as noting the gate-well potential used in this model features 24 bound energy eigenstates, the degeneracy described above exists between gate-well energy eigenstates ψ0(x,t) and ψ1(x,t) with ψ22(x,t) and ψ23(x,t).
With the degeneracy between lowest and highest lying pairs of eigenstates present, we can write that
E1−E0=E23−E22=ℏωP (108)(7.5)
where each En is the energy eigenvalue for the ψn gate well eigenstate and ωp is the fundamental oscillation frequency of the gate well potential defined as E0=ℏωP/2. As a consequence of this degeneracy (equation 108 (7.5)), given a gate-well BEC occupying ψ0 with an excited portion occupying ψ1 and undergoing dipole oscillation at frequency ωP, the BEC couples the highest lying pair of energy eigenstates, ψ22 and ψ23. Therefore, as atoms occupying ψ22 and ψ23 collide with the oscillating BEC, stimulated emission and absorption of phonons with energy ℏωp results in continuous Rabi flopping between the upper two energy eigenstates. It is from this continuous atom-BEC interaction (as a result of the BEC coupling the highest lying pair of states) that leads to a gain mechanism in the transistor oscillator model. The presence and effects of such gain mechanism emerges directly from the gate-well many-bodied Hamiltonian.
Since the dipole oscillating BEC couples only the highest lying pair of energy eigenstates and is decoupled from any other transition, the quantized energy structure of the gate-well can be simplified by ignoring any intermediate energy eigenstates located in-between the lowest and highest lying pairs of states. That is, the model only considers states 0, 1, 22 and 23. Using this approximation, the gate well is reduced from 24 states to a four-level system. Accordingly, the pertinent energy eigenstates can be relabeled to |0, |1, |2 and |3, where the first two states (|0 and |1) are the ground and first excited states and the latter two states (|2 and 3) pertain to the highest lying pair of bound energy eigenstates in the gate well as shown in
E1−E0=E3−E2=ℏωP (109)(7.6)
In this reduced representation, as atoms transition between upper states |2 and |3 by colliding with the BEC, the excited condensate occupying |1 consequently gains (or loses) a phonon of energy ℏωP by by absorbing (or depositing) a quanta of energy ℏωP from (or into) |0. This is completely analogous to the Rahi model for two-level atoms interacting with a quantized electromagnetic field. Where the Rabi model features transitions between excited and ground electronic states |e and |g by absorbing (and emitting) photons from an oscillating electric field, the matter-wave transistor features similar transitions between upper harmonic oscillator states |2 and |3 by absorbing (and emitting) phonons from an oscillating BEC.
The intermediate eigenstates (3-21) of the gate-well may be ignored since omitting such states does not change the final results of the quantum model. Regardless of how many energy eigenstates make up the transistor oscillator gate-well, the model only requires a dipole oscillating BEC in the ground and first excited states to couple the highest lying pair of bound eigenstates and no other states. Whether the gate well has 24 states, 100 states, or 4 states, as long as this constraint is fulfilled, the gate-well many-bodied Hamiltonian, as well as the results of the model, will remain unchanged.
To explicitly show the emergence of a transistor oscillator gain mechanism in this model, the many-bodied Hamiltonian of the reduced gate-well confining a BEC undergoing dipole oscillation with frequency ωp can be formulated. In
where ℏωi is the single particle energy of gate well eigenstate |1 and where ai† and ai are the eigenstate's corresponding creation and annihilation operators. The factor of ½ in equation 7.7 prevents double counting of atom-BEC interactions. Since the quantum model uses harmonic oscillator eigenfunctions, the operators are defined in the number state basis as
ai|Ni=√{square root over (Ni)}|Ni−1
ai†|Ni=√{square root over (Ni+1)}|Ni+1
[ai,ai†]=1
[ai,aj†]=0. (111)(7.8)
Equation 7.8 shows the ground state annihilation and creation operators are a0|N0=√{square root over (N0)}|N0−1) and a0†|N0=√{square root over (N0+1)}|N0+1. However, if the BEC contains a large number of atoms, then N0≈N0−1 and reasonably approximate a0|N0≈√{square root over (N0)}|N0. Therefore, the result of this Bogoliubov approximation is to replace the ground state operators a0 and a0† with the real √{square root over (N0)}.
The factor gin equation 110 (7.7) is the interaction energy between two low energy, long wavelength atoms with equal masses m, given as
where a is the s-wave scattering length of the colliding atoms. Finally, the inner product (ij|V|kl) is the overlap integral over four harmonic oscillator eigenfunctions (of equation 107 (7.4)).
ij|V|kl=Vijkl=∫ψi(x)ψj*(x)ψk(x)ψl*(x)dx (113)(7.10)
giving the relative amplitude of a transition ai†aj†ak al between the normalized reduced gate well eigenstates |I, |j, |k and |I. Using the definitions in equations 111 (7.8), 112 (7.9), and 113 (7.10), the first summation of equation 111(7.7) physically represents all energy contributions from individual gate atoms occupying eigenstates |I, |j, |k and |I, while the second summation is energy contributions from atom-BEC collisions that result in transitions ai†aj†ak al between gate-well eigenstates.
Since the gate-well eigenstates are orthogonal, expanding the single-particle contributions to equation 111 (7.7) (H0) yields
where δij is the Kronecker delta. The model treats the BEC occupying state |0 of the gate-well as a macroscopic occupation of the gate-well ground state and accordingly uses the Bogoliubov approximation to state:
a0≈√{square root over (N0)} (115)(7.12)
where N0 is the number of atoms in the ground state BEC. Using this approximation, equation 114 (7.11) can be expressed as
H0=ℏω0N0+ℏω1a1†a1+ℏω2a2†a2+ℏω3a3†a3. (116)(7.13)
From equation 110 (7.7), the summation for energy contributions due to atom-BEC collisions, given
is considerably larger than the summation for single particle energies (equation 7.11) and contains 256 terms when evaluated exactly. Thankfully, some helpful approximations can eliminate all but four of the terms without altering the results of the model. Justification for excluding certain collective sets of terms in the summation in equation 7.14 is provided as follows.
First, mean-field energy shifts are neglected. That is, terms in the summation containing factors of ai†ai†aiai or ai†ai†akai are excluded from the model. Since mean-field shifts increase all gate-well energy levels equally, the shifts do not change the overall eigenstate layout of the gate-well and, more importantly, preserve the coupling between the highest pair of bound eigenstates via the oscillating condensate. Whether the mean-field energy shifts are included or not, the physics of the model remains unchanged. Accordingly, the model ignores the mean-field terms in equation 117 (7.14). Second, nonsensical transitions with the form ai†ai†aiaj are ignored by the model. These represent meaningless transitions that are either unphysical such as a3†a3†a3a0 or transitions that leave the system unchanged, such as a3†a3a3a3†a1†a1a1 a01. Third, the rotating wave approximation is used to eliminate the terms a0†a1†a2a3. These contributions in equation 117 (7.14) oscillate rapidly compared to the other terms and thus average approximately to zero. Fourth, the Bogoliubov approximation (equation 115 (7.12)) is once again used since the number of atoms occupying the ground state BEC is considered to be large.
Using the first two bullet points in the above described set of approximations, the Hamiltonian for atom-BEC collisions, which is now labeled the interaction Hamiltonian Hint, is reduced to
Using the commutation relations for creation and annihilation operators [ai†aj†]=[al, aj†]=δij, equation 7.15 is further reduced to only six unique terms and is expressed as
Finally, using the third and fourth of approximations, the rotating wave and Bogoliubov approximations reduce the interaction Hamiltonian to its most condensed form in the harmonic oscillator basis:
The Hamiltonians H0 (equation 116 (7.13)) and Hit (equation 120 (7.17)) for single-particle and atom-BEC contributions to the gate well many-bodied Hamiltonian are expressed in the harmonic oscillator (or number) basis. Considering that the dipole oscillating BEC couples upper states |2 and |3, the reduced gate-well (shown in
where |+ and |− are the symmetric and antisymmetric normal modes, respectively. (The normal mode basis, with eigenstates |+ and |−, is simply a 2D vector space rotated θ=π/4 from the original |2 and |3 basis states.) Additionally, in the normal mode basis, the oscillating, excited portion of the gate-well condensate that occupies |1 is described as a coherent state |α with magnitude |α|.
When using matrix notation, the symmetric and antisymmetric modes given in equations 7.18a and 7.18b are defined as
To convert H0 (equation 116 (7.13)) and Hint (equation 120 (7.17)) into the normal mode basis, the following symmetric and antisymmetric operator transformations are introduced:
Applying the operator transformations 123a, 123b, and 123c (7.20a, 7.20b, and 7.20c) to equation 116 (7.13) outputs
and applying the operator transforms to equation 120 (7.17) produces
Combining equations 124 (7.21) and 125 (7.22) gives the many-bodied Hamiltonian for the reduced gate-well in the normal mode basis:
Equation 7.23 can be simplified further by using the following identities:
a+†a+=|++| (127a)(7.24a)
a−†a−=|−−| (127b)(7.24b)
(a+†±ia+†)(a+∓ia−)=1±σ2 (127c)(7.24c)
(a+a†)(a+†a+−a−†a−)=(a=a†)σ3 (127d)(7.24d)
where 1 is the identity matrix, σ2 is the Pauli y-matrix, and σ3 is the Pauli z-matrix defined as
Plugging equations 127a-127d (7.24a-7.24d) into equation 126 (7.23) and simplifying gives the many-bodied Hamiltonian for the reduced gate-well in the normal mode basis as
Since the normal modes |+i and |− are orthogonal and complete, matter-waves can be treated as occupying the reduced, coupled gate-well as a linear combination of normal modes:
where A± are the probability amplitudes for the |± modes, there exists matter-wave interference between the constituent (and simultaneously occupying) |+ and |− modes. As a result, the total wave function in the gate-well has a beat note that oscillates with a relative phase ϕξ. This phase described here is a relative phase defined with respect to some reference in the gate well, such as the ground state BEC occupying |0. An example of the superposition of gate well normal modes, as well as the relative phase, is shown in
The amplitude of the atom-BEC interaction Hamiltonian, gV√{square root over (N0)}, gains a phase dependence and can be expressed as a complex valued function ξ as
ξ=gV√{square root over (N0)}eiϕξ (131)(7.28)
Using equation 131 (7.28) in equation 129 (7.26), the entire gate-well many-bodied Hamiltonian can be expressed in the normal mode basis as a sum of a free space term H0, a condensate term Hc, and the atom-BEC interaction term Hint as
where
Other than a change of basis and the fixed condensate energy, equation 132 (7.29) is identical to the Rabi Hamiltonian for a two-level atom interacting with an electromagnetic field. Similar to the Rabi Hamiltonian, in the gate-well free space Hamiltonian component H0 derived in equation 132 (7.29), the σ2 term represents spontaneous emission: Hse=(ω3−ω2) σ2 ℏ/2 and couples states of opposite symmetric or anti-symmetric parity. Where the Rabi Hamiltonian has electronic states interacting with photons, the gate-well Hamiltonian HG involves center-of-mass states (meaning the symmetric and antisymmetric normal modes) interacting with phonons.
As shown in equation 7.29, the atom-BEC interaction Hamiltonian expressed in the normal mode basis is
Hint(ξ*a+ξa†)σ3 (133)(7.30)
Eigenstates of the interaction Hamiltonian are |α, ±=|α⊗|±, where |αis a coherent state with magnitude |α|, and |± are the symmetric and antisymmetric normal modes of the gate-well. The corresponding eigenvalues Eint of the interaction Hamiltonian are
where φ is the relative phase difference between ϕξ, the phase of the matter-wave interference of the gate-well normal modes, and ϕα, the phase of the gate-well dipole oscillating BEC. Therefore, the physical interpretation of φ=ϕξ−ϕα is that it represents the phase of a difference of phases. The oscillating condensate and matter-wave interference are shown relative to each other in
=2|ξα|≤Eint≤+2|ξα| (135)(7.32)
The interaction energy eigenvalue in equation 7.31 can be written in a more complete form by taking the expressions for g (equation 48 (2.41)) and V (equation 113 (7.10)) and plugging them into the expression for ξ (equation 131 (7.28)), giving
Additionally, since the dipole oscillating BEC occupying |1 is described as a coherent state |α the definition of the coherent state magnitude can be used to write
|α|=√{square root over (N0)} (137)(7.34)
Taking the results of equation 136 (7.33) and 137 (7.34) and inserting them into equation 134 (7.31) yields a complete expression for the eigenvalue of the interaction Hamiltonian Hint in the normal mode basis:
The gain mechanism of the matter-wave transistor oscillator can be understood by comparing the flow of a matter-wave current resonantly tunneling from the source-well, through the gate-well, and into the drain-well in two different cases: 1) when no dipole oscillating BEC occupies the transistor gate-well; that is when the gate-well is “empty” and Hint=0; and 2) when a dipole oscillating BEC occupies the gate-well and couples high lying gate-well eigenstates; that is when Hint≠0.
A greater flux of source-well matter-waves flows via resonant tunneling through the gate-well and into the drain-well when a dipole oscillating BEC occupies the gate (when Hint≠0) as opposed to when the gate has no such BEC (when Hint=0). As a result, a greater matter-wave probability amplitude exists in the transistor drain-well when the interaction potential Hint is present compared to when it is missing from the system, thus the interaction potential Hint qualifies as a matterwave gain mechanism.
Consider the triple-well transistor potential with an “empty” gate, that is, no oscillating BEC present in the gate-well. As a source-well matterwaves flows through the gate-well, it's probability amplitude in the drain-well is a function of the transmission coefficients of the gate-well. Such transmission coefficients are calculated for our specific gate-well potential by utilizing a numerical technique called the general impedance method, which cleverly formulates an equivalence between matterwaves tunneling through the gate-well (whose solutions are governed by the 1D time-independent Schrödinger equation) and electromagnetic waves propagating through transmission line junctions (whose solutions are governed by transmission line theory from classical electrodynamics). The general impedance method is capable of calculating transmission coefficients for arbitrary potential barriers in quantum systems, making it useful for determining transmission properties of the asymmetric, double Gaussian gate-well for the atomtronic transistor.
Examination of the transmission spectrum,
Using the reduced gate-well nomenclature, when incident source-well matterwaves match a gate-well resonance boundary condition by having energy E2 or E3, the matterwaves couple into the gate-well and resonate in the gate-well, followed by subsequent tunneling into the drain-well with probability τ2≈10−4 or τ3≈10−4, respectively. Slight deviations from these transmission resonances (by having an energy E just off-resonance from E2 or E3) results in a substantial decrease in matterwave probability amplitudes in the drain well. Consequently, for the case of no BEC in the gate well: the resulting observable effect is that an appreciable atom flux Φ flows into the drain-well only when source-well atoms are initially prepared to match a gate-well eigenstate.
With the presence of the BEC, the atom-BEC interaction Hamiltonian Hint (equation 133 (7.30)) must be included in the gate-well many bodied Hamiltonian. By including Hint, the matterwave current (an atom flux) flowing from source to drain increases. To understand why the current is increased, consider a single matterwave ψs(x,t) with frequency ωs originating in the source-well and whose energy does not coincide with a gate-well eigenstate; that is ωs≠ω2≠ω3. The model treats the source-well matterwaves as plane wave solutions (equation 106 (7.3)). The model operates in the normal mode basis and additionally, treats the source-well matterwave as a linear combination of the gate-well symmetric s+(x,t) and antisymmetric ψs−(x,t) normal modes:
ψs(x,t)=A+ψs+(x,t)+A−ψs−(x,t) (139)(7.36)
with amplitudes A+ and A (and s denotes a source-well quantity). Using the definitions for the symmetric and anti-symmetric modes (equation 7.18a and 7.18b) produces expressions for ψs+(x,t) and ωs−(x,t) as
with wavenumber ks and frequency ωs. Since the symmetric (equation 140a (7.37a)) and antisymmetric (equation 140b (7.37b)) components of ψs(x,t) are sinusoids that are offset by π/2 radians, they are characterized as the in-phase and quadrature components of the source-well matterwave, respectively.
Using the above source-well matterwave definitions, when ψs(x,t) penetrates the SG barrier and enters the gate-well, the atom-BEC interaction, described by Hine, occurs. For the incoming wave ψs(x,t) to simultaneously match the gate-well boundary conditions and conserve energy, the wave must deposit energy into, or absorb energy from the atom-BEC interaction energy. The energy shift ΔE imparted by the atom-BEC interaction is calculated in equation 134 (7.31) as ΔE=±2|ξα|cos φ where the sign of the energy shift is positive for the symmetric component ψs+(x,t) and negative for the antisymmetric component ψs−(x,t). It follows that for an initially off-resonant matterwave ψs(x,t) with frequency ωs to match a gate-well eigenstate with frequency ωg and satisfy gate-well resonance boundary conditions, the atom-BEC interaction must change the incoming wave energy by Δ=ℏ ωs−ωg. Using equation 134 (7.31), this energy shift defines an in-phase condition
whereby, if satisfied, the in-phase, symmetric component ψs+(x,t) (equation 140a (7.37a)) is shifted onto a gate-well eigenstate and couples into the gate-well, while the quadrature, antisymmetric component ψs−(x,t) (equation 7.37b) is completely reflected since ψs−(x,t) is shifted even further from the gate-well resonance. Note, one could, without changing the model results, define instead a quadrature condition of
whereby, if satisfied, the antisymmetric component ψs−(x,t) couples into the gate with the in-phase, symmetric component being reflected.
The domain of equation 141 (7.38) reveals that the range of source-well matterwave frequencies ψs that may be shifted onto a gate-well resonance with frequency ωg via the atom-BEC interaction is given by
which indicates the interaction Hamiltonian broadens the gate-well resonant tunneling bandwidths A from the very narrow peaks shown at the left in in
where subscripts 2 and 3 correspond to the resonant frequencies ω2 or ω3 of gate-well eigenstates |2 and |3, respectively.
The result of equations 144a and 144b are shown pictorially in
From equations 144a (7.40a) and 144b (7.40b), it is apparent that the resonant tunneling bandwidth is controlled by the magnitude |ξα|, which is a tunable parameter in the sense that the larger the gate-well condensate, the greater the magnitude of |ξα|, and hence a larger interaction potential Hint, which results in broader resonant tunneling bandwidths Λ2 and Λ3.
By comparing the matterwave flux with and without a gate-well BEC, it is apparent that there is a matterwave gain mechanism for the transistor oscillator. Compared to when no BEC is present in the gate-well, the atom-BEC interaction Hint increases the range of matterwaves that may couple into the gate well with subsequent resonant tunneling out to the drain well. The gain mechanism thus amplifies the matterwave probability density in the drain-well, which results in the observable effect of intensifying the atom flux tunneling into the drain-well when a gate well BEC is present compared to when it is absent.
Since each of the symmetric and antisymmetric components ψS+(x,t) and ψs−(x,t) constitute fractions |A+|2 and |A−|2 of the total wavefunction ψs(x,t), thus, given a gate-well dipole oscillating BEC described by |α with magnitude |α|2 that couples gate-well eigenstates |2 and |3, if a source-well matterwave with frequency ψs is within either frequency bands Λ2 or Λ3, then the symmetric component ψs+(x,t) is transmitted through the gate-well 100×|A+|2 percent of the time with probability τ2 or τ3, respectively. The remaining fraction is reflected with certainty.
The physics in the transistor oscillator model draws very similar comparisons to the cold atom micromaser, where a slow, excited-state atom approaching a microwave cavity is randomly reflected from or transmitted through with 50% probability. This behavior occurs due to a coherent interaction between a quantized field of microwave photons and two-level atoms, which creates a potential barrier that reflects one component of the atom wave function with certainty and transmits the other component with probability close to unity. In the micromaser case, since each of the components constitutes half of the wave function, about half of the atoms are transmitted through the maser cavity with the other half reflected. The matterwave transistor oscillator is analogous, except that it additionally demands resonant tunneling for transmission.
Finally, as an example to characterize the gain mechanism, the general impedance method is used to numerically calculate the gate-well transmission spectrum when a dipole oscillating BEC (described as a coherent state |0) with N0=20,000 atoms occupies the ground state |0 and N1=|α|2=5 atoms occupies the first excited state |1 of the gate-well. The transmission spectrum for this specific case is shown in
With the DOBEC coupling upper states |2 and |3, the atom-BEC interaction potential described by Hint (equation 133 (7.30)) broadens the resonant transmission peaks by Λ2=Λ3=150 Hz (equations 144A and 144B (7.40a and 7.40b). To approximate the gain of the transistor oscillator in this example, the ratio of matterwave transmission curves when the BEC is present to when the BEC is absent is plotted in
A finite temperature 7Rb Bose-Einstein condensate of approximately 7*103 atoms at a temperature of 50 nanoKelvin (nK) serves as the atom ensemble used to source the transistor oscillator system. The BEC is formed using an atom chip in the Atomtronics Double MOT UHV of
A summary 9000 of all laser cooling, repump, optical pumping, and probe transitions are shown in
Cooling setup 9100 is shown in
The 780 nm cooling master ECDL is locked to the 85Rb |F=3→|F′=2/|F′=4 crossover transition by picking off a small amount of light from the output of the ECDL using Doppler free saturated absorption spectroscopy. An additional 5 mW of cooling light is picked off from the main ECDL output to be used as a frequency reference which is stabilized using a phase locked loop. Using a frequency offset locking setup, the light is blue-shifted 1.2 GHz such that the laser is red-detuned by 2Γ from the 87Rb|F=2 →|F=3 cooling transition. Additionally, 40 mW of cooling light is used to seed a Toptica BoosTA tapered amplifier, which then outputs 1000 mW of 780 nm laser light. Using a polarizing beam splitter 9132, the TA output is split into two paths, where 500 mW of 2D MOT cooling light and 70 mW of 3D MOT cooling light are coupled with 50% efficiency into their individual polarization maintaining fibers and sent to the apparatus.
A repump setup 9200 is shown in
For repumping light, a 780 nm repump master ECDL 9202 is locked to the 87Rb |F=1 →|F=1/|F′=2 crossover transition by picking off a small amount of light from the output of the ECDL and using Doppler free saturated absorption spectroscopy. Approximately 1 mW of repump light is picked off from the main beam line and used to seed a repump slave laser 9220. The output of both the master and slave lasers are each sent through a respective AOM 9210, 9224 driven at 87 MHz to blue-shift the beams onto resonance with the 87Rb|F=1 →|F′=2 transition. Each beam is then coupled into a polarization maintaining fiber 9212, 9226, where the master light is sent to the 2D MOT and the slave light is sent to the 3D MOT. With the described configuration, the master repump setup provides 1 mW of 3D repump light and the slave provides 12 mW of 2D repump light.
A pump and probe laser setup 9300 is shown in
The system used for producing a BEC is the optical pumping and probe setup 9300 shown in
An in-situ laser setup 9400 is shown in
The production of a BEC in the Atomtronics Double MOT UHV cell begins by generating a 3D MOT of 87Rb inside the 3D chamber of the UHV cell (
This configuration provides 3D cooling along the orthogonal {circumflex over (x)}, ŷ, and {circumflex over (z)} directions. Anti-Helmholtz coils oriented around the 3D chamber as in
After loading the 3D MOT, the 2D MOT lasers are shuttered and the 3D MOT is compressed into a CMOT by ramping the anti-Helmholtz field gradient up an additional 10 G/cm while simultaneously detuning the cooling lasers by an additional 40 MHz and attenuating the 3D repump to 50 μW of optical power. Typically, the CMOT stage lasts 15 ms and the MOT atom temperature after compression is approximately T≈130 to 160 μK. After CMOT, sub-Doppler cooling is performed. The MOT anti-Helmoltz coils are switched off while x-bias, y-bias, and z-bias coils are turned on in order to eliminate any stray magnetic fields at the location of the atoms such that a PGC stage can properly be performed. After eliminating stray magnetic fields at the atoms, w the 3D cooling lasers are red-detuned an additional 105 MHz for 3 ms, which cools the atoms to temperatures of approximately T≈20 μK. After PGC, all cooling light is extinguished while optical pumping light illuminates the atoms for 800 μs to pump as many atoms as possible to the |F=2, mF=2 Zeeman magnetic sublevel. As shown in equations 24 and 25 (2.18 and 2.19), since mF>0, this permits the optically pumped atoms to be trapped in a minimum of a spatially varying magnetic field. At the conclusion of the optical pumping stage, all lasers are extinguished and no more laser cooling stages are applied.
Given that the atoms are prepared in the |F=2, mF=2 Zeeman sublevel, if current through the 3D MOT coils is linearly ramped in the Helmholtz configuration over a span of 35 ms, the optically pumped atoms are then captured in a magnetic quadrupole field with a 60 G/cm field gradient. To vertically transport the atoms from the location of the 3D MOT up to the atom chip, current is linearly run through a pair of magnetic transport coils, also in the Helmholtz configuration to produce a transport field gradient of 110 G/cm while the quadrupole field is linearly ramped to zero over the same time period. This operation vertically shifts the zero point of the quadrupole field that the atoms are trapped in. As long as the movement is done adiabatically, the atoms remain trapped in the magnetic quadrupole field as it is transported vertically. Upon reaching the height of the atom chip, the transported atoms are transferred into a single-guide wire Ioffe-Pritchard trap produced on the atom chip.
To transfer atoms into this chip trap, the transport field is first compressed over a time span of 65 ms such that the magnetic quadrupole field gradient matches the (anticipated) chip trap gradients. After such compression, the chip guide-wire and H-wire currents, as well as the y-bias current are linearly ramped on over a span of 10 ms and a magnetic potential, as shown in section 3.3.1, is formed. Approximately 20*106 atoms at a temperature of 17 to 20 μK remain in the Ioffe-Pritchard trap after this transfer sequence and are located approximately 150 μm below the bottom surface of the atom chip.
To condense the atoms into a BEC, a forced RF evaporation stage is applied to the trapped atom ensemble to further cool the atoms below the BEC critical temperature (equation 2.29) necessary for condensation to occur. In order to perform efficient RF evaporation, the chip magnetic trap is compressed along the radial and longitudinal directions in 500 ms in order to increase the atom collision and re-thermalization rate during the evaporative cooling stage. Using an RF antenna located 3 mm above the top surface of the atom chip, RF radiation with a frequency of νRF=40 MHz is illuminated onto the atoms. The RF frequency is then linearly swept from νRF=40 MHz to 26 MHz in 150 ms, which expels the absolute warmest atoms from the chip trap. Afterwards, the RF frequency is swept over a span of 1800 ms from
νRF=26 MHz→11 MHz→6.5 MHz→4.8 MHz→2.62 MHz (145)(8.2)
where the final RF value is approximately 100-200 kHz above the trap bottom, leaving on the coldest of the cold atoms remaining in the magnetic trap.
After the above described RF evaporation sequence, a Bose-Einstein condensate of about 7*103 atoms is produced in the chip guide-wire trap. This BEC is produced under the chip guide wire and thus blocks optical access to the condensate. To move the BEC out from under the wire and align it with the center of the atom chip window (as shown in
ω(x,y,z)=2π×(67,1500,1500)Hz. (146)(8.3)
Absorption images 9711, 9712, and 9713, shown in
The atom density profile of the condensate shown in
The triple-well transistor oscillator potential is formed by overlapping a harmonic magnetic potential with two repulsive Gaussian barriers. The harmonic magnetic potential is the same cigar shaped Ioffe-Pritchard trap with radial and longitudinal trapping frequencies of ωrad=2π*67 Hz and ωrad=2π*1600 Hz, respectively that confines the 50 nK BEC. The magnetic potential is located approximately 150 μm below the bottom surface of the atom chip window. Next, the repulsive source-gate (SG) and gate-drain (GD) barriers are formed by focusing two blue-detuned, 760 nm Gaussian laser beams with equal beam waists of ω0=2:0 μm and a separation of d=4.8 m onto the center of this harmonic magnetic potential.
The 760 nm laser setup for producing the SG and GD barriers is shown in
A schematic for projecting the SG and GD barriers onto the magnetic potential is shown in
The same lens 9812 hat focuses the barriers also serves to for in-situ absorption imaging of atoms, e.g., such as image 9900 in the upper portion of
A total of 35 mW of 760 nm laser light enters the barrier projection system from a polarization maintaining, single mode fiber and is collimated to a beam waist of ω=4 mm with an f=15 mm aspherical lens. The 760 nm light then propagates into an IntraAction Corp DTD-6010RH29 2D acousto-optic detector (AOD) 9804, which can dynamically control the shape and positions of the SG and GD barriers. (The 0 order beams are blocked leaving only the +1, +1 diffracted order beams.)
To generate the longitudinal (or x) position of the barriers, two RF tones, one for each barrier, drive the AOD and diffracts each beam into its +1 order along the {circumflex over (x)} direction. The RF tone for the GD barrier is generated by a VCO while the tone for the SG barrier is produced by a Lecroy ArbStudio Arbitrary Waveform Generator (AWG). Because the RF driving voltages are easily adjustable, this scheme provides a mechanism to dynamically change the longitudinal positions of the SG and GD barriers. The radial (or y) position of the both SG and GD barriers is generated by using a second channel of the AWG to continuously raster the SG and GD barriers back-and-forth in the y direction at a frequency of 20 kHz.
Using an Infinity Probe TS-160 objective in conjunction with a Zeiss LD Plan Neouar 0.6 NA, 40× magnification microscope objective, the optical barrier potential is projected onto the harmonic magnetic potential (
This barrier projection system also doubles as an imaging system to observe the flow of ultra-cold atoms through the triple-well transistor potential. To observe atoms in the transistor potential, in-situ absorption imaging is performed (as opposed to time-of-flight imaging) while the atoms are still confined to the triple well potential. A probe beam resonant with the 87Rb D2|F=2 →|F′=3μ transition is propagated vertically through the vacuum cell (entering through the base of the 2D MOT chamber and exiting though the atom chip window, see
A Bose-Einstein condensate of approximately 7*103 atoms with a temperature of T=50 nK (
In
Ultra-cold atoms are loaded into the transistor source-well by turning on the SG barrier and adiabatically sweeping it from the right-most longitudinal boundary of the harmonic magnetic potential to near the center of the magnetic potential in 40 ms (
The drain terminator, which is a 10−13 W, 780 nm beam focused to a waist of 16 μm, is projected onto the center of the transistor drain-well using the same optical projection system for the SG and GD barriers and can remove all drain-well atoms within 5 ms without disturbing any atoms confined to the source or gate well. After clearing any residual atoms not swept into the source-well, the SG barrier height is ramped down to its desired height of VSG=30 kHz while simultaneously, the GD barrier is ramped on to a height of VGD=31 kHz and is displaced from the SG barrier by a distance d=4:8 μm. After this sequence, the matterwave transistor is prepared in a state shown in
With a current energy of approximately 1 kHz, the source-well atom ensemble does not flow through the transistor-oscillator potential. As shown in
This is accomplished by rapidly sweeping the SG barrier back-and-forth at a frequency of 19 kHz, extending over a majority of the source well, while ramping up the optical power of laser beam forming the SG barrier. Specifically, an arbitrary waveform generator (AWG) is used to ramp the driving voltage of the acousto-optic deflector (AOD) that controls the longitudinal position of the SG barrier.
The total sweeping frequency of the SG barrier is 19 kHz, which is considerably larger than both the 1.6 kHz and 125 Hz radial and longitudinal trapping frequencies of the transistor source well potential. As a result, by continuously repeating this waveform, the atoms experience the time averaged intensity of the sweeping optical field, which is shown in
Time-of-flight imaging performed on the atom ensemble in the raised source-well shows that the atom temperature increases to approximately 78 nK and the ensemble is no longer condensed into a BEC. The heating is attributed to technical noise from position fluctuations of the magnetic potential due to current noise of 35 nA/√{square root over (Hz)} in the current drivers that supply the atom chip currents Ix and 120 nA/√{square root over (Hz)} in the current drivers that form the z-bias field βz account for a 1 nK/ms heating rate. No measureable heating has been observed from sweeping the SG barrier back and forth.
With an ultracold atom ensemble prepared in the raised source-well, the flow of an atom current through the transistor oscillator is observed by probing the system in-situ. The raised source-well can be loaded with approximately N=7*103 atoms at a temperature of approximately Ts=78 nK and a chemical potential of μs=27 kHz. The atom current begins to flow at time t=0, which is defined as the time when the source-well chemical potential is raised to μs=27 kHz and the SG and GD barriers are set to their desired heights of VSG=30 kHz and VGD=31 kHz.
The transistor oscillator system is allowed to evolve for sometime Δt, during which an atom current flows throughout the transistor potential landscape. As the system evolves, the drain-well is constantly illuminated with the drain terminator beam to remove any atoms that flow into the drain-well. The terminator ultimately prevents the transistor oscillator system from reaching an equilibrium. Note however, that the terminator beam is extinguished approximately 6 ms prior to taking an absorption image in order to take an instantaneous snapshot of the atom current accumulating in the drain-well. The 6 ms corresponds to approximately ¼νx where νx is the longitudinal trapping frequency of the harmonic magnetic potential.
After an ultracold atom current flows for sometime Δt, the occupancy of the source, gate, and drain wells of the transistor oscillator is observed by performing in-situ absorption imaging.
Time-of-flight measurements show that during the first 10 ms of the atom current flow time, the source-well ensemble temperature increases from 78 nK to approximately 109 nK. However, once the current begins to flow into the drain-well, the source-well ensemble continuously cools. This shows that as the atom current flows into the gate and drain wells, it carries away energy from the source-well. This source-well heating and cooling behavior can be expected since the in-situ absorption images show that an atom current isn't flowing during the first 10 ms, indicating that there is no flow of energy out of the source well. Once current begins to flow into the gate and drain wells, there is a noticeable drop in source-well temperature.
The ultracold atom current flowing through the transistor oscillator in
The longitudinal position of the atoms collecting in the drain-well is indicative of their energy. The greater total energy the drain-well atoms have, the greater distance they flow into the drain-well until reaching the classical turning point. A corollary of this is that the longitudinal position of the drain-well atoms is also an indication of their momentum. It follows that the longitudinal width of the drain-well ensemble is thus a measurement of the distribution of the momenta of the drain-well current. For example, if every single atom entering the drain-well had the same momentum, every single atom would travel to the same longitudinal turning point in the drain-well. However, if the drain-well current contains atoms with a spread of momenta Δp, then there will be a corresponding spread of the drain-well ensemble at the location of the longitudinal turning point.
By comparing the drain-well ensembles in which the current either does or does not form a BEC, it is clear that when a significant atom current flows into the drain-well, the BEC drain atom flux of
The absorption images in
Φ(p)∝P;P=constant (147)(8.4)
where Φ(p) is centered about the longitudinal turning point of the drain-well and where the width of Φ(p) is proportional to the spread of momenta of the atom flux flowing into the drain-well, resulting in the step function A710 of some constant height and width shown in
Consider an ensemble of atoms confined to a harmonic oscillator potential. As the ensemble evolves in the harmonic oscillator potential, the atoms have a minimum momentum and a maximum displacement from the center of the potential when they reach the longitudinal turning points. Consequently, when located at the center of the harmonic oscillator potential, the atom displacement is at a minimum and the atom momentum is at a maximum. That is, position and momentum are conjugate variables and are Fourier transform duals. The drain-well momentum wavefunction Φ(p) and the gate-well position wavefunction Ψ(x) can be treated as canonically conjugate variables that are Fourier transform duals of each other. Since the total momentum space wavefunction Φ(p) of the atom current in the drain-well is a step function (or for lack of better term, a “box function”), the envelope of the total position space wavefunction Ψ(x) in the gate-well can be obtained by performing the Fourier transform of Φ(p). That is
Ψ(x)∝sinc(x) (148)(8.5)
as shown in
When an atom current begins to flow out of the source-well and atoms become trapped in the gate-well, the wavefunctions of all of the atoms confined in the gate-well interfere in a superposition and the resulting envelope of the total position space wavefunction Ψ(x) is a sinc function (equation 148 (8.)5). Source-well matterwaves that couple into the gate and flow into the drain also interfere in the drain-well. The matterwave probability amplitudes suggest that destructive interference occurs in the drain-well between the GD barrier and the drain-well longitudinal turning point. On the contrary, constructive interference occurs at the turning point. The resulting superposition creates an envelope of nearly constant probability amplitude (equation 147 (8.4)) at the drain-well turning point whose width is directly proportional to the range of unique matterwave momenta able to couple into the transistor gate-well.
Recalling the gain mechanism theorized above, matterwave currents that feature a greater coupling strength to the gate-well, as a result of the gain mechanism, will have a greater range of energies (and by extension, a greater range of momentum) permitted to flow into the drain-well. This gain mechanism theory can be descried using the above definitions for Ψ(x) and Φ(p). Using a wavefunction description, since Ψ(x) and Φ(p) are Fourier duals, matterwave currents that feature a greater coupling to the gate-well, as a result of the gain mechanism, have a broader momentum distribution in the drain-well and corresponding tighter spatial distribution in the gate-well. This exact behavior was observed in
Additionally, depending on whether the prepared source-well atom ensemble is a purely thermal ensemble or a condensate, different spreads of momenta occur in the drain-well current. This indicates that the gate-well coupling strengths for a matterwave current depends on the temperature and coherence of the ultracold atoms flowing through the transistor. When preparing the source-well with a purely thermal atom ensemble with initial temperature of T=78 nK, the imaging results in
The observed increase in the spread of the drain-well momentum wavefunction Φ(p), which is also an indication of a greater matterwave coupling to the drain well, when using colder atoms as bottom image A804 of
The results of the matterwave current in
This fact that a drain-well atom flux only appears at the longitudinal turning point of the drain-well is attributed to destructive interference occurring within the matterwave superposition at longitudinal positions between the GD barrier and the drain-well turning point. Moreover, at the drain-well turning point, the matterwave superposition constructively interferes, which results in the observed drain-well atom flux. This result suggests that the gain mechanism is a coherent gain mechanism. This coherence can be credited to the gain mechanism only transmitting a symmetric (or antisymmetric) matterwave while always reflecting the mode of opposite parity. The transmitted modes, all of which have the same symmetric (or antisymmetric) parity, add up coherently in the drain-well, resulting in destructive interference at longitudinal positions after the GD barrier and constructive interference at the drain-well turning point.
If the atom current was incoherent, then matterwave destructive interference in the drain-well cannot occur and atoms will appear in all locations throughout the entire drain-well as shown in
Finally, the empirical observation of the canonically conjugate relationship between the envelopes of the gate-well position space wavefunction T(x) and the drain-well momentum space wavefunction Φ(p) is also suggestive that the matterwaves flowing into the drain-well have coherence. Without coherence, this canonically conjugate relationship would not be observed in the transistor oscillator results.
The quantum model predicts dynamics different from the semi-classical results only when the atoms flowing through the transistor oscillator have sufficiently low enough temperatures such that the motional state of a dipole oscillating BEC, placed in the transistor itself, couples atom transitions between high lying transistor energy eigenstates. In this ultracold temperature regime, interesting physics that were not observed in the semi-classical theory or experiments occurs, specifically the manifestation of a coherent gain mechanism that broadens the range of energy (and by extension, momenta) of a matterwave current flowing from the source-well, through the gate-well, and into the drain-well.
Time-of-flight temperature measurements of the transistor source-well ensemble confirm that as an atomic current flows through the transistor, the current constantly carries energy away from the source-well ensemble, which provides continuous cooling to the source-well atoms. Interestingly, our experimental results showed that a wider momentum spread in the atom flux flowing into the drain-well is detected when the source-well ensemble contains a condensate, compared to when the source-well ensemble is a purely thermal collection of atoms. Thus, the atom current sourced with a condensate has a greater coupling to the transistor gate-well. Additionally, this result confirms a canonically conjugate relationship between the spatial wavefunction of the gate-well and the momentum wavefunction of the drain-well, which suggests coherence in the transistor gain mechanism and the drain-well atom current. Evidence of a coherent drain-well atom current can also be observed in the complete absence of an atom flux at all positions in the drain-well with the exception of an atom flux located at the longitudinal turning point. This observation is due to a superposition of coherent matterwaves where destructive interference occurs throughout the entire drain-well except at the longitudinal turning point, where the matterwave superposition constructively interferes.
In order to measure the differential matterwave gain due to the atom-BEC interaction (equation 133 7.30) change in atom flux owing into the drain-well as a function of the change of occupancy of the gate-well BEC must be measured. In an experimental regime where a matterwave current only flows via tunneling through the SG and GD barriers, this measurement provides a more robust method of quantifying matterwave gain, as it removes any impact of stray source-well atoms that may randomly enter the gate and drain-wells (as noise) by classically traversing the SG and GD barriers. Changes in drain-well flux as a function of the change in gate-well occupancy, rather than strictly measuring the drain-well flux as a function of gate-well atom number could be measured. This method would essentially “ignore” noise in the atom flux.
A challenge of measuring differential gain is changing the occupancy of the gate-well BEC without altering the gate-well potential. From the semi-classical models and experiments, the size of the gate-well BEC is controlled by the ratio of SG and GD barrier heights, VSG and VGD. While we can readily control the BEC occupancy by adjusting VSG and VGD, the quantum model shows that doing so changes the degeneracy (and subsequent coupling) between the dipole oscillating BEC and the highest lying pair of gate-well eigenstates. This degeneracy and coupling, being an essential precondition to matterwave gain, cannot be altered. Instead, we can adjust the flux of source-well atoms that enter the gate-well by increasing or decreasing the source-well chemical potential, μs. Extending the semiclassical formalism, an atom current flowing from the source-well to the gate-well ISG is proportional to the fugacity factor
ISG∝exp(μs/kTs) (149)(9.1)
where k is the Boltzmann constant, Ts is the source-well temperature, and μs is the source-well chemical potential. Since the source-well chemical potential can be freely adjusted, the parameter μs becomes a tunable characteristic of the transistor without impacting the SG and GD barrier heights or the gate-well eigenstates. From this, there is an experimentally adjustable “knob” that one can use to control the amount of atoms that enter the gate-well, allowing us to measure the transistor differential gain.
Explicitly measuring the transmission of the gate-well symmetric or antisymmetric modes ψ+(x) and ψ−(x) can be accomplished by performing a Hermite-Gaussian transform of the drain-well atom flux in-situ absorption images. A Hermite-Gaussian transform is a mathematical operation that decomposes the drain-well probability amplitude signal (which note is just the drain-well absorption image), such as in
{Hm(x)}=∫−∞∞e−x
This operation is completely analogous to the Fourier transform which decomposes a function of time into the frequencies that make it up.
To give an example comparing the Fourier and Hermite-Gaussian transforms, consider some sine wave function, such as ƒ(x)=sin(2×)+sin(5×). The Fourier transform of ƒ(x) is expressed in the frequency basis and would feature delta functions F(ω)=δ(ω±2)+δ (ω±5) corresponding to the frequencies ω=2 and ω=5 present in the original function ƒ(x). Now, given some harmonic oscillator wavefunction, such as ψ(ω)=2(x)+5(x), the Hermite-Gaussian transform of Ψ(ω) decomposes the function to the harmonic oscillator basis and would feature delta functions (n)=δ(n−2)+δ(n−5) corresponding to the Hermite-Gaussian modes n=2 and n=5 present in Ψ(ω).
In the transistor-oscillator, the observable signal (i.e. the drain-well in-situ absorption image) is comprised of the square of Hermite-Gaussian functions (since the absorption image really represents the matterwave probability amplitude squared). Thus, the transform of the square of the wavefunction is taken. The corresponding Hermite-Gaussian transform for this signal can be found by taking the definition of the Hermite-Gaussian transform (equation 150 9.2) and applying the convolution theorem to arrive at
{Hm(x)}=∫−∞∞e−x
It then follows that the Hermite-Gaussian transform of a product of Hermite-Gaussian functions is
where the “⊙” operation denotes the convolution. This formalism can be applied to determine the expected Hermite-Gaussian transform of the drain-well signal. Given the symmetric mode that transmits into the drain-well is
ψ+(x)=[ψ22(x)+ψ23(x)]/√{square root over (2)} (153)(9.5)
the observable signal is
|ψ+*(x)ψ+(x)|=½[|ψ22(x)|2+|ψ23(x)|2+2|ψ22(x)ψ23(x)|]. (154)(9.6)
Applying equation 151 (9.3) to take the Hermite-Gaussian transform of each product, yields:
{|ψ+(x)|2}=½[{|ψ22ψ22(x)|}+{|ψ23ψ23(x)|}+{|ψ22ψ23(x)|}]
∝δ(n−44)+δ(n−45)+δ(n−46) (155)(9.7)
Therefore, the Hermite-Gaussian transform of the observed symmetric mode transmitted into the drain-well takes the form of three delta functions in the harmonic oscillator basis: a first at twice the mode number corresponding to the highest lying gate-well eigenstate (corresponding to |3 in the reduced model); a second at twice the mode number corresponding to the second highest lying gate-well eigenstate (corresponding to |2 in the reduced model); and a third at the sum of both gate-well eigenstate modes making up the symmetric normal mode.
A plot showing the result of a Hermite-Gaussian transform applied to a hypothetical transmitted symmetric mode is shown in
This explicitly measures what mode has transmitted into the transistor gate-well, which identifies which modes are having their probability amplified as a result the gate-well gain mechanism. As a generalized extension of equation 154 (9.6) and 155 (9.7), if the constructed gate-well potential has N bound energy eigenstates, then the Hermite-Gaussian transform of the drain-well signal in the presence of the ultracold gain mechanism is
{|ψ+(x)|2}∝δ(n−(2N−2))+δ(n−(2N−1)+δ(n−2N) (156)(9.8)
As is apparent from equation 134 (7.31), the atom-BEC interaction responsible for transmitting symmetric (or antisymmetric) modes into the drain well has a phase dependence φ. Consequently, the phase of the transmitted symmetric mode that enters the drain-well may not be the same for each successive run of the experiment. This complicates the problem of averaging drain-well images (such as averaging in-situ absorption images to filter random noise) due to the symmetric mode having a different phase between shot-to-shot. Because the Hermite-Gaussian transform of the drain-well signal produces the same result irrespective of the phase of the transmitted symmetric mode, one must perform many data runs of the experiment and average all of the Hermite-Gaussian transforms of the absorption images rather than the raw absorption image themselves in order to collect reliable data to verify the presence of a specific symmetric of antisymmetric normal mode owing into the drain-well.
The foregoing demonstrates a coherent matterwave gain mechanism as the operating principle of a matter-wave transistor that can serve as a building block for complex multi-transistor integrated atomtronic circuit, e.g., permitting atomtronic equivalents of logic gates. For example, an atomtronic controlled NOT (CNOT) gate can be implemented with two independently controlled matter-wave transistors. The CNOT gate is a “universal” gate in that any other logic gate can be constructed from sets of CNOT gates. The atomtronic CNOT gate can thus be used to build atomtronic computers.
In addition, the matter-wave transistor can be used to implement a continuous, coherent matterwave source, i.e., an atom laser. The continuous coherent matterwaves can be used to source integrated matterwave circuits. These integrated circuits leverage the physics of coherent matterwaves launched into waveguides where the matterwave current can switched, divided, and recombined as the matterwaves flow throughout the system. The matter-wave transistor can be used as a continuous, rather than pulsed, coherent matterwave source.
As explained above, a source-to-drain matterwave flux can be increased by switching from a first matter-wave transistor state to a second matter-wave transistor state, where, in the second state, the gate well of the matter-wave transistor includes a DOBEC, and, in the first state, the gate well does not include a DOBEC. For example, in the first state, the gate well: 1) may not contain any matter-wave particles; 2) may include matter-wave particles that have not formed a BEC; or 3) may include a BEC that is not without an oscillating dipole. In case where a gate well includes a DOBEC, the source-to-drain matter-wave flux can be increased by increasing the population of the DOBEC. Herein, “ultra-cold” means less than 100 nanoKelvin; in some embodiments, the matter-wave flux is characterized by a temperature below 60 nanoKelvin. Likewise, the matter-wave flux may or may not include a BEC.
Herein, “resonant tunneling conditions” refers to conditions, e.g., of a gate well and/or its barriers with a source well and a drain well, that affect the likelihood of a particle tunneling through a source-gate barrier and/or a gate-drain barrier. The “breadth” of resonant tunneling conditions refers to a range or a combination of ranges of resonant conditions affecting the likelihoods of tunneling through a barrier.
Herein, all art labeled “prior art”, if any, is admitted prior art; all art not labelled “prior art”, if any, is not admitted prior art. The illustrated embodiments, variations thereupon and modifications thereto are provided for by the present invention, the scope of which is defined by the following claims.
Number | Name | Date | Kind |
---|---|---|---|
20180107938 | Morello | Apr 2018 | A1 |
20200395448 | Conradson | Dec 2020 | A1 |
Number | Date | Country | |
---|---|---|---|
20200161446 A1 | May 2020 | US |
Number | Date | Country | |
---|---|---|---|
62769868 | Nov 2018 | US |