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
0answers
45 views

Numerical problem with Hartree-Fock equations for dilute Bose gas

I have to solve the following set of equations self-consistently: $$\begin{align} n_c(\mathbf{r}) & = \frac{1}{g}\left[\mu - V_{\rm ext}(\mathbf{r}) - 2 g n_{T}(\mathbf{r}) \right] \\[3mm] ...
2
votes
0answers
36 views

Is the finite element method easily extrapolable to new disciplines? [closed]

I am a structural engineer with a good solid mechanics background. I am reasonably familiar with finite element methods as applied to linear/nonlinear solids and structures, I would like to know how ...
2
votes
0answers
35 views

Comparing two versions of the same hydrodynamic code and their error

So I have two versions of a hydrodynamic code that has the same underlying physics. Lets call them code A and B. However code B is more optimized and more object oriented. I was trying to compare the ...
2
votes
0answers
45 views

How to let Pythia to output the complete event record in LHEF or StDHep format?

I am currently using Pythia 6.4 to simulate some processes (I do not have the intention to upgrade to 8). Now, I need Pythia to output its event record in the format of LHEF or StDHeP. I accidently ...
2
votes
0answers
243 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: ...
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: $$ ...
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 ...
2
votes
0answers
83 views

Modelling gravitational potential of a galaxy

I am interested in modelling the gravitational potential of a disc-shaped galaxy with radius $R$, i.e. solving the 2D Poisson equation numerically by Gauss-Seidel relaxation: $$\nabla^2 \phi = 4\pi G ...
2
votes
0answers
210 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 ...
2
votes
0answers
37 views

Coefficients and Parameters for contracted Gaussian basis sets

This is a repost from Chemistry.stackexchange in the hopes that someone here will be able to help me. Any help at all would be greatly appreciated. As far as I understand, an STO-NG contracted ...
2
votes
2answers
210 views

How do you measure numerically the central charge of a system?

Let's say that you are doing some Monte-Carlo simulations of a statistical system on a lattice and you observe scale invariance, meaning that you are at a conformal point. Can you get a numerical ...
2
votes
1answer
345 views

Unexpected eccentricity in moon orbit simulation

I've used Mathematica's NDSolve function to calculate the orbit of the Moon around the Earth. I used the following initial positions (perigee): ...
2
votes
0answers
106 views

Doing numerics in physics [closed]

Soon, I am going to write my master thesis in theoretical physics. I assume there, and later on in my career, I will have to do more serious numerics than I did up to this point. That's why I want to ...
2
votes
0answers
63 views

Euler equation with single state variables [duplicate]

I want to simulate a flow using the Euler equations. For this reason I'm wondering if there's a modified version of the Euler equations. At the moment my formula looks like this: $$ ...
2
votes
1answer
262 views

Exact diagonalization to resolve ground state degeneracies

I am studying a perturbed Toric Code model that is not analytically solvable. On a torus the ground state degeneracy of the unperturbed model is 4. Once we turn on the perturbation there is a change ...
2
votes
1answer
175 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, ...
2
votes
0answers
400 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 ...
2
votes
0answers
435 views

simple MHD simulation code for (self) education and play with [closed]

I would like some super simple computational code for solving magnetohydrodynamics problems. High accuracy nor performance is not my concern. I wan't it just to visually explore qualitative behavior ...
2
votes
0answers
63 views

Trotter splitting and entanglement entropy

I have heard that a numerical solution to the Schrodinger equation using the Trotter splitting formula for a many-body Hamiltonian can cause an artificial increase in the entanglement entropy. I was ...
2
votes
2answers
268 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 ...
2
votes
0answers
381 views

fun computational physics projects [closed]

I want to learn how to code in Java and am in need of some cool computational physics projects that would provide motivation to do so. I was wondering if you guys could list some medium-level ...
1
vote
1answer
1k views

Why do people use C++ in Computational Physics? [closed]

Once day, I talked to my Prof. about Computational Physics and he showed me his work. I'm quite surprised that why Prof. use C++ for computation because C++ quite complicated. I heard he talked about ...
1
vote
3answers
147 views

Why does the Metropolis algorithm allow changes even for ∆E > 0?

In the Metropolis Monte Carlo algorithm, why can you accept changes even for ∆E > 0 (provided that a random number is less than a given probability ratio, e.g. exp(-β∆E))?
1
vote
2answers
196 views

What is a good model for computing water dropping on a surface?

Before introducing my question I would like to outline the fact that I'm a coder, so I can be wrong when using some kind of terminology. What is a good model for computing the flow, the shape and all ...
1
vote
4answers
1k views

Non-computable numbers in Physics

It is known that there are real numbers that can't be calculated (non-computable numbers). Quite probably that some physical phenomena (it is possible still undetected) depend on this numbers. Whether ...
1
vote
1answer
75 views

Difference between Enzo & Gadget astronomy simulation codes

Enzo and Gadget are simulation codes used in astronomy. What are the largest differences between them both in terms of physics they simulate and in their implementations?
1
vote
1answer
388 views

I'm getting weird autocorrelations when simulating an Ising model below the critical temperature

So I'm simulating an Ising model using Monte Carlo and the Metropolis algorithm. After letting it reach equilibrium, I try to calculate the autocorrelation of the magnetization. As long as the system ...
1
vote
2answers
2k 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
1answer
26 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 ...
1
vote
2answers
168 views

Do electrons actually reside in orbitals?

(Following this answer on Chemistry.SE) Let's define "dwell"/"reside" as the area in which the probability to find an electron is more than 90%, for example. Calculated orbitals comprise a basis ...
1
vote
1answer
72 views

SDSS spectra FITS format: what do the rows mean?

This might be in the FITS Header, but I couldn't find it. I am looking at the FITS files for SDSS 1D spectro images from DR7 using the SDSS Data Archive Server at das.sdss.org (there might be an ...
1
vote
2answers
98 views

Langevin Equation - Stochastic Differential Equation. What are the subtleties?

I am trying to find out the motion of a particle in 3D governed by the Langevin equation, numerically. Anyway, the Langevin equation is given by $$m \ddot{x} = -(6\pi a\nu) \dot{x} + F_b $$ where ...
1
vote
2answers
112 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, ...
1
vote
2answers
763 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
372 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
1k 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
46 views

Gravitational field of thin 2D ring - numerical simulation

I'm aware of Newton's Shell Theorem, which states that inside of a thin ring of uniform density, the gravitational force exerted on a point mass should be zero. I wrote a quick field simulation to ...
1
vote
1answer
134 views

Connection between bond-dimension of a matrix product state and entanglement

The bond dimension is the dimension of the truncated matrix product state (MPS). Let us assume that I am simulating some many-body system with high entanglement via the density matrix renormalization ...
1
vote
1answer
260 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
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 ...
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
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 ...
1
vote
1answer
388 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
172 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
146 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
2answers
71 views

Can lens designs be brute-force computed to create optically precise equipment? [closed]

It seems lens design hasn't changed much over the last 50 years or so with many older manual lenses still very much being used today. Can the order, distance between and type of lens elements be brute ...
1
vote
1answer
28 views

Does *advection* describe the change of density of massless infinitesimal tiny *thingies* injected into a fluid?

I'm considering an incompressible Newtonian fluid with uniform density and try to figure out what's meant by the term advection. Let $\Omega_0\subseteq\mathbb R^d$ be an (infinitesimal small) bounded ...
1
vote
1answer
68 views

Numerical relativity in causally pathological spacetimes

To perform numerical relativity simulations one almost universally adopts the so called "3+1" approach: spacetime is divided up into spacelike slices, each representing a "moment in time". After some ...
1
vote
1answer
36 views

Reading ensdf files

I understand that this is more of a computing question, but I still think the appropriate place to ask it is here. The Evaluated Nuclear Structure Data Files are on the internet, available for ...
1
vote
2answers
66 views

Guidelines for $n$-body problem inital conditions

I'm trying to solve a 5-body problem (for the sake of the argument it can be any given $n$) numerically, yet every initial condition I choose seems to eventually drive the bodies apart to infinity or ...