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)

3
votes
1answer
147 views

Calculating a two-dimensional orbital path with infinite granularity (non-Euler integration)

For a game I am making, I am trying to calculate the position of an orbiting object around one or more bodies. I have successfully implemented this gravity simulation by calculating the force, then ...
3
votes
1answer
159 views

Looking for a good introductory-level review of pseudopotential methods

I'm looking for a good introductory-level review of pseudopotential methods. In particular, I'd like to understand how the self-consistent pseudopotential methods work.
3
votes
1answer
186 views

Calculation of long-range forces in Molecular Dynamics - Ewald summation

I am trying to write a code to calculate the potential and forces, for the same using ewald summation.For this purpose, the formula for potential and force I have used is : $$ U = U^{(r)} + U^{(k)} ...
3
votes
1answer
101 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 ...
3
votes
1answer
311 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
1answer
149 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 ...
3
votes
1answer
250 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 ...
3
votes
1answer
131 views

How to determine predicted CP violation for a given SUSY point?

I'm currently studying at the spectra of some supersymmetric models, and would like to know whether the parameter points I'm looking at are ruled out due to excessive CP violation. I am using SPheno, ...
3
votes
1answer
246 views

Using Relaxation Method to Model Negative Dielectrics in an electric field?

How can you use the relaxation method to model negative dielectrics? The relaxation method is usually used to model electrostatics problems but negative dielectrics are only see in dynamic systems. ...
3
votes
0answers
73 views

Are there resources for simulating and/or theoretically describing solitons?

Recently there are striking new ideas emerging on "lower level" dynamics with respect to quantum mechanics involving fluid mechanics principles, including hints of soliton-like aspects to particle ...
3
votes
0answers
139 views

black body simulation

black body radiation is typically understood from Planck's argument of light resonance in a box, from which the density of states is computed. Now, suppose I want to simulate a black body ...
3
votes
0answers
71 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 ...
3
votes
0answers
92 views

Non-linear Wave Equation - Numerical Methods

Motivation: I'm working with a highly non-linear spherical wave-like equation (second order PDE). The equation can be written on the form $$\ddot{u} = f(t, \dot{u},\dot{u}',u',u'')$$ where ...
3
votes
1answer
196 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
0answers
149 views

Difference between a Fixed Point and a Limit Point in implementations of the Renormalization Group (RNG) in Large Eddy Simulation (LES) model

In the introduction of this paper, it is explained that and how the application of a dynamic subrid scale model for turbulence into a large eddy simulation (LES) model corresponds to doing one ...
3
votes
0answers
968 views

Comparison of different ab-initio codes

One may find on the web a lot of different computational packages to perform "ab-initio" calculations of electron structure of the solids. Usually, the documentation is not quite transparent about the ...
3
votes
2answers
123 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 ...
2
votes
4answers
319 views

Computer Science Modeling of Physical World

I am curious what efforts have been made to date to define virtual computer worlds based on the physics we know in the real world? I think it would be awesome to say start off with an atom defining ...
2
votes
3answers
349 views

Solving the two body problem numerically

I'm trying to solve the two body problem numerically, setting up $G$, $m1$ and $m2$ to be equal to 1. then I located each mass on positions -5 and 5 respectively along the $x$ axis and gave them both ...
2
votes
3answers
94 views

Is it possible to reconstruct the wavefunction of a molecule from a collection of spectra?

Spectra of a molecule can be calculated if the wavefunction is known. Is it possible to do the opposite?
2
votes
3answers
172 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 ...
2
votes
1answer
52 views

Quantification of uncertainty, or information loss, when modeling a physical system?

Say I have a physical circuit (say a battery, wires and an incandescent lamp) that I want to model. To model it I could perhaps use: traditional circuit theory some kind of finite element analysis ...
2
votes
1answer
33 views

Testing the accuracy of a solar system simulation

I've written a java program to model the solar system, and I need to test its accuracy. What would be some good ways to do that? I've checked for conservation of energy and angular momentum. I'm using ...
2
votes
2answers
251 views

Hamiltonian mechanics really useful for numerical integration? Lagrangian can become 1st-order

(I'm talking about the classical mechanics.) Many texts say that Euler-Lagrange equations are difficult to treat numerically because they are second-order ODEs, ${f_i(\boldsymbol{q, \dot{q}, ...
2
votes
2answers
213 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 ...
2
votes
2answers
140 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$. ...
2
votes
2answers
205 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 ...
2
votes
2answers
70 views

time step in kinetic montecarlo simulation

I have a doubt about how the time step is calculated in kinetic montecarlo simulations. One state with index $i$ is connected to other $N$ states, indexed by $j=1...N$, by transitions that happen ...
2
votes
2answers
117 views

Exact Diagonalization of a BdG Hamiltonian on a Finite Lattice

I would like to numerically find the edge modes of a $p_x$ + $i p_y$ BdG Hamiltonian. The lattice version is given by H = $\sum\left[-t \left(c_{m+1,n}^{\dagger} c_{m,n} + \text{h.c} \right) - ...
2
votes
3answers
207 views

Textbook recommendation for computational physics

Can anyone recommend some good textbooks for undergraduate level computational physics course? I am using numerical recipe but find it not a very good textbook.
2
votes
2answers
437 views

Why can't we set the lattice spacing 'a' in lattice QCD?

My question is to do with lattice QCD. In the lattice action there is a parameter, 'a', the lattice spacing in physical units. However, if we want to generate a configuration with a certain lattice ...
2
votes
1answer
165 views

Finding the Riemann tensor for the surface of a sphere with sympy.diffgeom

I have implemented a SymPy program that can calculate the Riemann curvature tensor for a given curve element. However, I am encountering problems solving for the case when the curve element is the ...
2
votes
1answer
127 views

9-point stencil “equivalent” for advection equation

So I inherited from some people a code that solves the advection-diffusion-reaction equation for a particular system. The original code was first implemented in 1D which worked fine in cartesian ...
2
votes
1answer
95 views

Damping a spring force

I'm modelling particles in a system using a spring and damper force. $$F= kx -cv$$ $x=x_i -x_f$, where $x_i$=centre of spring and $x_f$= displaced position. Above $x$ is the displacement and $v$ is ...
2
votes
1answer
690 views

Tight binding model in a magnetic field

The standard way to treat a tight binding method in a magnetic is to replace the hopping matrix element: $t_{i,j}\rightarrow e^{i\int_i^j\mathbf{A(x)}.d\mathbf{x}}$ the so called "Peierls ...
2
votes
2answers
1k views

Gas viscosity at high pressure, high temperature

EDIT 1 PER COMMENTS I am wanting to model nitrogen gas viscosity as a function of pressure and temperature OR learn of an existing equation that models nitrogen viscosity for the pressure and ...
2
votes
1answer
82 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
2answers
1k views

Lennard Jones liquid vapour phase transition

I am trying to generate a phase diagram of pressure versus temperature in a NPT ensemble using molecular dynamics method. I am using the Lennard Jones potential with 108 particles and Berendsen ...
2
votes
2answers
419 views

Computational Science involve programming? [closed]

I read what is computational science in Wikipedia but the explanation and understanding are not very clear. So, I could you please give a simple example computational science project and what all ...
2
votes
1answer
720 views

How is acceleration computed in the leap-frog algorithm in computer simulations of particles?

I am reading the book Computer Simulation of Liquids by Allen and Tildesley (here is another link). On page 80, the authors describe the leap-frog algorithm, which is used extensively in molecular ...
2
votes
1answer
160 views

If the graph isomorphism problem can't be solved in polynomial time, do spin-networks in loop quantum gravity violate the Church-Turing thesis?

If the graph isomorphism problem can't be solved in polynomial time, do spin-networks in loop quantum gravity violate the Church-Turing thesis? To determine whether or not two graphs are isomorphic ...
2
votes
1answer
44 views

Current Status of the Monte Carlo Sign Problem

I've been reading about the Monte Carlo sign problem, and I am a little confused about its current status. Specifically, after reading this post When is the "minus sign problem" in quantum ...
2
votes
1answer
56 views

Difference between finite volume, characteristic method and plug flow models of a pipe

I have to model pipes (of a district heating network) with ODE's. My background is computer-science, so it is not that easy for me to understand different approaches. Finite volume approach. Method ...
2
votes
1answer
175 views

Problems implementing Lattice Boltzmann method

I want to make a basic 2D implementation of the Lattice Boltzmann method (in javascript) to simulate gas/fluid flow, but I have run into some problems and I can't seem to find the cause of it. ...
2
votes
2answers
227 views

Why does GW-DFT give higher bandgaps in semiconductors

Usually the GW Density Functional Theory (DFT) gives larger band gaps in semiconductors compared to the LDA and GGA methods. This seems to be related to the screened potential in GW, but it is not ...
2
votes
2answers
94 views

Is there any untamed physics in the CPU?

In a x86/amd64 CPU could I find an instruction or a combination of instructions that leak inside the system? Could some other physics properties (not accounted by the design) be measured or interacted ...
2
votes
1answer
163 views

Density Matrix Renormalization Group (DMRG) Simulation of a String-Net Model

In the following paper, Dr. Xiao Gang-Wen et. al. introduce the idea that string-net condensed states can be represented in terms of tensor product states: http://arxiv.org/pdf/0809.2821.pdf The ...
2
votes
1answer
717 views

How to measure the spin-spin correlation in a Monte Carlo simulation of the Ising model?

I'm simulating the Ising Model in 2D up to 5D and I want to calculate the spin-spin correlation, correlation length, and critical exponent of the system. What is a good way to go about doing this? ...
2
votes
1answer
74 views

Estimating divergence of set of vectors

I have a set of points where directions and intensities of a flow are given (in 3D). Is it possible to estimate the divergence of the flow defined by those vectors? I only need a rough estimate and I ...
2
votes
1answer
301 views

Simulation of fluid flow using Euler equation

I have been looking on Euler's equations for a while and can't grasp one thing. Suppose we have initial system state with volumes of fluid "hanging" in air (time is frozen and equal to zero), each of ...