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)

0
votes
0answers
10 views

Swendsen Wang different implementations

I have recently a confusion about the different adding probabilities of bonds in the Swendsen Wang cluster algorithms with e.g. application to the Ising or Potts model. In literature there are ...
0
votes
0answers
43 views

Two Particles in a Harmonic Oscillator with repulsive short-range potential

Do bear with me, I am attempting to learn to write some simulations on the computer and learn some simple MD, so I defined sort of a toy problem. I have two particles confined in a Harmonic Potential ...
1
vote
2answers
51 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 ...
-4
votes
0answers
27 views

Testing Bell's inequality violated by quantum mechanics [closed]

I made this code to test Bell's inequality that violated by quantum mechanics that said the max value is 2*sqrt(2), but when i run, it gave me the max value only 2. so could you please tell me what i ...
1
vote
1answer
20 views

Adding fluctuations to a hydrodynamic simulation to trigger instabilities

I'm doing a 2D combustion hydrodynamic simulation and there's a hydrodynamic instability that should be triggered because of the particular physical properties of my system. The key to the instability ...
0
votes
0answers
20 views

Calculating atomic forces from total force and torque on molecule in molecular dynamics (MD)

I guess this problem applies to any body that is composed of a small number of particles, but for the sake of simplicity and concreteness, let's consider a classical model of a molecule consisting of ...
0
votes
0answers
35 views

MEMS Accelometer and Gyroscope data fusion for counting cycles

I will simplify my scenario as follows: a 6-DOF MEMS gyroscope and accelerometer unit is placed on the edge of a rotating wheel and the goal is to count the number of cycles using any of the units of ...
0
votes
0answers
9 views

Arbitrary waveform generator (How does the input of AWG affects the output)

How should one choose the sampled input for AWG so that it gives a waveform similar to the ideal waveform? I have a dataset for voltage signal sampled over time of 1 ns, which I am entering as an ...
1
vote
0answers
25 views

Minimal possible energy to make calculation [duplicate]

Are there any restrictions to energy that must be consumed to make calculations? Can calculations theoretically be done without consuming it?
2
votes
3answers
94 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 ...
3
votes
2answers
103 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
2answers
139 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 ...
5
votes
0answers
81 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 ...
2
votes
1answer
62 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
33 views

Lennard Jones Total system energy [closed]

i am trying to implement/extend an implementation of Lennard-Jones potential simulation regarding Xenon molecules (for curious ones, the code can be found here functions (force_naive->lj_force)). ...
0
votes
1answer
64 views

Classical Computation without NOT [closed]

Is it possible to do universal classical computation using bits and 2-bit gates when you cannot perform a NOT operation on a single bit (hence cant do CNOT and so on). If yes, what are the possible ...
0
votes
1answer
37 views

How to produce a 3D density map of a time-depenent system of particles?

I have a time-dependent system of varying number of particles (~100k particles). In fact, each particle represents an interaction in a 3D space with a particular strength. Thus, each particle has ...
0
votes
1answer
43 views

The meaning of $\vec {\dot r}$ in the solution of a particle's movement in the space

Suppose I have a charged particle $q$ with mass $m$, and an infinite long wire that lies on the $z$ axis, in which flows a constant current $I$. My ODE that describes its motion is $$\vec {\ddot ...
0
votes
2answers
37 views

How to find time when temperature crossed certain threshold?

I have a data from experimental system concerning temperature of molten lead. In this system, temperature is increased in one place and then is measured when this increase of temperature travels to a ...
0
votes
0answers
17 views

Can we use Variational Monte Carlo for degenerate cases?

Consider Simple Example of Bose-Hubbard model $$H=-J\sum\limits_{<i,j>}b_i^{\dagger}b_j+h.c.+\frac{U}{2}\sum\limits_{i}n_i(n_i-1) . \tag{1}$$ We can solve this Hamiltonian by Variational ...
1
vote
0answers
23 views

Computational method for finding edge states?

I am actually interested to learn how to calculate edge states in 1D topological systems using computational methods, Q. can anyone tell me which method is best suited and easy to calculate edge ...
1
vote
1answer
69 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
0answers
67 views

What are the differences between the Jetphox, Pythia and Herwig event generators?

I know Jetphox is a parton-level event NLO generator program. But I want to know more about other generator programs such as Pythia and Herwig. What are the differences? I am undergraduate student so ...
1
vote
1answer
41 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?
0
votes
0answers
49 views

Importance sampling for Coulomb potential

The integral I have to solve is: $$I=\int\int d \mathbf{r}d \mathbf{r}' \frac{\Phi(\mathbf{r})\Phi(\mathbf{r}')}{|\mathbf{r}-\mathbf{r}'|}$$ It is a six-dimensional integral which I am going to ...
2
votes
2answers
108 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 ...
1
vote
1answer
64 views

What is the velocity of Sun due to Earth alone?

If Earth and Sun were in a isolated system, will the Earth's motion around Sun will be similar? What will be Sun's and Earth's velocity when Earth is at its aphelion? Please note that it's not a ...
0
votes
0answers
19 views

Modelling Fourier Transform Profilometry

Basically I want to simulate a surface profilometry technique through Matlab. For that I want to create a GUI with controls for generating a grating pattern of light at a particular angle with respect ...
3
votes
2answers
77 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 ...
2
votes
1answer
59 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 ...
1
vote
0answers
67 views

Selecting physical solutions in numerical eigenvalue problems

I try to solve a certain time-independent Schrodinger equation numerically, using the method of finite differences. My boundary conditions are such that the finite difference method gives me an ...
0
votes
0answers
66 views

Limitations of Quantum Simulations on a Classical Computer

My question is on the simulation of a quantum computer on a classical machine. I understand that a classical computer to simulate any quantum algorithm--the problem is that the quantum computer does ...
0
votes
1answer
40 views

Sampling a 2-dimensional surface: How many samples along X & Y axes?

I have a set of first 25 Zernike polynomials. Below are shown few in Cartesin co-ordinate system. z2 = 2*x z3 = 2*y z4 = sqrt(3)*(2*x^2+2*y^2-1) : : ...
2
votes
1answer
111 views

What is this nested bracket notation?

The following is an excerpt from K. Varga's paper, Precise solution of few-body problems with stochastic variational method on correlated Gaussian basis: ...The function $θ_{LM_L}(\mathbf{x})$ in ...
0
votes
1answer
139 views

How to calculate impulse force on collision?

I've made a program which is supposed to simulate free fall of a ball on a ground. For collision with the ground I relied on calculating the impulse force with the answer in this Phys.SE question ...
1
vote
2answers
113 views

Incompressible Navier-Stokes pressure solve in simulations

I am a complete newcomer when it comes to fluid simulations. I'm currently working through some tutorials to understand the idea of of the discretized Navier-Stokes equations for numerical ...
0
votes
0answers
23 views

Is there any way I can tell whether increase in total energy in molecular dynamics simulation is due to large time step or not?

Does increase due to large time step has any type of signature that will tell for certain (or almost certainly) that that is the reason?
0
votes
2answers
93 views

Which units should I use in molecular dynamics simulation?

I have written a simple molecular dynamics simulation program. The simulation runs fine but the physical properties (in particular, I have calculated temperature) are off by many scales. I understand ...
2
votes
2answers
90 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
143 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: ...
5
votes
5answers
330 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 ...
19
votes
3answers
2k 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
48 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
38 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
49 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
73 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
92 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
84 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
62 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
82 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 ...