I need to solve, as mentioned, Schrodinger's equation in a Bunimovich stadium-shaped infinite potential well with Dirichlet BC Numerically (this isn't possible analytically). In order to do so, I need to define the discrete Hamiltonian matrix of the system. Since the shape isn't defined well in any coordinates system, I have no idea how to do so.

Normally (as I did with a square well), I'd define the size of the well and then solve 1D TISE for each dimension using MATLAB eig function. Can it be solved similarly?

Any help will be greatly appreciated.

  • $\begingroup$ Can you describe the equations and the form of the potential you want to solve for? This may allow people to understand the question who are not already familiar with the Bunimovich Stadium. $\endgroup$ Apr 28, 2016 at 14:32
  • 1
    $\begingroup$ I don't see how is different from this question. You can use Finite Differences or Finite Methods for that purpose. BTW, the drawing you are referring to is really confusing. $\endgroup$
    – nicoguaro
    Apr 28, 2016 at 19:50
  • 1
    $\begingroup$ Not different at all. In the approaches I mentioned (and in the codes) the potentials are turn into a discrete representation. $\endgroup$
    – nicoguaro
    Apr 28, 2016 at 21:06
  • 1
    $\begingroup$ You still discretize the 2D domain, and the row and column numbers correspond to your discretization process. It depends if you use Finite Differences or Finite Elements but in essence is the same. $\endgroup$
    – nicoguaro
    Apr 28, 2016 at 23:15
  • 2
    $\begingroup$ That's right. You replace your domain by a set of triangles or quadrilaterals, and then you just do the regular finite element approach on it (or finite differences, if you're so inclined). $\endgroup$ Apr 29, 2016 at 2:45


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.