def run_dwave(): # example usage as follows : print('2D Migdal with superconductivity') basedir = '../test_dwave2/' if not os.path.exists(basedir): os.makedirs(basedir) params = {} params['nw'] = 128 params['nk'] = 16 params['t'] = 1.0 params['tp'] = -0.3 params['omega'] = 0.17 params['dens'] = 0.8 params['renormalized'] = True params['sc'] = True params['band'] = band_square_lattice params['beta'] = 40.0 params['fixed_mu'] = -1.11 lamb = 1 / 6 W = 8.0 params['g0'] = mylamb2g0(lamb, params['omega'], W) print('g0 is ', params['g0']) ks = np.arange(-np.pi, np.pi, 2 * np.pi / params['nk']) params['gk2'] = 2 * (0.5 * (np.cos(ks[:, None]) - np.cos(ks[None, :])))**2 migdal = Migdal(params, basedir) sc_iter = 2000 savedir, mu, G, D, S, GG = migdal.selfconsistency(sc_iter, S0=None, PI0=None, frac=0.2) print('------------------------------------------') print('simulation took', time.time() - time0, 's')
params['band'] = band_square_lattice params['dim'] = 2 #ibeta = int(sys.argv[1]) #betas = np.linspace(1.0, 50.0, 50) #params['beta'] = betas[ibeta] params['beta'] = float(sys.argv[1]) renorm = True if sys.argv[2] == 'True' else False params['renormalized'] = renorm #lamb = 0.2 #lamb = 0.3 # try this one as well...... lamb = float(sys.argv[3]) W = 8.0 params['g0'] = mylamb2g0(lamb, params['omega'], W) params['gq2'] = 1 params['dw'] = 0.001 params['wmin'] = -4.1 params['wmax'] = +4.1 params['idelta'] = 0.010j basedir = '/home/groups/simes/bln/data/elph/production/' if not os.path.exists(basedir): os.makedirs(basedir) migdal = RealAxisMigdal(params, basedir) sc_iter, S0, PI0 = 2000, None, None migdal.selfconsistency(sc_iter, S0=S0,
def run_dwave(l, beta=16.0, S0=None, PI0=None): time0 = time.time() basedir = '../test_dwave/' if not os.path.exists(basedir): os.makedirs(basedir) print('2D Renormalized Migdal') params = {} params['nw'] = 128 params['nk'] = 16 params['t'] = 1 params['tp'] = -0.3 params['omega'] = 0.17 params['dens'] = 0.8 params['renormalized'] = True params['sc'] = False params['band'] = band_square_lattice params['beta'] = beta params['dim'] = 2 params['fixed_mu'] = -1.11 lamb = l W = 8.0 params['g0'] = mylamb2g0(lamb, params['omega'], W) print('g0 is ', params['g0']) ks = np.arange(-np.pi, np.pi, 2 * np.pi / params['nk']) params['gk2'] = (np.cos(ks[:, None]) - np.cos(ks[None, :]))**2 migdal = Migdal(params, basedir) sc_iter = 1000 savedir, mu, G, D, S, GG = migdal.selfconsistency(sc_iter, S0=S0, PI0=PI0, frac=0.2) if G is None: return False, None, None PI = params['g0']**2 * GG save(savedir + 'S.npy', S) save(savedir + 'PI.npy', PI) save(savedir + 'G.npy', G) save(savedir + 'D.npy', D) sc_iter = 400 Xsc, Xcdw = migdal.susceptibilities(savedir, sc_iter, G, D, GG, frac=1.0) save(savedir + 'Xsc.npy', [Xsc]) save(savedir + 'Xcdw.npy', [Xcdw]) if Xcdw is None: return False, None, None print('------------------------------------------') print('simulation took', time.time() - time0, 's') return True, G, S
renorm = int(sys.argv[1]) params = {} params['nw'] = 128 params['nk'] = 120 params['t'] = 1.0 params['tp'] = -0.3 params['dens'] = 0.8 #params['fixed_mu'] = -1.11 params['omega'] = 0.17 params['renormalized'] = True if renorm == 1 else False params['sc'] = False params['band'] = band_square_lattice params['beta'] = 16.0 params['dim'] = 2 params['g0'] = mylamb2g0(lamb=1 / 6, omega=params['omega'], W=8.0) params['Q'] = None ''' params['dw'] = 0.001 params['wmin'] = -4.2 params['wmax'] = +4.2 params['idelta'] = 0.002j migdal = Migdal(params, basedir) ''' params['dw'] = 0.01 params['wmin'] = -7.2 params['wmax'] = +7.2 params['idelta'] = 0.05j
def run(): # example usage as follows : print('2D Migdal with superconductivity') params = {} params['nw'] = 512 params['nk'] = 40 params['t'] = 1.0 params['tp'] = -0.3 params['omega'] = 0.1 params['dens'] = 0.8 params['renormalized'] = True params['sc'] = True params['band'] = band_square_lattice #params['beta'] = 1.0 params['Q'] = None #params['g0'] = 0.525 #params['dim'] = 2 #basedir = '/home/groups/simes/bln/data/elph/imagaxis/example/' basedir = '/scratch/users/bln/elph/data/sc2d_200223/' if not os.path.exists(basedir): os.makedirs(basedir) #lamb = 0.3 #lamb = float(sys.argv[1]) lamb = 0.2 print('lamb = ', lamb) W = 8.0 params['g0'] = mylamb2g0(lamb, params['omega'], W) print('g0 is ', params['g0']) params['beta'] = 5.0 path = sys.argv[1] if path: S0 = load(path + 'S.npy') PI0 = load(path + 'PI.npy') mu0 = load(path + 'mu.npy')[0] else: S0, PI0, mu0 = None, None, None migdal = Migdal(params, basedir) sc_iter = 2000 savedir, mu, G, D, S, GG = migdal.selfconsistency(sc_iter, S0=S0, PI0=PI0, mu0=mu0, frac=0.1) if not params['renormalized']: PI = None else: PI = params['g0']**2 * GG save(savedir + 'S.npy', S) save(savedir + 'PI.npy', PI) save(savedir + 'G.npy', G) save(savedir + 'D.npy', D) print('------------------------------------------') print('simulation took', time.time() - time0, 's')