Giovanna Morigi

Ion Coulomb crystals are ordered structures formed by laser-cooled ions in traps that are characterized by interparticle distances of several micrometers and energy scales on the order of $μ$eV. Their crystalline structure emerges from the interplay between Coulomb repulsion and the external confining potential, which can be readily tuned. Moreover, individual ions can be precisely manipulated with lasers and imaged via resonance fluorescence. These unusual and unique properties make ion crystals a powerful platform for studying phases of matter in the strongly correlated regime and at low temperatures where their dynamics is manifestly quantum mechanical. This review examines the theoretical framework and experimental characterization of ion Coulomb crystals from a condensed-matter perspective. We discuss their dynamical and thermodynamic properties in one, two, and three dimensions, and review recent investigations into their out-of-equilibrium behavior. We provide outlooks on future directions for exploring novel condensed matter phenomena with trapped ion crystals, as well as for exploiting these features for scientific and technical applications.

Trapped ions in a periodic potential are a paradigm of a frustrated Wigner crystal. The dynamics is captured by a long-range Frenkel-Kontorova model. The classical ground state can be mapped to the one of an antiferromagnetic spin chain with long-range interactions in a magnetic field, whose strength is determined by the mismatch between chain's and substrate lattice's periodicity. The mapping is exact when the substrate potential is a piecewise harmonic potential and holds for any two-body interaction decaying as $1/r^α$ with the distance $r$. The ground state is a devil's staircase of regular, periodic structures as a function of the mismatch, whose range of stability depends also on the coefficient $α$. While the staircase is well defined in the thermodynamic limit for $α>1$, for Coulomb interactions, $α=1$, it disappears and the sliding-to-pinned transitions becomes crossovers. However, due to the logarithmic convergence to the thermodynamic limit characteristic of the Coulomb potential, the staircase is found for any finite number of ions. We discuss the experimental parameters as well as the features that allow one to observe and reveal our predictions in experimental platforms. These dynamics are a showcase of the versatility of trapped ion platforms for exploring the interplay between frustration and interactions.

The ability for real-time control of topological defects can open up prospects for dynamical manipulation of macroscopic properties of solids. A sub-category of these defects, formed by particle dislocations, can be effectively described using the Frenkel-Kontorova chain, which characterizes the dynamics of these particles in a periodic lattice potential. This model is known to host solitons, which are the topological defects of the system and are linked to structural transitions in the chain. This work addresses three key questions: Firstly, we investigate how imperfections present in concrete implementations of the model affect the properties of topological defects. Secondly, we explore how solitons can be injected after the rapid change in lattice potential or nucleated due to quantum fluctuations. Finally, we analyze the propagation and scattering of solitons, examining the role of quantum fluctuations and imperfections in influencing these processes. Furthermore, we address the experimental implementation of the Frenkel-Kontorova model. Focusing on the trapped ion quantum simulator, we set the stage for controllable dynamics of topological excitations and their observation in this platform.

A Wigner chain in a periodic potential is a paradigmatic example of geometric frustration with long-range interactions. The dynamics emulates the Frenkel-Kontorova model with Coulomb interactions. In the continuum approximation, dislocations are sine-Gordon solitons with power-law decaying tails. We show that their action is mapped into a massive, long-range (1+1) Thirring model, where the solitons are charged fermionic excitations over an effective Dirac sea. We identify the corresponding mean field theory and show that the Coulomb interactions destabilize structures commensurate with the periodic substrate, suppressing their onset and giving rise to {\it interaction-induced} lubrication. Our study identifies the role of long-range interactions on determining nanofriction. Our predictions can be probed in state-of-the-art trapped ion experiments.

Topological defects have strong impact on both elastic and inelastic properties of materials. In this article, we investigate the possibility to controllably inject topological defects in quantum simulators of solid state lattice structures. We investigate the quench dynamics of a Frenkel-Kontorova chain, which is used to model discommensurations of particles in cold atoms and trapped ionic crystals. The interplay between an external periodic potential and the inter-particle interaction makes lattice discommensurations, the topological defects of the model, energetically favorable and can tune a commensurate-incommensurate structural transition. Our key finding is that a quench from the commensurate to incommensurate phase causes a controllable injection of topological defects at periodic time intervals. We employ this mechanism to generate quantum states which are a superposition of lattice structures with and without topological defects. We conclude by presenting concrete perspectives for the observation and control of topological defects in trapped ion experiments.

We investigate the normal mode spectrum of a trapped ion chain at the symmetry-breaking linear to zigzag transition and at finite temperatures. For this purpose we modulate the amplitude of the Doppler cooling laser in order to excite and measure mode oscillations. The expected mode softening at the critical point, a signature of the second-order transition, is not observed. Numerical simulations show that this is mainly due to the finite temperature of the chain. Inspection of the trajectories suggest that the thermal shifts of the normal-mode spectrum can be understood by the ions collectively jumping between the two ground state configurations of the symmetry broken phase. We develop an effective analytical model, which allows us to reproduce the low-frequency spectrum as a function of the temperature and close to the transition point. In this model the frequency shift of the soft mode is due to the anharmonic coupling with the high frequency modes of the spectrum, acting as an averaged effective thermal environment. Our study could prove important for implementing ground-state laser cooling close to the critical point.

We theoretically analyse the equation of topological solitons in a chain of particles interacting via a repulsive power-law potential and confined by a periodic lattice. Starting from the discrete model, we perform a gradient expansion and obtain the kink equation in the continuum limit for a power law exponent $n \ge 1$. The power-law interaction modifies the sine-Gordon equation, giving rise to a rescaling of the coefficient multiplying the second derivative (the kink width) and to an additional integral term. We argue that the integral term does not affect the local properties of the kink, but it governs the behaviour at the asymptotics. The kink behaviour at the center is dominated by a sine-Gordon equation and its width tends to increase with the power law exponent. When the interaction is the Coulomb repulsion, in particular, the kink width depends logarithmically on the chain size. We define an appropriate thermodynamic limit and compare our results with existing studies performed for infinite chains. Our formalism allows one to systematically take into account the finite-size effects and also slowly varying external potentials, such as for instance the curvature in an ion trap.

Slow variations (quenches) of the magnetic field across the paramagnetic-ferromagnetic phase transition of spin systems produce heat. In systems with short-range interactions the heat exhibits universal power-law scaling as a function of the quench rate, known as Kibble-Zurek scaling. In this work we analyze slow quenches of the magnetic field in the Lipkin-Meshkov-Glick (LMG) model, which describes fully connected quantum spins. We analytically determine the quantum contribution to the residual heat as a function of the quench rate $δ$ by means of a Holstein-Primakoff expansion about the mean-field value. Unlike in the case of short-range interactions, scaling laws in the LMG model are only found for a ramp ending at the critical point. If instead the ramp is symmetric, as in the typical Kibble-Zurek scenario, after crossing the critical point the system tends to reabsorb the defects formed during the first part of the ramp: the number of excitations exhibits a crossover behavior as a function of $δ$ and tends to a constant in the thermodynamic limit. Previous, and seemingly contradictory, theoretical studies are identified as specific limits of this dynamics. Our results can be tested on several experimental platforms, including quantum gases and trapped ions.

Crystals of repulsively interacting ions in planar traps form hexagonal lattices, which undergo a buckling instability towards a multi-layer structure as the transverse trap frequency is reduced. Numerical and experimental results indicate that the new structure is composed of three planes, whose separation increases continuously from zero. We study the effects of thermal and quantum fluctuations by mapping this structural instability to the six-state clock model. A prominent implication of this mapping is that at finite temperature, fluctuations split the buckling instability into two thermal transitions, accompanied by the appearance of an intermediate critical phase. This phase is characterized by quasi-long-range order in the spatial tripartite pattern. It is manifested by broadened Bragg peaks at new wave vectors, whose line-shape provides a direct measurement of the temperature dependent exponent $η(T)$ characteristic of the power-law correlations in the critical phase. A quantum phase transition is found at the largest value of the critical transverse frequency: here the critical intermediate phase shrinks to zero. Moreover, within the ordered phase, we predict a crossover from classical to quantum behavior, signifying the emergence of an additional characteristic scale for clock order. We discuss experimental realizations with trapped ions and polarized dipolar gases, and propose that within accessible technology, such experiments can provide a direct probe of the rich phase diagram of the quantum clock model, not easily observable in condensed matter analogues. Therefore, this works highlights the potential for ionic and dipolar systems to serve as simulators for complex models in statistical mechanics and condensed matter physics.

A protocol is discussed for preparing a spin chain in a generic many-body state in the asymptotic limit of tailored non-unitary dynamics. The dynamics require the spectral resolution of the target state, optimized coherent pulses, engineered dissipation, and feedback. As an example, we discuss the preparation of an entangled antiferromagnetic state, and argue that the procedure can be applied to chains of trapped ions or Rydberg atoms.

The dynamics of cold trapped ions in a high-finesse resonator results from the interplay between the long-range Coulomb repulsion and the cavity-induced interactions. The latter are due to multiple scatterings of laser photons inside the cavity and become relevant when the laser pump is sufficiently strong to overcome photon decay. We study the stationary states of ions coupled with a mode of a standing-wave cavity as a function of the cavity and laser parameters, when the typical length scales of the two self-organizing processes, Coulomb crystallization and photon-mediated interactions, are incommensurate. The dynamics are frustrated and in specific limiting cases can be cast in terms of the Frenkel-Kontorova model, which reproduces features of friction in one dimension. We numerically recover the sliding and pinned phases. For strong cavity nonlinearities, they are in general separated by bistable regions where superlubric and stick-slip dynamics coexist. The cavity, moreover, acts as a thermal reservoir and can cool the chain vibrations to temperatures controlled by the cavity parameters and by the ions phase. These features are imprinted in the radiation emitted by the cavity, which is readily measurable in state-of-art setups of cavity quantum electrodynamics.

The efficient initialization of a quantum system is a prerequisite for quantum technological applications. Here we show that several classes of quantum states of a harmonic oscillator can be efficiently prepared by means of a Jaynes-Cummings interaction with a single two-level system. This is achieved by suitably tailoring external fields which drive the dipole and/or the oscillator. The time-dependent dynamics that leads to the target state is identified by means of Optimal Control Theory (OCT) based on Krotov's method. Infidelities below $10^{-4}$ can be reached for the parameters of the experiment of the ENS group in Paris, where the oscillator is a mode of a high-Q microwave cavity and the dipole is a Rydberg transition of an atom. For this specific situation we analyze the limitations on the fidelity due to parameter fluctuations and identify robust dynamics based on pulses found using ensemble OCT. Our analysis can be extended to quantum-state preparation of continuous-variable systems in other platforms, such as trapped ions and circuit QED.

We study the crossover from classical to quantum phase transitions at zero temperature within the framework of $φ^4$ theory. The classical transition at zero temperature can be described by the Landau theory, turning into a quantum Ising transition with the addition of quantum fluctuations. We perform a calculation of the transition line in the regime where the quantum fluctuations are weak. The calculation is based on a renormalization group analysis of the crossover between classical and quantum transitions, and is well controlled even for space-time dimensionality $D$ below 4. In particular, for $D=2$ we obtain an analytic expression for the transition line which is valid for a wide range of parameters, as confirmed by numerical calculations based on the Density Matrix Renormalization Group. This behavior could be tested by measuring the phase diagram of the linear-zigzag instability in systems of trapped ions or repulsively-interacting dipoles.

Ions of the same charge inside confining potentials can form crystalline structures which can be controlled by means of the ions density and of the external trap parameters. In particular, a linear chain of trapped ions exhibits a transition to a zigzag equilibrium configuration, which is controlled by the strength of the transverse confinement. Studying this phase transition in the quantum regime is a challenging problem, even when employing numerical methods to simulate microscopically quantum many-body systems. Here we present a compact analytical treatment to map the original long-range problem into a short-range quantum field theory on a lattice. We provide a complete numerical architecture, based on Density Matrix Renormalization Group, to address the effective quantum phi-four model. This technique is instrumental in giving a complete characterization of the phase diagram, as well as pinpoint the universality class of the criticality.

We show how to measure the structural witnesses proposed in [P. Krammer et al., Phys. Rev. Lett. 103, 100502 (2009)] for detecting entanglement in a spin chain using photon scattering. The procedure, moreover, allows one to measure the two-point correlation function of the spin array. This proposal could be performed in existing experimental platforms realizing ion chains in Paul traps or atomic arrays in optical lattices.

Trapped singly-charged ions can crystallize as a result of laser cooling. The emerging structure depends on the number of particles and on the geometry of the trapping potential. In linear multipole radiofrequency traps, the geometry of the radial potential can lead to the formation of single-ring structures. We analyse the conditions and stability of single rings as a function of the number of poles. For larger numbers of ions the rings form tubes in which the arrangement of the ions corresponds to a triangular lattice folded onto a cylinder. The stability of these tubular structures is numerically studied for different lattice constants and their normal mode spectrum is determined.

We analyze the dynamics of a chain of singly-charged ions confined in a linear Paul trap and which couple with the mode of a high-finesse optical resonator. In these settings the ions interact via the Coulomb repulsion and are subject to the mechanical forces due to scattering of cavity photons. We show that the interplay of these interactions can give rise to bistable equilibrium configurations, into which the chain can be cooled by cavity-enhanced photon scattering. We characterize the resulting equilibrium structures by determining the stationary state in the semiclassical limit for both cavity field and crystal motion. The mean occupation of the vibrational modes at steady state is evaluated, showing that the vibrational modes coupled to the cavity can be simultaneously cooled to low occupation numbers. It is also found that at steady state the vibrations are entangled with the cavity field fluctuations. The entanglement is quantified by means of the logarithmic negativity. The spectrum of the light at the cavity output is evaluated and the features signaling entanglement are here identified.

We theoretically show how two impurity defects in a crystalline structure can be entangled through coupling with the crystal. We demonstrate this with a harmonic chain of trapped ions in which two ions of a different species are embedded. Entanglement is found for sufficiently cold chains and for a certain class of initial, separable states of the defects. It results from the interplay between localized modes which involve the defects and the interposed ions, it is independent of the chain size, and decays slowly with the distance between the impurities. These dynamics can be observed in systems exhibiting spatial order, viable realizations are optical lattices, optomechanical systems, or cavity arrays in circuit QED.

We study the nonequilibrium dynamics of the linear to zigzag structural phase transition exhibited by an ion chain confined in a trap with periodic boundary conditions. The transition is driven by reducing the transverse confinement at a finite quench rate, which can be accurately controlled. This results in the formation of zigzag domains oriented along different transverse planes. The twists between different domains can be stabilized by the topology of the trap and under laser cooling the system has a chance to relax to a helical chain with nonzero winding number. Molecular dynamics simulations are used to obtain a large sample of possible trajectories for different quench rates. The scaling of the average winding number with different quench rates is compared to the prediction of the Kibble-Zurek theory, and a good quantitative agreement is found.

The generation of entanglement between two oscillators that interact via a common reservoir is theoretically studied. The reservoir is modeled by a one-dimensional harmonic crystal initially in thermal equilibrium. Starting from a separable state, the oscillators can become entangled after a transient time, that is of the order of the thermalization time scale. This behavior is observed at finite temperature even when the oscillators are at a distance significantly larger than the crystal's interparticle spacing. The underlying physical mechanisms can be explained by the dynamical properties of the collective variables of the two oscillators which may decouple from or be squeezed by the reservoir. Our predictions can be tested with an ion chain in a linear Paul trap.

A quasi one--dimensional system of trapped, repulsively interacting atoms (e.g., an ion chain) exhibits a structural phase transition from a linear chain to a zigzag structure, tuned by reducing the transverse trap potential or increasing the particle density. Since it is a one dimensional transition, it takes place at zero temperature and therefore quantum fluctuations dominate. In [Fishman, et al., Phys. Rev. B 77, 064111 (2008)] it was shown that the system close to the linear-zigzag instability is described by a $φ^4$ model. We propose a mapping of the $φ^4$ field theory to the well known Ising chain in a transverse field, which exhibits a quantum critical point. Based on this mapping, we estimate the quantum critical point in terms of the system parameters. This estimate gives the critical value of the transverse trap frequency for which the quantum phase transition occurs, and which has a finite, measurable deviation from the critical point evaluated within the classical theory. A measurement is suggested for atomic systems which can probe the critical trap frequency at sufficiently low temperatures T. We focus in particular on a trapped ion system, and estimate the implied limitations on T and on the interparticle distance. We conclude that the experimental observation of the quantum critical behavior is in principle accessible.

This work theoretically addresses the trapping an ionized atom with a single valence electron by means of lasers, analyzing qualitatively and quantitatively the consequences of the net charge of the particle. In our model, the coupling between the ion and the electromagnetic field includes the charge monopole and the internal dipole, within a multipolar expansion of the interaction Hamiltonian. Specifically, we perform a Power-Zienau-Woolley transformation, taking into account the motion of the center of mass. The net charge produces a correction in the atomic dipole which is of order $m_e/M$ with $m_e$ the electron mass and $M$ the total mass of the ion. With respect to neutral atoms, there is also an extra coupling to the laser field which can be approximated by that of the monopole located at the position of the center of mass. These additional effects, however, are shown to be very small compared to the dominant dipolar trapping term.

A string of trapped ions at zero temperature exhibits a structural phase transition to a zigzag structure, tuned by reducing the transverse trap potential or the interparticle distance. The transition is driven by transverse, short wavelength vibrational modes. We argue that this is a quantum phase transition, which can be experimentally realized and probed. Indeed, by means of a mapping to the Ising model in a transverse field, we estimate the quantum critical point in terms of the system parameters, and find a finite, measurable deviation from the critical point predicted by the classical theory. A measurement procedure is suggested which can probe the effects of quantum fluctuations at criticality. These results can be extended to describe the transverse instability of ultracold polar molecules in a one dimensional optical lattice.

Structural defects in ion crystals can be formed during a linear quench of the transverse trapping frequency across the mechanical instability from a linear chain to the zigzag structure. The density of defects after the sweep can be conveniently described by the Kibble-Zurek mechanism. In particular, the number of kinks in the zigzag ordering can be derived from a time-dependent Ginzburg-Landau equation for the order parameter, here the zigzag transverse size, under the assumption that the ions are continuously laser cooled. In a linear Paul trap the transition becomes inhomogeneous, being the charge density larger in the center and more rarefied at the edges. During the linear quench the mechanical instability is first crossed in the center of the chain, and a front, at which the mechanical instability is crossed during the quench, is identified which propagates along the chain from the center to the edges. If the velocity of this front is smaller than the sound velocity, the dynamics becomes adiabatic even in the thermodynamic limit and no defect is produced. Otherwise, the nucleation of kinks is reduced with respect to the case in which the charges are homogeneously distributed, leading to a new scaling of the density of kinks with the quenching rate. The analytical predictions are verified numerically by integrating the Langevin equations of motion of the ions, in presence of a time-dependent transverse confinement. We argue that the non-equilibrium dynamics of an ion chain in a Paul trap constitutes an ideal scenario to test the inhomogeneous extension of the Kibble-Zurek mechanism, which lacks experimental evidence to date.

Ultracold polar molecules, in highly anisotropic traps and interacting via a repulsive dipolar potential, may form one-dimensional chains at high densities. According to classical theory, at low temperatures there exists a critical value of the density at which a second order phase transition from a linear to a zigzag chain occurs. We study the effect of thermal and quantum fluctuations on these self-organized structures using classical and quantum Monte Carlo methods, by means of which we evaluate the pair correlation function and the static structure factor. Depending on the parameters, these functions exhibit properties typical of a crystalline or of a liquid system. We compare the thermal and the quantum results, identifying analogies and differences. Finally, we discuss experimental parameter regimes where the effects of quantum fluctuations on the linear - zigzag transition can be observed.

The light scattered by a cold trapped ion, which is in the stationary state of laser cooling, presents features due to the mechanical effects of atom-photon interaction. These features appear as additional peaks (sidebands) in the spectrum of resonance fluorescence. Among these sidebands the literature has discussed the Stokes and anti-Stokes components, namely the sidebands of the elastic peak. In this manuscript we show that the motion also gives rise to sidebands of the inelastic peaks. These are not always visible, but, as we show, can be measured in parameter regimes which are experimentally accessible.

In order to use a collection of trapped ions for experiments where a well defined preparation of vibrational states is necessary, all vibrational modes have to be cooled to ensure precise and repeatable manipulation of the ions' quantum states. A method for simultaneous sideband cooling of all axial vibrational modes is proposed. By application of a magnetic field gradient the absorption spectrum of each ion is modified such that sideband resonances of different vibrational modes coincide. The ion string is then irradiated with monochromatic electromagnetic radiation, in the optical or microwave regime, for sideband excitation. This cooling scheme is investigated in detailed numerical studies. Its application for initializing ion strings for quantum information processing is extensively discussed.

The density of ions trapped in a harmonic potential in one dimension is not uniform. Consequently the eigenmodes are not phonons. We calculate the long wavelength modes in the continuum limit, and evaluate the density of states in the short wavelength limit for chains of $N\gg 1$ ions. Remarkably, the results that are found analytically in the thermodynamic limit provide a good estimate of the spectrum of excitations of small chains down to few tens of ions. The spectra are used to compute the thermodynamic functions of the chain. Deviations from extensivity of the thermodynamic quantities are found. An analytic expression for the critical transverse frequency determining the stability of a linear chain is derived.

We show that the unitary evolution of a harmonic oscillator coupled to a two-level system can be undone by a suitable manipulation of the two-level system -- more specifically: by a quasi-instantaneous phase change. This enables us to isolate the dissipative evolution to which the oscillator may be exposed in addition. With this method we study the decoherence of a photon mode in cavity QED, and that of the quantized harmonic motion of trapped ions. We comment on the relation to spin echoes and multi-path interferometry.

We study theoretically Doppler laser-cooling of a cluster of 2-level atoms confined in a linear ion trap. Using several consecutive steps of averaging we derive, from the full quantum mechanical master equation, an equation for the total mechanical energy of the one dimensional crystal, defined on a coarse-grained energy scale whose grid size is smaller than the linewidth of the electronic transition. This equation describes the cooling dynamics for an arbitrary number of ions and in the quantum regime. We discuss the validity of the ergodic assumption (i.e. that the phase space distribution is only a function of energy). From our equation we derive the semiclassical limit (i.e. when the mechanical motion can be treated classically) and the Lamb-Dicke limit (i.e. when the size of the mechanical wave function is much smaller than the laser wavelength). We find a Fokker-Planck equation for the total mechanical energy of the system, whose solution is in agreement with previous analytical calculations which were based on different assumptions and valid only in their specific regimes. Finally, in the classical limit we derive an analytic expression for the average coupling, by light scattering, between motional states at different energies.