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
vote
1answer
61 views

Physical simulation in python [on hold]

Is there a "standard" python package used to aid in physical simulations? What is the most popular? edit: perhaps I should have worded this question differently. Something more to the effect: Could ...
1
vote
1answer
34 views

Position based dynamics constraint scaling factor

Reading through Müller et al., Position Based Dynamics, 2007 I got lost when passing from equation (5) $$\Delta p = \frac{C(p)}{|\nabla_pC(p)|^2}\nabla_pC(p)$$ to equation (6) (and applying the ...
1
vote
0answers
38 views

Does Planck's constant imply limits to computing *results*

... I don't mean quantum effects limiting hardware fabrication sizing. Such small scales have for some time exhibited issues. Rather, along the lines of imagining the smallest possible divisions of ...
1
vote
1answer
13 views

How do I set up the tridiagonal matrix for a heat diffusion with layers of different thermal diffusivity?

I have Scala code that recreates the Crank-Nicolson solutions for the diffusion equations, and matches 'Excel for Scientists and Engineers' (Joe Billo, Wiley). However, I would like to be able to ...
0
votes
0answers
17 views

Electron scattering with CSDA simulation

I am currently in the process of writing simulation code for a "computational physics lab course" and in the present case the task is to simulate the trajectories of 20 [MeV] electrons impinging on a ...
0
votes
0answers
28 views

Examples of application of detour matrices in physics?

Are there any good examples of application of detour matrices in physics?
1
vote
1answer
64 views

A puzzle of thermalization in simulating the 3D XY-model

I am learning the classical Monte Carlo simulation. When I simulate the 3D XY-model $$ \beta H = -\beta J \sum_{<i,j>} cos(\theta_i-\theta_j) $$ where $\beta$ is the inverse of the temperature ...
2
votes
0answers
47 views

Textbooks on algorithms for the perturbative calculation of High energy physics

For the perturbative calculation of High energy physics, I have known some packages such as FeynArts, FeynCalc, MadGraph, CompHEP, GiNaC, and so on. But I am wondering whether there exists a textbook ...
1
vote
0answers
24 views

Computer coding of Perdew Burke and Ernzerho (PBE) method

I have been assigned the task to write a computer code to implement the Perdew Burke and Ernzerhof (PBE) method. Does anyone know a good reference which can make the coding of this method easier?
0
votes
1answer
30 views

Broken Thin Lens Algorithm [closed]

We all know the thin lens equation. For $o$ being a horizontal object distance and $f$ being the focal length, the horizontal image distance $i$ is described by: $$\frac{1}{f} = \frac{1}{o} + ...
0
votes
0answers
38 views

Monte carlo simulation for continuous spin model (e.g. XY or Heisenberg model)

Unlike the Ising model, the XY model and the Heisenberg model have a continuous spectrum. So one need discretize them for a numerical simulation. But how to make sure the discretization procedure ...
1
vote
3answers
109 views

Why does the Metropolis algorithm allow changes even for ∆E > 0?

In the Metropolis Monte Carlo algorithm, why can you accept changes even for ∆E > 0 (provided that a random number is less than a given probability ratio, e.g. exp(-β∆E))?
1
vote
1answer
64 views

Question about Metropolis Monte Carlo in the case of equal energies

If configuration A is equal to configuration B in a Metropolis Monte Carlo method, do you still do the attempted update?
2
votes
0answers
14 views

Coefficients and Parameters for contracted Gaussian basis sets

This is a repost from Chemistry.stackexchange in the hopes that someone here will be able to help me. Any help at all would be greatly appreciated. As far as I understand, an STO-NG contracted ...
2
votes
1answer
49 views

Unexpected eccentricity in moon orbit simulation

I've used Mathematica's NDSolve function to calculate the orbit of the Moon around the Earth. I used the following initial positions (perigee): ...
3
votes
1answer
55 views

Complex semi-definite programming

I'm doing some calculations and I want to simulate them in python or matlab (or whatever). However I use hermitian matrices and I don't really manage to find a library which enables me to calculate ...
2
votes
2answers
102 views

Modeling a potential well

I attempted to simulate the interaction of a moving particle and a potential well in Mathematica. The particle should experience a force of -$1/r^2$, if the equation for the potential well is -$1/r$. ...
1
vote
1answer
55 views

Is there a normalized form of the Euler equation discretized with finite volumes?

I want to calculate a flux on my fpga using the Euler equations with the finite volume method. Unfortunately the values of the state variables differ a lot. For example the pressure has a value of ...
1
vote
0answers
31 views

XCrysDen Structure file [on hold]

Does anyone know how to directly convert a .cif file to a Xcryden structure file(.xsf) ? I know how to extract the lattice vectors and the atom positions from a .cif file, but don't whether the .cif ...
3
votes
2answers
147 views

Boundary Element Method or Boundary Integral Method Computational Aspects

I have to solve a Helmholtz equation inside a simply connected domain. I know that in general the boundary integral can be written as, $$\phi(x)=\int_V G(x,x') \rho(x')\ d^3x'+\int_S ...
4
votes
0answers
95 views

Nature as a calculator [closed]

I wondered whether physical experiments can find the outcome of otherwise intractable mathematical problems. For example, solitons or other non-perturbative effects in QFT cannot be seen at any ...
0
votes
0answers
34 views

Poisson equation solver with specific boundary condition

I want to solve 2d Poisson equation with this Boundary conditions below $$u(-5,y) = 0 , \\ \frac{\partial u(x,y)}{ \partial x} = 0 \,\,{\rm at}\,\, x = -5\\ u(x,-5) = u(x,5)$$ Now my question is ...
2
votes
0answers
62 views

Doing numerics in physics [closed]

Soon, I am going to write my master thesis in theoretical physics. I assume there, and later on in my career, I will have to do more serious numerics than I did up to this point. That's why I want to ...
1
vote
1answer
64 views

Ground state Phase Diagram of Bose-Hubbard Model

The Hamiltonian of Bose-Hubbard model reads as $$H=-J\sum\limits_{<i,j>}b_i^{\dagger}b_j+h.c.+\frac{U}{2}\sum\limits_{i}n_i(n_i-1)-\mu\sum\limits_in_i~~~~~~~~~(1)$$ For this we plot phase ...
0
votes
0answers
19 views

Verlet timestep

I am using velocity Verlet in molecular dynamics. Is just a simple question, if a simulation using time-step femtosecond, in velocity Verlet is just necessary $dt = 10^{-15}$ to use femtosecond?
1
vote
0answers
57 views

Which are the good universities for doing PHD in computational physics? [closed]

Currently i am doing MSC(physics) from Pune university, India. And I want to do PHD in computational physics as it interests me the most. So I want to know which universities will be a good option for ...
1
vote
1answer
50 views

Calculate Euler equations of fluid dynamics without division?

I'm working on the calculation of the Euler equations with the finite volume method. Unfortunately I'm not allowed to do a division. So I'm wondering if there's a form which does not need a division. ...
2
votes
1answer
62 views

Numerical Ising Model: Swendsen–Wang algorithm, Percolation theory?

When you look at the original paper of Swendsen and Wang in 1987: "Nonuniversal critical dynamics in Monte Carlo simulations" it is somewhat mentioned that the proposed algorithm uses percolation ...
1
vote
2answers
63 views

Difference in calculated and simulated ellipsies

My task here is to determine orbit parameters, using current values: $\mu=GM$ - standard gravitational parameter $r$ - distance to the object with Mass $M$ $v$ - speed of the object in the point $r$ ...
0
votes
0answers
25 views

Computational package to find the ground state of a particle in 3D domain

I am developing a numerical algorithm to find the ground state of a Hermitian matrix. Obvious applications are quantum many-body systems and particles in various potentials. I am a little stuck with ...
3
votes
1answer
81 views

Euler Equations, Sod shock tube & conservation

Conservation of momentum? I am considering the Euler equations in conservative form and solving the Sod shock tube problem I have written a Godunov finite volume type solver. It solves for density ρ, ...
2
votes
0answers
55 views

Euler equation with single state variables

I want to simulate a flow using the Euler equations. For this reason I'm wondering if there's a modified version of the Euler equations. At the moment my formula looks like this: $$ ...
0
votes
0answers
61 views

Need help simulating a fan

I'm trying to write some code to approximately simulate a number fans in an arbitrarily shaped room. I'm not real strong in physics, and could use some help. If it helps I'm not interested in fine ...
8
votes
2answers
265 views

How do I calculate the Reynolds number in multiphase flows?

I am modeling a gas flowing through a liquid. How do I calculate the Reynolds number in multiphase flows? And, at what Reynolds number should I consider the flow to be turbulent? The problem is of a ...
3
votes
1answer
60 views

Update velocity or position first in computation?

I am trying to make a simulation of a vibrating string. The string is divided into $n$ points, and each point along the string is acted upon by a force due to the positioning of its neighbors. I ...
5
votes
1answer
213 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
0answers
35 views

How to easily calculate lengths (or relative lengths) of paths between symmetry points in BZ

I am trying to easily calculate the length between special kpoints within the BZ of the 32 point groups in a crystal system. I am calculating the lengths in order to scale k point sampling along these ...
7
votes
3answers
187 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} ...
4
votes
1answer
89 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 ...
13
votes
3answers
455 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
80 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 ...
6
votes
0answers
78 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
221 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
118 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 ...
1
vote
0answers
131 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
1answer
130 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 ...
3
votes
0answers
97 views

What is the probability of quantum tunneling occurring in this CPU?

You may have noticed over the last few years that Moore's law is no longer applying to the real world. This observation states that over the history of computing hardware, the number of transistors on ...
2
votes
1answer
97 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
52 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 ...
0
votes
1answer
60 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 ...