https://github.com/caley/OpenEPSim
Raw File
Tip revision: 7e09553e6bd7750f808546745e6a4e819193313b authored by Caley Finn on 30 September 2018, 09:31:54 UTC
Added JOSS status badge
Tip revision: 7e09553
asep.json
// Sample configuration for the ASEP in the high density phase
// corresponding to parameters
//
//    p = 1, q = 0.3, a = 2, b = 3, c = -0.2, d = -0.2
//
// See https://arxiv.org/abs/cond-mat/0312457 for definitions of the
// parameters and the ASEP phase diagram, and (for example)
// https://arxiv.org/abs/cond-mat/0609645 for the definition of the
// transition matrices.
//
// Expected density rho = 0.75, current J = 0.13
{
    // Lattice length
    "L": 100,

    // Simulation length (time)
    "tMax": 1000000,

    // Start recording statistics from t = 0
    "tSkip": 0,

    // Number of states each site can be in (empty or occupied for the ASEP)
    "nstates": 2,

    // Specify the transitions
    "transitions": [

        // Left boundary: particles enter with rate 0.29 and exit with
        // rate 0.12 (exit rates are calculated automatically).
        //
        // This local transition matrix affects one site (window = 1), and is
        // applied to site 1 (from = to = 1).
        {
            "window": 1,
            "from": 1,
            "to": 1,
            "count": true,
            "rates": [
                [0, 0.12],
                [0.29, 0]
             ]
        },

        // Bulk hopping: particles hop right into empty sites with
        // rate 1, and left with rate 0.3.
        //
        // This local transition matrix affects two sites (window = 2), and is
        // applied to sites 1 to L - 1(from 1, to = -2, negative indices count
        // from the right starting at -1 for site L).
        {
            "window": 2,
            "from": 1,
            "to": -2,
            "rates": [
                [0, 0, 0, 0],
                [0, 0, 1, 0],
                [0, 0.3, 0, 0],
                [0, 0, 0, 0]
             ]
        },

        // Right boundary: particles exit with rate 0.22 and exit with
        // rate 0.13
        //
        // This local transition matrix affects one site (window = 1), and is
        // applied to site L (from = to = -1).
        {
            "window": 1,
            "from": -1,
            "to": -1,
            "rates": [
                [0, 0.22],
                [0.13, 0]
             ]
        }
    ]
}
back to top