Benchmarking a Novel Efficient Numerical Method for Localized 1D Fermi-Hubbard Systems on a Quantum Simulator

Quantum simulators have made a remarkable progress towards exploring the dynamics of many-body systems, many of which offer a formidable challenge to both theoretical and numerical methods. While state-of-the-art quantum simulators are in principle able to simulate quantum dynamics well outside the domain of classical computers, they are noisy and limited in the variability of the initial state of the dynamics and the observables that can be measured. Despite these limitations, here we show that such a quantum simulator can be used to in-effect solve for the dynamics of a many-body system. We develop an efficient numerical technique that facilitates classical simulations in regimes not accessible to exact calculations or other established numerical techniques. The method is based on approximations that are well suited to describe localized one-dimensional Fermi-Hubbard systems. Since this new method does not have an error estimate and the approximations do not hold in general, we use a neutral-atom Fermi-Hubbard quantum simulator with $L_{\text{exp}}\simeq290$ lattice sites to benchmark its performance in terms of accuracy and convergence for evolution times up to $700$ tunnelling times. We then use these approximations in order to derive a simple prediction of the behaviour of interacting Bloch oscillations for spin-imbalanced Fermi-Hubbard systems, which we show to be in quantitative agreement with experimental results. Finally, we demonstrate that the convergence of our method is the slowest when the entanglement depth developed in the many-body system we consider is neither too small nor too large. This represents a promising regime for near-term applications of quantum simulators.Comment: 24 pages, 10 figure

On the uniqueness of the surface sources of evoked potentials

The uniqueness of a surface density of sources localized inside a spatial region $R$ and producing a given electric potential distribution in its boundary $B_0$ is revisited. The situation in which $R$ is filled with various metallic subregions, each one having a definite constant value for the electric conductivity is considered. It is argued that the knowledge of the potential in all $B_0$ fully determines the surface density of sources over a wide class of surfaces supporting them. The class can be defined as a union of an arbitrary but finite number of open or closed surfaces. The only restriction upon them is that no one of the closed surfaces contains inside it another (nesting) of the closed or open surfaces.Comment: 16 pages, 5 figure

Observing non-ergodicity due to kinetic constraints in tilted Fermi-Hubbard chains

The thermalization of isolated quantum many-body systems is deeply related to fundamental questions of quantum information theory. While integrable or many-body localized systems display non-ergodic behavior due to extensively many conserved quantities, recent theoretical studies have identified a rich variety of more exotic phenomena in between these two extreme limits. The tilted one-dimensional Fermi-Hubbard model, which is readily accessible in experiments with ultracold atoms, emerged as an intriguing playground to study non-ergodic behavior in a clean disorder-free system. While non-ergodic behavior was established theoretically in certain limiting cases, there is no complete understanding of the complex thermalization properties of this model. In this work, we experimentally study the relaxation of an initial charge-density wave and find a remarkably long-lived initial-state memory over a wide range of parameters. Our observations are well reproduced by numerical simulations of a clean system. Using analytical calculations we further provide a detailed microscopic understanding of this behavior, which can be attributed to emergent kinetic constraints.Comment: accepted in Nature Communication

Fast Optimal Transport Averaging of Neuroimaging Data

Knowing how the Human brain is anatomically and functionally organized at the level of a group of healthy individuals or patients is the primary goal of neuroimaging research. Yet computing an average of brain imaging data defined over a voxel grid or a triangulation remains a challenge. Data are large, the geometry of the brain is complex and the between subjects variability leads to spatially or temporally non-overlapping effects of interest. To address the problem of variability, data are commonly smoothed before group linear averaging. In this work we build on ideas originally introduced by Kantorovich to propose a new algorithm that can average efficiently non-normalized data defined over arbitrary discrete domains using transportation metrics. We show how Kantorovich means can be linked to Wasserstein barycenters in order to take advantage of an entropic smoothing approach. It leads to a smooth convex optimization problem and an algorithm with strong convergence guarantees. We illustrate the versatility of this tool and its empirical behavior on functional neuroimaging data, functional MRI and magnetoencephalography (MEG) source estimates, defined on voxel grids and triangulations of the folded cortical surface.Comment: Information Processing in Medical Imaging (IPMI), Jun 2015, Isle of Skye, United Kingdom. Springer, 201

The Unique Determination of Neuronal Currents in the Brain via Magnetoencephalography

The problem of determining the neuronal current inside the brain from measurements of the induced magnetic field outside the head is discussed under the assumption that the space occupied by the brain is approximately spherical. By inverting the Geselowitz equation, the part of the current which can be reconstructed from the measurements is precisely determined. This actually consists of only certain moments of one of the two functions specifying the tangential part of the current. The other function specifying the tangential part of the current as well as the radial part of the current are completely arbitrary. However, it is also shown that with the assumption of energy minimization, the current can be reconstructed uniquely. A numerical implementation of this unique reconstruction is also presented

Non-Equilibrium Mass Transport in the 1D Fermi-Hubbard Model

We experimentally and numerically investigate the sudden expansion of fermions in a homogeneous one-dimensional optical lattice. For initial states with an appreciable amount of doublons, we observe a dynamical phase separation between rapidly expanding singlons and slow doublons remaining in the trap center, realizing the key aspect of fermionic quantum distillation in the strongly-interacting limit. For initial states without doublons, we find a reduced interaction dependence of the asymptotic expansion speed compared to bosons, which is explained by the interaction energy produced in the quench

Nonequilibrium Mass Transport in the 1D Fermi-Hubbard Model.

We experimentally and numerically investigate the sudden expansion of fermions in a homogeneous one-dimensional optical lattice. For initial states with an appreciable amount of doublons, we observe a dynamical phase separation between rapidly expanding singlons and slow doublons remaining in the trap center, realizing the key aspect of fermionic quantum distillation in the strongly interacting limit. For initial states without doublons, we find a reduced interaction dependence of the asymptotic expansion speed compared to bosons, which is explained by the interaction energy produced in the quench

Dipole source analysis of auditory P300 response in depressive and anxiety disorders

This paper is to study auditory event-related potential P300 in patients with anxiety and depressive disorders using dipole source analysis. Auditory P300 using 2-stimulus oddball paradigm was collected from 35 patients with anxiety disorder, 32 patients with depressive disorder, and 30 healthy controls. P300 dipole sources and peak amplitude of dipole activities were analyzed. The source analysis resulted in a 4-dipole configuration, where temporal dipoles displayed greater P300 amplitude than that of frontal dipoles. In addition, a right-greater-than-left hemispheric asymmetry of dipole magnitude was found in patients with anxiety disorder, whereas a left-greater-than-right hemispheric asymmetry of dipole magnitude was observed in depressed patients. Results indicated that the asymmetry was more prominent over the temporal dipole than that of frontal dipoles in patients. Patients with anxiety disorder may increase their efforts to enhance temporal dipole activity to compensate for a deficit in frontal cortex processing, while depressed patients show dominating reduction of right temporal activity. The opposite nature of results observed with hemispheric asymmetry in depressive and anxiety disorders could serve to be valuable information for psychiatric studies

Neural Signatures of Stimulus Features in Visual Working Memory—A Spatiotemporal Approach

We examined the neural signatures of stimulus features in visual working memory (WM) by integrating functional magnetic resonance imaging (fMRI) and event-related potential data recorded during mental manipulation of colors, rotation angles, and color–angle conjunctions. The N200, negative slow wave, and P3b were modulated by the information content of WM, and an fMRI-constrained source model revealed a progression in neural activity from posterior visual areas to higher order areas in the ventral and dorsal processing streams. Color processing was associated with activity in inferior frontal gyrus during encoding and retrieval, whereas angle processing involved right parietal regions during the delay interval. WM for color–angle conjunctions did not involve any additional neural processes. The finding that different patterns of brain activity underlie WM for color and spatial information is consistent with ideas that the ventral/dorsal “what/where” segregation of perceptual processing influences WM organization. The absence of characteristic signatures of conjunction-related brain activity, which was generally intermediate between the 2 single conditions, suggests that conjunction judgments are based on the coordinated activity of these 2 streams