Researcher profile

Kaspar Sakmann

Kaspar Sakmann contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2015arXiv

Single shot simulations of dynamic quantum many-body systems

The single-particle density is the most basic quantity that can be calculated from a given many-body wave function. It provides the probability to find a particle at a given position when the average over many realizations of an experiment is taken. However, the outcome of single experimental shots of ultracold atom experiments is determined by the $N$-particle probability density. This difference can lead to surprising results. For example, independent Bose-Einstein condensates (BECs) with definite particle numbers form interference fringes even though no fringes would be expected based on the single-particle density [1-4]. By drawing random deviates from the $N$-particle probability density single experimental shots can be simulated from first principles [1, 3, 5]. However, obtaining expressions for the $N$-particle probability density of realistic time-dependent many-body systems has so far been elusive. Here, we show how single experimental shots of general ultracold bosonic systems can be simulated based on numerical solutions of the many-body Schrödinger equation. We show how full counting distributions of observables involving any number of particles can be obtained and how correlation functions of any order can be evaluated. As examples we show the appearance of interference fringes in interacting independent BECs, fluctuations in the collisions of strongly attractive BECs, the appearance of randomly fluctuating vortices in rotating systems and the center of mass fluctuations of attractive BECs in a harmonic trap. The method described is broadly applicable to bosonic many-body systems whose phenomenology is driven by information beyond what is typically available in low-order correlation functions.

preprint2014arXiv

Universality of Fragmentation in the Schrödinger Dynamics of Bosonic Josephson Junctions

The many-body Schrödinger dynamics of a one-dimensional bosonic Josephson junction is investigated for up to ten thousand bosons and long times. The initial states are fully condensed and the interaction strength is weak. We report on a universal fragmentation dynamics on the many-body level: systems consisting of different numbers of particles fragment to the same value at constant mean-field interaction strength. The phenomenon manifests itself in observables such as the correlation functions of the system. We explain this universal fragmentation dynamics analytically based on the Bose-Hubbard model. We thereby show that the extent to which many-body effects become important at later times depends crucially on the initial state. Even for arbitrarily large particle numbers and arbitrarily weak interaction strength the dynamics is many-body in nature and the fragmentation universal. There is no weakly interacting limit where the Gross-Piatevskii mean-field is valid for long times.

preprint2013arXiv

Excitation spectra of many-body systems by linear response: General theory and applications to trapped condensates

We derive a general linear-response many-body theory capable of computing excitation spectra of trapped interacting bosonic systems, e.g., depleted and fragmented Bose-Einstein condensates (BECs). To obtain the linear-response equations we linearize the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method, which provides a self-consistent description of many-boson systems in terms of orbitals and a state vector (configurations), and is in principle numerically-exact. The derived linear-response many-body theory, which we term LR-MCTDHB, is applicable to systems with interaction potentials of general form. From the numerical implementation of the LR-MCTDHB equations and solution of the underlying eigenvalue problem, we obtain excitations beyond available theories of excitation spectra, such as the Bogoliubov-de Gennes (BdG) equations. The derived theory is first applied to study BECs in a one-dimensional harmonic potential. The LR-MCTDHB method contains the BdG excitations and, also, predicts a plethora of additional many-body excitations which are out of the realm of standard linear response. In particular, our theory describes the exact energy of the higher harmonic of the first (dipole) excitation not contained in the BdG theory. We next study a BEC in a very shallow one-dimensional double-well potential. We find with LR-MCTDHB low-lying excitations which are not accounted for by BdG, even though the BEC has only little fragmentation and, hence, the BdG theory is expected to be valid. The convergence of the LR-MCTDHB theory is assessed by systematically comparing the excitation spectra computed at several different levels of theory.

preprint2012arXiv

Dynamics and symmetries of a repulsively bound atom pair in an infinite optical lattice

We investigate the dynamics of two bosons trapped in an infinite one-dimensional optical lattice potential within the framework of the Bose-Hubbard model and derive an exact expression for the wavefunction at finite time. As initial condition we chose localized atoms that are separated by a distance of $d$ lattice sites and carry a center of mass quasi-momentum. An initially localized pair ($d=0$) is found to be more stable as quantified by the pair probability (probability to find two atoms at the same lattice site) when the interaction and/or the center of mass quasi-momentum is increased. For initially separated atoms ($d \neq 0$) there exists an optimal interaction strength for pair formation. Simple expressions for the wavefunction, the pair probability and the optimal interaction strength for pair formation are computed in the limit of infinite time. Whereas the time-dependent wavefunction differs for values of the interaction strength that differ only by the sign, important observables like the density and the pair probability do not. With a symmetry analysis this behavior is shown to extend to the $N$-particle level and to fermionic systems. Our results provide a complementary understanding of the recently observed [Winkler \textit{et al.}, Nature (London) \textbf{441}, 853 (2006)] dynamical stability of atom pairs in a repulsively interacting lattice gas.

preprint2012arXiv

Exact quantum dynamics of bosons with finite-range time-dependent interactions of harmonic type

The exactly solvable quantum many-particle model with harmonic one- and two-particle interaction terms is extended to include time-dependency. We show that when the external trap potential and finite-range interparticle interaction have a time-dependency the exact solutions of the corresponding time-dependent many-boson Schrödinger equation are still available. We use these exact solutions to benchmark the recently developed multiconfigurational time-dependent Hartree method for bosons (MCTDHB) [Phys. Rev. Lett. {\bf 99}, 030402 (2007), Phys. Rev. A {\bf 77}, 033613 (2008)]. In particular, we benchmark the MCTDHB method for: (i) the ground state; (ii) the breathing many-body dynamics activated by a quench scenario where the interparticle interaction strength is suddenly turned on to a finite value; (iii) the non-equilibrium dynamic for driven scenarios where both the trap- and interparticle-interaction potentials are {\it time-dependent}. Excellent convergence of the ground state and dynamics is demonstrated. The great relevance of the self-consistency and time-adaptivity, which are the intrinsic features of the MCTDHB method, is demonstrated by contrasting the MCTDHB predictions and those obtained within the standard full configuration interaction method spanning the Fock space of the same size, but utilizing as one-particle basis set the fixed-shape eigenstates of the one-particle potential. Connections of the model's results to ultra-cold Bose-Einstein condensed systems are addressed.

preprint2011arXiv

Number fluctuations of cold spatially split bosonic objects

We investigate the number fluctuations of spatially split many-boson systems employing a theorem about the maximally and minimally attainable variances of an observable. The number fluctuations of many-boson systems are given for different numbers of lattice sites and both mean-field and many-body wave functions. It is shown which states maximize the particle number fluctuations, both in lattices and double-wells. The fragmentation of the states is discussed, and it is shown that the number fluctuations of some fragmented states are identical to those of fully condensed states.

preprint2011arXiv

Optimal time-dependent lattice models for nonequilibrium dynamics

Lattice models are abundant in theoretical and condensed-matter physics. Generally, lattice models contain time-independent hopping and interaction parameters that are derived from the Wannier functions of the noninteracting problem. Here, we present a new concept based on time-dependent Wannier functions and the variational principle that leads to optimal time-dependent lattice models. As an application, we use the Bose-Hubbard model with time-dependent Wannier functions to study a quench scenario involving higher bands. We find a separation of times scales in the dynamics and show that under some circumstances the multi-band nonequilibrium dynamics of a quantum system can be obtained essentially at the cost of a single-band model.

preprint2011arXiv

Recursive formulation of the multiconfigurational time-dependent Hartree method for fermions, bosons and mixtures thereof in terms of one-body density operators

The multiconfigurational time-dependent Hartree method (MCTDH) [Chem. Phys. Lett. {\bf 165}, 73 (1990); J. Chem. Phys. {\bf 97}, 3199 (1992)] is celebrating nowadays entering its third decade of tackling numerically-exactly a broad range of correlated multi-dimensional non-equilibrium quantum dynamical systems. Taking in recent years particles' statistics explicitly into account, within the MCTDH for fermions (MCTDHF) and for bosons (MCTDHB), has opened up further opportunities to treat larger systems of interacting identical particles, primarily in laser-atom and cold-atom physics. With the increase of experimental capabilities to simultaneously trap mixtures of two, three, and possibly even multiple kinds of interacting composite identical particles together, we set up the stage in the present work and specify the MCTDH method for such cases. Explicitly, the MCTDH method for systems with three kinds of identical particles interacting via all combinations of two- and three-body forces is presented, and the resulting equations-of-motion are briefly discussed. All four possible mixtures of fermions and bosons are presented in a unified manner. Particular attention is paid to represent the coefficients' part of the equations-of-motion in a compact recursive form in terms of one-body density operators only. The recursion utilizes the recently proposed Combinadic-based mapping for fermionic and bosonic operators in Fock space [Phys. Rev. A {\bf 81}, 022124 (2010)] and successfully applied and implemented within MCTDHB. Our work sheds new light on the representation of the coefficients' part in MCTDHF and MCTDHB without resorting to the matrix elements of the many-body Hamiltonian with respect to the time-dependent configurations. It suggests a recipe for efficient implementation of the schemes derived here for mixtures which is suitable for parallelization.

preprint2009arXiv

Accurate multi-boson long-time dynamics in triple-well periodic traps

To solve the many-boson Schrödinger equation we utilize the Multiconfigurational time-dependent Hartree method for bosons (MCTDHB). To be able to attack larger systems and/or to propagate the solution for longer times, we implement a parallel version of the MCTDHB method thereby realizing the recently proposed [Streltsov {\it et al.} arXiv:0910.2577v1] novel idea how to construct efficiently the result of the action of the Hamiltonian on a bosonic state vector. We study the real-space dynamics of repulsive bosonic systems made of N=12, 51 and 3003 bosons in triple-well periodic potentials. The ground state of this system is three-fold fragmented. By suddenly strongly distorting the trap potential, the system performs complex many-body quantum dynamics. At long times it reveals a tendency to an oscillatory behavior around a threefold fragmented state. These oscillations are strongly suppressed and damped by quantum depletions. In spite of the richness of the observed dynamics, the three time-adaptive orbitals of MCTDHB(M=3) are capable to describe the many-boson quantum dynamics of the system for short and intermediate times. For longer times, however, more self-consistent time-adaptive orbitals are needed to correctly describe the non-equilibrium many-body physics. The convergence of the MCTDHB($M$) method with the number $M$ of self-consistent time-dependent orbitals used is demonstrated.

preprint2009arXiv

Quantum dynamics of attractive versus repulsive bosonic Josephson junctions: Bose-Hubbard and full-Hamiltonian results

The quantum dynamics of one-dimensional bosonic Josephson junctions with attractive and repulsive interparticle interactions is studied using the Bose-Hubbard model and by numerically-exact computations of the full many-body Hamiltonian. A symmetry present in the Bose-Hubbard Hamiltonian dictates an equivalence between the evolution in time of attractive and repulsive Josephson junctions with attractive and repulsive interactions of equal magnitude. The full many-body Hamiltonian does not possess this symmetry and consequently the dynamics of the attractive and repulsive junctions are different.