supersolids.helper.cut_1d

Animation for the numerical solver for the non-linear time-dependent Schrodinger equation.

Module Contents

Functions

cut_1d(System, slice_indices = [0, 0, 0], psi_sol_3d_cut_x = None, psi_sol_3d_cut_y = None, psi_sol_3d_cut_z = None, dir_path = Path(__file__).parent.parent.joinpath('results'), y_lim = (0.0, 1.0))

Creates 1D plots of the probability function of the System \(|\psi|^2\)

prepare_cuts(func, N, alpha_z, e_dd, a_s_l_ho_ratio)

Helper function to get \(R_r\) and \(R_z\) and set it for the given func.
supersolids.helper.cut_1d.cut_1d(System, slice_indices=[0, 0, 0], psi_sol_3d_cut_x=None, psi_sol_3d_cut_y=None, psi_sol_3d_cut_z=None, dir_path=Path(__file__).parent.parent.joinpath('results'), y_lim=(0.0, 1.0))

Creates 1D plots of the probability function of the System \(|\psi|^2\) and if given of the solution.

Parameters

  • System (supersolids.Schroedinger.Schroedinger) – Schrödinger equations for the specified system

  • slice_indices (numpy.ndarray) – Numpy array with indices of grid points in the directions x, y, z (in terms of System.x, System.y, System.z) to produce a slice/plane in mayavi, where \(\psi_{prob}\) = \(|\psi|^2\) is used for the slice Max values is for e.g. System.Res.x - 1.

  • psi_sol_3d_cut_x (Optional[Callable]) – 1D function after cut in x direction.

  • psi_sol_3d_cut_y (Optional[Callable]) – 1D function after cut in y direction.

  • psi_sol_3d_cut_z (Optional[Callable]) – 1D function after cut in z direction.

  • dir_path (pathlib.Path) – Path where to save 1d cut plots

  • y_lim (Tuple[float, float]) – Limit of y for plotting the 1D cut

Return type

None

supersolids.helper.cut_1d.prepare_cuts(func, N, alpha_z, e_dd, a_s_l_ho_ratio)

Helper function to get \(R_r\) and \(R_z\) and set it for the given func.

Parameters

  • func (Callable) – Function to take cuts from

  • N (int) – Number of particles

  • alpha_z (float) – Ratio between z and x frequencies of the trap \(w_{z} / w_{x}\)

  • e_dd (float) – Factor \(\epsilon_{dd} = a_{dd} / a_{s}\)

  • a_s_l_ho_ratio (float) – \(a_s\) in units of \(l_{HO}\)

Returns

func with fixed \(R_r\) and \(R_z\) (zeros of \(func_{125}\)), if no singularity occurs, else None.

Return type

Optional[Callable]