Subdiscipline of theoretical physics which consists in the study of physics problems using numerical algorithms. PLEASE NOTE that questions about computational methods and/or programming are OFF-TOPIC on Physics.SE.

learn more… | top users | synonyms (1)

6
votes
1answer
801 views

Numerical analytic continuation for Green's function

Recently, I happened to hear about the possibility of doing analytic continuation numerically. That sounds attractive for the ubiquitous $\mathrm{i}\omega_n\rightarrow\omega+\mathrm{i}0^+$ procedure, ...
0
votes
1answer
65 views

Numerical Computation of Linbald Equation [closed]

Can anybody suggest me a good algorithm for the time evolution of the reduced density matrix using Linbald equation. My Hamiltonian is time dependent. I am aware about Qotoolbox and Qutip. I have ...
0
votes
2answers
434 views

Numerical solving of 2D and 3D Schrodinger equations

I am studying 2D quantum scattering models for my Bachelor's thesis. Somewhat like these: ,with Dirichlet ($\psi \mid_\Gamma = 0$) boundary conditions on the "walls" of the waveguide and the ...
7
votes
3answers
324 views

Runge Kutta Method for a Lindblad Equation

I am solving a Lindblad equation for a dissipative Harmonic Oscillator. My Hamiltonian is time dependent, My Lindblad Equation can be written as \begin{equation} ...
5
votes
1answer
319 views

Discretizing the Wave Equation in polar coordinates

I want to discretize the wave equation $$\frac{1}{c^2}\frac{\partial^2\psi\left(\vec{r},t\right)}{\partial t^2}=\triangle\psi\left(\vec{r},t\right)$$ in polar coordinates. I find the following ...
14
votes
3answers
798 views

What limitations are there in measuring physical properties accurately?

In a StackOverflow answer, I attempted to explain why a 32-bit float was perfectly adequate for representing the questioner's weight measurement: Physical properties are inaccurately measured ...
1
vote
1answer
387 views

I'm getting weird autocorrelations when simulating an Ising model below the critical temperature

So I'm simulating an Ising model using Monte Carlo and the Metropolis algorithm. After letting it reach equilibrium, I try to calculate the autocorrelation of the magnetization. As long as the system ...
7
votes
2answers
208 views

Feynman's infinite amount of logic for one tiny bit of space

Watching one of Feynman's lectures, I came across something that puzzled me. What was Feynman referring to when he said the following? What goes on in no matter how tiny a region of space and no ...
7
votes
2answers
233 views

Monte-Carlo and $O(n)$ models for non-integer n

$O(n)$ lattice statistical models can be generalized to non integer values of n, starting from their (expanded and resumed in graphs) partition function: $$Z = \sum_{\text{loop configurations}} n^{\# ...
3
votes
2answers
201 views

Simulate the universe?

Alright, Lets assume that I have a computer with limited calculation speed (1-4GHz) but unlimited parallel processing capability and unlimited memory capacity to go with it. Under this assumption ...
2
votes
1answer
262 views

Exact diagonalization to resolve ground state degeneracies

I am studying a perturbed Toric Code model that is not analytically solvable. On a torus the ground state degeneracy of the unperturbed model is 4. Once we turn on the perturbation there is a change ...
1
vote
0answers
579 views

Hysteresis Curve and how to implement it using Preisach model or other models

As a homework I need to draw a hysteresis curve (preferably an interactive one) using Matlab or any other programming language. The problem is I have trouble finding a good algorithm to do so. I need ...
3
votes
2answers
115 views

Numerical Tools to find Braiding Statistics of Quasiparticles

While certain classes of systems that exhibit topological order can be solved exactly (such as the Toric Code, Abelian FQH Edges, etc.) there also exist systems (think of perturbed versions of the ...
3
votes
1answer
277 views

Epsilon Tensor in FeynCalc

A few days ago I started to use the Mathematica package FeynCalc and one thing confuses me: Assume we have a four-vector $p_\mu$ and we contract it with the epsilon tensor. FeynCalc produces ...
4
votes
1answer
427 views

How to solve the heat equation for compound materials with different heat conductivities numerically?

I'm solving the heat equation with time dependent boundary conditions numerically in a 2D system using the ADI scheme. For the purpose of this question, let's assume a constant heat conductivity and ...
2
votes
1answer
175 views

How much computer power would we need to model every physical aspect of the universe [closed]

So I know there are equations that give you the total energy in the Universe(Friedman). So about if we we to simulate every single interaction from neurons to physical forces, strong and weak, ...
3
votes
1answer
132 views

Why is the objective function used in Nudged Elastic Band method reasonable?

In Nudged Elastic Band (NEB) method, which is used to find reaction pathways when both initial and final states are know, an objective function is first constructed and then minimized to find reaction ...
11
votes
3answers
1k views

Correct way to write the eigenvector of a diagonalized hamiltonian in second quantization

I am studying diagonalization of a quadratic bosonic Hamiltonian of the type: $$ H = \displaystyle\sum_{<i,j>} A_{ij} a_i^\dagger a_j + \frac{1}{2}\displaystyle\sum_{<i,j>} [B_{ij} ...
0
votes
1answer
152 views

How to Extend Relaxation Methods for 2D Laplace Equation given in Jackson E&M to 3D?

In Jackson (3 ed) chapter 1.13 an outline is given for using relaxation to solve laplace equation in 2D. The general procedure in 3D involves minimizing the quantity $$I[\Psi]=\frac{1}{2}\int_V ...
3
votes
0answers
303 views

Comparison between Cadabra and other Symbolic Computer Algebra software [closed]

Does anyone has some experience about working with Cadabra and it's (dis)advantage in comparison to other Symbolic Computer Algebra software such as Maple and Mathematica (physics package) in the ...
1
vote
1answer
154 views

Rigid body collision, 3 circles in contact

I'm working on a 2d physics simulation. It's a continuous time simulation, that is, it uses swept shapes over the time-frame and geometrical/vector 'analysis' to determine most immediate time of ...
1
vote
0answers
59 views

How to get discretization coefficients of matrix $A$ in Finite Volume Method (FVM)? [closed]

First we have Discretization of the Transport Equation $$ \frac{\partial \rho \phi}{\partial t} + \nabla(\rho U \phi) - \nabla (\rho \Gamma_\phi \nabla \phi) = S_\phi (\phi) $$ In Finite Volume ...
0
votes
1answer
87 views

Failure of a simple stat mech simulation

so I did a simple simulation that I thought would yield a Boltzmann distribution, but it failed to, and I was wondering if anyone has insight into why it failed. Ok, so I had a simple discrete system ...
1
vote
1answer
79 views

Getting nonphysical results when solving for the index of refraction of a slab?

I'm trying to computationally find the refractive indices (real and imaginary) for a thin slab suspended in air (so the only indices to deal with are air and my material's). I've experimentally taken ...
4
votes
0answers
1k views

Numerical problem in solving the Bogoliubov de Gennes equations- methods to solve?

I am trying to solve an assignment on solving the Bogoliubov de Gennes equations self-consistently in Matlab. BdG equations in 1-Dimension are as follows:- $$\left(\begin{array}{cc} ...
1
vote
0answers
79 views

Application of wavelets in computational electromagnetics

I am a first-year graduate student (electrical eng.); wavelets is one of the courses that I am taking presently. As part of the course, I plan to explore application of wavelets in computational ...
0
votes
2answers
198 views

Modelling multiple heat sources?

I am working on a project to model heat transfer in a data center. I have basic models for a one dimensional heat transfer using conduction and convection. Such as, taking a server and modelling the ...
3
votes
1answer
241 views

confusion in discrete transform to solve kronig penney matrix equation in fourier space

I have a periodic potential $$V(x) =\sum_{K}e^{iKx}V_{K} =\sum_{n}e^{\iota2\pi nx/a}V_{n} $$ where $K =\frac{2\pi n}a$ is the reciprocal lattice vector and $a$ is the lattice constant and $n =\pm ...
1
vote
0answers
577 views

Numerical solution of two coupled second order differential equations of motion [closed]

Is there a numerical algorithm for solving a pair of coupled second order differential equations? This question arises from a homework problem that I have that involves two dimensional projectile ...
1
vote
0answers
125 views

How does the Lennard Jones Potential changes for interaction between particles of different sizes?

I am interested in incorporating a Lennard-Jones potential in a simulation. When the interaction only involves the same type of particle, with same characteristics, we can use reduced units, scaling ...
6
votes
2answers
274 views

Eigenvalue problem for differential equations in QM

I have a very simple question with regard to numerical methods in physics. I want to solve the eigenvalue problem for a particle moving in an arbitrary potential. Let's take 1D to be concrete. I.e. I ...
1
vote
2answers
196 views

What is a good model for computing water dropping on a surface?

Before introducing my question I would like to outline the fact that I'm a coder, so I can be wrong when using some kind of terminology. What is a good model for computing the flow, the shape and all ...
1
vote
0answers
120 views

Gerstner Wave Formula to Vertex information $(x,y,z)$

I have created a program, that draws a plane mesh, and allows to me edit or animate all the separate vertices along it. The idea behind it, is for me to be able to create a gridMesh based wave ...
1
vote
1answer
117 views

Motion of rigid body system in absense of work

In the absence of work on the system, is there a closed form equation for the motion of a set of constrained rigid bodies (let's say, using Revolute (ie: simple pivot) constraints)? If the bodies are ...
2
votes
2answers
183 views

Sampling a distribution (from a galaxy model)

I am reading the following article: http://www.kof.zcu.cz/st/dis/schwarzmeier/galaxy_models.html and am currently at section 5.6 (positions of bodies in a galaxy). I am trying to redo the simulations ...
1
vote
0answers
605 views

Is Bremermann's limit redundant with Landauer's limit for all practical purposes?

Trying to understand the physical limits to computation, I notice that among these we have two types of limits that constrain the minimum allowable energy for a computation. Limits that constrain ...
1
vote
1answer
207 views

Smoothed Particle Hydrodynamics, Chemical Reactions and Stiffness [closed]

Dear people from stackexchange physics. I have been using SPH for quite a while to simulate free boundary flow, and just recently we tried to include some kind of (simple) chemical reaction in our ...
4
votes
1answer
153 views

How to find positions of $n$ masses in Newton mechanics?

I ran into a problem while doing research. The problem can be described as: consider the original $n$-body problem, and if we fix the position of them(unknowns), no interaction among them, they don't ...
1
vote
1answer
70 views

Find out ground sates for large 2D classical spin model

Reaching the ground state of a large 2D classical spin model (e.g. classical Heisenberg model) might be a relatively difficult task while using conventional "flip/reject" Monte Carlo method. The ...
2
votes
3answers
170 views

Do algorithms have an intrinsic time direction?

This article says There is no intrinsic time direction in Newton's mechanics nor in the differential equations of the new physics. My question is, do other types of mathematics, say a cellular ...
4
votes
2answers
282 views

Can computers accurately model all of the details (to the subatomic level) of macro objects in collisions?

Frequently when trying to solve cosmology questions physicists turn to computer simulations of the universe (albeit massively simplified) in order to verify or disprove their hypotheses. This got me ...
10
votes
3answers
2k views

How are physics and computer science getting united? [closed]

How is theoretical computer science getting united with physics? Phenomena like Quantum Computing uses Quantum Mechanics to be able to compute things, how are computers helping not just to model our ...
-2
votes
1answer
150 views

I am looking for 'Programming+Electronics+ Physics' field [closed]

I wish to pursue a career that somehow involves programming, electronics, and physics. What are such careers? Also, I have heard of some 1 year post-graduate diplomas/courses for specialization in ...
1
vote
2answers
137 views

Discrete approximation of charge density

Given the electric potential $\Phi(r)$ and the Poisson's equation: $$ \nabla^2 \Phi(r) = - 4\pi \rho(r)$$ Consider the 2-dimensional case and let's say that I want to discretize this using a square ...
1
vote
1answer
95 views

what is two species Bose Hubbard Model? what chage it takes on Fock basis representation?

The Hamiltonian will look like \begin{equation} \ H= -t^a\sum_{<i,j>}^{}(a_i^\dagger a_i + h.c.)-t^b\sum_{<i,j>}^{}(b_i^\dagger b_i + ...
0
votes
0answers
51 views

Should particles in Smoothed Particle Hydrodynamics (SPH) always move during a simulation?

Or can they just be used as an interpolation points and use some other "transported property" which are just evolved and propagated from boundary conditions like for eg. heat conduction through a ...
0
votes
2answers
143 views

Grand canonical ensemble with interaction, simulation doubts

First of all, this is a follow-up of my first question. The idea is the same, every state would consist of some particles in a line (with an energy associated) and the particles can not be nearer than ...
0
votes
1answer
132 views

Plane Wave expansion method

I really don't know if this is the right forum to ask the question...but please help me if you can!! I was going through the Plane Wave Expansion Method Derivation...But to be honest I could not find ...
1
vote
2answers
763 views

Acceptance probability 2D Ising Model

Disclaimer: I just found a possible solution - eventhough i don't really understand, whats wrong with my prior approach. Edit: I just tried to calculate it from scratch and found the following: $E ...
0
votes
0answers
70 views

Can peridynamics give accurate quantitative results?

I'm applying for a position that is interested in studying crack growth using peridynamics. I have experience with MD on polymeric systems, but I had never knew of this method before. It seems very ...