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)

2
votes
2answers
96 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
0answers
352 views

Precession of Mercury (Python simulation)

I was trying to simulate the precession of Mercury based on the perturbed solution, and my questions about its implementation in python can be seen here: http://scicomp.stackexchange.com/questions/...
6
votes
5answers
464 views

Resolving General relativity and Newtonian mechanics to a computer [closed]

I know this is considered an old subject long ridiculed by many as the folly of layman. But I work in the field of computer simulation, specifically in producing fully functional 3D interactive ...
21
votes
3answers
3k views

Programming in physics [closed]

How important is programming in physics? I am studying physics at university and these first years there is actually no approach (as it is understandable) to what working in physics is like. Now, I ...
0
votes
1answer
72 views

What exactly is closed orbit theory and what assumptions go into it?

I am just beginning a research project on the study of closed orbits, specifically as related to hydrogen, and I wanted a little bit more information on what exactly makes something a "closed orbit". ...
1
vote
0answers
47 views

Is there a numerical calculation for black hole - neutron star merges?

I think it might be even a more frequent event as the black hole merges. It would be similar to the black hole merges in the gravitational wave spectrum, but it could have a very clear neutrino and ...
1
vote
1answer
93 views

Meaning of imaginary eigenvalues in advection equation

I'm reading the book Fluid Simulation for Computer Graphics (Amazon link), and stumped by the following sentence in Chapter 3: what's happening is that the eigenvalues of the Jacobian generated ...
0
votes
2answers
835 views

Correlation energy - is it the difference between the Hartree-Fock energy and exact energy, or Hartree-Fock PE and exact PE?

For some reason I can't find anything stating it either way explicitly. What I'm talking about is this. Is this difference referring to potential energies or just energies in general? I assume it's ...
0
votes
2answers
183 views

What is the fastest algorithm to solve a Nonlinear second order differential equation numerically

I am trying to solve a second order non linear differential equation in one variable. Using RK4 I am getting good accuracy and is working fine. But the problem is my range is very high so it will ...
3
votes
2answers
114 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 ...
0
votes
1answer
112 views

Force Needed To Maintain Velocity, Increasing Mass

Hoping somebody could give me a hand with this. Unfortunately I haven't taken physics yet, my background is in math/computer science. I'm studying computational science and the problem is essentially ...
2
votes
1answer
184 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 ...
1
vote
0answers
101 views

Mixed spin Ising Model

As we know ferrimagnets can be modeled by the Ising model. I came across this equation in "Compensation Temperature of the Mixed-Spin Ising Model on the Hexagonal Lattice" by W. Figueiredo, M. Godoy, ...
1
vote
1answer
318 views

Determing Velocity of Moons

I have a question that I believe is relatively easy to answer, I am working on an $N$-body simulation of a fictional star system and am having trouble finding the velocity of moons so that they will ...
3
votes
1answer
508 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, ...
2
votes
1answer
801 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 ...
0
votes
1answer
97 views

How to solve highly oscillating differential equation [closed]

The equation looks like: $$x''(t)+bx'(t)+c x(t)+dx^3(t)=0.$$ This is the motion of a particle in a potential $cx^2/2+dx^4/4$ with friction force $bx'$. In my case, the friction term is very small and ...
0
votes
0answers
68 views

What is the largest atom with a reliable configuration-interaction (CI) calculation?

The simplest approximation for calculating the ground state of an atom is the Hartree-Fock approximation. To get accurate result for the ground state energy, one has to do configuration-interaction (...
0
votes
0answers
51 views

Q: Any standard terminology for “gold standard” Monte Carlo solutions

The "standard" way in our field to compute scattering statistics in wave propagation problems is via a "deterministic" wave prop code operating in a Monte Carlo context where we iterate over many ...
0
votes
1answer
121 views

Connections between Density Matrix Renormalization Group and Conformal Field Theory

Can we use the density matrix renormalization group (DMRG) method to understand problems in conformal field theory? I have been trying to find some connections, but nothing is coming up when I search. ...
1
vote
2answers
114 views

How to verify/falsify the existence of localised edge states numerically?

I have to consider a Hamiltonian given in second quantized form in real space $$H = \sum c_k^\dagger h_{kl} c_l \, ,$$ describing fermions on a 2d hypercubic lattice. The concrete form of the matrix ...
0
votes
0answers
159 views

How to solve the Dirac equation numerically?

The effective Hamiltonian for my system is: $$ H=\nu_{F} \sigma\cdot\left(q-By\hat x\right) $$ where $\sigma$ and $q$ are the Pauli matrices and the momentum operator respectively and $\nu_{F}$ and $B$...
1
vote
2answers
146 views

Anyonic braiding statistics from density matrix renormalization group (DMRG) simulations

How does the ground state energy of the system change when we braid two anyons? Can the braiding of anyons be simulated with a computational method such as the density matrix renormalization group, ...
0
votes
1answer
53 views

Amount of unknown parameters in compressible Euler equations

I'm looking at this page for the compressible Euler equations. To me it seems, in the 1-dimensional case, there should be 3 unknowns: density, velocity, and pressure. This is because the energy $E$ ...
-6
votes
1answer
187 views

Computational physics using mathematica [closed]

So I was confused about this question on how to exactly begin to answer it. I am a novice in mathematica and I am teaching myself thus I require help in this question. From what I think I should do, ...
1
vote
1answer
369 views

Finding the initial velocity vector of an orbiting body

I'm writing a program that simulates Newton's law of universal gravitation by simply calculating the force and applying it on the objects. The simulation works very well, but now I want to simulate ...
1
vote
1answer
68 views

Formulating a symplectic integrator for a non-local Hamiltonian

I recently asked two questions, Q. [1] and Q. [2], regarding reformulating non-local Lagrangians as Hamiltonians. In these questions, the Hamiltonian is formulated as an integral because of it's non-...
0
votes
1answer
120 views

Using the Fourier transform to find the natural frequencies of coupled oscillators

How can I find the natural frequencies of a system consisting of a pair of coupled oscillators using Fourier transforms? The System consists of two masses and three springs. One of the springs ...
1
vote
0answers
43 views

Introduction material for lattice field simulation

I am looking for introductory material of lattice field theory simulation. It is better start from simple example (e.g. \phi^4) and include some source code. Is there any on-line or book which is ...
4
votes
2answers
3k 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 ...
1
vote
1answer
59 views

Numerical construction of phase space for a dynamical system

Suppose I have a standard, deterministic dynamical system. For concreteness I'll assume it's a two variable system of the form, $$ \dot x_1 = f(x_1,x_2; \theta_1)\\ \dot x_2 = g(x_1,x_2; \theta_2) $$ ...
2
votes
1answer
818 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? ...
1
vote
0answers
120 views

Ising model Monte Carlo simulations in 4D and 5D

I'm going to be simulating the Ising Model in 4D and above to calculate spin-spin correlations and critical exponents and am wondering how to tackle this algorithmically. For example, in 1D, use an ...
0
votes
1answer
29 views

Is it realistic for soundwaves under water to “sink” or “float”?

I'm studying soundwaves under water and I had a numerical problem that I was asking about. http://stackoverflow.com/questions/28904017/are-my-matlab-iterations-working Now I wonder if you can tell ...
1
vote
0answers
144 views

Numerical Solution of 1D Boltzmann Transport Equation

I need to solve the one-dimensional Boltzmann transport equation in a semiconductor numerically, and I want to take a deterministic approach toward the problem (i.e. not use Monte-Carlo or similar non-...
2
votes
1answer
78 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 ...
0
votes
0answers
39 views

What does invariant exactly mean and how does it get the invariant?

I have read many journal about simulation of regularized long wave. In numerical test section,many researcher use invariant of mass,momentum and energy to check accuracy of their method but i found ...
2
votes
2answers
270 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}, \ddot{q}}...
1
vote
1answer
138 views

Hartree Fock equations

I don't understand how the Hartree Fock equations define an iterative method! For this discussion, I am referring to the HF equations as described here: click me! Basically if you guess a bunch of ...
0
votes
1answer
229 views

Probability density function of a particle for computation [closed]

I'm writing a program, part of which relies on a particle being able to change location similar to a how a real particle would behave (pardon my physics). For example, on a grid of 100x100, a ...
0
votes
0answers
52 views

Why are functional representations of systems important in physics or computational physics?

This was an addendum to a previous question I asked, but I figured I should make it it's own discussion. Assuming I am able derive a functional representation for any dynamical system (dissipative, ...
1
vote
2answers
312 views

Leapfrog method in Particle-in-cell

Recently I wanted to write a 3D electromagnetic Particle-in-cell code with c++.I know that I should use leapfrog method.For example,when I calculates the position and velocity of particles,i should ...