Questions tagged [quantum-mechanics]
Numerical methods and problems involving the solution of the Schrodinger equation and related subatomic models.
0
votes
1answer
34 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 ...
-2
votes
1answer
79 views
How to plot 1/energy units in the energy (E) vs density of states (DOS) plot?
I have calculated the eigenvalues of Hamiltonian by exact diagonalization. Now I want to plot density of states (DOS) on y-axis and E on x-axis. DOS counts the number of energy levels in unit interval ...
0
votes
0answers
33 views
Boundary conditions for solving the time-independent SE for the hydrogen atom
I am trying to solve the schrodinger equation for the hydrogen atom numerically, using finite elements, with matlab's solvepdeeig(). I have a hard time getting the solution to be right, and it seems ...
0
votes
0answers
37 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}\, .$...
1
vote
0answers
33 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 ...
3
votes
1answer
201 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 ...
2
votes
1answer
57 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 ...
1
vote
0answers
27 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 ...
1
vote
1answer
375 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 ...
2
votes
1answer
69 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 = \...
1
vote
1answer
151 views
Calculate partial trace of an outer product in Python?
I have a python implementation of calculating the partial trace over select dimensions.
...
1
vote
0answers
108 views
Split step Fourier method to solve Schrödinger equation for moving potential
I'm trying to use the excellent Schrodinger Python class by Jake VanderPlas (https://jakevdp.github.io/blog/2012/09/05/quantum-python/) to simulate a wave packet within a moving Gaussian potential. I ...
-2
votes
1answer
202 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) = - \...
2
votes
1answer
76 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 ...
3
votes
0answers
143 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
1answer
86 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 ...
2
votes
1answer
58 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
...
7
votes
1answer
173 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
votes
0answers
146 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(...
2
votes
1answer
74 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 $...
3
votes
1answer
101 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 ...
1
vote
0answers
53 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 ...
1
vote
2answers
1k 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 ...
3
votes
1answer
481 views
Numerical propagation of a density matrix using Liouville von Neumann equation
I want to look at time evolution of the density matrices of some, very simple, spin systems, but I am having trouble with my approach.
I want to use a simple for-...
1
vote
1answer
114 views
Matrix exponential by eigenvectors - implementation issues
I posted a similar question yesterday but I deleted it since I think that I had to reformulate it after some insights.
I want to calculate
$$
\exp(-i\Delta t\,\mathcal{H}) = V\,\mathrm{diag}(\{\exp(-...
3
votes
2answers
2k views
Calculating partial trace of array in NumPy
A simulation I'm doing requires me to calculate the partial trace of a large density matrix. I am trying to calculate it using tools from numpy, but my code seems to be having some problems. For ...
1
vote
0answers
125 views
Second Quantization in Matlab
This question may be more suited for physics.stackexchange, but I saw this post was recommended for StackOverflow or Computational Science, so I'm asking my question here.
I am trying to write a ...
2
votes
2answers
548 views
Schrödinger equation with time dependent Hamiltonian
I need to solve the Schrödinger equation with a time dependent Hamiltonian
$$i\hbar \frac{\partial}{\partial t} \Psi = \left[-\frac{\hbar^2}{2m}\nabla^2 +\frac{1}{2} k(t)(x^2+y^2) + V(r)\right]\Psi $...
1
vote
0answers
92 views
lightweight implementation of semiempirical quantum chemistry (e.g. MNDO,AM-1,PM3)
I'm searching for semi-empirical quantum chemistry solver which would be easy to integrate into my own software. I found a few implementations which can be in principle used e.g. MOPAC, ORCA, SQM some ...
1
vote
0answers
47 views
Parallelizing molecular simulation with full configuration energy
First, let just me stress that I'm not a an expert in computation chemistry, so now the problem:
We have GCMC molecular simulation, in the Grand Canonical ensemble, using the standard metropolis ...
2
votes
1answer
2k views
Simple open-source Quantum chemistry or DFT code in C/C++
I know lot of density functional packages in fortran, including one which we are developing in our group (Fireball-DFT) but I don't like fortran very much and I would like something which is easier ...
0
votes
1answer
93 views
Stability of dark solitons in a harmonic trap
This question is based upon a research article which I am trying to reproduce. One of the main result of this paper is the condition on transverse confinement of the Bose-Einstein Condensate(BEC) to ...
1
vote
0answers
71 views
constrained quadratic binary problems and quantum adiabatic evolution
I'm going through an article with title "Solving constrained quadratic binary problems via quantum adiabatic evolution" (reference 1). And there are several points confusing me a lot.
This article is ...
0
votes
1answer
567 views
Best books and notes for beginners in the parallel FORTRAN 90 programming
I am a PhD student and my field of study is Quantum Information and Computation in the theoretical aspect. Actually I write the computer codes/program in FORTRAN 90 which are serial codes but now I ...
2
votes
0answers
289 views
Numerical solution of Dirac equation (eigenvalue problem)
Suppose we have equation of the form:
$$H \Psi = E \Psi $$
where $H$ is Dirac Hamiltonian (also my question can be answered by people who are not familiar with Dirac Hamiltonian but familiar with ...
0
votes
1answer
146 views
How to compute matrix representation of $\hat{y}\frac{\partial}{\partial x}$?
I have a 2-dimensional system which I would like to solve numerically (I'm using finite difference method right now), and its an eigenvalue problem. I have a term that looks like $H\psi(x,y) = [-\frac{...
1
vote
0answers
89 views
How to get the eigenvalues of Hamiltonian in an over complete basis
Let $|\psi_i\rangle$, $i=1...N+m$, be a set of overcomplete basis vector in a $N$-dim Hilbert space. The following are known: (Einstein's summation convention assumed)
$$\hat{H}|\psi_i\rangle=H_{ji}|\...
6
votes
2answers
726 views
Numerical Solution to Schrödinger Equation--Multiple Wells
I am trying to solve for the allowed wavefunctions and energies for a 1D quartic potential well.
To do this I am using the patching method (https://engineering.dartmouth.edu/microeng/otherweb/...
2
votes
2answers
516 views
Zero Eigenvalues in Lanczos Algorithm
I need to find the smallest few eigenvalues of a Hamiltonian (exact diagonalization) I use Python, and SciPy's built-in sparse eigenvalue solver. I notice, however, that for my small system (only a ...
0
votes
1answer
1k views
Solving Time dependent Schrodinger equation using MATLAB ode45 [closed]
The Schrodinger equation for time-dependent Hamiltonian is
$$i\hbar\frac{d}{dt}\psi(t) = H(t)\psi(t) \, .$$
I try to implement solve the Schrodinger equation for time-dependent Hamiltonian in ODE ...
2
votes
0answers
117 views
Solving Schrodinger's Equation Numerically in a Bunimovich Stadium
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 ...
0
votes
3answers
108 views
Quantum Chemical Calculations is there a book for which method to use with what problem?
Does anyone know of a book that will outline which quantum chemical methods are appropriate for what problems? I am trying to make informed choices before I start using computational resources. It is ...
1
vote
1answer
1k views
imaginary time propagation to find ground state wavefunction
I understand the basic idea of imaginary time propagation method:
The wavefunction $\psi(x,t)$ as a superposition of energy eigenstates $\phi_m(x)$:
$$
\psi(x,t)=\sum_m \phi_m(x)e^{-iE_mt/\hbar}
$$
...
3
votes
1answer
370 views
ground state from the Schroedinger equation with a central potential what happens to the origin
I have code that attempts to implement a solution to the Schrödinger equation where there is a central potential (more or less im thinking of hydrogen),
in 1-D using the numerov method to construct ...
-1
votes
1answer
572 views
Gaussian and external basis set. How to write the input correctly [closed]
I have to use the LANL2DZpd basis set. It can be obtained from EMSL Basis Set Exchange.
But trying to calculate simple Me2Se molecule :
...
1
vote
1answer
2k views
Transfer Matrix Method in a rectangular potential well
I am trying to follow an algorithm that is described in Elementary Quantum Mechanics in 1D.
I want to compute eigen-energies and functions in bound states in the basic case in rectangular potential ...
3
votes
0answers
331 views
C60 orbital calculation
I am currently trying to reproduce the results published by Hands et al.1 using MATLAB. They calculated the bases of the C60 wave functions of HOMO, LUMO and LUMO+1.
I did the following:
I ...
7
votes
1answer
190 views
Forcing an ODE solver to preserve the norm
I have an ODE of the form
$$ \frac{dy}{dt} = -i H y \enspace .$$
where $y$ is a complex vector and $H$ is a time dependent Hermitian matrix.
The norm of the solution $y(t)$ at any point in time ...
6
votes
1answer
2k views
Computational methods for finding the energy eigenvalues of the time-independent Schrodinger equation with arbitrary potential
I have seen in some papers that the energy levels in some arbitrary potential are specified. How can one find the energy levels in such arbitrary potentials. For example, $V(x)=\sin^2(x/2)$ with $x\in[...
1
vote
2answers
760 views
Parallelization of element-wise matrix multiplication
I use Armadillo as an interface to OpenBLAS. In my current program, I have a loop, in which I do multiplications of the form
...
