karthik/Calculations
1
0
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity

Script to plot ansatzes for different phases.

Browse Source
This commit is contained in:
Karthik 2025-03-20 15:32:47 +01:00
parent f6394cd6a0
commit 22ae799879
1 changed files with 90 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

90
Estimations/CosineModulatedAnsatz.m Normal file
View File

@ -0,0 +1,90 @@
%% STRIPES
% 2-D
% Parameters
c = 1; % Fourier coeffecient
k = 2; % wavenumber
n = 2; % order
% Range for x and y
x = linspace(-2, 2, 500);
y = linspace(-2, 2, 500);
% Create a meshgrid for 2D
[~, Y] = meshgrid(x, y);
% Define the 2D function
f_xy = (1 + (c * cos(n * k * Y))) / (1 + (0.5 * c^2));
% Plot the 2D image using imagesc
figure(1);
imagesc(x, y, f_xy);
axis xy; % Make sure the y-axis is oriented correctly
colorbar;
xlabel('$x l_o$', 'Interpreter', 'latex', 'FontSize', 14)
ylabel('$y l_o$', 'Interpreter', 'latex', 'FontSize', 14)
title('Stripe Lattice ansatz for $|\Psi(x,y)|^2$', 'Interpreter', 'latex', 'FontSize', 16);
xlabel('x');
ylabel('y');
colormap parula;
%% TRIANGULAR LATTICE
% 2-D
% Parameters
c1 = 0.2;
c2 = 0.2;
k = 3;
n = 1;
% Range for x and y
x = linspace(-2, 2, 500);
y = linspace(-2, 2, 500);
% Create a meshgrid for 2D
[X, Y] = meshgrid(x, y);
% Define the 2D function for a triangular lattice
f_xy = 1 + (c1 * cos(n * k * (2/sqrt(3)) * Y)) + (2 * c2 * cos(n * k * (1/sqrt(3)) * Y) .* cos(n * k * X));
% Plot the 2D image using imagesc
figure(2);
clf
imagesc(x, y, f_xy);
axis xy; % Make sure the y-axis is oriented correctly
colorbar;
xlabel('$x l_o$', 'Interpreter', 'latex', 'FontSize', 14)
ylabel('$y l_o$', 'Interpreter', 'latex', 'FontSize', 14)
title('Triangular Lattice ansatz for $|\Psi(x,y)|^2$', 'Interpreter', 'latex', 'FontSize', 16);
xlabel('x');
ylabel('y');
colormap parula;
%% HONEYCOMB LATTICE
% 2-D
% Parameters
c1 = 0.2;
c2 = 0.2;
k = 3;
n = 1;
% Range for x and y
x = linspace(-2, 2, 500);
y = linspace(-2, 2, 500);
% Create a meshgrid for 2D
[X, Y] = meshgrid(x, y);
% Define the 2D function for a honeycomb lattice
f_xy = 1 - (c1 * cos(n * k * (2/sqrt(3)) * X)) - (2 * c2 * cos(n * k * (1/sqrt(3)) * X) .* cos(n * k * Y));
% Plot the 2D image using imagesc
figure(3);
clf
imagesc(x, y, f_xy);
axis xy; % Make sure the y-axis is oriented correctly
colorbar;
xlabel('$x l_o$', 'Interpreter', 'latex', 'FontSize', 14)
ylabel('$y l_o$', 'Interpreter', 'latex', 'FontSize', 14)
title('Honeycomb Lattice ansatz for $|\Psi(x,y)|^2$', 'Interpreter', 'latex', 'FontSize', 16);
xlabel('x');
ylabel('y');
colormap parula;