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
2answers
1k views

How do I integrate the Poisson equation to determine the electric potential along a particular direction (e.g., $z$)?

This question is a sequel of sorts to my earlier (resolved) question about a recent paper. In the paper, the authors performed molecular dynamics (MD) simulations of parallel-plate supercapacitors, ...
1
vote
2answers
334 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 ...
1
vote
2answers
247 views

How can one determine at which distance the Lennard-Jones potential reaches a given value?

My question is fairly simple, but I do need clarification on how to get the inverse of the Lennard-Jones potential V(x). I am working with the following expression: $$ V(x) = e\times[(R/x)^{12} ...
1
vote
3answers
724 views

Calculate kinematics of body movement from the set of spatial coordinates

Short intro I have a set of 3D (x,y,z) spatial coordinates of arm movement obtained using motion capture system. The example set of such coordinates looks like this (rounded up): ...
1
vote
1answer
70 views

Turning a finite difference equation into code (2d Schrodinger equation)

I am trying to convert the following finite difference equations into code (taken from the bottom of page 12 of this thesis by Maike Schulte Numerical Solution of the Schrodinger Equation on Unbounded ...
1
vote
1answer
59 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
1answer
116 views

Computational Fluid Dynamics methods

I have read some articles about the finite difference method on a cartesian orthogonal grid. I understand how it works when Dirichlet boundary conditions are used, or when Neumann boundary conditions ...
1
vote
1answer
142 views

Strange Behavior in Spring Computer Model

To learn more about oscillatory motion which I am learning about in my high school physics class, I have created a computer model of a damped spring where the damping force is proportional to ...
1
vote
1answer
29 views

Kernel normalization in Smoothed Particle Hydrodynamcs

I am trying to write some code for a fluid simulation. I have read about Smoothed Particle Hydrodynamics (SPH) and my question is related to the properties of the smoothing kernel. How do I set the ...
1
vote
1answer
68 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?
1
vote
2answers
65 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$ ...
1
vote
1answer
60 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 ...
1
vote
1answer
122 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 ...
1
vote
2answers
89 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
119 views

Deriving the change in the Helmholtz free energy in the context of the free energy perturbation method

I am reading Free Energy Calculations: Theory and Applications in Chemistry and Biology by Chipot and Pohorille. At the beginning of the text (page 19, for example), the authors define the Helmholtz ...
1
vote
1answer
114 views

reversible cellular automata

Let's suppose a cellular automaton has a value $b(r,t)$ belongs to $Q$ at site $r$ and time $t$, where $Q$ is the set of possible states at each site. Let $N(r, t)$ be the values of the states of all ...
1
vote
2answers
165 views

Why are Markov chains better than naive algorithm for Random picking?

I have a question about the problem of coming up with a random distribution of non-overlapping coins in a given two dimensional space. I don't understand why using the random sequential adsorption ...
1
vote
2answers
259 views

How to decide convergence of the ensemble average in a Monte Carlo process

I am simulating electromagnetic scattering off a rough surface. The usual process is to do a Monte Carlo simulation, which is briefly described as follows. Generate a randomly rough surface, and ...
1
vote
1answer
348 views

Does the universe follow a general physical and computational principle of evolution?

Evolution is a principle in biology, whereby organisms evolve their ability to replicate and multiply in numbers over successive generations. From a computational point of view, the organisms employ a ...
1
vote
1answer
541 views

Calculating Kramers-Kronig using Mathematica

First of all, I know there is a Mathematica group in beta, but I don't think the problem of the following is directly a Mathematica issue. I am trying to calculate the change of the refractive index ...
1
vote
1answer
72 views

Physical simulation in python [closed]

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
69 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 ...
1
vote
1answer
66 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 ...
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. ...
1
vote
1answer
79 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
1answer
76 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 ...
1
vote
1answer
72 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 + ...
1
vote
1answer
409 views

Magnetic field of permanent magnet [closed]

Hi I am looking for some code to (computationally) calculate the magnetic field of a rectangular permanent magnet in 3D space. Preferably I could edit it to add in my own dimensions. Does anyone know ...
1
vote
1answer
114 views

Parameter determining argon phase

Currently I am working a molecular simulation to determine phases of an argon NPT ensemble using Lennard Jones potential. Mainly I use the radial distribution function to determine solid, liquid, or ...
1
vote
1answer
279 views

Help: 3D visualization of magnetic field around moving point charge

The diagrams I'm able to find online only show the concentric field lines in the particle's plane, perpendicular to the motion, which of course generalizes to the cylindrical shape around a conductor. ...
1
vote
2answers
263 views

Where does Computer Science background students fit in Theoretical Physics [closed]

I am basically an Electronics student - background in computer science (that's where I want to work). I applied for an internship in USA in a research institute where the group is focused in ...
1
vote
1answer
251 views

Numerical algorithms to generate a random wavefunction from a thermal ensemble

I am seeking an algorithm to generate a random wavefunction = $\sum {c_i |\varphi _i\rangle }$ from a thermal ensemble, whose density matrix $\rho \sim e^{-\beta H}$, without the need to diagonalize ...
1
vote
1answer
179 views

What software can I use to simulate supersonic gas particle flow? [closed]

Please please help me out with this. I am trying to find a software/program that I could use the simulate the scattering effect of light when it strikes gas particles moving at supersonic speeds? ...
1
vote
0answers
21 views

How to compute matsubara frequency summation over computer?

Matsubara frequency sum usually takes the following form: $S_\eta = \frac{1}{\beta}\sum_{i\omega_n} g(i\omega_n).$ But in my problem, $g(i\omega_n)$ is a lengthy expression, which can not be ...
1
vote
1answer
43 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
40 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
17 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 ...
1
vote
0answers
28 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?
1
vote
0answers
38 views

XCrysDen Structure file [closed]

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 ...
1
vote
0answers
60 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
0answers
169 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 ...
1
vote
0answers
37 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 ...
1
vote
0answers
61 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 ...
1
vote
1answer
40 views

2nd - 3rd Year Physics Labs [closed]

I'm looking to propose a new physics laboratory experiment, does anyone have any particularly fond memories of an experiment they carried out whilst at University? I'd like to include a section that ...
1
vote
0answers
380 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
66 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 ...
1
vote
0answers
87 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
0answers
368 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
0answers
50 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 ...
1
vote
0answers
188 views

What is the drag coefficient of an open wedge?

To check my two dimensional CFD calculation I am looking for reference data on the drag coefficient of an open wedge. The geometry is shown below, together with the flow direction. I have found ...