Basis.m

Basis set control. This is the second mandatory function (after create.m) that must be called in every calculation to build the spin_system data structure.

Syntax

    spin_system=basis(spin_system,bas)

Parameters

    spin_system   - primary Spinach data structure, the output
                    of create.m function

    bas           - basis set specification structure described
                    in detail in the online manual

Outputs

    spin_system   - primary Spinach data structure, updated with
                    the basis set and related information

Examples

A simple example of a spin system and basis set specification is given below. There are hundreds of further examples in the example set.

    % Set the spin system
    sys.isotopes={'1H','1H'};
    inter.zeeman.scalar={0.0 0.0};
    inter.coordinates={[0.00 0.00 0.00];
                       [0.00 0.00 2.00]};
    
    % Magnet field
    sys.magnet=14.1;
    
    % Basis set
    bas.formalism='sphten-liouv';
    bas.approximation='none';
    
    % Spinach housekeeping
    spin_system=create(sys,inter);
    spin_system=basis(spin_system,bas);

Notes

It is important to understand the factors that influence basis set selection in spin dynamics simulations - see our paper (http://link.aip.org/link/doi/10.1063/1.3624564) for further information on this subject.

Bosonic mode couplings declared in inter.modes are added to the coupling graph from which reduced (IK-n) basis sets are built, so that spin-mode correlations survive the approximation. Complete basis sets are unaffected.

See also

create.m, assume.m, operator.m, state.m, hamiltonian.m, evolution.m, reduce.m, thermalize.m, device.m, Kernel_functions

Version 2.4, authors: Ilya Kuprov