Questions tagged [quantum-mechanics]
Numerical methods and problems involving the solution of the Schrodinger equation and related subatomic models.
114
questions
0
votes
1
answer
57
views
How to get a Wannier function of tight-binding model numerically (i.e. Python)
I have a question on construction of a Wannier function for tight-binding model. Let's say we consider the tight-binding model of 1D chain with two atoms( site A and B in a unit cell). In k-space we ...
1
vote
0
answers
30
views
Solving the non-linear Hamiltonian using Scipy's root finding method
I am a complete novice to computational physics and am finding difficulty in implementing a code to iteratively solve for a $2\times2$ nonlinear Hamiltonian using Scipy's root solver. I can't seem to ...
0
votes
0
answers
57
views
Transparent Boundary Conditions for Finite Difference ADI PR 2D TDSE solution
I want to put (non-dirichlet) boundary conditions inside the code I wrote to solve the 2dim TDSE using the alternating direction implicit Peaceman - Rachford method.
$$
(1 + iB\Delta t/2 ) \psi^{n+1/2}...
- 131
1
vote
1
answer
181
views
Schrodinger's Equation differential
I am working on a modified version of Schrodinger's equation (time-independent) where $\frac{d^2ψ}{dx^2}=-2(E-V)ψ$, where I have to consider $V = 0$ at all times.
I have been asked to use Python in ...
- 13
0
votes
0
answers
31
views
Solving a time-independent 2D Schrödinger equation using PDE Toolbox in MATLAB
I want to solve the following PDE eigenvalue problem,
$$-\left(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}\right)\psi_n(x,y)+x^2y^2\psi_n(x,y)=E_n\psi_n(x,y)$$
inside a square with ...
- 129
2
votes
0
answers
27
views
Numerical calculation of out-of-time order correlators (OTOCs)
I want to numerically calculate OTOCs for different quantum mechanical systems. Consider the following Hamiltonian
$$H=p_x^2+p_y^2+x^2y^2$$
and I want to calculate the following OTOC
$$C_T(t)=-\left&...
- 129
1
vote
1
answer
56
views
Roothaan equations and Galerkin method
When we do Hartree-Fock computations by solving the Roothaan equations $FC = SC \varepsilon$ is it a Galerkin method?
- 331
1
vote
1
answer
209
views
Numerical Solution of the Schrödinger equation for hydrogen
I'm trying to solve the Schrödinger equation for the hydrogen atom in the following form numerically:
$$\left[-\frac{\hbar^2}{2m}\frac{d^2}{dr^2}+V(r)+\frac{\hbar^2l(l+1)}{2mr^2}\right]R(r)=ER(r).$$
...
- 41
0
votes
0
answers
62
views
Numerov method for 1D potential barrier problem: numerical errors
I wanted to solve a 1D potential barrier problem. It seemed easy enough to do: set a travelling wave coming from the left, part of it gets reflected by the barrier, and part of it gets gets ...
- 101
0
votes
1
answer
130
views
Computing eigenvalues of Schrodinger equation with spin
I want to solve a 2-dimensional particle in box problem with two electrons in the quantum well.I would like to take into account spin of electrons and Coulomb interactions to compute singlet and ...
- 21
1
vote
0
answers
228
views
How to compute the Eigenvalue and Eigenstates of Quantum well with Effective mass using finite difference method in Python?
I want to compute the eigenvalues and eigenstates of a quantum well with different effective masses of electron in the barrier and in the quantum well. As can be seen [1]: https://github.com/mholtrop/...
- 21
0
votes
0
answers
43
views
Numerically solving schrodinger equation [duplicate]
Consider the potential given above:
$$U(r) = \frac{U_0}{\exp[(r - r_0)/\epsilon] + 1}\, .$$
How to solve the Schrodinger equation with this potential numerically and find the eigenvalues?
- 1
1
vote
0
answers
47
views
Weird behavior in for solving TISE in harmonic oscillator potential using the shooting method
I was solving the time independent Schrödinger equation using the shooting method for harmonic oscillator potential. This is the code that I wrote for that with the results (code is written in julia):
...
- 11
1
vote
1
answer
143
views
0
votes
1
answer
79
views
Problems with the time-dependent Schrödinger equation solutions
I have this Mathematica problem that solves numerically the time-dependent Schrödinger equation in a box:
...
- 23
0
votes
0
answers
373
views
Numerov method for solving Schrödinger equation
I have just begun learning computer science to apply it to Physics and I am trying to write a code for solving Schrödinger's equation of the harmonic oscillator (setting $V=\frac{x^2}{2}$) in one ...
- 27
2
votes
0
answers
85
views
Finding the extrema of a transition probability function for a quantum walker on a graph
The goal
Implement some Python code to find the extrema points of a function that is strongly oscillating.
The background
Let $G$ be a connected graph with $n$ points with Laplacian matrix $L(G)$. We ...
3
votes
0
answers
106
views
Python routine to calculate shape resonances of H2
I am currently doing a project in which my aim is to write a program that can be used to calculate single and multi-channel shape resonances. So I'm looking at bound states and quasi-bound states.
...
- 31
2
votes
1
answer
168
views
A problem with Poisson equation
I'm computing the Hartree potentials of atoms by solving the Poisson equation and I use hydrogen atom as a test case. The Poisson equation for hydrogen atom in atomic units is given by
$$\nabla^2 V_H =...
- 331
2
votes
1
answer
160
views
Why is Time evolving block decimation so efficient?
I have a short question about Time evolving block decimation (TEBD). During a lecture I was told that this method is very efficient in evolving 1D quantum spin systems with only nearest neighbor ...
- 131
2
votes
1
answer
902
views
How can I evaluate more accurate energy eigenvalues from Schrodinger equation using shooting method?
I am trying to use the "shooting method" for solving Schrodinger's equation for a reasonably arbitrary potential in 1D. But the eigenvalues so evaluated in the case of potentials that do not have hard ...
- 193
-1
votes
2
answers
299
views
2d Schrodinger Equation via matrix diagonalization in C
I am trying to solve the time-independent Schrodinger equation in two dimensions via discrete matrix diagonalization. I want the energy eigenvalues and the corresponding eigenfunctions for a given ...
- 95
-1
votes
1
answer
104
views
If not MATLAB, what software/programming language should I use to simulate/animate wave functions in various potentials + more? (example given)
I want to integrate programming into my learning in math and science in a very specific way. I want to create visualizations and simulations of concepts I am learning. When I learn a numerical method ...
- 11
1
vote
1
answer
207
views
Operator splitting to solve time dependent Schrödinger equation
I encountered the split operator method to solve the time dependent Schrödinger equation during a lecture. I understand the method on a theoretical basis (I think at least), but I'm struggling to ...
- 131
4
votes
2
answers
517
views
Numerical solution of zero-potential time-dependent Schrödinger equation in 1D
I want to solve numerically the one-dimensional time-dependent Schrödinger equation $$i \psi_t(x,t)=-\frac{\hbar}{2m} \psi''(x,t)$$
My issue is that I don't have the physical background to understand ...
- 520
4
votes
1
answer
285
views
Recommendation for a fixed-step ODE solver?
My problem involves the solution of a second-order ODE with a fixed-step (input and output). Specifically, this ODE is the radial part of Dirac and Schrödinger equation for a spherical symmetric ...
- 143
2
votes
0
answers
105
views
Numerical integration of SDE: choice of $dt$ and algorithm
I am working on the following Stochastic Differential Equation (SDE) in the Quantum Mechanics context:
$$dX_{t} = a X_{t} dt + b X_{t} dW$$
where $X_{t}$ is my stochastic varible, $dt$ is my ...
- 21
2
votes
1
answer
270
views
Numerical solution to the Landau-Zener problem
I tried to use a midpoint method and numerically solve the Schrödinger equation for the original Landau-Zener (LZ) problem: a $2\times 2$ Hamiltonian
$$\left(\begin{array}{c} \alpha t\\ \delta \end{...
- 21
2
votes
1
answer
289
views
Numerov method for Schrodinger equation
While learning about numerical methods for solving the Schrödinger equation I came across Numerov's method.
I want to get the solution for the harmonic oscillator by alreading giving the eigenvalues. ...
- 21
1
vote
1
answer
598
views
What will be the impact of quantum computing on existing numerical techniques (e.g. CFD)?
Quantum computing seems to be a very active and promising development area in computer science. However, I am curious as to what impact (if any) quantum computing will have on existing classical ...
- 121
1
vote
1
answer
240
views
How to do Weierstrass-transform in MATLAB?
I have a diagonalization problem. I have the eigenstates correctly, and I want to do a Gaussian-smearing (Weierstrass-transform) on them. So I have the wave functions ($\Psi$), and the continuous ...
- 13
0
votes
0
answers
329
views
Finding second excited state of Schrödinger equation with secant Runge Kutta method
In our assignment, we are required to find the energies of the ground state and the first two excited states of the Schrödinger equation in a harmonic potential:
$$V = \frac{50 x^2}{(10^{-11})^2}\, .$...
- 13
1
vote
0
answers
89
views
How to numerically calculate the transition dipole integral in periodic systems?
Now I have wave functions $\psi_a$ and $\psi_b$ of two states in Gaussian CUBE format. I'd like to evaluate the transition dipole moment integral $\pmb\mu$ between these two states. As my simulation ...
- 111
3
votes
1
answer
289
views
MATLAB: Faber approximation of the matrix exponential to solve Liouville-von-Neumann equation
EDIT: I moved the full code to my github page so the post can be read more easily.
I am writing a script to take the Faber approximation approach outlined in Hassan Fahs paper (free access)
and apply ...
- 73
2
votes
1
answer
491
views
Incorporating a potential barrier in a wave-packet simulation (Fourier Transform method)
I'm trying to simulate the scattering of a wave-packet at a potential barrier in Python. I'm using a Fourier Transform method (not sure if its the same as the Split-Step method), where I apply Fourier ...
- 21
1
vote
0
answers
47
views
Hydrogen-like wavefunction as starting guess for atomic solver?
I've been looking into radial solvers for quantum wave equations (Schroedinger and Dirac). In both cases, the suggestion seems to be to go with the "shooting method", with integration schemes of ...
- 191
1
vote
1
answer
4k
views
Double potential well with Python
I'm trying to understand the Schrödinger equation and solving it a bit better, and I'm running into some doubts while coding, even though I am adapting the code to this situation. Also I tried asking ...
- 113
2
votes
1
answer
156
views
Discrepancies between numerical and analytical solution for particle in a finite potential well?
Analytical
Inside the box, the wavefunction is:
\begin{equation}
\frac{\hbar^2}{2m} \frac{d^2 \psi(x)}{dx^2} = E \psi(x) \iff
\frac{d^2 \psi(x)}{dx^2} = k^2 \psi(x)
\end{equation}
where $k = \...
- 125
2
votes
1
answer
731
views
Calculate partial trace of an outer product in Python?
I have a python implementation of calculating the partial trace over select dimensions.
...
- 131
-1
votes
1
answer
852
views
Solving Schrodinger equation numerically
Here is the Schrodinger equation that is to be solved:
A 1D hard wall potential in $[0, 1]$. The potential within the potential well is given by a linear combination of Gaussian dips
$$v(x) = - \...
- 101
2
votes
1
answer
155
views
Magnetization Vector from XY Model for an AntiFerromagnetic System
I am working on an XY model and I'm trying to calculate the magnetization and direction for an anti-ferromagnetic (AF) system. So I have a collection of spins in the $XY$ plane represented as vectors ...
- 217
5
votes
0
answers
404
views
Avoiding divergent solutions with `odeint`? shooting method
I am trying to solve an equation in Python. Basically what I want to do is to solve the equation:
$$
\frac{1}{x^2}\frac{d}{dx}\left(Gam \frac{dL}{dx}\right)+L\left(\frac{a^2x^2}{Gam}-m^2\right)=0
$$
...
4
votes
1
answer
357
views
Matrix exponential of hermitian matrix with eigenvectors from generalized eigenvalue problem
I want to calculate the following expression
$$
\exp(-i\Delta t\mathbf{H})
$$
where $\mathbf{H}\in\mathbb{C}^{n\times n}$ is a hermitian matrix. Since I have a highly optimized eigensolver in the code ...
- 83
2
votes
1
answer
159
views
Two variables integration matlab
I'm trying to solve physical problem in quantum mechanics of helium atoms, the solution require numerical integration over 2 variables. However when i'm trying to run the next code
...
- 45
7
votes
1
answer
276
views
Non-hermitian discretizations in quantum mechanics
Consider the Schroedinger equation
$$\left(-\frac12\frac{\partial^2}{\partial x^2} + V(x) \right) \psi(x) = E \psi(x)$$
The usual way to solve it is to introduce a discretization of $\psi(x)$. This ...
- 2,274
2
votes
0
answers
773
views
Numerical solution to Time-dependent Schrodinger equation with time dependent hamiltonian
Currently I am facing the problem to solve numerically the following equation for a double well harmonic potential:
$iℏ\frac{\partial}{\partial t}ψ(x,t)=−\frac{ℏ}{2m}\frac{\partial ^2}{∂x^2}ψ(x,t)+V(...
- 21
2
votes
1
answer
171
views
Jump-Diffusion process: practical solver beyond Euler method?
A jump-diffusion process is a stochastic process where both continuous noise (in my case complex Wiener noise $dZ,dZ^*$ such that $dZ^2=dZ^{*2}=0,|dZ|^2=dt$) and discrete Jumps (in my case Poissonian $...
- 123
3
votes
1
answer
127
views
Normalizing a density matrix at each iteration
I need to numerically evolve a density matrix using this formula(Actually I have more terms but right nows I am starting with this and facing problems):
$$\dot\rho(t) = -i[H(t), \rho(t)]$$
$H(t)$ is ...
- 133
1
vote
0
answers
63
views
Is there a software package for Quantum Monte Carlo estimation of the exchange correlation functional?
Quantum Monte Carlo (QMC) calculations historically have been used to parameterise the exchange-correlation functionals for Density Functional Theory (DFT). For example, this article explains a way to ...
- 191
1
vote
2
answers
4k
views
Solving the 1D Particle-in-a-Box using C++
I've just finished learning the physics behind the problem and would like to write a program in C++ than can solve the problem. I'm actually stuck at the start. I've quite a bit of research, the ...
- 111