Homework Solution: The Schlogl model is famous for its bistable steady-state distribution. The reactions are B_1 2X rightwardsharpo…

    3. Schlögl model. The Schlögl model is famous for its bistable steady-state distri- bution. The reactions are where Bi and B2 denote buffered species whose respective molecular populations Ni and N2 are assumed to remain essentially constant over the time interval of interest. There is only one time-varying species, X; the state change vectors are v, = V3 = 1, ½ = Va =-1; and the propensity functions are al (x) a2(x) 우Nix(x-1) 똥x(x-1)(x-2). = = a4(x) = car For some values of the parameters this model has two stable states, and that is the case for the parameter values we have chosen here: 3 × 10-7. C3 = 10-3, №=2x105 c; C2 = 10-4, c-3.5. Ni=1x105, In the following questions, we assume the initial condition is at X(0) 250. (a) Write down the Chemical Master Equation for the model. (b) Write down the reaction rate equations and simulate the equation in matlab. (c) Write a program to run the SSA simulation and compare the result with the reaction rate equations result. You need to show at least five different trajectories.

    The Schlogl model is famous for its bistable steady-state distribution. The reactions are B_1 2X rightwardsharpoonoverleftwardsharpoon 3X, B_2 rightwardsharpoonoverleftwardsharpoon X, where B_1 and B_2 denote buffered species whose respective molecular populations N_1 and N_2 are assumed to remain essentially constant over the time interval of interest. There is only one time-varying species, X: the state change vectors are v_1 = V_3 = 1, v_2 = V_4 = -1: and the propensity functions are a_1(x) = c_1/2 N_1x(x - 1), a_2(x) = c_2/6x(x - 1)(x - 2) a_3(x) = c_3N_2, a_4(x) = c_4x. For some values of the parameters this model has two stable states, and that is the case for the parameter values we have chosen here: c_1 = 3 times 10^-7, c_2 = 10^-4, c_3 = 10^-3, c_4 = 3.5, N_1 = 1 times 10^5, N_2 = 2 times 10^5. In the following questions, we assume the initial condition is at X(0) = 250. (a) Write down the Chemical Master Equation for the model. (b) Write down the reaction rate equations and simulate the equation in matlab. (c) Write a program to run the SSA simulation and compare the result with the reaction rate equations result. You need to show at least five different trajectories.

    Expert Answer

    3. Schlögl standard. The Schlögl standard is glorious control its biunwavering steady-specify distri- bution. The reactions are where Bi and B2 portray buffered office whose appertaining molecular populations Ni and N2 are conducive to accrue essentially firm aggravate the opportunity intermission of share. There is simply individual opportunity-varying office, X; the specify vary vectors are v, = V3 = 1, ½ = Va =-1; and the proclivity functions are al (x) a2(x) 우Nix(x-1) 똥x(x-1)(x-2). = = a4(x) = car Control some values of the parameters this standard has couple unwavering specifys, and that is the fact control the parameter values we accept separated here: 3 × 10-7. C3 = 10-3, №=2x105 c; C2 = 10-4, c-3.5. Ni=1x105, In the subjoined questions, we pretend the primal term is at X(0) 250. (a) Write down the Chemical Master Equation control the standard. (b) Write down the reaction blame equations and impersonate the equation in matlab. (c) Write a program to course the SSA assumption and assimilate the development with the reaction blame equations development. You scarcity to likeness at meanest five divergent trajectories.

    The Schlogl standard is glorious control its biunwavering steady-specify dispensation. The reactions are B_1 2X rightwardsharpoonoverleftwardsharpoon 3X, B_2 rightwardsharpoonoverleftwardsharpoon X, where B_1 and B_2 portray buffered office whose appertaining molecular populations N_1 and N_2 are conducive to accrue essentially firm aggravate the opportunity intermission of share. There is simply individual opportunity-varying office, X: the specify vary vectors are v_1 = V_3 = 1, v_2 = V_4 = -1: and the proclivity functions are a_1(x) = c_1/2 N_1x(x – 1), a_2(x) = c_2/6x(x – 1)(x – 2) a_3(x) = c_3N_2, a_4(x) = c_4x. Control some values of the parameters this standard has couple unwavering specifys, and that is the fact control the parameter values we accept separated here: c_1 = 3 opportunitys 10^-7, c_2 = 10^-4, c_3 = 10^-3, c_4 = 3.5, N_1 = 1 opportunitys 10^5, N_2 = 2 opportunitys 10^5. In the subjoined questions, we pretend the primal term is at X(0) = 250. (a) Write down the Chemical Master Equation control the standard. (b) Write down the reaction blame equations and impersonate the equation in matlab. (c) Write a program to course the SSA assumption and assimilate the development with the reaction blame equations development. You scarcity to likeness at meanest five divergent trajectories.

    Expert Reply

     

    ** Plea