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)

5
votes
1answer
111 views

In a Monte Carlo NVT simulation How do I determine equilibration

I'm running an NVT (constant number of particles, volume and temperature) Monte Carlo simulation (Metropolis algorithm) of particles in two dimensions interacting via Lennard-Jonse potential ($U = ...
5
votes
0answers
142 views

Numerical Solution of the Propagation-Dispersion equation

I have asked this question on Computational Science and also on Mathoverflow, but no satisfactory answers so far. I thought maybe the physics community could shed some insight on the issue. I am ...
5
votes
0answers
658 views

Interpretation of Stiffness Matrix and Mass Matrix in Finite Element Method

I would like to have a general interpretation of the coefficients of the stiffness matrix that appears in FEM. For instance if we are solving a linear elasticity problem and we modelize the relation ...
5
votes
1answer
312 views

Ray tracing in General Relativity

I would like to find out what one would see at the Schwarzschild radius of a massive non-rotating black hole, if the black hole is surrounded by a bright ring. For that, I would place the observer at ...
4
votes
3answers
368 views

Future Computer Performance

Moore's law has succesfully predicted up to now that integrated circuit transister density doubles every two years. However, computer performance is dependent on additional factors like architecture, ...
4
votes
3answers
1k views

Nanorobots. What stops us from producing them yet?

If we can already predicts accuratelly motion on molecular levels, what stops us from developing small robots to, for instance, navigate through our blood vessels looking for cancerous cells and ...
4
votes
1answer
301 views

What states are satisfying an entropic area law and why do they satisfy it? More specificly why do matrix product states satisfy it?

I am currently reading some papers concerning the question why the density matrix renormalization group (DMRG) method is working well for simulating one dimensional systems and bad for higher ...
4
votes
1answer
828 views

Industry application of computational quantum mechanics?

I was wondering if anybody knew of an industry application of computational quantum mechanics. For example, the efficient placement of circuit elements on a PCB is in part motivated by classical FDTD ...
4
votes
1answer
105 views

buiding circuits from color superconductors

caveat: the sort of exotic matter engineering in here is currently beyond the reach of our technology, but, that having been said: Has their been any research on building models of these sorts of ...
4
votes
2answers
1k views

Radial Schrodinger equation with inverse power law potential

Recently I read a paper about solving radial Schrodinger equation with inverse power law potential. Consider the radial Schrodinger equation(simply set $\mu=\hbar=1$): ...
4
votes
3answers
327 views

Good source for numerical simulations of Wigner function?

I'm interested in simulating the time evolution of a Wigner function for a harmonic oscillator (and possibly some other potentials) and I can't seem to find a good resource for that. My background in ...
4
votes
1answer
188 views

Using the monte carlo method to compute the magnetic field of a curent carrying loop

I have written a program in cpp that computes the magnetic field at a point from a current carrying loop. It uses the biot savart law and the monte carlo technique to carry out the integral. The ...
4
votes
2answers
931 views

Can I use imaginary time propagation for many-body problems?

There are various ways to numerically find the ground state energy and wavefunction of a many-body Hamiltonian. You can diagonalize the Hamiltonian and pick out the lowest eigenstate, or you use ...
4
votes
2answers
2k views

On numerically solving the Schrödinger equation

I just read a paper 'A pocket calculator determination of energy eigenvalues' by J Killingbeck (1979). Link: http://iopscience.iop.org/0305-4470/10/6/001 I have some questions about section 2. Why ...
4
votes
1answer
473 views

Integral of Sedov's self-similar solution to the spherical blast wave problem

I'm studying the Taylor-Sedov self-similar solution to the problem of a strong explosion in a homogenoeus atmosphere. The problem is discussed in Landau & Lifschitz VI (in the 2nd edition it's ...
4
votes
1answer
528 views

Numeric method to calculate the charge distribution on a conducting surface?

If I have an arbitrary (closed?) conducting surface and a nearby charge density, is there a simple numeric way of computing the induced charge distribution on the surface?
4
votes
1answer
488 views

Philae lander simulation off by factor of ~3

I'm trying to simulate the Philae landing by writing a program to compute the position of the lander vs time. According to various mission websites, the orbiter will match its orbit to the rotation of ...
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 ...
4
votes
1answer
420 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 ...
4
votes
0answers
36 views

How to calculate dispersion relation from a Finite Difference (FD) wave simulation

I have a python code that calculates the solution of the inhomogeneous acoustic wave equation for a 2D medium with any velocity and source configuration. It was implemented using Finite Differences ...
4
votes
0answers
106 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 ...
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} ...
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 ...
3
votes
3answers
408 views

Can a wavefunction be solved to any arbitrary precision, given enough computer time?

I learned that the wavefunction for the hydrogen atom can be solved analytically (we did the derivation in class), but that for more complicated atoms it is "impossible" to solve and that only ...
3
votes
2answers
281 views

Rotation matrix of Euler's equations of rotation relative to inertial reference frame

I was playing with simulation of Euler's equations of rotation in MATLAB, $$ I_1\dot{\omega}_1 + (I_3 - I_2)\omega_2\omega_3 = M_1, $$ $$ I_2\dot{\omega}_2 + (I_1 - I_3)\omega_3\omega_1 = M_2, $$ ...
3
votes
3answers
898 views

What's the difference between “numerical methods” & “mathematical analysis” as said by Feynman in his lectures?

While reading his lectures, I came to these lines: On the basis of Newton's second law of motion,which gives the relation between the acceleration of any body & the force acting on it,any ...
3
votes
2answers
147 views

Square root of a matrix appears in massive gravity. How to solve $\sqrt{A+B}$ perturbatively

$A=\text{diag}\{\lambda_1,...,\lambda_n\}$, where $\lambda_i$ can be any number and not necessarily a small number, $\lambda_i>0$, $B$ is a positive definite symmetric matrix, and ...
3
votes
2answers
93 views

Why is my two body model not working? [closed]

I have created a simple 1 Dimensional two body particle model in C++. In the model, particle 1 starts at position (0,0) and particle 2 starts at position (1,0). The particles are accelerated towards ...
3
votes
2answers
355 views

Mixed conductive and convective heat transfer

I want to make a finite element analysis of a cold airflow through warmer pipes. In particular I want to see how the pipes cool down and the air heats up, as it travels through the pipes. Wich are the ...
3
votes
2answers
250 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 ...
3
votes
1answer
255 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 ρ, ...
3
votes
2answers
200 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 ...
3
votes
2answers
185 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 ...
3
votes
1answer
52k views

How does force relate to velocity

I had originally asked this question on math overflow and it was suggested that I ask it here. So I know that a force will change the magnitude of velocity if it is at an angle other that 90 degrees. ...
3
votes
2answers
99 views

Numerical modelling of a step function in time in a hydrodynamic system. (Runge Kutta fourth order)

So I'm trying to model a hydrodynamic system that introduces a sudden "jump" in the value of a function at a specific time. The system is solved with a Runge-Kutta fourth order method. I have a ...
3
votes
1answer
245 views

Is there an algorithm for N body simulations in General Relativity [duplicate]

I am new to general relativity but have a background in computer science. Why is it so hard to do n-body simulations in GR? For example, there could be a program which takes the properties (mass, ...
3
votes
1answer
116 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 ...
3
votes
3answers
528 views

Adding rotation to internal coordinates

I'm optimizing the geometry of a system composed of several interacting masses (a molecule). The energy of the system depends on the relative position of the masses, and all velocities are zero. In a ...
3
votes
1answer
980 views

Some hints/introductions/textbooks for LBM (Lattice Boltzmann methods) fluid simulation?

I come from a computer science background and I want to learn details about LBM in fluid simulation. I have been searching around for a concise tutorial that tells me how LBM is used to solve a ...
3
votes
1answer
702 views

Numerical computation of the Rayleigh-Lamb curves

The Rayleigh-Lamb equations: $$\frac{\tan (pd)}{\tan (qd)}=-\left[\frac{4k^2pq}{\left(k^2-q^2\right)^2}\right]^{\pm 1}$$ (two equations, one with the +1 exponent and the other with the -1 exponent) ...
3
votes
1answer
75 views

Numerically solving a simple Schrodinger equation with fast Fourier transforms

While trying to solve a stochastic Gross-Piaevskii equation I have found a problem that can be tracked down to something buggy occurring in the simplest Schrodinger equation possible: $$\partial_t ...
3
votes
2answers
2k views

Numerical solution to Schrödinger equation - eigenvalues

This is my first question on here. I'm trying to numerically solve the Schrödinger equation for the Woods-Saxon Potential and find the energy eigenvalues and eigenfunctions but I am confused about how ...
3
votes
1answer
125 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
154 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
123 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
93 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
276 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
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 ...
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 ...
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, ...