Trace preserving quantum dynamics using a novel reparametrization-neutral summation-by-parts difference operator

Ålund, Oskar; Akamatsu, Yukinao; Laurén, Fredrik; Miura, Takahiro; Nordström, Jan; Rothkopf, Alexander

doi:10.1016/j.jcp.2020.109917

Trace preserving quantum dynamics using a novel reparametrization-neutral summation-by-parts difference operator

Oskar Ålund
(
- Linkoping U., Math. Inst.
)
,
Yukinao Akamatsu
(
- Osaka U., Inst. Phys.
)
,
Fredrik Laurén
(
- Linkoping U., Math. Inst.
)
,
Takahiro Miura
(
- Osaka U., Inst. Phys.
)
,
Jan Nordström
(
- Linkoping U., Math. Inst. and
- Johannesburg U.
)

Apr 9, 2020

34 pages

Published in:

J.Comput.Phys. 425 (2021) 109917,
J.Comput.Phys. 425 (2021) 109917

Published: Oct 13, 2020
and
Published: Jan 15, 2021

e-Print:

2004.04406 [physics.comp-ph]

DOI:

10.1016/j.jcp.2020.109917 (publication)

View in:

pdf

reference search13 citations

Citations per year

Abstract: (Elsevier Inc.)

We develop a novel numerical scheme for the simulation of dissipative quantum dynamics, following from two-body Lindblad master equations. It exactly preserves the trace of the density matrix and shows only mild deviations from hermiticity and positivity, which are the defining properties of the continuum Lindblad dynamics. The central ingredient is a new spatial difference operator, which not only fulfills the summation by parts (SBP) property, but also implements a continuum reparametrization property. Using the time evolution of a heavy-quark anti-quark bound state in a hot thermal medium as an explicit example, we show how the reparametrization neutral summation-by-parts (RN-SBP) operator enables an accurate simulation of the full dissipative dynamics of this open quantum system. •New difference operator with continuum property: neutrality under coordinate transformation

(x, y) \to (z = x - y, z^{'} = x + y)

.•First direct simulation of the high temperature Lindblad equation for quarkonium particles with exact trace preservation.

Note:

34 pages, 7 figures, open-access code available via https://doi.org/10.5281/zenodo.3744460

Time integration
Initial boundary value problems
Dissipative systems
Open quantum systems
Summation-by-parts operators
Mimetic operator
reparametrization
dissipation
master equation
density matrix

References(46)

Figures(8)

[1]

A. Lampo
,
S.H. Lim
,
M.A. Garcá-March
,
M. Lewenstein

- Quantum 1 (2017) 30
•
e-Print:
- 1704.07623
•
DOI:
- 10.22331/q-2017-09-27-30

[2]

Heavy Quarkonium in Extreme Conditions

Alexander Rothkopf
(
- Stavanger U.
)

- Phys.Rept. 858 (2020) 1-117
•
e-Print:
- 1912.02253
•
DOI:
- 10.1016/j.physrep.2020.02.006

[3]

H.P. Breuer
,
F. Petruccione

[4]

J.J. Sakurai
,
J. Napolitano

[5]

Completely Positive Dynamical Semigroups of N Level Systems

Vittorio Gorini
(
- Texas U.
)
,
Andrzej Kossakowski
(
- Texas U.
)
,
E.C.G. Sudarshan
(
- Texas U. and
- Bangalore, Indian Inst. Sci.
)

- J.Math.Phys. 17 (1976) 821
•
DOI:
- 10.1063/1.522979

[6]

On the Generators of Quantum Dynamical Semigroups

Goran Lindblad
(
- Royal Inst. Tech., Stockholm
)

- Commun.Math.Phys. 48 (1976) 119,
- Commun.Math.Phys.48 (1976) 119-130 and Stockholm Royal Tech Inst - TRITA-TFY-75-1 (75,REC.MAR) 34p
•
DOI:
- 10.1007/BF01608499

[7]

M. Svärd
,
J. Nordström

- J.Comput.Phys. 268 (2014) 17-38

[8]

D.C.D.R. Fernández
,
J.E. Hicken
,
D.W. Zingg

- Comput. Fluids 95 (2014) 171-196

[9]

T. Lundquist
,
J. Nordström

- J.Comput.Phys. 270 (2014) 86-104

[10]

J. Nordström
,
S. Eriksson
,
P. Eliasson

- J.Comput.Phys. 231 (2012) 14, 4867-4884

[11]

J. Nordström
,
K. Forsberg
,
C. Adamsson
,
P. Eliasson

- Appl.Numer.Math. 45 (2003) 4, 453-473

[12]

M.H. Carpenter
,
T.C. Fisher
,
E.J. Nielsen
,
S.H. Frankel

- SIAM J.Sci.Comput. 36 (2014) 5, B835-B867

[13]

M.H. Carpenter
,
D. Gottlieb

- J.Comput.Phys. 129 (1996) 1, 74-86

[14]

P. Castonguay
,
D.M. Williams
,
P.E. Vincent
,
A. Jameson

- Comput.Methods Appl.Mech.Eng. 267 (2013) 400-417

[15]

H.T. Huynh

[16]

H. Ranocha
,
P. Öffner
,
T. Sonar

- J.Comput.Phys. 311 (2016) 299-328

[17]

G.J. Gassner

- SIAM J.Sci.Comput. 35 (2013) 3, A1233-A1253

[18]

J.S. Hesthaven
,
D. Gottlieb

- SIAM J.Sci.Comput. 17 (1996) 3, 579-612

[19]

J. Nordström

- J.Sci.Comput. 71 (2017) 1, 365-385

[20]

K. Mattsson
,
J. Nordström

- J.Comput.Phys. 199 (2004) 2, 503-540

[21]

K. Mattsson

- J.Comput.Phys. 274 (2014) 432-454

[22]

J. Nordström
,
T. Lundquist

- J.Comput.Phys. 251 (2013) 487-499

[23]

J. Nordström
,
T. Lundquist

- SIAM J.Sci.Comput. 38 (2016) 3, A1561-A1586

[24]

M.H. Carpenter
,
D. Gottlieb
,
S. Abarbanel

- J.Comput.Phys. 111 (1994) 2, 220-236

[25]

Real-time quantum dynamics of heavy quark systems at high temperature

Yukinao Akamatsu
(
- KMI, Nagoya and
- Nagoya U. (main)
)

- Phys.Rev.D 87 (2013) 4, 045016
•
e-Print:
- 1209.5068
•
DOI:
- 10.1103/PhysRevD.87.045016

1-25 of 46
1
2
25 / page