Creating your own MCTDH-X simulations

After redoing the simulations provided in this tutorial, the readers are encouraged to create their own simulations tailored to their own specific problems at hand. In this case, the following steps should be taken:

- Create a working folder, copy the exemplary input file
`MCTDHX.inp`into the working folder using command`inpcp`and modify it according to the following steps:- Determine the type of the problem, i.e., whether it is static (ground state) or dynamic (time-evolving state). Related input variables in
`MCTDHX.inp`:`Job_Prefactor`. - Determine the Hamiltonian of the problem, i.e., determine the number and statistics of the particles, the one-body potential and the two-body interaction. Related input variables in
`MCTDHX.inp`:`JOB_TYPE`,`Npar`,`whichpot`,`parameter1`,`parameter2`, ,`xlambda_0`,`which_interaction`,`Interaction_Type`,`Interaction_Parameter1`,`Interaction_Parameter2`... Customized potentials and interactions can be created using the source files`Get_1bodyPotential.F`and`Get_InterParticle_Potential.F`, respectively. - Determine the other physical parameters used in the simulation, including the number of orbitals, the dimension of the system, the spatial extent, the spatial grid, etc. Related input variables in
`MCTDHX.inp`:`Morb`,`DIM_MCTDH`,`NDVR_X`,`x_initial`,`x_final`, etc. - Determine the initial state of the problem. This is usually done by using a randomly generated initial state or preparing the initial state with another simulation. In the latter case, the user should
set up the ground state or a time-evolving state of a particular Hamiltonian as initial state
with another simulation directory. Then the user should copy the generated files
`Header`,`PSI_bin`and`CIc_bin`into the working folder. Related input variables in`MCTDHX.inp`:`GUESS`,`Binary_Start_Time`. - Determine the integration parameters. Related input variables in
`MCTDHX.inp`:`Time_Begin`,`Time_Final`,`Output_TimeStep`,`Integration_Stepsize`,`Coefficients_Integrator`.

- Determine the type of the problem, i.e., whether it is static (ground state) or dynamic (time-evolving state). Related input variables in
- Copy the executables
`MCTDHX_gcc`/`MCTDHX_intel`and`analysis_gcc`/`analysis_intel`into the working folder using command`bincp`. Run the executable using command`MCTDHX`. - Sanity test of the result
- Check the one-body potential in
`***orbs.dat`to see if they are correctly generated as specified by the Hamiltonian. Check the evolution of the energy and orbital occupations as functions of real/imaginary time in the generated file`NO_PR.out`and the evolution of the real-space density distribution in`***orbs.dat`to see if they are physically sensible. If not, confirm that there is no mistake/typo in the input file. - Check if the real-space density distribution vanishes at the boundary. If not, use a wider spatial extent.
- Test the convergence in different parameters. Use a different spatial grid and orbital number and check if the results remain unchanged. A detailed discussion on the convergence in orbital number can be found in Sec. S3 and in Ref. [1,2].

- Check the one-body potential in
- Analysis of the result
- Specify the desired time interval of analysis. Related input variables in
`analysis.inp`:`Time_From`,`Time_To`. - Specify the desired quantities of interest. Refer to the manual at http://ultracold.org/data/manual
for related input variables in
`analysis.inp`. - Visualize the generated quantities of interest.

- Specify the desired time interval of analysis. Related input variables in