READ ME

This text describes the data presented in the paper: 
"Non-equilibrium atomic condensates and mixtures: collective modes, condensate growth and thermalisation"
J. Phys. B: At. Mol. Opt. Phys. 49, 214003 (2016), arxiv:1607.06939

========================
Introductory information
========================
Files included in the data deposit (within the compressed file raw_data.zip): 

1) directory "dipole": data generated from dynamical simulations of the ZNG model, showing the time-evolution of the centers-of-mass of the atoms. (Fig. 1 and Fig. 2).

Fig. 1: com_1compo.dat

Fig. 2: com_2compo.dat

2) directory "thermalization_1compo": data generated from dynamical simulations of the ZNG model for a single-component Bose gas of Rb87, showing:

Fig 3: init_tpvel.dat, final_tpvel.dat, the Fortran binary files that store the velocities of test particles, at the beginning and end of simulations. Used together with histo.py to create the figure.

Fig. 4 & 6: fraction_Ec**.dat, showing the condensate fraction as a function of time, where ** stands for the cutoff energy.

Fig. 5, 6, 8: ave_ke_Ec**.dat, showing the average kinetic energy (or T_eff) or the thermal atoms as a function time, where ** stands for the cutoff energy.

Fig. 6: width_Ec2400.dat, showing the width of the cloud as a function of time, at cutoff energy=2400nK.

Fig. 7: conz_Ec2400_final.dat and thz_Ec2400_final.dat, showing the condensate and thermal cloud density profles along the z-axis at the end of simulation, for cutoff energy=2400nK

conz_Ec2400_eq.dat and thz_Ec2400_eq.dat, showing the equilibrium condensate and thermal cloud density profles along the z-axis to match with those at the end of simulation, for cutoff energy=2400nK

3) directory "thermalization_2compo": data generated from dynamical simulations of the ZNG model for a two-component Bose-Bose mixture of Rb87-K41, showing:

Fig. 8 and 10: *_fraction*.dat, showing the condensate fraction of Rb87 and K41 as a function of time.
        scaled_ave_ke.dat, showing the average kinetic energy (or T_eff) of the thermal atoms as a function time.

Fig. 9: conz_*.dat and thz_*.dat - condensate and thermal cloud density profiles along the z-axis, at the end of simulation

Key words used to describe the data:
Centers-of-mass, widths of atomic clouds, condensate numbers, density profiles

========================== 
Methodological information
==========================
A brief method description – what the data is, how and why it was collected or created, and how it was processed:
The purpose of this data set is to show that the ZNG model can capture essential physics, and it can be used to investigate the interesting physics of thermalization. We simulate the dipole oscillation, to show that the center-of-mass of the whole cloud oscillates undamped in a harmonic trap, despite of the interactions and temperature. We also show the interesting monopole oscillation emerges from a quenched cooling process.


Instruments, hardware and software used: 
Data is generated with ZNG model using code written by Kean Loon Lee (Newcastle University, email: leekeanloon@gmail.com). The code is written in Fortran, parallelised with OpenMP. Data is simulated on SAgE cluster (topsy), and cross-checked with simulations on N8 computing cluster. 

Details on the numerical implementation of the ZNG model can be found in the followings:
1) B. Jackson and E. Zaremba, Phys. Rev. A 66, 033606 (2002)

2) A. Griffin, T. Nikuni and E. Zaremba, "Bose-condensed Gases at Finite Temperatures", Cambridge University Press (2009)

3) A. J. Allen, PhD Thesis (Newcastle University, 2012), "Non-Equilibrium and Finite Temperature Trapped Bose Gases: Interactions and Decay of Macroscopic Excitations"

4) K.L. Lee and N.P. Proukakis, "Non-equilibrium atomic condensates and mixtures: collective modes, condensate growth and thermalisation", J. Phys. B: At. Mol. Opt. Phys. 49, 214003 (2016)

Date(s) of data collection: 
September 2015 -- June 2016


Geographic coverage of data:
na


Data validation (how was the data checked, proofed and cleaned):
Data from the two clusters have been cross-checked. For equilibrium profiles, we check the density profiles with those published in journal articles by other researchers. For ZNG simulations, we verify the correctness of our code by looking at other dynamical situations, e.g. Kohn mode, monopole oscillations, etc. that have been investigated and published by Jackson and Zaremba.



Overview of secondary data, if used:
na


=========================
Data-specific information
=========================
Definitions of names, labels, acronyms or specialist terminology uses for variables, records and their values: 

1) "dipole" directory:

Center-of-mass (COM) of a single-component Rb87 gas (com_1compo.dat)
column 1 - omega*t
column 2 to 4 - COM of condensates along x,y,z directions, in units of ell_Rb
column 5 to 7 - COM of thermal cloud along x,y,z directions, in units of ell_Rb
column 8 to 10 - COM of whole cloud (condensate + thermal cloud) along x,y,z directions, in units of ell_Rb

Center-of-mass (COM) of a two-component Rb87-K41 Bose-Bose mixture (com_2compo.dat)
column 1 - omega*t
column 2 to 4 - COM of Rb87 condensates along x,y,z directions, in units of ell_Rb
column 5 to 7 - COM of K41 condensates along x,y,z directions, in units of ell_Rb
column 8 to 10 - COM of Rb87 thermal cloud along x,y,z directions, in units of ell_Rb
column 11 to 13 - COM of K41 thermal cloud along x,y,z directions, in units of ell_Rb
column 14 to 16 - COM of Rb87 whole cloud (condensate + thermal cloud) along x,y,z directions, in units of ell_Rb
column 17 to 19 - COM of K41 whole cloud (condensate + thermal cloud) along x,y,z directions, in units of ell_Rb

2) "thermalization_1compo" directory

fraction_Ec**.dat
column 1 -- w_ref*t, w_ref = 2pi * 200 Hz, in unit of 1.e-4
column 2 -- condensate fraction
ignore the rest

ave_ke_Ec**.dat
column 1 -- w_ref*t, w_ref = 2pi * 200 Hz, in unit of 1.e-4
column 2 -- average kinetic energy of the thermal atoms, in units of hbar*w_ref

width_Ec2400.dat
column 1 - omega*t
column 2 to 4 - average width of condensates along x,y,z directions, in units of ell_Rb^2
column 5 to 7 - average width of thermal cloud along x,y,z directions, in units of ell_Rb^2
column 8 to 10 - average width of whole cloud (condensate + thermal cloud) along x,y,z directions, in units of ell_Rb^2

conz_Ec2400_**.dat
column 1 - position along z axis, in unit of ell_Rb
column 3 - condensate density along z axis, in units of per ell_Rb^3
ignore the rest

thz_Ec2400_**.dat
column 1 - position along z axis, in unit of ell_Rb
column 2 - thermal cloud density along z axis, in units of per ell_Rb^3
ignore the rest


3) "thermalization_2compo" directory
***_fraction_**.dat
*** - scaled --> scaled collisions
    - full   --> full model
    - nocxc  --> no new term (exchange C12)
column 1 -- w_ref*t, w_ref = 2pi * 200 Hz, in unit of 1.e-4
column 2 -- condensate fraction
ignore the rest

conz_*.dat 
column 1 - position along z axis, in unit of ell_Rb
column 3 - condensate density along z axis, in units of per ell_Rb^3
ignore the rest

thz_*.dat 
column 1 - position along z axis, in unit of ell_Rb
column 2 - thermal cloud density along z axis, in units of per ell_Rb^3
ignore the rest


Explanation of weighting and grossing variables: 
na


Outline any missing data:
na


=======
Contact
=======
Please contact rdm@ncl.ac.uk for further information