autorelaxed_twistronic_graphenes/Layer.m at master · stcarr/autorelaxed_twistronic_graphenes · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
classdef Layer
    % Basic datastructure to store geometric info for a 2D layer.

    properties
        layer_index     % index in the heterostructure
        A               % global lattice
        G               % global recip lattice
        A_local         % local lattice (theta = 0)
        G_local         % local recip lattice (theta = 0)
        theta           % twist/orientation angle
        orbPos
    end

    methods
        function obj = Layer(index,theta)
            %Layer Construct an instance of this class
            %   Detailed explanation goes here
            obj.layer_index = index;
            obj.theta = theta;

            alpha = 1.42*sqrt(3);

            obj.A_local = alpha*[   1       1/2
                                    0 sqrt(3)/2];
            obj.G_local = getRecip(obj.A_local);

            R_base = [cos(theta) -sin(theta);
                       sin(theta) cos(theta)];

            obj.A = R_base*obj.A_local;
            obj.G = getRecip(obj.A);

            obj.orbPos(1,:) = [0 0]; % A orbital
            obj.orbPos(2,:) = 1/3*(obj.A(:,1) + obj.A(:,2)); % B orbital

        end
    end
end