This Readme provides an introduction to the contents in the directory

1. VT rate coefficients
These are the VT rate coefficients calcalated and used in this study. They are in binary format. A sample python code to read these files is as follows

import numpy as np
binfname='./VTrate_i_L15_v5k_high.bin'
binf = open(binfname, "r")
coefarr= np.fromfile(binf, dtype=np.double).reshape(15,15,order='F')


2. N2 vibrational distribution
These are the N2 vibrational distribution versus different Ne and Te values calcalated in this study. They are in binary format. A sample python code to read these files is as follows

import numpy as np
binfname='./nv_L15_v5k_low.bin'
binf = open(binfname, "r")
arr= np.fromfile(binf, dtype=np.double)
num_Ne=16
num_Te=17
Nearr=arr[0:num_Ne]
Tearr=arr[num_Ne:num_Ne+num_Te]
N2varr=arr[(num_Ne+num_Te):(num_Ne+num_Te+15*num_Ne*num_Te)].reshape(15,num_Ne,num_Te,order='F')


3. A demonstrative IDL program to calculate transition probability
The actual calculation of VT transition conducted in this study is based on cuSOLVER (documents can be found in https://docs.nvidia.com/cuda/cusolver/index.html). However, cuSOLVER is a closed-source package and requires pre-installation. It is also hardware-dependent (GPU) and not entirely cross-platform. The program built for our calculaion is speific to hardware/platform and not designed for general use.
Here we offer a standalone IDL program to demonstrate the numerical procedures involved in the quantum calculation of transition probability, to facilitate readers to understand and quickly check the result.
The core procedure here is the linear-equation system solver. In this standalone IDL program we have used the built-in LU-decomposion functions in IDL, it can be replaced by cuSOLVER API (once installed) or any other linear equation solver.