Designing Arbitrary One-dimensional Potentials on an Atom Chip

TL;DR

This paper demonstrates the creation of customizable one-dimensional potentials for ultracold atoms on an atom chip using laser shaping, enabling advanced studies in quantum dynamics and simulations.

Contribution

It introduces a method to generate arbitrary 1D optical potentials on an atom chip by laser shaping with a digital micro-mirror device, combining optical and magnetic trapping advantages.

Findings

Realized various 1D potentials including box-like and sinusoidally modulated traps.

Enabled tunable barriers within the potentials for quantum control.

Provided new tools for studying 1D quantum systems and quantum thermodynamics.

Abstract

We use laser light shaped by a digital micro-mirror device to realize arbitrary optical dipole potentials for one-dimensional (1D) degenerate Bose gases of 87Rb trapped on an atom chip. Superposing optical and magnetic potentials combines the high flexibility of optical dipole traps with the advantages of magnetic trapping, such as effective evaporative cooling and the application of radio-frequency dressed state potentials. As applications, we present a 160 ${\mu}$m long box-like potential with a central tuneable barrier, a box-like potential with a sinusoidally modulated bottom and a linear confining potential. These potentials provide new tools to investigate the dynamics of 1D quantum systems and will allow us to address exciting questions in quantum thermodynamics and quantum simulations.

Equations

$${ϵ}_{\mathrm{RMS}}^i=\sqrt{\frac{1}{j-l}\sum _{k=j}^l{\left(\frac{{\rho }_{\mathrm{T}}^i(z_k)-{\bar{\rho }}^i(z_k)}{{\rho }_{\mathrm{T}}^i(z_k)}\right)}^2},$$

$$V(z)=\mu -\frac{g}{2\pi a_{\perp }^2}\frac{\rho (z)}{\sqrt{1+2a_{\mathrm{s}}\rho (z)}}-\frac{\hslash {\omega }_{\perp }}{2}\left(\sqrt{1+2a_{\mathrm{s}}\rho (z)}+\frac{1}{\sqrt{1+2a_{\mathrm{s}}\rho (z)}}\right).$$

Full text

††journal: oe

1 Introduction

Ensembles of ultra cold atoms are an ideal and versatile tool to investigate a plethora of physical phenomena. Their near-perfect isolation from the environment, paired with the availability of powerful techniques to probe and manipulate them, allows for the observation of coherent quantum evolution over long times with unprecedented detail [1, 2]. When building such model systems, flexible and precise control over the confining potential often represents a key ingredient to access new phenomena. The realization of homogeneous systems, for example, can unveil effects otherwise masked by an inhomogeneous density distribution [3, 4, 5], which allowed for the observation of recurrences in our isolated quantum many-body system [6]. Full local and dynamic control of the potential landscape can allow for the controlled excitation of collective modes [7], the realization of sonic black holes [8, 9] or exploring highly entangled topological states[10]. This flexibility is often achieved through optical dipole potentials [11] of shaped light fields realized by spatial light modulators (SLMs) [12, 13, 14, 15].

Another versatile technique for trapping and manipulating cold neutral atoms are micro-traps generated by atom chips [16, 17]. In a number of experiments a large variety of potential configurations have been realized using combinations of static magnetic [18, 19], electric [20], optical [21] and radio-frequency (RF) fields [22]. Atom chips are well suited to create precisely controlled, steep and robust confinements, in particular highly elongated trapping configurations. In such traps, the atoms are tightly confined along two spatial dimensions and loosely confined along the third dimension, realizing a quasi one-dimensional (1D) configuration (see Fig. 1). In our setup, the tight transverse confinement is controlled by the current flowing in a single micro-wire and a homogeneous external bias field. The weak longitudinal confinement, on the other hand, is generated by additional current carrying wires on the chip. The detailed shape of the longitudinal confinement is influenced by imperfections in the current flow [23] which can be measured by magnetic field microscopy [24, 25, 26, 27].

In this paper we demonstrate that the spatial structure of the longitudinal trapping potential and its imperfections can be ‘corrected’ by applying an additional structured dipole potential to the existing atom chip trap, which maintains the transverse trapping and manipulation. A far blue detuned light field will create a conservative repulsive potential for the atoms. Applying the structured light field from a direction orthogonal to the 1D chip trap allows nearly arbitrary modification of the potential landscape along the 1D trap.

This setup preserves the advantages of the tight magnetic trap in the transverse direction, such as the ability to deform the trap into a double-well potential using RF dressed state potentials, while greatly improving the control over the longitudinal 1D confinement. We use a digital micro-mirror device (DMD) to shape arbitrary intensity profiles necessary for realizing a large variety of 1D optical dipole potentials.

2 Basic design

The potential created by a light field interacting with induced dipole moments of neutral atoms is proportional to the intensity of the light field and the sign of the detuning, $V_{\mathrm{dip}}(\bm{r})\propto\mathrm{sgn}(\omega-\omega_{0})I(\bm{r})$  [11]. Here $\omega_{0}$ is the optical transition frequency, $\omega$ and $I(\bm{r})$ are driving frequency and intensity of the light field, respectively. The sign of the detuning $\Delta=\omega-\omega_{0}$ determines whether the potential is repulsive ($\Delta>0$, blue-detuned) or attractive ($\Delta<0$, red-detuned). Figure 2 shows an example of how the right intensity distribution can be used to create a box-like potential by correcting a rough harmonic confinement through the addition of blue- or red-detuned dipole potentials.

The spatial structure of optical dipole potentials is determined by the spatial structure of the light beam. As such, beam-shaping techniques based on different types of SLMs can be exploited to modulate not only the intensity but also the phase of light fields. We use a DMD which is a reflective SLM offering fast and flexible intensity modulation of light fields.

DMDs are powerful devices in terms of controllability and speed of performance, but ultimately, they can only produce pixelated binary patterns. To create gray-scale images from binary DMD patterns, time and/or spatial averaging has to be applied.

Imaging the 2D binary pattern of a DMD results in the 2D intensity distribution $I(y,z)$ on the image plane (plane of atoms). Consider a 1D Bose gas extending along the $z$-axis in a highly anisotropic harmonic trap ($\omega_{\perp}/\omega_{z}\gg 1$). For typical experimental values of $\omega_{\perp}$, the transverse width of such an atomic cloud is much smaller than the resolution of commonly used optical imaging systems. For example, in our experiment with $\omega_{\perp}\sim 2\pi\times 2\,$\mathrm{kHz}, the width of the Gaussian ground state wave function is $a_{\perp}=\sqrt{\hbar/(M\omega_{\perp})}\sim 250\,$\mathrm{nm}, with $M$ being the mass of a \ce^87_Rb atom. Thus, assuming that the cloud is located at $y_{\mathrm{cloud}}$, only a tiny slice from the image, $I(z)\equiv I(y=y_{\mathrm{cloud}},z)$, is relevant for shaping the longitudinal potential. In the rest of the manuscript, $I(z)$ refers to a 1D intensity profile in the longitudinal direction at $(y=y_{\mathrm{cloud}},z)$.

To achieve a fine intensity resolution in the projected light pattern we must implement spatial averaging of the binary DMD pattern. This can be accomplished by designing the imaging system in such a way, that through large demagnification, many binary pixels contribute to a single diffraction limited spot in the image plane (see Fig. 3). In our case of controlling the longitudinal potential of a 1D system, the demagnification has an additional advantage: the DMD pixels contribute to the intensity at $y_{\mathrm{cloud}}$ depending on their vertical distance from the central axis, and therefore we obtain an improved intensity resolution at $y=y_{\mathrm{cloud}}$.

An additional knob of control is spatial filtering: manipulating the Fourier components of light affects the optical resolution of the system. Blocking high frequency components worsens the optical resolution, which ultimately extends the contribution range of a single pixel in the binary pattern of the DMD.

In the case of 1D systems, only the spatial resolution along the longitudinal direction is relevant. Reducing the numerical aperture along the perpendicular direction will leave the resolution in the longitudinal direction unchanged. But more pixels from the perpendicular direction will contribute to the 1D intensity. This improves the tuneability of the potential at $y=y_{\mathrm{cloud}}$. The price to pay for performing spatial filtering is a loss in overall intensity.

Beyond spatial averaging fast intensity modulation by the DMD can be exploited to extend the dynamical range of the created potential. This modulation has to be fast compared to the dynamical timescales of the trapped atoms, such that only a time-averaged potential is seen by the atoms [28, 29].

An alternative to imaging the DMD directly is to employ it in the Fourier plane of an imaging system as a binary filter to modulate both phase and amplitude of the light field [30]. However, if phase modulation is not necessary, direct imaging turns out to be more straightforward.

3 Experimental implementation

In our experiment, we realize 1D degenerate Bose gases of \ce^87_Rb atoms with an atom chip [19]. Current carrying micro-fabricated wires on the chip surface along with an external bias field create a highly anisotropic Ioffe-Pritchard type trap with approximately 3D harmonic confinement: $V_{\mathrm{mag}}(\bm{r})=M\omega_{\perp}^{2}\,(x^{2}+y^{2})/2+M\omega_{z}^{2}\,z^{2}/2$, with $M$ being the mass of a single atom, $\omega_{\perp}$ the transverse frequency and $\omega_{z}$ the longitudinal frequency of the trap. Typical parameters are $\omega_{z}\sim 2\pi\times 10\,$\mathrm{Hz} and $\omega_{\perp}\sim 2\pi\times 2\,$\mathrm{kHz}. Furthermore, the cigar shaped cloud can be transversally split into two by applying RF fields [22].

Absorption imaging is used to obtain projected 2D atomic density distributions, $\rho(y_{k},z_{k})$, in time of flight (TOF) after switching off the trap [32]. The spatial grid $(y_{k},z_{k})$ correspond to camera pixels in the plane of atoms. Density profiles along the longitudinal axis of the 1D system, $\rho(z_{k})$, are extracted by integrating over the $y$-direction. The scheme of our absorption imaging system (AIS) is shown in Fig. 4. The projected size of the camera pixel in the object plane of the AIS (plane of atoms) is $1.05\,$\mathrm{\SIUnitSymbolMicro m}. The influence of the imaging resolution of the AIS onto the measured 1D profiles can be approximately modeled by a convolution with a Gaussian point spread function with $\sigma_{\mathrm{AIS}}=2.1\,$\mathrm{\SIUnitSymbolMicro m} [33].

The small transverse width of the cloud renders direct in situ imaging not practical. The weak longitudinal confinement suppresses the longitudinal expansion and allows us to measure the longitudinal density profile at short TOFs (between $1$ and $3\,$\mathrm{ms}$$) with negligible influence of the expansion. For cold samples, in situ thermal fluctuations are also negligible at this short TOFs [34, 35, 36].

Although a variety of trapping geometries can be produced using RF dressed state potentials, once the atom chip is designed, adjusting the longitudinal confinement is greatly limited. Moreover, there exists potential roughness due to imperfections of micro-wires on the chip [23]. Compensating for this potential roughness without an additional potential is possible in principle, but requires special chip designs [37]. Adding an optical dipole potential allows to correct the potential roughness and shape the longitudinal trapping potential in a very flexible way.

To create the dipole potentials we use a $\lambda_{\mathrm{dip}}=660\,$\mathrm{nm} laser (opus660, Laser Quantum) which provides $1.5\,$\mathrm{W} of spatially single mode light, blue-detuned with respect to the D line of \ce^87_Rb ($\lambda_{\mathrm{D_{2}}}=780\,$\mathrm{nm}, $\lambda_{\mathrm{D_{1}}}=795\,$\mathrm{nm}). For our purpose, which is primarily correcting an existing potential and creating potentials with steep walls in 1D, a blue-detuned light field is advantageous [11]. The choice of $\lambda_{\mathrm{dip}}$ is mainly motivated by the anti-reflection coating of the existing optical elements in the experiment (see Fig. 4) which favours wavelengths ranging from $650$ to $1150\,$\mathrm{nm}. Being far-detuned from the resonance guarantees that for $U_{\mathrm{dip}}/\hbar\sim 2\pi\times 1\,$\mathrm{kHz} the scattering rate is $\sim 2\pi\times 0.1\,$\mathrm{mHz}$$. This leads to only one spontaneous scattering per second for an ensemble of $10^{4}$ atoms, causing negligible atom loss and heating.

We apply the concepts discussed in Sec. 2 by imaging DMD patterns directly onto the plane of atoms. Our DMD chip is a DLP9500 from Texas Instruments with $1080\times 1920$ micro-mirrors, each $10.8\times 10.8\,$\mathrm{\SIUnitSymbolMicro m}. The refresh rate for the whole array is $17\,$\mathrm{kHz}, which increases to $50\,$\mathrm{kHz}$$ when the number of active rows is limited to $50$. The DMD chip is a part of the V-9501 module by ViALUX. This module also contains a V4395 main board which includes a completely configured high-speed FPGA logic and firmware.

The optical setup imaging DMD patterns onto the plane of atoms is shown in Figure 4. It consists of two demagnification stages: first, a 4f system is used to perform the horizontal-pass spatial filtering using an adjustable rectangular slit. Secondly, a plano-convex lens along with the objective of the AIS images the DMD pattern onto the plane of atoms.

The DMD chip is placed at the focal point of the first lens and the distance between lens 1 and lens 2 is the sum of the two focal lengths, $f_{1}+f_{2}=30\,$\mathrm{cm}$+10\,$\mathrm{cm}$=40\,$\mathrm{cm}. The width of the rectangular slit sitting in the Fourier plane of the system can be adjusted independently in vertical and horizontal direction from [math] to $12\,$\mathrm{mm}. The image formed by the first two lenses is further demagnified by lens 3 and the objective of the AIS. The focus of the whole system is adjusted by lens 3 which is mounted on a motorized stage. This is necessary such that the AIS and the dipole light can be focused independently. The total magnification of the DMD imaging system is $1/25.9$ and it has a minimum spot size (Airy diameter) of $4\,$\mathrm{\SIUnitSymbolMicro m}$$ in the plane of atoms. This implies that more than 70 DMD pixels contribute to the Airy disk at a single point. The DMD chip is illuminated by a Gaussian beam whose horizontal extension is elongated by a cylindrical telescope. This is done to match the available intensity to the size of the cloud and thereby efficiently use the available power. The polarization of the dipole potential beam is also set such that the beam gets reflected in the PBS cube in the second demagnification stage with minimal loss.

4 Real-time pattern optimization

With the optical setup described in the previous section, an area of interest (AOI) on the DMD chip with $n_{y}\times n_{z}$ pixels (typically $20\times 500$ pixels for a $200\,$\mathrm{\SIUnitSymbolMicro m}$$ long 1D trap) will be imaged onto the plane of atoms. With these pixels, $2^{(n_{y}\times\,n_{z})}$ different patterns can be created. The challenge is to find a pattern which produces a 1D optical dipole potential which results in the desired target potential $V_{\mathrm{T}}(z)$ when superposed with the magnetic confinement of the atom chip. In principle, we can calculate the 2D image of DMD pattern numerically and design the pattern to obtain desired potentials. However, due to interferences the actual dipole potential is sensitive to precise alignment of optical elements and imperfections in the setup and is difficult to predict precisely in a real world setting. Hence, an iterative pattern optimization is used to obtain the best possible reconstruction of the target potential.

In the experiment, we measure density profiles $\rho(z_{k})$. Thus, the ultimate goal of our pattern optimization process is to find a DMD pattern which creates a 1D density profile $\rho(z_{k})$ as close as possible to a simulated target density $\rho_{\mathrm{T}}(z_{k})$. The target density is obtained by calculating the ground state of a 1D non-polynomial Schrödinger equation (NPSE) [38] with the external target potential $V_{\mathrm{T}}(z_{k})$.

Before starting the pattern optimization, a calibration is performed to infer the corresponding position in the plane of atoms for the different DMD pixels in horizontal and vertical directions separately. For horizontal (longitudinal) calibration, a number of DMD columns are turned on at different positions. Their position is then detected on the 1D density and used to derive a linear map (Fig. 5(a)). To perform vertical calibration, a simple barrier is created by a thin super pixel in the center of the atomic cloud (typically $3$ rows with $3$ to $50$ pixels each). Shifting this super pixel vertically changes the barrier height. At the point where this height is maximum, the center of the super pixel is selected as the central row of the DMD pattern (Fig. 5(b)). This row corresponds to the position of the elongated atomic cloud ($y=y_{\mathrm{cloud}}$) in the plane of atoms and has the maximum contribution to $I(z)$. The vertical calibration also defines the vertical AOI, the area outside where the DMD pixels do not contribute to $I(z)$.

The process of optimizing the DMD pattern starts by choosing an initial pattern $p^{0}$. This pattern can be either totally black (all DMD pixels are off), a pattern that was calculated to correct the actual potential $V(z_{k})$ estimated from a measurement of $\rho(z_{k})$, or a pattern from a previous optimization processes. In each iteration $i$, with a DMD pattern $p^{i}$, the average 1D density profile, ${\bar{\rho}}^{i}(z_{k})$, and the mean atom number, $\bar{N}^{i}$, are obtained from $n$ absorption pictures. The mean atom number $\bar{N}^{i}$ is then used to calculate $\rho_{\mathrm{T}}^{i}(z_{k})$.

To estimate how well the measured density ${\bar{\rho}}^{i}(z_{k})$ and the calculated target density $\rho_{\mathrm{T}}^{i}(z_{k})$ agree, we calculate the normalized root-mean-square deviation,

[TABLE]

in each iteration for a specific region from $z_{j}$ to $z_{l}$ and compare it with a target value, $\epsilon_{\mathrm{T}}$. Ultimately, the quality required for a specific target potential, depends on the application. Note that however, there is a lower limit for $\epsilon_{\mathrm{T}}$ which is given by the noise in the density measurement: $\epsilon_{\mathrm{T}}$ should not be smaller than the average of the standard error of the mean density, $\sigma_{\bar{\rho}}$, in the same region.

We terminate the optimization process when $\epsilon_{\mathrm{RMS}}^{i}<\epsilon_{\mathrm{T}}$. In that case, the final DMD pattern is $p^{i}$. If $\epsilon_{\mathrm{RMS}}^{i}\geq\epsilon_{\mathrm{T}}$ a new pattern $p^{i+1}$ is generated using an update procedure discussed below and then the process is repeated.

In the DMD pattern update procedure, the states of individual pixels are altered based on the difference between $\bar{\rho}^{i}(z_{k})$ and $\rho_{T}^{i}(z_{k})$. In regions where the measured density is higher than the target density more pixels are switched on, thus yielding a locally higher potential. Conversely, by switching pixels off or shifting active pixels away from the center, the optical dipole potential in the corresponding region is decreased resulting in higher densities. The logic behind the optimization is based on the one-to-one correspondence between the trapping potential and the (Thomas-Fermi) density distribution (see Eq. (2)). In a trapped degenerate ensemble with constant number of atoms, locally increasing (decreasing) the potential decreases (increases) the local density distribution.

In Fig. 6, the image of seventeen rows of the DMD chip is shown in the plane of atoms, where the zeroth row is the central row discussed before. The orange dashed lines represent the camera pixels mapped in the plane of atoms, defining the spatial grid points $z_{k}$. The length of each grid point ($1.05\,$\mu\mathrm{m}) is equivalent to the length of almost $2.5$ DMD pixels mapped in the plane of atoms ($2.5\times 10.8\,$\mu\mathrm{m}$/25.9=1.04\,$\mu\mathrm{m}).

For each grid point $z_{k}$, the corresponding region on the DMD remains unchanged if $|\bar{\rho}^{i}(z_{k})-\rho_{T}^{i}(z_{k})|<\sigma^{i}_{\bar{\rho}}(z_{k})$ where $\sigma^{i}_{\bar{\rho}}(z_{k})$ is the standard error of the mean density, represented by the error bars in this figure. In our example here, this is the case for $z_{1}$, $z_{6}$, $z_{7}$, $z_{8}$, $z_{11}$, $z_{12}$ and $z_{13}$. In the other cases where $|\bar{\rho}^{i}(z_{k})-\rho_{T}^{i}(z_{k})|>\sigma^{i}_{\bar{\rho}}(z_{k})$, the magnitude of the total deviation determines how many pixels corresponding to grid point $z_{k}$ will be changed.

If $\bar{\rho}^{i}(z_{k})-\sigma^{i}_{\bar{\rho}}(z_{k})>\rho_{T}^{i}(z_{k})$, the potential has to be increased by turning additional pixels on. To chose these, first the distances of the outermost active pixels for the columns involved are compared together to select the column for which this distance is the smallest. Consider three columns involved in the region defined by $z_{10}$: the maximum distance of an active pixel from the central row is only 4 for the column in the left. However, for the central column and for the column in the right this number is $5$. So in this case the left column will be chosen. In the selected column, if the outermost pixel is located below the central row, the pixel above the central row which is one row further away will be turned on, if the outermost pixel is above the central row, then the pixel with the same distance but below the central row will be turned on. For $z_{10}$, the pixel marked with the black dot will be turned on. For $z_{9}$, the pixel marked by the orange dot will be turned on.

If $\bar{\rho}^{i}(z_{k})+\sigma^{i}_{\bar{\rho}}(z_{k})<\rho_{T}^{i}(z_{k})$, the potential has to be decreased by shifting the outermost pixel away from the central row in each region. For example in the region corresponding to $z_{3}$ ($z_{5}$) the pixel in the row number $-5$ ($-6$) will be shifted to the row number $6$ ($7$) in the same column, as marked in the Fig. 6. Note that in principle, the pixel $-5$ ($-6$) could also be shifted to the row $-6$ ($-7$). Although due to the interference results may be different, choosing one over the other is arbitrary. If the selected pixel is in the border of AOI it will be simply turned off.

5 Results

We used the technique described in Sec. 4 to realize a variety of 1D trapping potentials including a box-like potential, two box-like potentials separated with barriers with different heights, a box-like potential with sinusoidal modulations and a V-shaped potential. Results are shown in Figs. 7 and 8.

In these examples, the optimization process started with different initial DMD patterns. In the first and the fourth example (Figs. 7(a) and 7(d)), the optimization started with all DMD pixels turned off. Thus the initial density profiles (potential) plotted in Figs. 7(a) and 7(d) (orange) represent the 1D magnetic potential. In the examples presented in Figs. 7(b) and 7(c) we started with a previously optimized DMD pattern.

The initial and final potentials, plotted in Fig. 7, are obtained using the NPSE for 1D condensates  [38]. In the limit where the kinetic energy of motion in $z$-direction is neglected (Thomas-Fermi approximation), the solution for the stationary state reads:

[TABLE]

Here, $\rho(z)$ is the measured 1D density, $a_{\mathrm{s}}$ is the s-wave scattering length, $g=4\pi\hbar^{2}a_{\mathrm{s}}/M$ is the 3D scattering amplitude, $a_{\perp}$ is the oscillator length mentioned in Sec. 2, and $\mu$ is the chemical potential which constitutes only an offset. Although this measured potential is not used in the optimization process, a qualitative comparison between the measured potential and the target potential is instructive.

For the examples presented here, the number of experimental realizations was varied in each iteration depending on the normalized root-mean-squared deviation of the previous step. This helps speed up the optimization process in the beginning where only a rough measurement of the density is needed. At the initial stages of the process, only $n=2$ shots were used to calculate the averaged 1D density and the mean atom number. At the more advanced stages of the optimization this number was increased to $n=5$ and $n=10$ accordingly. For all measurements, the optical dipole trap was turned on during the evaporative cooling in a transverse double-well. Absorption pictures are taken after $2\,$\mathrm{ms}$$ TOF.

In Fig. 9, the normalized root-mean-squared deviation in each iteration, $\epsilon^{i}_{\mathrm{RMS}}$, is plotted for different examples. The number of iterations needed to fulfill the termination criterion was different depending on initial and target potentials. For example, a local modification to an existing box-like potential (adding a barrier) was done faster than achieving the box-like potential starting with a harmonic potential. Given the fact that each experimental realization takes approximately $25\,$\mathrm{s}$$, the total optimization time for these examples was between $60$ to $135$ minutes.

To check the stability of achieved potentials after optimization, $\epsilon^{i}_{\mathrm{RMS}}$ is plotted over a time span of more than six hours in Fig. 10(a). In this case, the optimization stopped after $40$ iterations and the achieved pattern remained unchanged through the rest of the measurement. In Fig. 10(b), the orange curve is the normalized density measured immediately after the optimization stopped. Fig. 10(b) shows a good agreement between normalized densities after $200$ (green), $400$ minutes (blue) and the orange curve.

6 Conclusion

We have implemented a robust and versatile way to control the longitudinal confinement of 1D Bose gases created on an atom chip. This is accomplished by modifying the longitudinal potential landscape defined by the magnetic fields by adding an optical dipole potential that is locally controlled by a DMD. We have furthermore demonstrated an autonomous optimization procedure to achieve any given longitudinal potential.

Our system opens up many new possibilities for experiments with 1D quantum systems. A box-like potential allows to create homogeneous systems, additional barriers will allow longitudinal splitting and tunneling. Switching between patterns will allow quench experiments. For example, imprinting a density modulation corresponding to the density quadrature of a specific mode will allow to excite this mode precisely and look at its evolution and decay. Furthermore, DMD patterns can be changed faster than the typical time scales of the 1D quantum gas. This will allow us to implement time varying potentials. One application will be to build and explore thermal quantum machines with 1D quantum fields [39].

Funding

Fonds zur Förderung der wissenschaftlichen Forschung (FWF) (I 3010-N27, 3980-N27, W 1210-N25); John Templeton foundation (JTF) (60478); Wiener Wissenschafts- und Technologie Fonds (WWTF) (MA16-066); Fundação para a Ciência e Tecnologia, Portugal (PD/BD/128641/2017); Österreichische Akademie der Wissenshaften (ÖAW) (ESQ Fellowship 801110)

Acknowledgment

The authors would like to thank Kaspar Sakmann who contributed to the initial phase of this project. We also thank Andrew N. Kanagin for helping proof read the manuscript.