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)

1
vote
0answers
42 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 ...
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 ...
0
votes
0answers
73 views

Ising Model 2D Correlation Length

I'm using Metropolis and Wolff Clustering algorithms to estimate the spin-spin correlation function $$<s_os_r>$$ I know that this is related to the correlation length but how do we determine ...
1
vote
1answer
51 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
425 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
97 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
28 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
124 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 ...
2
votes
1answer
73 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
38 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
215 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}, ...
1
vote
1answer
129 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
0answers
31 views

What is the transfer function in fft beam propagation for unpolarized light?

What is the transfer function in fft beam propagation for unpolarized light ? How to construct the fft beam propagation ? This is for homework. For coherent light the beam propagation is E(x,z) ...
0
votes
1answer
174 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
50 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
245 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 ...
3
votes
0answers
70 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 ...
0
votes
0answers
15 views

What is the phase-amplitude numerical method?

What is the phase-amplitude numerical method? I heard its used to calculate long range interactions numerically, but I cannot find any papers discussing its method of implementation.
0
votes
0answers
125 views

Stability of Mathieu's equation and parameteric resonance

I am given the following equation (Mathieu's equation) in my subject of Numerical Analysis : $$ \frac{d^2 x}{dt^2}=-\omega^2(1+\epsilon\cos(t))x $$ I am supposed to find those frequencies $\omega$ ...
2
votes
3answers
91 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?
1
vote
0answers
274 views

Problems while numerically computing band structure using k.p theory

I want to use k.p theory to numerically compute the band structure of a bulk semiconductor. The band I like to include are the lowest conduction band (cb), the heavy-hole (hh), the light-hole (lh) and ...
0
votes
0answers
91 views

Differential Equations for Two-Dimensional N-Body Simulation

So I recently asked a question about improving the stability of a two-dimensional orbital simulation (as before I was using Euler integration). I was told to use Runge-Kutta 4 for a more stable ...
2
votes
1answer
113 views

Hard processes in particle collisions

This is in reference to a MC generator called Pythia, but has wider applicability because the phrases are used widely in particle phyiscs, even though I have not managed to find a clear, ...
2
votes
1answer
52 views

Cellular automata and time T

Reading around about various CA models I saw that there are even those who are following the track to provide a single mechanicist model of the universe. For example Gerard van 't Hooft claims that ...
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 ...
4
votes
1answer
476 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 ...
0
votes
1answer
90 views

How do we define capillary number in 2D (two dimensions)

In 3D the capillary number $Ca$ is defined as: $$Ca=\frac{\nu \rho U}{\gamma}$$ where $\nu$ is the kinematic viscosity ($m^2/s$), $U$ is the velocity ($m/s$)and $\gamma$ is the interfacial tension ...
0
votes
0answers
24 views

Why breakdown voltage decreases with increasing FWHM of an asperity?

I am running a simulation for a dielectric breakdown. I am getting breakdown voltage decreases with FWHM, which does not make sense as field enhancement is higher for a small width. Simulations solve ...
1
vote
1answer
37 views

Effects of terrain on outdoor sound propagation?

I'm looking for a computationally fast model for capturing the effects of varying terrain topology in far-field outdoor acoustic propagation, and was wondering if anyone had any insight or could point ...
2
votes
3answers
183 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.
0
votes
1answer
51 views

2D room with electron-point-charges (stability)

Is there a fast way to find out stable #"point charges" & constellation for this kind of Simulation. The 2D repeating itselve in each direction till infinity, so i've a kind of symetrie here. ...
2
votes
0answers
59 views

How can the fictitious mass in the Car-Parrinello method reproduce the “real” dynamics?

In the Car-Parrinello method, to solve simultaneously the classical equations of motion for the atoms and the Kohn-Sham equations for the electrons, the following effective Lagrangian is used: $$ ...
1
vote
1answer
86 views

Molecular dynamics and detailed balance

In developing methods to perform Monte Carlo simulations one sufficient condition to preserve the stationarity of the target probability distribution is to impose detailed balance i.e. [Gardiner page ...
0
votes
2answers
514 views

How to determine velocity vector direction with respect to acceleration.

I'm currently writing a program that attempts to simulate particle movement in a gravitational field with more than one object exerting a force on it. I decided that I'd have the particle move by ...
-1
votes
1answer
82 views

Possibility of simulating quantum physics and relativity [closed]

I must be clear; this is not a 'specific' question, but rather a hypothetical one. As of now, computer physics simulations are becoming more and more advanced, and I am programming one now. I was ...
3
votes
0answers
133 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 ...
1
vote
1answer
240 views

Numerical solution for Friedmann equations

My problem today is to solve the Friedmann equations, for those who aren't familiar with them, here they are (in my specific case): $$ \left ( \frac{\dot{a}}{a^2} \right )^2 = \frac{\rho_1}{a^4} - ...
1
vote
1answer
265 views

Can I use Runge-Kutta to solve these equations?

Edit: I'm going to give some more background and derivation to show how I got to these equations. I am basically following the derivation that is found in the appendix of the following paper: ...
6
votes
1answer
306 views

The use of Artificial Intelligence in physics resarch [closed]

I have been reading about a machine that observed a double pendulum and created equations that both described its motion and associated conservation laws. The authors claims: We have developed a ...
0
votes
0answers
154 views

Boundary conditions of stream function

I have to do an problem about solving numerically the flow that goes under an airfoil. The airfoil has a flap deployed downwards and I need to solve the mesh that it's under the airfoil. I have drawn ...
3
votes
3answers
914 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 ...
0
votes
1answer
419 views

My python trajectory calculator for a small mass in a gravitational field - why is energy conservation violated?

I wrote this thing to calculate the next position of a small mass, given its current position and velocity, and a vector field determines acceleration. It sort of works - except I can never get the ...
2
votes
1answer
49 views

Mass scaling in molecular dynamics

I've noticed some authors scale the mass of particles in molecular dynamics simulations while leaving the force field parameters the same in order to achieve materials of different densities. Does ...
1
vote
1answer
120 views

Heat equation with heat radiation and heat transfer

If I want to calculate steady temperature distribution on a one-dimensional stick, and I need to consider both the heat radiation and heat transfer, then my equation will be in the form: $$ ...
1
vote
1answer
139 views

Thermal conduction through 1D rod with changing thermal diffusivity?

I've been working with solving partial differential equations in 1D numerically using the finite difference scheme. My case is for example this: I've got a rod that's insulated around the sides such ...
8
votes
2answers
354 views

Arrhenius Fit: Linear or exponential form?

I have a seemingly easy question about performing an Arrhenius fit to the equation $$y = A \times \exp \left( -\frac{E_A}{RT} \right)$$ I can either fit this in the exponential form using a ...
1
vote
1answer
672 views

Why is the ground state energy of the Heisenberg XXZ Model unbounded for some values of $J$?

At the moment, I'm looking at numerically studying the Heisenberg XXZ model. The Hamiltonian is given below: $$ H=\sum_{j=1}^{N-1}\left(J S_j^z ...
2
votes
1answer
69 views

FDM versus fdtd for computational electromagnetics

Every article I've seen on computational electromagnetics uses the Yee fdtd algorithm, and I'm wondering if would be possible to do the same analyses using out of the box FDM with more data points, ...
0
votes
1answer
154 views

Ideas for physics-related simulations to code up as weekend (or monthly) projects [closed]

What are some good ideas for simulations that I could try to implement/code in my spare time just for educational purposes? I'm thinking that it would be a good way for me to delve into the underlying ...
5
votes
1answer
152 views

How close to the critical point is sufficient close for measuring critical exponents?

I am learning Monte Carlo and just manage to simulate a phase transition by computing the heat capacity or the susceptibility. I wish I can also compute critical exponents.To this purpose, I have read ...