Charge waves and dynamical signatures of topological phases in Su–Schrieffer–Heeger chains

The file main.f90 contains the program code written in Fortran90 for numerical calculations of stationary LDOS (Local Density of States) and charge distribution along atomic SSH chains placed on a mettalic surface. The calculations are performed within the tight-binding Hamiltonian framework and the Green function technique.The file main2.f90 contains the program code written in Fortran90 for numerical calculations of time dependent LDOS and dynamics of charge distribution along atomic SSH chains placed on a mettalic surface. The calculations are performed within the tight-binding Hamiltonian framework and the evolution operator method.The attached files contain the data used to produce the figures presented in the paper: T.Kwapinski, M.Kurzyna, L. E. F. Foa Torres, Charge waves and dynamical signatures of topological phases in SSH chains, accepted for publication in Physical Review B. Data for figure fig01.eps of this repository (fig1 in the main text of the paper).The files 'fig01A_0lad' (data for panel A), 'fig01B_0lad' (data for panel B) and 'fig01C_0lad' (data for panel C) include charge occupancies along the normal chain (panel A), SSH0 chain (panel B) and SSH1 chain (panel C) of the length N=50 as a function of the site number for differnt values of the on-site energies ranging from e_0=-5 up to +5:column1: site numbercolumn2: charge at each site of the chaincolumn3: normalization of LDOS functioncolumn4: on-site energy Data for figure fig02.jpg of this repository (fig2 in the main text of the paper).The files 'fig02A_vlad1_n3', 'fig02A_vlad11_n3' (data for panel A)'fig02B_vlad1_n5', 'fig02B_vlad11_n5' (data for panel B)'fig02C_vlad1_n7', 'fig02C_vlad11_n7' (data for panel C) and'fig02C_vlad1_n10', 'fig02C_vlad11_n10' (data for panel D) include charge occupancies along the SSH1 chain of length N=30 (panels A and B) and N=50 (panels C and D) as a function of the site number for the average filling of the chain equals: 1/3, 1/5, 1/7 and 1/10:column1: site numbercolumn2: charge at each site of the chaincolumn3: normalization of LDOS function Data for figure fig03.jpg of this repository (fig10 in the main text of the paper).The files 'fig03A_vlad1_n5', 'fig03AA_vlad1_n5' and fig03AAA_vlad1_n5 (data for panel A) include charge occupancies along the SSH1 chain of length N=30 (panels A) as a function of the site number for different geometry of the system as described in the text:column1: site numbercolumn2: charge at each site of the chaincolumn3: normalization of LDOS functionThe files 'fig03B_vdos1_n5', 'fig03BB_vdos1_n5' and fig03BBB_vdos1_n5 (data for panel B) include normalized Local Density of States as a function of energy for a chain of N=30 linear atoms corresponding to data from panel A:column1: energycolumn2: LDOS function at each site of the chaincolumn3: site number Data for figure fig04.eps of this repository (fig3 in the main text of the paper).The file 'fig04A_0dos14' includes normalized LDOS as a function of energy for a chian of N=30 linear atoms coupled directly with the surface with the couplig parameters t_1=1 and t_2=4 for th on-site energies e_1=0 and e_2=0:column1: energycolumn2: Local Density of States for e1 sitecolumn3: Local Density of States for e2 sitecolumn4: Local Density of States for e15 sitecolumn5: Local Density of States for e3 sitecolumn6: Local Density of States for e4 sitecolumn7: Local Density of States for e5 sitecolumn8: Local Density of States for e7 siteThe file 'fig04A_0dos14e2' includes normalized LDOS as a function of energy for a chian of N=30 linear atoms coupled directly with the surface with the couplig parameters t_1=1 and t_2=4 for th on-site energies e_1=0 and e_2=2:column1: energycolumn2: Local Density of States for e1 sitecolumn3: Local Density of States for e2 sitecolumn4: Local Density of States for e15 sitecolumn5: Local Density of States for e3 sitecolumn6: Local Density of States for e4 sitecolumn7: Local Density of States for e5 sitecolumn8: Local Density of States for e7 siteThe file 'fig04A_0dos14e4' includes normalized LDOS as a function of energy for a chian of N=30 linear atoms coupled directly with the surface with the couplig parameters t_1=1 and t_2=4 for th on-site energies e_1=0 and e_2=4:column1: energycolumn2: Local Density of States for e1 sitecolumn3: Local Density of States for e2 sitecolumn4: Local Density of States for e15 sitecolumn5: Local Density of States for e3 sitecolumn6: Local Density of States for e4 sitecolumn7: Local Density of States for e5 sitecolumn8: Local Density of States for e7 site Data for figure fig05.eps of this repository (fig4 in the main text of the paper).The files 'fig05A_xlad1_n7', 'fig05A_xlad1e4_n7' and 'fig05A_xlad1e4_n10' include charge occupancies along the SSH1 chain of length N=50 as a function of the site number for the average filling of the chain equals: 1/7 or 1/10:column1: site numbercolumn2: charge at each site of the chaincolumn3: normalization of LDOS function Data for figure fig06.jpg of this repository (fig5 in the main text of the paper).The file 'fig06ABC_dos14t'includes local DOS for a chain of N = 30 sites, evaluated at sites i = 1, 2 and 15 (panel a, b and c) as a function of energy and time, following a sudden change in the coupling parameters at t = 30:column1: timecolumn2: energycolumn3: LDOS function at 1-st sitecolumn4: LDOS function at 2-st sitecolumn5: LDOS function at 15-st site Data for figure fig08.eps of this repository (fig7 in the main text of the paper).The files 'fig08A_lad14' (data for panel A) and 'fig08B_lad14' (data for panel B) include time evolution of the charge occupancies, n_i(t) at sites i = 1, 2, 3, 4 in a chain of length N = 30 following a sudden change in the coupling parameters at time t = 30. Initially, the system is a uniform chain with couplings t_i = 2 (panel a) and t_i = 4 (panel b), and after the quench, the couplings are changed to those of the SSH1 chain with (t1, t2) = (1, 4):column1: timecolumn2: charge at each chain sitecolumn3: normalized chargecolumn4: site number of the chain Data for figure fig09.eps of this repository (fig8 in the main text of the paper).The files 'fig09_lad14', 'fig09_lad14x1', 'fig09_lad14x2', 'fig09_lad14x3', 'fig09_lad14x4' and 'fig09_lad14x5' include time evolution of the charge occupancy at the end atom, n_1(t), following a sudden change in the coupling parameters at time t = 30, in a chain of length N = 30, for different on-site energies e_i = −2.5, 0, 2.5, and 4.4 and for the transition from a normal chain to the SSH1 (SSH0) topological phase, i.e., the couplings change from t_i = 4 to (t1, t2) = (1, 4) or (4, 1):column1: timecolumn2: charge at each chain sitecolumn3: normalized chargecolumn4: site number of the chain Data for figure fig10.eps of this repository (fig9 in the main text of the paper).The files 'fig10A_lad14x' (data for panel A) and 'fig10_lad14x' (data for panel B) include time-dependent occupancies ni(t) at the first four sites (i = 1, 2, 3, 4) in a chain of length N = 30, following a sudden change in the coupling parameter t_23 (between sites 2 and 3), which was set to zero in the time interval from t = 30 up to t = 40. Initially, the system is in the topologically trivial (panel a) and nontrivial (panel b) phase, with couplings(t1, t2) = (4, 1) and (t1, t2) = (1, 4), respectively, and the on-site energies are e_i = 4.4:column1: timecolumn2: charge at each chain sitecolumn3: normalized chargecolumn4: site number of the chain **************************** END OF FILE