Code: Select all
spinSite = SpinSite(S=0.5, conserve='Sz')
lattice = Chain(L=N, site=spinSite)
model_params = {...}
model = SpinModel(model_params)
[...]
for i in range(8):
psi = MPS.from_product_state(sites, State, "finite")
mixer_params = {}
dmrg_params = {"trunc_params": {"chi_max": chi_max, "svd_min": svd_min}, "mixer": mixer,"mixer_params":mixer_params,
'combine': True,"orthogonal_to": excited_states, 'verbose':False}
info = dmrg.run(psi, model, dmrg_params)
excited_states.append(psi)
# Energy
energy = info['E']
bond_energies = model.bond_energies(psi)
[...]
Code: Select all
energy: 0.5
bond_energies: [0.25 0.25]
sum(bond_energies): 0.5
energy: -1.0
bond_energies: [-0.5 -0.5]
sum(bond_energies): -1.0000000000000004
energy: -1.3877787807814466e-17
bond_energies: [ 2.91433544e-16 -3.19189120e-16]
sum(bond_energies): -2.775557561562896e-17
energy: -1.1148336360798884e-16 # <---
bond_energies: [ 0.24669963 -0.44867985] # <---
sum(bond_energies): -0.20198021966036506 # <---
energy: -1.0
bond_energies: [-0.5 -0.5]
sum(bond_energies): -0.9999999999999998
energy: -2.4903616914531226e-18 # <---
bond_energies: [ 0.15147753 -0.70883147] # <---
sum(bond_energies): -0.5573539334676407 # <---
energy: -7.320883228133898e-18 # <---
bond_energies: [-0.08633658 0.07732684] # <---
sum(bond_energies): -0.009009746969017227 # <---
energy: 0.5
bond_energies: [0.25 0.25]
sum(bond_energies): 0.5