How to solve 2d problems numerically.

  • Context: Graduate 
  • Thread starter Thread starter ehj
  • Start date Start date
  • Tags Tags
    2d
Click For Summary
SUMMARY

This discussion focuses on solving two-dimensional problems in quantum mechanics numerically, specifically using the Hamiltonian \(\hat{H} = \frac{\widehat{p}_{x}^{2}+\widehat{p}_{y}^{2}}{2m}+x^{2}y^{2}\). Two effective methods are highlighted: first, employing a basis of harmonic oscillator eigenfunctions combined with a diagonalization routine for symmetric matrices, such as LAPACK; second, utilizing a grid of points with finite difference approximations for derivatives, followed by matrix diagonalization. Both methods are essential for finding eigenvalues and eigenfunctions in complex, non-separable systems.

PREREQUISITES
  • Understanding of quantum mechanics principles, particularly Hamiltonians
  • Familiarity with numerical methods for solving differential equations
  • Knowledge of linear algebra, specifically matrix diagonalization
  • Experience with LAPACK for numerical linear algebra operations
NEXT STEPS
  • Research the implementation of LAPACK for diagonalizing symmetric matrices
  • Learn about finite difference methods for numerical differentiation
  • Explore the use of harmonic oscillator eigenfunctions in quantum mechanics
  • Study grid-based numerical methods for solving partial differential equations
USEFUL FOR

Quantum mechanics students, physicists working on numerical simulations, and computational scientists interested in advanced numerical methods for solving two-dimensional problems.

ehj
Messages
79
Reaction score
0
I havn't had much classes on numerical methods in quantum mechanics and I'm wondering how one would solve a general problem involving 2d motion. With general, I mean a problem that cannot be separated. Consider for instance the hamiltonian

[itex]\hat{H} = \frac{\widehat{p}_{x}^{2}+\widehat{p}_{y}^{2}}{2m}+x^{2}y^{2}[/itex]

How does one find the eigenvalues and eigen functions numerically?
 
Last edited:
Physics news on Phys.org
1. Use a basis of e.g. harmonic oscillator eigenfunctions and a diagonalization routine for symmetric matrices (e.g. Lapack).
2. Use a grid of points and finite difference approximation for the derivatives. Then diagonalize the matrix like in 1.
 

Similar threads

  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 14 ·
Replies
14
Views
2K
Replies
1
Views
2K
  • · Replies 13 ·
Replies
13
Views
1K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 4 ·
Replies
4
Views
8K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 7 ·
Replies
7
Views
3K