5.4.12 Dust: dustscheme=dust_bins_01

In this section the parameters for dustscheme=dust_bins_01 (multi-size-bin dust scheme for Forsterite dust) are described:

• integer n_dustgrainradius:
  The number of dust grain radius bins (including one for the monomers) is specified with
  integer n_dustgrainradius f=I8 b=4 n='Number of dust grain radius bins' &
    c1='Includes: bin for monomers'
  8
  
  The default value is 1 (to avoid an empty array). The value should be the same as the upper dimension in real ar_dustgrainradius.

• real ar_dustgrainradius:
  The radii of the dust grains are specified with e.g.
  real ar_dustgrainradius f=E15.8 b=4 p=1 d=(1:8) n='Dust grain radii' u='cm'
  0.0
  0.00390625E-04
  0.0156250E-04
  0.06250E-04
  0.25000E-04
  1.00000E-04
  4.00000E-04
  16.0000E-04
  
  The default value is 0.0.

• real c_dust01:
  The density of the grain material is specified with e.g.
  real c_dust01 f=E15.8 b=4 n='Density of grain material' u=g/cm^3
  3.3E+00
  

• real c_dust02:
  The atomic weight of the "dust monomer" is specified with e.g.
  real c_dust02 f=E15.8 b=4 n='Atomic weight of dust monomer' u=u
  140.71
  
  The value 140.71 should be appropriate for Forsterite (Mg$_2$SiO$_4$). The unit is the atomic mass unit. The default value is 0.0.

• real c_dust03:
  The number fraction of the rarest component is specified with e.g.
  real c_dust03 f=E15.8 b=4 n='Number fraction of rarest component' u=1 &
    c0='For Mg2SiO4, Mg is the rarest, its fraction is 3.1187E-05' &
    c1='Mg2SiO4 contains 2 Mg -> * 1/2 -> 1.55935e-05'
  1.55935e-05
  
  For Mg$_2$SiO$_4$ the rarest component is (usually) Mg$_2$. So, half the magnesium abundance is required as value. The default value is 0.0.

• real c_dust04:
  The lower dust limit fraction is specified with e.g.
  real c_dust04 f=E15.8 b=4 n='Lower dust limit fraction' u=1
  1.0E-10
  
  If the density in a bin relative to the gas density falls below this value, the bin is completely emptied and the material is added to the bin with the largest dust or monomer density. Typical values are 1.0E-09 or so. The default value is 0.0.

• real c_dust06:
  The sticking coefficient for the condensation rate is specified with e.g.
  real c_dust06 f=E15.8 b=4 n='Sticking coefficient' u=1 &
    c0='For C-rich dust: 0.37; can be between 0 and 1'
  1.0
  
  Typical values are between 0.0 (deactivating this process) and 1.0. The default value is 0.0. The standard choice is 1.0.

• real c_dust07:
  The sticking coefficient for the coagulation rate is specified with e.g.
  real c_dust07 f=E15.8 b=4 n='Sticking coefficient for coagulation' u=1 &
    c0='Typically between 0 and 1'
  1.0
  
  Typical values are between 0.0 (deactivating this process) and 1.0. The default value is 0.0.

• real c_dust08:
  The sticking coefficient for the coalescence rate is specified with e.g.
  real c_dust08 f=E15.8 b=4 n='Sticking coefficient for coalescence' u=1 &
    c0='Typically between 0 and 1'
  1.0
  
  Typical values are between 0.0 (deactivating this process) and 1.0. The default value is 0.0.

• real c_dust09:
  The sticking coefficient for the nucleation rate is specified with e.g.
  real c_dust09 f=E15.8 b=4 n='Sticking coefficient for nucleation' u=1 &
    c0='Typically between 0 and 1'
  1.0
  
  Typical values are between 0.0 (deactivating this process) and 1.0. The default value is 0.0. The standard choice is 1.0.