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
2answers
56 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
1answer
22 views

Expression for spontaneous emission in solids

I have yet to find a reliable source for the expression of the rate of spontaneous emission in a solid. Can anyone confirm if the following is correct? The basic ingredients of the calculation are ...
0
votes
1answer
43 views

An impossible relationship between variables

This book chapter defines $CV$ as follows: $$ CV^2 = 2\pi\nu^2\int^\frac{V_{th}-V_{ss}}{\sigma_V}_\frac{V_{r}-V_{ss}}{\sigma_V}dxe^{x^2}\int^x_{-\inf}dye^{y^2}(1 + erf(y)) $$ However, Figure 15.2 in ...
0
votes
1answer
31 views

Solving for the firing rate of a model neuron

I'm trying to decipher Figure 15.2A on page 442 of this book chapter. It plots the firing rate $\nu$ of a model neuron against its mean current $\mu_c$. The equation used to calculate the firing ...
2
votes
0answers
23 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 ...
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
44 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
54 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
22 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
36 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
96 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
111 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
82 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
64 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
38 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
42 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
71 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
42 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
110 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
70 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
67 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
154 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
116 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
26 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
97 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
332 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
50 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 ...