$Device
{

  Region well
  {
    material = GaAs

    doping = 1e15
    doping_type = donor
    doping_level = 0.025

  }

  Region  buffer1
  {
    mesh_regions  = ( barrier1_a,barrier1_q)
    material = AlGaAs
    x = 0.3

    doping = 1e17
    doping_type = donor
    doping_level = 0.025
    
   
  }
  
  

  Region  buffer2
  {
    mesh_regions  = (barrier2_a, barrier2_q  )
    material = AlGaAs
    x = 0.3

    doping = 1e17
    doping_type = donor
    doping_level = 0.025
   
  }


  Cluster Quantum_1
  {
    mesh_regions = (well ,barrier1_q, barrier2_q   )
  }


}

$Models
{

  model macrostrain
  {

    options
    {          
      simulation_name = str

      physical_regions = all
    }

  
  
    
    BC_Regions  
    {
      BC_Region substrate
      {
	type = substrate
	material = GaAs
      }
    }

  } 


  model  driftdiffusion
  { 

    options
    {          
      simulation_name = dd
      
      physical_regions = all
    }

    

  }


 


  model efaschroedinger

  {

   
    options
    {
     simulation_name = quantum_el
     physical_regions = Quantum_1
    }


  }


 
  

  model quantumdensity
  {

  
    options 
    {
      simulation_name = dens_el
      physical_regions =  Quantum_1
    }


  }

  



}

$Solver
{

 macrostrain
  {
    substrate = substrate

  }

 
  driftdiffusion
  {
    coupling = poisson
  }



  quantum_el
  {
    Dirichlet_bc_everywhere = false
    particle = el
    number_of_eigenstates = 22
    poisson_model_name = dd
    strain_model_name = macrostrain
    convergent_density = true

  }



  dens_el
  {
    k_space_dimension = 2
    k-space_basis = true
    k1 = (0, 0, 0.1)
    k2 = (0, 0.1, 0)
    

    number_of_nodes = (4,4)
    wedge = quarter

    quantum_simulation = quantum_el

    degeneracy = 2
    refine_fraction = 0.20
    relative_accuracy = 0.01
    refine_k_space = true
    output_density_in_k_space = true
    uniform_refinement = false
    mesh_order = FIRST

    initial_eigenstates_number = 10 #20
    analytic = true # false # true # 
#   --------  if  false ->  k integration  performed,  if true,  analitical formula for 1D 
  }


  Selfconsistent
  {
    poisson_schroedinger
    {
      simulations = (dens_el, dd)
      max_iterations = 10
      abs_tolerance = 1e-4
      rel_tolerance = 1e-6
    }
  }

}


$Physics
{

  driftdiffusion
  {
    electron_quantum_density = dens_el
    statistics = FD
    use_density_predictor = true
  }

  quantum_el 
  {
    model = conduction_band    
  }


}

$Simulation
{
  #searchpath = ../../materials       
  meshfile = quantum_poisson.msh
  mesh_units = 1e-9
  dimension = 1
  temperature = 300
  solve = (str, dd, poisson_schroedinger)
  resultpath = output
  output_format = grace
  plot = (strain, polarization, Ec, Ev, QFermi_e, QFermi_h, EField , eDensity, hDensity, EigenFunctions, EigenEnergy, EnergyLevels,quantum_density)
}