March 2020 - Comp-phys-alps-users

extracting outputs
by Negin Mohseni 24 Mar '20

24 Mar '20

Dear All, My question is related to the calculated results of the below code. This entails MPS implantation to reach energy gap of a Heisenberg chain in the presence of a magnetic field. We wrote: %%%%%%%%%%%%%%% Start %%%%%%%%%%%%%%%%%%%%%%% import pyalps import numpy as np import matplotlib.pyplot as plt import pyalps.plot parms = [] for hz in [0.0, 0.2, 0.4, 0.6, 0.8, 1.0]: parms.append( { 'LATTICE' : "open chain lattice", 'MODEL' : "spin", 'CONSERVED_QUANTUMNUMBERS' : 'N,Sz', 'Sz_total' : 0, 'J' : 1, 'h' : hz, 'SWEEPS' : 4, 'L' : 32, 'MAXSTATES' : 100, 'NUMBER_EIGENVALUES' : 2 } ) input_file = pyalps.writeInputFiles('parm_spin_one_half_gap',parms) res = pyalps.runApplication('mps_optim',input_file,writexml=True) data = pyalps.loadEigenstateMeasurements(pyalps.getResultFiles(prefix='parm_spin_one_half_gap')) energies = np.empty(0) for s in data[0]: if s.props['observable'] == 'Energy': energies = s.y else: print(s.props['observable'], ':', s.y[0]) energies.sort() print('Energies:', end=' ') for e in energies: print(e, end=' ') print('\nGap:', abs(energies[1]-energies[0])) %%%%%%%%%%%%%%% END %%%%%%%%%%%%%%%%%%%%%%% Although, in the file "parm_spin_one_half_gap.in.xml" there are expensive information such as energies for any "hz", we need to collect all results in a *.dat or *.txt file. For instance, a "file.txt" as: hz gap 0.0 0.11 0.2 0.12 0.4 0.14 0.6 0.15 How can we generate such a file?!

1 0

ALPS installation via Conda and SVN Trunk package
by Manpreet Singh 10 Mar '20

10 Mar '20

Dear all, I have installed ALPS via Conda on Linux Ubuntu 16.04 LTS and it seems to work fine. However, there is an ALPS Trunk version available via SVN at alps.comp-phys.org/svn/alps1/trunk/alps. I would like to know which one of these versions has the latest and updated libraries. Also, if ALPS Trunk is the latest version, how do I update the existing libraries which were installed via Conda? Any help would be appreciated. Thanks in advance! Best regards- Manpreet

2 1

drawing an arbitrary lattice
by Negin Mohseni 01 Mar '20

01 Mar '20

So sorry there was a mistake in the last email. I has corrected recently and asked it again. Dear all Firstly, thanks to Matthias for his previous answers. We are looking for a way to map our lattice. With guidance presented in the web site we wrote our lattice. Although, the unit cell of what we wanted has been brought by default in the lattice collections of the ALPS, we are trying to apply and see what is in our mind explicitly. We wrote: <LATTICES> <LATTICE name="stripm" dimension="2"> <PARAMETER name="a" default="1.73205"/> <PARAMETER name="b" default="1.73205"/> <PARAMETER name="phi" default="Pi/3"/> <BASIS> <VECTOR> a 0 </VECTOR> <VECTOR> b*cos(phi) b*sin(phi) </VECTOR> </BASIS> </LATTICE> <FINITELATTICE name="honeyladder" dimension="2"> <LATTICE ref="stripm"/> <PARAMETER name="L" default="10"/> <PARAMETER name="W" default="2" /> <EXTENT dimension="1" size="L" /> <EXTENT dimension="2" size="W" /> <BOUNDARY dimension="1" type="periodic" /> <BOUNDARY dimension="2" type="open" /> </FINITELATTICE> <UNITCELL name="honeyl" dimension="2" vertices="2"> <VERTEX id="1" ><COORDINATE> cos(phi/2) cos(phi/2) </COORDINATE></VERTEX> <VERTEX id="2"><COORDINATE> 2*cos(phi/2) 2* cos(phi/2) </COORDINATE></VERTEX> <EDGE><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="1" offset="0 1"/></EDGE> </UNITCELL> <LATTICEGRAPH> <FINITELATTICE ref="honeyladder"/> <UNITCELL ref="honeyl"/> </LATTICEGRAPH> </LATTICES> *Question*: How to see the output of the .xml file? Thanks in advance

2 1

draw an arbitrary lattice
by Negin Mohseni 01 Mar '20

01 Mar '20

Dear all Firstly, thanks to Matthias for his previous answers. We are looking for a way to map our lattice. With guidance presented in the web site we wrote our lattice. Although, the unit cell of what we wanted has been brought by default in the lattice collections of the ALPS, we are trying to apply and see what is in our mind explicitly. We wrote: <LATTICES> <LATTICE name="stripm" dimension="2"> <PARAMETER name="a" default="1.73205"/> <PARAMETER name="b" default="1.73205"/> <PARAMETER name="phi" default="Pi/3"/> <BASIS> <VECTOR> a 0 </VECTOR> <VECTOR> b*sin(phi) b*cos(phi) </VECTOR> </BASIS> </LATTICE> <FINITELATTICE name="honeyladder" dimension="2"> <LATTICE ref="stripm"/> <PARAMETER name="L" default="10"/> <PARAMETER name="W" default="2" /> <EXTENT dimension="1" size="L" /> <EXTENT dimension="2" size="W" /> <BOUNDARY dimension="1" type="periodic" /> <BOUNDARY dimension="2" type="open" /> </FINITELATTICE> <UNITCELL name="honeyl" dimension="2" vertices="2"> <VERTEX id="1" ><COORDINATE> cos(phi) cos(phi) </COORDINATE></VERTEX> <VERTEX id="2"><COORDINATE> 2*cos(phi) 2* cos(phi) </COORDINATE></VERTEX> <EDGE><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="1" offset="1 0"/></EDGE> <EDGE><SOURCE vertex="2"/><TARGET vertex="1" offset="0 1"/></EDGE> </UNITCELL> <LATTICEGRAPH> <FINITELATTICE ref="honeyladder"/> <UNITCELL ref="honeyl"/> </LATTICEGRAPH> </LATTICES> *Question*: How to see the output of the .xml file? Thanks in advance

1 0