The simple mode of DSQSS/DLA ISSP, UTokyo Kazuyoshi Yoshimi - PowerPoint PPT Presentation

The simple mode of DSQSS/DLA ISSP, UTokyo Kazuyoshi Yoshimi 2019-06-06 @ ISSP

2 The simple mode of DSQSS/dla Users can simulate of a predefined model 
 • on a predefined lattice from one text file (std.toml file). Schematic calculation flow of dla Input files for dla Input file Output file dla_pre: Input files 
 dla: Main engine generator for dla Automatically generated

Hamiltonian Jz = -1.0 # coupling constant, negative for AF Specify information of [hamiltonian] section 
 Specify information of lattice [lattice] section 
 bc = false # open boundary L = 2 # number of sites along each direction dim = 1 # dimension lattice = "hypercubic" # hypercubic, periodic [lattice] h = 0.0 # magnetic field Jxy = -1.0 # coupling constant, negative for AF M = 1 # S=1/2 model = "spin" [hamiltonian] <latexit sha1_base64="8ame4KnhU16fJze/+fjCaN0DsAM=">AChXichVG7SgNBFD2uj2h8RW0Em2BQbAx3RVRsFG3EyiTGCBpkd53EJftidxOIwR8QW7WwUrAQP8APsPEHLPIJYhnBxsKbzYJoUO8wM2fO3HNnDld1DN3zieodUmdXd0+kty/aPzA4NBwbGd3x7LKriaxmG7a7qyqeMHRLZH3dN8Su4wrFVA2RU0vrzftcRbieblvbftUReVMpWnpB1xSfqfTs5kEsQUkKIt4O5BAkEMaWHXvAPg5hQ0MZJgQs+IwNKPB47EGwWEujxpzLiM9uBc4QZS1Zc4SnKEwW+K1yKe9kLX43KzpBWqNXzF4uqyMY4qe6Y4a9ET39EIfv9aqBTWaf6nyra0wjkYPh3PvP+rMn3cfSl+kOhcvbfnwUsBR40dmbEzBNl1qrfuX4spFZTk/VpumGXtnfNdXpkR1alTftNiXSV4hyg+Sf7WgHO3NJmZJyaj6xuha2qhcTmMQM92MRq9jAFrL8bgFnOMeFJFmpXlpoZUqdYSaMXwLaeUT50eP+A=</latexit> <latexit sha1_base64="8ame4KnhU16fJze/+fjCaN0DsAM=">AChXichVG7SgNBFD2uj2h8RW0Em2BQbAx3RVRsFG3EyiTGCBpkd53EJftidxOIwR8QW7WwUrAQP8APsPEHLPIJYhnBxsKbzYJoUO8wM2fO3HNnDld1DN3zieodUmdXd0+kty/aPzA4NBwbGd3x7LKriaxmG7a7qyqeMHRLZH3dN8Su4wrFVA2RU0vrzftcRbieblvbftUReVMpWnpB1xSfqfTs5kEsQUkKIt4O5BAkEMaWHXvAPg5hQ0MZJgQs+IwNKPB47EGwWEujxpzLiM9uBc4QZS1Zc4SnKEwW+K1yKe9kLX43KzpBWqNXzF4uqyMY4qe6Y4a9ET39EIfv9aqBTWaf6nyra0wjkYPh3PvP+rMn3cfSl+kOhcvbfnwUsBR40dmbEzBNl1qrfuX4spFZTk/VpumGXtnfNdXpkR1alTftNiXSV4hyg+Sf7WgHO3NJmZJyaj6xuha2qhcTmMQM92MRq9jAFrL8bgFnOMeFJFmpXlpoZUqdYSaMXwLaeUT50eP+A=</latexit> <latexit sha1_base64="8ame4KnhU16fJze/+fjCaN0DsAM=">AChXichVG7SgNBFD2uj2h8RW0Em2BQbAx3RVRsFG3EyiTGCBpkd53EJftidxOIwR8QW7WwUrAQP8APsPEHLPIJYhnBxsKbzYJoUO8wM2fO3HNnDld1DN3zieodUmdXd0+kty/aPzA4NBwbGd3x7LKriaxmG7a7qyqeMHRLZH3dN8Su4wrFVA2RU0vrzftcRbieblvbftUReVMpWnpB1xSfqfTs5kEsQUkKIt4O5BAkEMaWHXvAPg5hQ0MZJgQs+IwNKPB47EGwWEujxpzLiM9uBc4QZS1Zc4SnKEwW+K1yKe9kLX43KzpBWqNXzF4uqyMY4qe6Y4a9ET39EIfv9aqBTWaf6nyra0wjkYPh3PvP+rMn3cfSl+kOhcvbfnwUsBR40dmbEzBNl1qrfuX4spFZTk/VpumGXtnfNdXpkR1alTftNiXSV4hyg+Sf7WgHO3NJmZJyaj6xuha2qhcTmMQM92MRq9jAFrL8bgFnOMeFJFmpXlpoZUqdYSaMXwLaeUT50eP+A=</latexit> <latexit sha1_base64="8ame4KnhU16fJze/+fjCaN0DsAM=">AChXichVG7SgNBFD2uj2h8RW0Em2BQbAx3RVRsFG3EyiTGCBpkd53EJftidxOIwR8QW7WwUrAQP8APsPEHLPIJYhnBxsKbzYJoUO8wM2fO3HNnDld1DN3zieodUmdXd0+kty/aPzA4NBwbGd3x7LKriaxmG7a7qyqeMHRLZH3dN8Su4wrFVA2RU0vrzftcRbieblvbftUReVMpWnpB1xSfqfTs5kEsQUkKIt4O5BAkEMaWHXvAPg5hQ0MZJgQs+IwNKPB47EGwWEujxpzLiM9uBc4QZS1Zc4SnKEwW+K1yKe9kLX43KzpBWqNXzF4uqyMY4qe6Y4a9ET39EIfv9aqBTWaf6nyra0wjkYPh3PvP+rMn3cfSl+kOhcvbfnwUsBR40dmbEzBNl1qrfuX4spFZTk/VpumGXtnfNdXpkR1alTftNiXSV4hyg+Sf7WgHO3NJmZJyaj6xuha2qhcTmMQM92MRq9jAFrL8bgFnOMeFJFmpXlpoZUqdYSaMXwLaeUT50eP+A=</latexit> 2 1 <latexit sha1_base64="UvOZW/4t5LzxWzoHwQCgzm7qQ=">ACrXichVE9SxRFD2OidlsEl21CaR5ZDHYZHNHAogLNpIKnWzKji6zDyf6+B8MfN2QYf9Axa2FlYKFiFtIJU2afwDFv6EkNJAmhS5OzsoiWju8Oade94973DdSLPTRVZ/R/+jxwJPC0+Kz5y8Gh0rDI8tJ2IqlqsvQC+NVx06U5waqrl3tqdUoVrbveGrF2Znrnq+0VZy4YfBR70Zq3bebgbvlSlsz1SiNp5a0PTHfETPi7QeRWo4vap2GKSy5GeqbfEI0SmWqUBbiLjBzUEYeC2HpKyxsIoRECz4UAmjGHmwk/K3BCFibh0pczEjNztX6KDI2hZXKa6wmd3hf5OztZwNO/2TDK15Fs8XjErBcbokj7RNV3QZ/pOv+/tlWY9um/Z5d3paVXUGNp/Wfv1X5XPu8b2reoBhcPVD3vS2MJU5sVlb1HGdF3KXv/23uF1bXpLH1DJ/SD/R3TFX1jh0H7pzxdVEtHKPKAzH/HcRcsT1RMqpiL78vV2XxUBbzCa4zPCZRxTwWUOd7D/AFZzg3hl1wzI2eqVGX64ZxV9hNP8ABRud4Q=</latexit> <latexit sha1_base64="UvOZW/4t5LzxWzoHwQCgzm7qQ=">ACrXichVE9SxRFD2OidlsEl21CaR5ZDHYZHNHAogLNpIKnWzKji6zDyf6+B8MfN2QYf9Axa2FlYKFiFtIJU2afwDFv6EkNJAmhS5OzsoiWju8Oade94973DdSLPTRVZ/R/+jxwJPC0+Kz5y8Gh0rDI8tJ2IqlqsvQC+NVx06U5waqrl3tqdUoVrbveGrF2Znrnq+0VZy4YfBR70Zq3bebgbvlSlsz1SiNp5a0PTHfETPi7QeRWo4vap2GKSy5GeqbfEI0SmWqUBbiLjBzUEYeC2HpKyxsIoRECz4UAmjGHmwk/K3BCFibh0pczEjNztX6KDI2hZXKa6wmd3hf5OztZwNO/2TDK15Fs8XjErBcbokj7RNV3QZ/pOv+/tlWY9um/Z5d3paVXUGNp/Wfv1X5XPu8b2reoBhcPVD3vS2MJU5sVlb1HGdF3KXv/23uF1bXpLH1DJ/SD/R3TFX1jh0H7pzxdVEtHKPKAzH/HcRcsT1RMqpiL78vV2XxUBbzCa4zPCZRxTwWUOd7D/AFZzg3hl1wzI2eqVGX64ZxV9hNP8ABRud4Q=</latexit> <latexit sha1_base64="UvOZW/4t5LzxWzoHwQCgzm7qQ=">ACrXichVE9SxRFD2OidlsEl21CaR5ZDHYZHNHAogLNpIKnWzKji6zDyf6+B8MfN2QYf9Axa2FlYKFiFtIJU2afwDFv6EkNJAmhS5OzsoiWju8Oade94973DdSLPTRVZ/R/+jxwJPC0+Kz5y8Gh0rDI8tJ2IqlqsvQC+NVx06U5waqrl3tqdUoVrbveGrF2Znrnq+0VZy4YfBR70Zq3bebgbvlSlsz1SiNp5a0PTHfETPi7QeRWo4vap2GKSy5GeqbfEI0SmWqUBbiLjBzUEYeC2HpKyxsIoRECz4UAmjGHmwk/K3BCFibh0pczEjNztX6KDI2hZXKa6wmd3hf5OztZwNO/2TDK15Fs8XjErBcbokj7RNV3QZ/pOv+/tlWY9um/Z5d3paVXUGNp/Wfv1X5XPu8b2reoBhcPVD3vS2MJU5sVlb1HGdF3KXv/23uF1bXpLH1DJ/SD/R3TFX1jh0H7pzxdVEtHKPKAzH/HcRcsT1RMqpiL78vV2XxUBbzCa4zPCZRxTwWUOd7D/AFZzg3hl1wzI2eqVGX64ZxV9hNP8ABRud4Q=</latexit> <latexit sha1_base64="UvOZW/4t5LzxWzoHwQCgzm7qQ=">ACrXichVE9SxRFD2OidlsEl21CaR5ZDHYZHNHAogLNpIKnWzKji6zDyf6+B8MfN2QYf9Axa2FlYKFiFtIJU2afwDFv6EkNJAmhS5OzsoiWju8Oade94973DdSLPTRVZ/R/+jxwJPC0+Kz5y8Gh0rDI8tJ2IqlqsvQC+NVx06U5waqrl3tqdUoVrbveGrF2Znrnq+0VZy4YfBR70Zq3bebgbvlSlsz1SiNp5a0PTHfETPi7QeRWo4vap2GKSy5GeqbfEI0SmWqUBbiLjBzUEYeC2HpKyxsIoRECz4UAmjGHmwk/K3BCFibh0pczEjNztX6KDI2hZXKa6wmd3hf5OztZwNO/2TDK15Fs8XjErBcbokj7RNV3QZ/pOv+/tlWY9um/Z5d3paVXUGNp/Wfv1X5XPu8b2reoBhcPVD3vS2MJU5sVlb1HGdF3KXv/23uF1bXpLH1DJ/SD/R3TFX1jh0H7pzxdVEtHKPKAzH/HcRcsT1RMqpiL78vV2XxUBbzCa4zPCZRxTwWUOd7D/AFZzg3hl1wzI2eqVGX64ZxV9hNP8ABRud4Q=</latexit> Tutorial: Heisenberg dimer (1) Sample of simple mode file (sample/dla/01_spindimer/std.toml) • S=1/2 antiferromagnetic Heisenberg dimer • − J H = − J S 1 · S 2

Specify simulation parameters seed = 31415 # seed of RNG [parameter] section 
 [parameter] beta = 100 # inverse temperature nset = 5 # set of Monte Carlo sweeps npre = 10 # MCSteps to estimate hyperparameter ntherm = 10 # MCSweeps for thermalization nmcs = 100 # MCSweeps for measurement Tutorial: Heisenberg dimer (2) Sample of simple mode file (sample/dla/01_spindimer/std.toml) • ����� ����� �������� ���� ������ ���� ������ ���� ������ ����

Tutorial: Heisenberg dimer (3) Run the dla_pre • param.in, lattice.xml, 
 $ dla_pre std.toml algorithm.xml output Run the dla • Random number parallelization $ dla param.in $ mpiexec -np 4 dla param.in Standard output (log) >>> The program is being run with MPI mode.( N_PROC = 1 ) +++++++++ input data +++++++++ RUNTYPE = 0 … +++++++++ input data +++++++++ Determining hyperparameter NCYC : 51 Start main calculation. 1 / 5 done. [Elapsed: 0.014687 sec. ETR: 0.058748 sec.] 2 / 5 done. [Elapsed: 0.195772 sec. ETR: 0.293658 sec.] 3 / 5 done. [Elapsed: 0.333755 sec. ETR: 0.222503 sec.] 4 / 5 done. [Elapsed: 0.471529 sec. ETR: 0.117882 sec.] ETR: Estimated Time Remaining 5 / 5 done. [Elapsed: 0.616055 sec. ETR: 0 sec.]

-3.75 (exact) P <name> = <value> Ground state energy 
 Other information. I <text> = <value> <error>： the statistical error of <mean>. <mean> ：the expected value Results of observables. R <name> = <mean> <error> Parameters read from the input files. Tutorial: Heisenberg dimer (4) Output file: sample.log I N_PROC = 1 P D = 1 P L = 2 … P SIMULATIONTIME = 0.000000 R sign = 1.00000000e+00 0.00000000e+00 R anv = 2.49300000e+01 2.96344394e-01 R ene = -3.74300000e-01 2.96344394e-03 R spe = 8.16400000e-02 1.47017825e+00 R som = 8.16400000e+00 1.47017825e+02 R len = 4.00127485e+00 2.92088648e-02 R xmx = 1.00031871e-02 7.30221620e-05 … R time = 1.37698000e-06 1.20496230e-07 I [the maximum number of segments] = 165 I [the maximum number of vertices] = 84 I [the maximum number of reg. vertex info.] = 1 Main Results are written in [R] section. Analyze the result • $ grep ene sample.log 
 R ene = -3.74300000e-01 2.96344394e-03

Tutorial: Heisenberg dimer (5) Main Results 1/3

Tutorial: Heisenberg dimer (6) Main Results 2/3

Tutorial: Heisenberg dimer (7) Main Results 3/3