trap_annealing_steps = [0.5] * 10
Vres = 0.50
Vtrap = 0.00
Vrg = -0.50
Vcg = 0.00
Vtg = -0.50
simulation_name = "Combined_Model"
save_path = r"/Users/gkoolstra/Desktop/Electron optimization/Realistic potential/Combined Model"
sub_dir = time.strftime("%y%m%d_%H%M%S_{}".format(simulation_name))
save = False
# Evaluate all files in the following range.
xeval = np.linspace(-5.0, 12.5, 501)
yeval = anneal.construct_symmetric_y(-5.0, 151)
master_path = r"/Volumes/slab/Gerwin/Electron on helium/Maxwell/M018 Yggdrasil/Greater Trap Area/V1"
#master_path = r"/Users/gkoolstra/Desktop/Electron optimization/Realistic potential/Potentials/V6 - Greater Trap Area"
x_eval, y_eval, output = anneal.load_data(master_path, xeval=xeval, yeval=yeval, mirror_y=True,
extend_resonator=True, do_plot=False)
# Note: x_eval and y_eval are 2D arrays that contain the x and y coordinates at which the potentials are evaluated
# Use resonator potential to solve for equilibrium position of resonator electrons
U_resonator = -Vres * output[0]['V'][-1, :]
# Solve resonator potential
RS = anneal.ResonatorSolver(y_eval[:, -1], U_resonator, efield_data=None,
box_length=resonator_box_length, spline_order_x=3, smoothing=smoothing)
res_initial_condition = anneal.setup_initial_condition(N_resonator_electrons,
#save_path = r"/Volumes/slab/Gerwin/Electron on helium/Electron optimization/Realistic potential/Resonator"
sub_dir = time.strftime("%y%m%d_%H%M%S_{}".format(simulation_name))
save = True
create_movie = True
# Load the data from the dsp file:
#path = r'/Volumes/slab/Gerwin/Electron on helium/Maxwell/M018 Yggdrasil/Greater Trap Area/V1big/DCBiasPotential.dsp'
path = r'S:\Gerwin\Electron on helium\Maxwell\M018 Yggdrasil\All simulation data M018V3\DCBiasPotential.dsp'
elements, nodes, elem_solution, bounding_box = import_data.load_dsp(path)
xdata, ydata, Udata = interpolate_slow.prepare_for_interpolation(
elements, nodes, elem_solution)
x0 = -2.0 # Starting point for y
k = 251 # This defines the sampling
xeval = anneal.construct_symmetric_y(x0, k)
fig0 = plt.figure(figsize=(5., 3.))
xinterp, yinterp, Uinterp = interpolate_slow.evaluate_on_grid(
xdata,
ydata,
Udata,
xeval=xeval,
yeval=np.linspace(-1, 1, 151),
clim=(0.00, 1.00),
plot_axes='xy',
linestyle='None',
cmap=plt.cm.viridis,
plot_data=True,
**common.plot_opt("darkorange", msize=6))
trap_annealing_steps = [1.0] * 10
Vres = 0.30
Vtrap = 0.00
Vrg = 0.00
Vcg = 0.00
Vtg = -0.30
simulation_name = "Combined_Model"
save_path = r"/Users/gkoolstra/Desktop/Electron optimization/Realistic potential/Combined Model"
sub_dir = time.strftime("%y%m%d_%H%M%S_{}".format(simulation_name))
save = False
# Evaluate all files in the following range.
xeval = np.linspace(-4, 9.5, 501)
yeval = anneal.construct_symmetric_y(-3.0, 121)
master_path = r"/Volumes/slab/Gerwin/Electron on helium/Maxwell/M018 Yggdrasil/Greater Trap Area/V1"
#master_path = r"/Users/gkoolstra/Desktop/Electron optimization/Realistic potential/Potentials/V6 - Greater Trap Area"
x_eval, y_eval, output = anneal.load_data(master_path, xeval=xeval, yeval=yeval, mirror_y=True,
extend_resonator=True, do_plot=False)
# Note: x_eval and y_eval are 2D arrays that contain the x and y coordinates at which the potentials are evaluated
# Use resonator potential to solve for equilibrium position of resonator electrons
U_resonator = -Vres * output[0]['V'][-1, :]
# Solve resonator potential
RS = anneal.ResonatorSolver(y_eval[:, -1], U_resonator, efield_data=None,
box_length=resonator_box_length, spline_order_x=3, smoothing=0)
res_initial_condition = anneal.setup_initial_condition(N_resonator_electrons,
x0=(inserted_res_length*1E-6 + box_length)/2., y0=0.0E-6, dx=0.40E-6)
x_trap_init, y_trap_init = anneal.r2xy(electron_initial_positions)
if platform.system() == 'Windows':
save_path = r"S:\Gerwin\Electron on helium\Electron optimization\Realistic potential\Single electron loading"
elif platform.system() == 'Linux':
save_path = r"/mnt/slab/Gerwin/Electron on helium/Electron optimization/Realistic potential/Single electron loading"
else:
save_path = r"/Users/gkoolstra/Desktop/Single electron loading/Data"
sub_dir = time.strftime("%y%m%d_%H%M%S_{}".format(simulation_name))
save = True
# Evaluate all files in the following range.
xeval = np.linspace(-4.0, box_length*1E6, 751)
yeval = anneal.construct_symmetric_y(-4.0, 201)
dx = np.diff(xeval)[0]*1E-6
dy = np.diff(yeval)[0]*1E-6
if platform.system() == 'Windows':
master_path = r"S:\Gerwin\Electron on helium\Maxwell\M018 Yggdrasil\M018V6\V6.2"
elif platform.system() == 'Linux':
master_path = r"/mnt/slab/Gerwin/Electron on helium/Maxwell/M018 Yggdrasil/M018V6/V6.2"
else:
master_path = r"/Users/gkoolstra/Desktop/Single electron loading/Potentials/M018V6/V6.2"
x_eval, y_eval, output = anneal.load_data(master_path, xeval=xeval, yeval=yeval, mirror_y=True,
extend_resonator=False, insert_resonator=True, do_plot=inspect_potentials,
inserted_res_length=inserted_res_length, smoothen_xy=(0.40E-6, 2*dy))