Rothkopf, Alexander
2020-04-08
<p>//*********************************************************************************<br>
//<br>
// THIS PROGRAM IMPLEMENTS A 1D VERSION OF THE MASTER EQUATION FOR HEAVY<br>
// QUARKONIUM DERIVED BY Y. AKAMATSU IN<br>
//<br>
// ``Heavy quark master equations in the Lindblad form at high temperatures,''<br>
// Phys. Rev. D 91, no. 5, 056002 (2015)<br>
// doi:10.1103/PhysRevD.91.056002<br>
// [arXiv:1403.5783 [hep-ph]].<br>
//<br>
// ON THE LEVEL OF THE DENSITY MATRIX AND INCLUDES FULL DISSIPATIVE EFFECTS.<br>
// THE SAME EQUATION OF MOTION HAS BEEN APPROXIMATELY UNRAVELLED VIA QUANTUM<br>
// STATE DIFFUSION IN<br>
//<br>
// T. Miura, Y. Akamatsu, M. Asakawa and A. Rothkopf,<br>
// ``Quantum Brownian motion of a heavy quark pair in the quark-gluon plasma,''<br>
// Phys. Rev. D 101, no.3, 034011 (2020)<br>
// doi:10.1103/PhysRevD.101.034011<br>
// [arXiv:1908.06293 [nucl-th]].<br>
//<br>
// This code deploys a novel discrete finite difference operator, specifically<br>
// designed to guarantee the trace conservation of the Lindblad dynamics<br>
// called a reparametrization neutral summation-by-parts operator. Its derivation<br>
// is outlined in an upcoming joint publication<br>
//<br>
// O.Ålund, Y. Akamatsu, F. Laurén, T. Miura, J. Nordström and A. Rothkopf<br>
// ``Trace preserving quantum dynamics using a novel reparametrization-neutral<br>
// summation-by-parts difference operator.''<br>
//<br>
// CODE AUTHOR: Dr. Alexander Karl Rothkopf DATE: April 4th 2020<br>
//<br>
//*********************************************************************************</p>
<p># ==================================================================<br>
# In order to compile the different source code files you need to<br>
# install the PETSC and SLEPC library<br>
# from https://www.mcs.anl.gov/petsc/ (tested on v.3.12)<br>
# and from https://slepc.upv.es (tested on v.3.12)<br>
# respectively. The Makefile requires that the environment variables<br>
# PETSC_DIR, PETSC_ARCH and SLEPC_DIR are set accordingly. The values<br>
# currently used in the Makefile thus need to be adjusted.<br>
#<br>
# The Lindblad dynamics code allows the user to compute the eigen-<br>
# functions of the Hamiltonian directly without resorting to a<br>
# distributed solver, which furthermore requires the EIGEN library to<br>
# be present, which can be found at<br>
# http://eigen.tuxfamily.org/ (tested on v.3.3.7)<br>
# Please adjust the EIGENPATH environmental variable accordingly.<br>
#<br>
# MasterPETSC1DxyNewD: Lindblad dynamics with RN-SBP operator<br>
#<br>
# run the program via example script file # FLB_1D_PETSC_XY_DEBYE_TM_T01_IC0.scr<br>
#<br>
# MasterPETSC1DxyNaiveD: Lindblad dynamics with naive operator<br>
#<br>
# run the program via example script file<br>
# FLB_1D_PETSC_XY_DEBYE_TM_T01_IC0_NAIVED.scr<br>
#<br>
# EigenStates1DNewD: Distributed computation of Eigenstates of<br>
# the Hamiltonian to be used in the Lindblad<br>
# dynamics as initial conditions and/or in<br>
# order to compute density matrix in the<br>
# state basis.<br>
# run the program as mpirun -n #nodes EigenStates1DNewD.run<br>
# ==================================================================</p>
The author of the code gladly acknowledges support by the Research Council of Norway under the FRIPRO Young Research Talent grant 286883. Developing and testing this code has utilized computing resources provided by
UNINETT Sigma2 - the National Infrastructure for High Performance Computing and Data Storage in Norway under project NN9578K-QCDrtX "Real-time dynamics of nuclear matter under extreme conditions"
https://doi.org/10.5281/zenodo.3744461
oai:zenodo.org:3744461
eng
Zenodo
https://doi.org/10.5281/zenodo.3744460
info:eu-repo/semantics/openAccess
Creative Commons Attribution 4.0 International
https://creativecommons.org/licenses/by/4.0/legalcode
Quarkonium
Lindblad equation
Crank-Nicolson
reparametrization neutral
summation-by-parts
open quantum system
1d High Temperature Quarkonium Lindblad Dynamics CN & RN-SBP Solver
info:eu-repo/semantics/other