Optical Detection and Manipulation of Pseudospin Orders in Wigner Crystals

OPTICAL DETECTION AND MANIPULATION OF PSEUDOSPIN ORDERS IN WIGNER CRYSTALSFIGURE DATA This folder contains 16 plain-text data files for Figures 1-7 of thearticle. These saved results can be plottedwithout running new simulations. File names identify the correspondingfigure and panel. File format and units The files contain space-separated columns. Lines beginningwith # are column headers. Frequencies are given in THz, as on the figure axes. Twist angles are indegrees. Carrier densities are quoted in units of 10^11 cm^-2, so a storedvalue of 1 corresponds to 10^11 cm^-2. Conductivity units are specified foreach file: e^2/h, arbitrary units (arb. u.), or dimensionless normalizedvalues. The reference momentum is q_M = 2*pi/(sqrt(3)*a). Momentum ratios, E_G/J and w/a aredimensionless. J denotes the exchange energy. Data files and columns 1. Fig1b.txt Phonon dispersion along the Gamma-K-M-Gamma path; 241 rows. Column 1: cumulative distance along the momentum path, multiplied by a (dimensionless). Column 2: q_x*a (dimensionless). Column 3: q_y*a (dimensionless). Column 4: transverse phonon frequency omega_T (THz). Column 5: longitudinal phonon frequency omega_L (THz). Plot column 1 against columns 4 and 5. 2. Fig1c.txt Optical conductivity for four pseudospin configurations; 200 rows. Column 1: frequency omega (THz). Column 2: Re sigma for the spin liquid (arb. u.). Column 3: Re sigma for the 120-degree antiferromagnet (arb. u.). Column 4: Re sigma for the stripe antiferromagnet (arb. u.). Column 5: Re sigma for the ferromagnet (arb. u.; zero in this dataset). Plot column 1 against columns 2-5. 3. Fig2b.txt Two-dimensional energy landscape; 100002 momentum-space samples. Column 1: q_x/q_M (dimensionless). Column 2: q_y/q_M (dimensionless). Column 3: energy E_G(q)/J (dimensionless). Columns 1 and 2 are the coordinates; column 3 gives the color value. 4. Fig2c.txt One-dimensional cuts of the energy landscape; 2001 rows. Column 1: q_y/q_M (dimensionless). Column 2: E_G(q)/J at pump field E = 0 (dimensionless). Column 3: E_G(q)/J at E = 0.65*10^3 V/cm (dimensionless). Column 4: E_G(q)/J at E = 1*10^3 V/cm (dimensionless). Plot column 1 against columns 2-4. 5. Fig2d.txt Phase diagram on a 501-by-501 frequency-field grid; 251001 rows. Column 1: driving frequency omega (THz). Column 2: electric field E in units of 10^3 V/cm. Column 3: ordering wave vector q_c/q_M (dimensionless). Columns 1 and 2 are the horizontal and vertical coordinates; column 3 gives the color value. A zero ordering wave vector denotes the ferromagnetic region; nonzero values denote stripe antiferromagnetic order. 6. Fig3a.txt, Fig3b.txt, Fig3c.txt and Fig3d.txt Twist-angle dependence of optical conductivity; 1000 rows per file. The twist angles are 1.0, 1.2, 1.4 and 1.6 degrees for panels a, b, c and d, respectively. Each file has the following columns: Column 1: frequency omega (THz). Column 2: Re sigma for the 120-degree antiferromagnet (e^2/h). Column 3: Re sigma for the stripe antiferromagnet (e^2/h). Plot column 1 against columns 2 and 3. 7. Fig4a.txt Carrier-density dependence of normalized optical conductivity for stripe antiferromagnetic order; 1000 rows. Column 1: frequency omega (THz). Columns 2-9: dimensionless, normalized conductivity profiles at the following densities(units of 10^11 cm^-2): 0.2000000000, 0.2628571429, 0.3257142857, 0.3885714286, 0.4514285714, 0.5142857143, 0.5771428571, 0.6400000000. The stored profiles include the displayed vertical offsets of 0, 0.7, 1.4, 2.1, 2.8, 3.5, 4.2 and 4.9, respectively. Plot column 1 against columns 2-9. 8. Fig4b.txt Pseudospin-order-induced peak frequency versus density; 8 rows. Column 1: carrier density n (10^11 cm^-2). Column 2: peak frequency omega_p (THz). Column 3: power-law fitted peak frequency (THz). Plot columns 1 and 2 for the data points, and columns 1 and 3 for the fitted line. The inset uses logarithmic axes. Column 3 contains fitted curve values at the densities in column 1. 9. Fig4c.txt Pseudospin-order-induced peak spectral weight versus density; 8 rows. Column 1: carrier density n (10^11 cm^-2). Column 2: spectral weight W_p in units of pi*n*e^2/(2*m). Column 3: power-law fitted W_p in the same units as column 2. Here m is the effective electron mass. Columns 2 and 3 are the dimensionless ratios W_p/[pi*n*e^2/(2*m)] and its fitted value. Plot columns 1 and 2 for the data points, and columns 1 and 3 for the fitted line. The inset uses logarithmic axes. 10. Fig5.txt Effective temperature and relative displacement width after pumping, as functions of twist angle; 5 rows. Column 1: twist angle theta (degrees). Column 2: effective temperature T_eff (K). Column 3: relative electron-displacement width w/a (dimensionless). Plot column 1 against column 2 on the left vertical axis, and against column 3 on the right vertical axis. 11. Fig6.txt Exchange energy versus carrier density; 10 rows. Column 1: carrier density n (10^11 cm^-2). Column 2: exchange energy J (ueV). Plot column 1 against column 2. 12. Fig7a.txt and Fig7b.txt Polarization dependence of optical conductivity; 200 rows per file. Fig7a.txt is for the 120-degree antiferromagnet; Fig7b.txt is for the stripe antiferromagnet. Each file has the following columns: Column 1: frequency omega (THz). Column 2: Re sigma for an x-polarized probe field E_x (arb. u.). Column 3: Re sigma for a y-polarized probe field E_y (arb. u.). Plot column 1 against columns 2 and 3.