Focus on Strongly Interacting Quantum Gases in One Dimension

We discuss the hydrodynamic approach to the study of the time evolution—induced by a quench—of local excitations in one dimension. We focus on interaction quenches: the considered protocol consists of creating a stable localized excitation propagating through the system, and then operating a sudden change of the interaction between the particles. To highlight the effect of the quench, we take the initial excitation to be a soliton. The quench splits the excitation into two packets moving in opposite directions, whose characteristics for short times can be expressed in a universal way. Our treatment allows for the description of the internal dynamics of these two packets in terms of the different velocities of their components. We confirm our analytical predictions through numerical simulations performed with the Gross–Pitaevskii equation and with the Calogero model (as an example of long range interactions and solvable with a parabolic confinement). Through the Calogero model we also discuss the effect of an external trapping on the protocol. The hydrodynamic approach shows that there is a difference between the bulk velocities of the propagating packets and the velocities of their peaks: it is possible to discriminate the two quantities, as we show through the comparison between numerical simulations and analytical estimates. We show that our analytical results capture with remarkable precision the findings of the numerical simulations also for intermediate times and we provide predictions for the time at which the two packets becomes distinguishable. In the realizations of the discussed quench protocol in a cold atom experiment, these different velocities are accessible through different measurement procedures.

We investigate the applicability of laser assisted tunneling in a strongly interacting one-dimensional (1D) Bose gas (the Tonks–Girardeau gas) in optical lattices. We find that the stroboscopic dynamics of the Tonks–Girardeau gas in a continuous Wannier–Stark-ladder potential, supplemented with laser assisted tunneling, effectively realizes the ground state of 1D hard-core bosons in a discrete lattice with nontrivial hopping phases. We compare observables that are affected by the interactions, such as the momentum distribution, natural orbitals and their occupancies, in the time-dependent continuous system, to those of the ground state of the discrete system. Stroboscopically, we find an excellent agreement, indicating that laser assisted tunneling is a viable technique for realizing novel ground states and phases with hard-core 1D Bose gases.

While the emergent field of quantum thermodynamics has the potential to impact energy science, the performance of thermal machines is often classical. We ask whether quantum effects can boost the performance of a thermal machine to reach quantum supremacy, i.e., surpassing both the efficiency and power achieved in classical thermodynamics. To this end, we introduce a nonadiabatic quantum heat engine operating an Otto cycle with a many-particle working medium, consisting of an interacting Bose gas confined in a time-dependent harmonic trap. It is shown that thanks to the interplay of nonadiabatic and many-particle quantum effects, this thermal machine can outperform an ensemble of single-particle heat engines with same resources, demonstrating the quantum supremacy of many-particle thermal machines.

We present a series of quantum states that are characterized by dark solitons of the nonlinear Schrödinger equation (i.e. the Gross–Pitaevskii equation) for the one-dimensional Bose gas interacting through the repulsive delta-function potentials. The classical solutions satisfy the periodic boundary conditions and we simply call them classical dark solitons. Through exact solutions we show corresponding aspects between the states and the solitons in the weak coupling case: the quantum and classical density profiles completely overlap with each other not only at an initial time but also at later times over a long period of time, and they move together with the same speed in time; the matrix element of the bosonic field operator between the quantum states has exactly the same profiles of the square amplitude and the phase as the classical complex scalar field of a classical dark soliton not only at the initial time but also at later times, and the corresponding profiles move together for a long period of time. We suggest that the corresponding properties hold rigorously in the weak coupling limit. Furthermore, we argue that the lifetime of the dark soliton-like density profile in the quantum state becomes infinitely long as the coupling constant approaches zero, by comparing it with the quantum speed limit time. Thus, we call the quantum states quantum dark soliton states.

Yrast states are the lowest energy states at given non-zero momentum and provide a natural extension of the concept of dark solitons to strongly interacting one-dimensional quantum gases. Here we study the yrast states of the balanced spin-$\tfrac{1}{2}$ Fermi gas with attractive delta-function interactions in one-dimension with the exactly solvable Yang–Gaudin model. The corresponding Bethe-ansatz equations are solved for finite particle number and in the thermodynamic limit. Properties corresponding to the soliton-like nature of the yrast excitations are calculated including the missing particle number, phase step, and inertial and physical masses. The inertial to physical mass ratio, which is related to the frequency of oscillations in a trapped gas, is found to be unity in the limits of strong and weak attraction and falls to $\approx 0.78$ in the crossover regime. This result is contrasted by one-dimensional mean field theory, which predicts a divergent mass ratio in the weakly attractive limit. By means of an exact mapping our results also predict the existence and properties of dark-soliton-like excitations in the super Tonks–Girardeau gas. The prospects for experimental observations are briefly discussed.

Inspired by recent experiments on Bose–Einstein condensates in ring traps, we investigate the topological properties of the phase of a one-dimensional Bose field in the presence of both thermal and quantum fluctuations—the latter ones being tuned by the depth of an optical lattice applied along the ring. In the regime of large filling of the lattice, quantum Monte Carlo simulations give direct access to the full statistics of fluctuations of the Bose-field phase, and of its winding number W along the ring. At zero temperature the winding-number (or topological-sector) fluctuations are driven by quantum phase slips localized around a Josephson link between two lattice wells, and their susceptibility is found to jump at the superfluid-Mott insulator transition. At finite (but low) temperature, on the other hand, the winding number fluctuations are driven by thermal activation of nearly uniform phase twists, whose activation rate is governed by the superfluid fraction. A quantum-to-thermal crossover in winding-number fluctuations is therefore exhibited by the system, and it is characterized by a conformational change in the topologically non-trivial configurations, from localized to uniform phase twists, which can be experimentally observed in ultracold Bose gases via matter–wave interference.

A one-dimensional world is very unusual as there is an interplay between quantum statistics and geometry, and a strong short-range repulsion between atoms mimics Fermi exclusion principle, fermionizing the system. Instead, a system with a large number of components with a single atom in each, on the opposite acquires many bosonic properties. We study the ground-state properties of a multicomponent repulsive Fermi gas trapped in a harmonic trap by a fixed-node diffusion Monte Carlo method. The interaction between all components is considered to be the same. We investigate how the energetic properties (energy, contact) and correlation functions (density profile and momentum distribution) evolve as the number of components is changed. It is shown that the system fermionizes in the limit of strong interactions. Analytical expressions are derived in the limit of weak interactions within the local density approximation for an arbitrary number of components and for one plus one particle using an exact solution.

We overview the main features of mobile impurities moving in one-dimensional superfluid backgrounds by modeling it as a mobile Josephson junction, which leads naturally to the periodic dispersion of the impurity. The dissipation processes, such as radiative friction and quantum viscosity, are shown to result from the interaction of the collective phase difference with the background phonons. We develop a more realistic depleton model of an impurity-hole bound state that provides a number of exact results interpolating between the semiclassical weakly interacting picture and the strongly interacting Tonks–Girardeau regime. We also discuss the physics of a trapped impurity, relevant to current experiments with ultra cold atoms.

A boson two-leg ladder in the presence of a synthetic magnetic flux is investigated by means of bosonization techniques and density matrix renormalization group (DMRG). We follow the quantum phase transition from the commensurate Meissner to the incommensurate vortex phase with increasing flux at different fillings. When the applied flux is ρπ and close to it, where ρ is the filling per rung, we find a second incommensuration in the vortex state that affects physical observables such as the momentum distribution, the rung–rung correlation function and the spin–spin and charge–charge static structure factors.

The correlation function is an important quantity in the physics of ultracold quantum gases because it provides information about the quantum many-body wave function beyond the simple density profile. In this paper we first study the M-body local correlation functions, g_M, of the one-dimensional (1D) strongly repulsive Bose gas within the Lieb–Liniger model using the analytical method proposed by Gangardt and Shlyapnikov (2003 Phys. Rev. Lett. 90 010401; 2003 New J. Phys. 5 79). In the strong repulsion regime the 1D Bose gas at low temperatures is equivalent to a gas of ideal particles obeying the non-mutual generalized exclusion statistics with a statistical parameter $\alpha =1-2/\gamma$, i.e. the quasimomenta of N strongly interacting bosons map to the momenta of N free fermions via ${k}_{i}\approx \alpha {k}_{i}^{F}$ with $i=1,\ldots ,N$. Here γ is the dimensionless interaction strength within the Lieb–Liniger model. We rigorously prove that such a statistical parameter α solely determines the sub-leading order contribution to the M-body local correlation function of the gas at strong but finite interaction strengths. We explicitly calculate the correlation functions g_M in terms of γ and α at zero, low, and intermediate temperatures. For M = 2 and 3 our results reproduce the known expressions for g₂ and g₃ with sub-leading terms (see for instance (Vadim et al 2006 Phys. Rev. A 73 051604(R); Kormos et al 2009 Phys. Rev. Lett. 103 210404; Wang et al 2013 Phys. Rev. A 87 043634). We also express the leading order of the short distance non-local correlation functions $\langle {{\rm{\Psi }}}^{\dagger }({x}_{1})\cdots {{\rm{\Psi }}}^{\dagger }({x}_{M}){\rm{\Psi }}({y}_{M})\cdots {\rm{\Psi }}({y}_{1})\rangle$ of the strongly repulsive Bose gas in terms of the wave function of M bosons at zero collision energy and zero total momentum. Here ${\rm{\Psi }}(x)$ is the boson annihilation operator. These general formulas of the higher-order local and non-local correlation functions of the 1D Bose gas provide new insights into the many-body physics.

The ground-state properties of one-dimensional electron-spin-polarized hydrogen ¹H, deuterium ²H, and tritium ³H are obtained by means of quantum Monte Carlo methods. The equations of state of the three isotopes are calculated for a wide range of linear densities. The pair correlation function and the static structure factor are obtained and interpreted within the framework of the Luttinger liquid theory. We report the density dependence of the Luttinger parameter and use it to identify different physical regimes: Bogoliubov Bose gas, super-Tonks–Girardeau gas, and quasi-crystal regimes for bosons; repulsive, attractive Fermi gas, and quasi-crystal regimes for fermions. We find that the tritium isotope is the one with the richest behavior. Our results show unambiguously the relevant role of the isotope mass in the properties of this quantum system.

We consider a mixture of one-dimensional strongly interacting Fermi gases with up to six components, subjected to a longitudinal harmonic confinement. In the limit of infinitely strong repulsions we provide an exact solution which generalizes the one for the two-component mixture. We show that an imbalanced mixture under harmonic confinement displays partial spatial separation among the components, with a structure which depends on the relative population of the various components. Furthermore, we provide a symmetry characterization of the ground and excited states of the mixture introducing and evaluating a suitable operator, namely the conjugacy class sum. We show that, even under external confinement, the gas has a definite symmetry which corresponds to the most symmetric one compatible with the imbalance among the components. This generalizes the predictions of the Lieb–Mattis theorem for a Fermionic mixture with more than two components.

We study the dynamics of quantum bosonic waves in a one-dimensional tilted optical lattice. An effective spatially localized nonlinear two-body potential barrier is set at the center of the lattice. This version of the Bose–Hubbard model can be realized in atomic Bose–Einstein condensates, with the help of localized optical Feshbach resonance, controlled by a focused laser beam, and in quantum optics, using an arrayed waveguide with selectively doped guiding cores. Our numerical analysis demonstrates that the central barrier induces anomalous quantum reflection of incident wave packets, which acts solely on bosonic components with multiple onsite occupancies, while single-occupancy components pass the barrier, allowing one to distill them in the interaction zone. As a consequence, in this region one finds a hard-core-like state, in which the multiple occupancy is forbidden. Our results demonstrate that this regime can be attained dynamically, using relatively weak interactions, irrespective of their sign. Physical parameters necessary for the experimental implementation of the setting in ultracold atomic gases are estimated.

We show that an adiabatic cycle excites Bose particles confined in a one-dimensional box. During the adiabatic cycle, a wall described by a δ-shaped potential is applied and its strength and position are slowly varied. When the system is initially prepared in the ground state, namely, in the zero-temperature equilibrium state, the adiabatic cycle brings all Bosons into the first excited one-particle state, leaving the system in a nonequilibrium state. The absorbed energy during the cycle is proportional to the number of Bosons.

We present an asymptotically exact renormalization-group theory of the superfluid–insulator transition in one-dimensional (1D) disordered systems, with emphasis on an accurate description of the interplay between the Giamarchi–Schulz (instanton–anti-instanton) and weak-link (scratched-XY) criticalities. Combining the theory with extensive quantum Monte Carlo simulations allows us to shed new light on the ground-state phase diagram of the 1D disordered Bose–Hubbard model at unit filling.

A density-dependent gauge field may induce density-induced geometric frustration, leading to a non-trivial interplay between density modulation and frustration, which we illustrate for the particular case of ultra-cold bosons in zig-zag optical lattices with a density-dependent hopping amplitude. We show that the density-induced frustration leads to a rich landscape of quantum phases, including Mott insulator, bond-order insulator, two-component superfluids, chiral superfluids, and partially paired superfluids. We show as well that the density-dependent hopping results in an effective repulsive or attractive interaction, and that for the latter case the vacuum may be destabilized leading to a strong compressibility. Finally, we discuss the characteristic momentum distribution of the predicted phases, which can be used to detect the phases in time-of-flight measurements.

We have developed an efficient computational method to treat long, one-dimensional systems of strongly interacting atoms forming self-assembled spin chains. Such systems can be used to realize many spin chain model Hamiltonians tunable by the external confining potential. As a concrete demonstration, we consider quantum state transfer in a Heisenberg spin chain and we show how to determine the confining potential in order to obtain nearly perfect state transfer.

We use the coordinate Bethe ansatz to exactly calculate matrix elements between eigenstates of the Lieb–Liniger model of one-dimensional bosons interacting via a two-body delta-potential. We investigate the static correlation functions of the zero-temperature ground state and their dependence on interaction strength, and analyze the effects of system size in the crossover from few-body to mesoscopic regimes for up to seven particles. We also obtain time-dependent nonequilibrium correlation functions for five particles following quenches of the interaction strength from two distinct initial states. One quench is from the noninteracting ground state and the other from a correlated ground state near the strongly interacting Tonks–Girardeau regime. The final interaction strength and conserved energy are chosen to be the same for both quenches. The integrability of the model highly constrains its dynamics, and we demonstrate that the time-averaged correlation functions following quenches from these two distinct initial conditions are both nonthermal and moreover distinct from one another.

We investigate a free one-dimensional spinless Fermi gas, and the Tonks–Girardeau gas interacting with a single impurity particle of equal mass. We obtain a Fredholm determinant representation for the time-dependent correlation function of the impurity particle. This representation is valid for an arbitrary temperature and an arbitrary repulsive or attractive impurity-gas δ-function interaction potential. It includes, as particular cases, the representations obtained for zero temperature and arbitrary repulsion in (Gamayun et al 2015 Nucl. Phys. B 892 83–104), and for arbitrary temperature and infinite repulsion in (Izergin and Pronko 1998 Nucl. Phys. B 520 594–632).

Adiabatic techniques can be used to control quantum states with high fidelity while exercising limited control over the parameters of a system. However, because these techniques are slow compared to other timescales in the system, they are usually not suitable for creating highly unstable states or performing time-critical processes. Both of these situations arise in quantum information processing, where entangled states may be isolated from the environment only for a short time and where quantum computers require high-fidelity operations to be performed quickly. Recently it has been shown that techniques like optimal control and shortcuts to adiabaticity can be used to prepare quantum states non-adiabatically with high fidelity. Here we present two examples of how these techniques can be used to create maximally entangled many-body NOON states in one-dimensional Tonks–Girardeau gases.

Synthetic ladders realized with one-dimensional alkaline-earth(-like) fermionic gases and subject to a gauge field represent a promising environment for the investigation of quantum Hall physics with ultracold atoms. Using density-matrix renormalization group calculations, we study how the quantum Hall-like chiral edge currents are affected by repulsive atom–atom interactions. We relate the properties of such currents to the asymmetry of the spin resolved momentum distribution function, a quantity which is easily addressable in state-of-art experiments. We show that repulsive interactions significantly enhance the chiral currents. Our numerical simulations are performed for atoms with two and three internal spin states.

We consider the Tonks–Girardeau gas subject to a random external potential. If the disorder is such that the underlying one-particle Hamiltonian displays localization (which is known to be generically the case), we show that there is exponential decay of correlations in the many-body eigenstates. Moreover, there is no Bose–Einstein condensation and no superfluidity, even at zero temperature.

We elaborate a model of the interferometer which, unlike previously studied ones, uses a local (δ-functional) nonlinear repulsive potential, embedded into a harmonic-oscillator trapping potential, as the splitter for the incident soliton. An estimate demonstrates that this setting may be implemented by means of the localized Feshbach resonance controlled by a focused laser beam. The same system may be realized as a nonlinear waveguide in optics. Subsequent analysis produces an exact solution for scattering of a plane wave in the linear medium on the δ -functional nonlinear repulsive potential, and an approximate solution for splitting of the incident soliton when the ambient medium is nonlinear. The most essential result, obtained by means of systematic simulations, is that the use of the nonlinear splitter provides the sensitivity of the soliton-based interferometer to the target, inserted into one of its arms, which is much higher than the sensitivity provided by the usual linear splitter.

By using a diffusion Monte Carlo (DMC) technique, we studied the behaviour of a mixture of spin-up and spin-down cold fermionic atoms confined in one-dimensional harmonic potentials. We considered only small balanced clusters (up to eight atoms) or arrangements in which the difference between the population of fermions with different spins is one (with a maximum total number of nine). The atom–atom interactions were modeled by a contact repulsive potential. We focused on the ordering of the confined fermions, finding that the probability of having a pure antiferromagnetic arrangement of atoms decreases strongly with the number of fermions, as to be almost negligible for our largest clusters. Exactly the same can be said of a pure ferromagnetic state.

We study one-dimensional (1D) strongly interacting Bose–Fermi mixtures by both the exact Bethe-ansatz method and variational perturbation theory within the degenerate ground state subspace of the system in the infinitely repulsive limit. Based on the exact solution of the 1D Bose–Fermi gas with equal boson–boson and boson–fermion interaction strengths, we demonstrate that the ground state energy is degenerate for different Bose–Fermi configurations and the degeneracy is lifted when the interaction deviates the infinitely interacting limit. We then show that the ground properties in the strongly interacting regime can be well characterized by using the variational perturbation method within the degenerate ground state subspace, which can be applied to deal with more general cases with anisotropic interactions and in external traps. Our results indicate that the total ground-state density profile in the strongly repulsive regime behaves like the polarized non-interacting fermions, whereas the density distributions of bosons and fermions display different properties for different Bose–Fermi configurations and are sensitive to the anisotropy of interactions.

We consider one-dimensional mixtures of an atomic Bose–Einstein condensate (BEC) and Tonks–Girardeau (TG) gas. The mixture is modeled by a coupled system of the Gross–Pitaevskii equation for the BEC and the quintic nonlinear Schrödinger equation for the TG component. An immiscibility condition for the binary system is derived in a general form. Under this condition, three types of BEC–TG interfaces are considered: domain walls (DWs) separating the two components; bubble-drops (BDs), in the form of a drop of one component immersed into the other (BDs may be considered as bound states of two DWs); and bound states of bright and dark solitons (BDSs). The same model applies to the copropagation of two optical waves in a colloidal medium. The results are obtained by means of systematic numerical analysis, in combination with analytical Thomas–Fermi approximations (TFAs). Using both methods, families of DW states are produced in a generic form. BD complexes exist solely in the form of a TG drop embedded into the BEC background. On the contrary, BDSs exist as bound states of TG bright and BEC dark components, and vice versa.

We study the equilibrium properties of the one-dimensional disordered Bose–Hubbard model by means of a gauge-adaptive tree tensor network variational method suitable for systems with periodic boundary conditions. We compute the superfluid stiffness and superfluid correlations close to the superfluid to glass transition line, obtaining accurate locations of the critical points. By studying the statistics of the exponent of the power-law decay of the correlation, we determine the boundary between the superfluid region and the Bose glass phase in the regime of strong disorder and in the weakly interacting region, not explored numerically before. In the former case our simulations are in agreement with previous Monte Carlo calculations.

The exact quantum decay of a one-dimensional Bose gas with inverse-square interactions is presented. The system is equivalent to a gas of particles obeying generalized exclusion statistics. We consider the expansion dynamics of a cloud initially confined in a harmonic trap that is suddenly switched off. The decay is characterized by analyzing the fidelity between the initial and the time-evolving states, also known as the survival probability. It exhibits early on a quadratic dependence on time that turns into a power-law decay, during the course of the evolution. It is shown that the particle number and the strength of interactions determine the power-law exponent in the latter regime, as recently conjectured. The nonexponential character of the decay is linked to the many-particle reconstruction of the initial state from the decaying products.

We analyze the quantum phases, correlation functions and edge modes for a class of spin-1/2 and fermionic models related to the one-dimensional Ising chain in the presence of a transverse field. These models are the Ising chain with anti-ferromagnetic long-range interactions that decay with distance r as $1/{r}^{\alpha }$, as well as a related class of fermionic Hamiltonians that generalize the Kitaev chain, where both the hopping and pairing terms are long-range and their relative strength can be varied. For these models, we provide the phase diagram for all exponents α, based on an analysis of the entanglement entropy, the decay of correlation functions, and the edge modes in the case of open chains. We demonstrate that violations of the area law can occur for $\alpha \lesssim 1$, while connected correlation functions can decay with a hybrid exponential and power-law behavior, with a power that is α-dependent. Interestingly, for the fermionic models we provide an exact analytical derivation for the decay of the correlation functions at every α. Along the critical lines, for all models breaking of conformal symmetry is argued at low enough α. For the fermionic models we show that the edge modes, massless for $\alpha \gtrsim 1$, can acquire a mass for $\alpha \lt 1$. The mass of these modes can be tuned by varying the relative strength of the kinetic and pairing terms in the Hamiltonian. Interestingly, for the Ising chain a similar edge localization appears for the first and second excited states on the paramagnetic side of the phase diagram, where edge modes are not expected. We argue that, at least for the fermionic chains, these massive states correspond to the appearance of new phases, notably approached via quantum phase transitions without mass gap closure. Finally, we discuss the possibility to detect some of these effects in experiments with cold trapped ions.

We employ a high-order perturbative expansion to characterize the ground state of the Mott phase of the one-dimensional Bose–Hubbard model. We compute for different integer filling factors the energy per lattice site, the two-point and density–density correlations, and expectation values of powers of the on-site number operator determining the local atom number fluctuations (variance, skewness, kurtosis). We compare these expansions to numerical simulations of the infinite-size system to determine their range of applicability. We also discuss a new sum rule for the density–density correlations that can be used in both equilibrium and non-equilibrium systems.

We introduce a lattice model of interacting spins and bosons that leads to Luttinger-liquid physics, and allows for quantitative tests of the theory of bosonization by means of trapped-ion or superconducting-circuit experiments. By using a variational bosonization ansatz, we calculate the power-law decay of spin and boson correlation functions, and study their dependence on a single tunable parameter, namely a bosonic driving. For small drivings, matrix-product-states (MPSs) numerical methods are shown to be efficient and validate our ansatz. Conversely, even static MPS become inefficient for large-driving regimes, such that the experiment can potentially outperform classical numerics, achieving one of the goals of quantum simulations.

We employ ab initio methods of quantum chemistry to investigate spin-1/2 fermions interacting via a two-body contact potential in a one-dimensional harmonic trap. The convergence of the total energy with the size of the one-particle basis set is analytically investigated for the two-body problem and the same form of the convergence formula is numerically confirmed to be valid for the many-body case. Benchmark calculations for two to six fermions with the full configuration interaction method equivalent to the exact diagonalization approach, and the coupled cluster (CC) method including single, double, triple, and quadruple excitations are presented. The convergence of the correlation energy with the level of excitations included in the CC model is analyzed. The range of the interaction strength for which single-reference CC methods work is examined. Next, the CC method restricted to single, double, and noniterative triple excitations, CCSD(T), is employed to study a two-component Fermi gas composed of 6–80 atoms in a one-dimensional harmonic trap. The density profiles of trapped atomic clouds are also reported. Finally, a comparison with experimental results for few-fermion systems is presented. Upcoming possible applications and extensions of the presented approach are discussed.

We numerically generate and then study the basic properties of dark soliton-like excitations in a dipolar gas confined in a quasi one-dimensional trap. These excitations, although very similar to dark solitons in a gas with contact interaction, interact with each other and can form bound states. During collisions these dipolar solitons emit phonons, losing energy but accelerating. Even after thousands of subsequent collisions they survive as gray solitons and finally reach dynamical equilibrium with background quasiparticles. Finally, in the frame of classical field approximation, we verified, that these solitons appear spontaneously in thermal samples, analogously to the type II excitations in a gas of atoms with contact interaction.

In this article, we show that eigenenergies and eigenstates of a system consisting of four one-dimensional hard-core particles with masses 6m, 2m, m, and 3m in a hard-wall box can be found exactly using Bethe Ansatz. The Ansatz is based on the exceptional affine reflection group ${\tilde{F}}_{4}$ associated with the symmetries and tiling properties of an octacube—a Platonic solid unique to four-dimensions, with no three-dimensional analogues. We also uncover the Liouville integrability structure of our problem: the four integrals of motion in involution are identified as invariant polynomials of the finite reflection group F₄, taken as functions of the components of momenta.