Questions tagged [computational-physics]

Questions with this tag should be about computational physics, which is the study of physical situations with the use of software (commercial or in-house). Please note that details of writing and/or debugging code is OFF-TOPIC and should be asked at either Computational Science, Code Review or Stack Overflow.

Filter by
Sorted by
Tagged with
4
votes
1answer
73 views

Why is FDTD derived directly from Maxwell's equations instead of the wave equation

I've been wondering why the Finite Difference Time Domain Method is derived directly from Maxwell's equations and not directly from the electromagnetic wave equation (that in theory is also derived ...
4
votes
2answers
548 views

Sample the Maxwell-Jüttner distribution

The Maxwell-Bolztmann distribution is easy to sample numerically, as it can be expressed directly with a unique variable $mv^2/k_B T$. One simply has to tabulate the cumulative distribution function ...
4
votes
1answer
121 views

Finding explicit unimodular transformations for Chern-Simons K-matrices

An invertible, symmetric matrix with integer entries, $K$, that encodes the braiding and statistics of an Abelian topologically ordered state, is equivalent to another such matrix, $K'$, if there ...
4
votes
1answer
110 views

How to conserve momentum when numerically-integrating the path of a charged particle through a magnetic field

I have the python script attached below, which is intended to track the trajectory of a charged particle in a static, uniform magnetic field. It is very simple to calculate the instantaneous force ...
4
votes
1answer
2k views

Is Machine Learning a useful tool for analyzing experimental data? [closed]

Apart from being a student of Physics I'm also a developer and recently in the context of cloud development I've been hearing a lot about cloud based Machine Learning. After looking at some examples ...
4
votes
1answer
557 views

Philae lander simulation off by factor of ~3

I'm trying to simulate the Philae landing by writing a program to compute the position of the lander vs time. According to various mission websites, the orbiter will match its orbit to the rotation of ...
4
votes
2answers
730 views

Epsilon Tensor in FeynCalc

A few days ago I started to use the Mathematica package FeynCalc and one thing confuses me: Assume we have a four-vector $p_\mu$ and we contract it with the epsilon tensor. FeynCalc produces $\...
4
votes
1answer
162 views

How to find positions of $n$ masses in Newton mechanics?

I ran into a problem while doing research. The problem can be described as: consider the original $n$-body problem, and if we fix the position of them(unknowns), no interaction among them, they don't ...
4
votes
1answer
110 views

Computronium and Time Dilation and Bremermann's Limit

I was contemplating computronium the other day and it occurred to me that there might be a point at which the computronium becomes so dense that gravitational time dilation effects would cause it to ...
4
votes
1answer
113 views

Physical meaning behind the choice of the proposal distribution in Markov Chain Monte Carlo (MCMC) methods

Let us consider the conventions on names used in the theoretical derivation of Metropolis-Hastings Monte Carlo as outlined here, for the sake of common nomenclature. What we are building is a step-...
4
votes
1answer
85 views

Finding $\pm 2 \pi$ defects in 2-D lattice nematic simulation

I'm working on a Monte Carlo simulation of a two-dimensional nematic system (XY-like model with even-order Legendre polynomial interactions, such that the director angle $\theta$ obeys $\theta \equiv \...
4
votes
1answer
925 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): $\vec{r}_{Earth}=\begin{pmatrix}-\frac{m_{Moon}}{m_{...
4
votes
1answer
2k views

How to solve the heat equation for compound materials with different heat conductivities numerically?

I'm solving the heat equation with time dependent boundary conditions numerically in a 2D system using the ADI scheme. For the purpose of this question, let's assume a constant heat conductivity and ...
4
votes
1answer
226 views

Difference between a Fixed Point and a Limit Point in implementations of the Renormalization Group (RNG) in Large Eddy Simulation (LES) model

In the introduction of this paper, it is explained that and how the application of a dynamic subrid scale model for turbulence into a large eddy simulation (LES) model corresponds to doing one ...
4
votes
1answer
9k views

Calculating displacement amplitude of ultrasonic power transducer

I am not a physicist, but my current project drives me to some physics-related computations, hence seeking help. I have some ultrasonic transducers, 5938D-25LBPZT, for which very limited information ...
4
votes
0answers
97 views

Method to build a polyhedral die with given probabilities [closed]

Let's define a die as a polyhedron that, if rolled over a perfect horizontal plane, ends up being in a physically stable unambiguous state labelled $n$. The die has $N$ states. Each state $n$ has a ...
4
votes
0answers
152 views

$N$-body gravity simulator: why does energy conservation break down when introducing an adaptive timestep?

I am playing with an N-body gravity simulator using the velocity-verlet algorithm. In the actual simulation I normalize everything so that it is nicely behaved numerically, but I'll express things in ...
4
votes
2answers
231 views

Dipole field issue in particle-mesh Ewald method with periodic boundary conditions

I am working on a thesis that makes a great use of molecular-dynamics simulations, and I am trying to understand how the particle-mesh Ewald method works. The problem is, I have difficulties ...
4
votes
0answers
106 views

Textbooks on algorithms for the perturbative calculation of High energy physics

For the perturbative calculation of High energy physics, I have known some packages such as FeynArts, FeynCalc, MadGraph, CompHEP, GiNaC, and so on. But I am wondering whether there exists a textbook ...
4
votes
0answers
117 views

Nature as a calculator [closed]

I wondered whether physical experiments can find the outcome of otherwise intractable mathematical problems. For example, solitons or other non-perturbative effects in QFT cannot be seen at any ...
4
votes
0answers
2k views

Numerical problem in solving the Bogoliubov de Gennes equations- methods to solve?

I am trying to solve an assignment on solving the Bogoliubov de Gennes equations self-consistently in Matlab. BdG equations in 1-Dimension are as follows:- $$\left(\begin{array}{cc} -\frac{\hbar^{2}}...
4
votes
2answers
473 views

Can computers accurately model all of the details (to the subatomic level) of macro objects in collisions?

Frequently when trying to solve cosmology questions physicists turn to computer simulations of the universe (albeit massively simplified) in order to verify or disprove their hypotheses. This got me ...
3
votes
2answers
1k views

Solving Schrödinger equation for an infinite potential well with an obstacle (Numerical/Analytic)

I'm a beginner in both quantum mechanics and computational physics, So please be as simple as possible. I know these topics can become so complicated. I've been trying to solve the Schrödinger ...
3
votes
2answers
594 views

System that is more efficient to simulate than to compute

I know some physics calculations require alot of computing to solve. Are there currently systems where it is more logical to do an experiment to "let the universe simulate the experiment for us" and ...
3
votes
4answers
2k views

Schrodinger equation for a Hamiltonian with explicit time-dependence?

Can I write a Schrodinger equation for time-dependent Hamiltonian like this: $$i\hbar\frac{d}{dt}\psi(t) = H(t)\psi(t)$$ and then perform Euler integration like this: $$\psi(t+\Delta t) = (1-\frac{...
3
votes
3answers
551 views

Can a wavefunction be solved to any arbitrary precision, given enough computer time?

I learned that the wavefunction for the hydrogen atom can be solved analytically (we did the derivation in class), but that for more complicated atoms it is "impossible" to solve and that only ...
3
votes
3answers
3k 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 ...
3
votes
2answers
138 views

How much of the standard model could have been discovered in a world without computers?

Computers were not needed for the discovery of quantum mechanics (1920s) or general relativity (1910s), and the standard model was proposed when computers were much less powerful than today (1970s). ...
3
votes
3answers
692 views

Computation theory and the simulation argument

Can physical states be treated as information (strings over some alphabet)? If (1) is true, isn't this a trivial conclusion that the universe can be simulated by a Turing machine or a cellular ...
3
votes
2answers
114 views

How to actually find a Hartree-Fock ground state?

I am interested in finding the Hartree-Fock ground state for a system of interacting fermions (with totally local scattering, so a delta-function interaction potential). I have read through some ...
3
votes
1answer
71 views

Laplace's heat equation equilibrium functions

I simulated the Laplace's heat equation on a grid (200x200) with boundary conditions of 0 (a cold source) on the top and left and 1 (of hot source) on the bottom and right. (see gif) I stopped ...
3
votes
1answer
68 views

How to identify binary stars in $N$-body simulation?

Binary stars constitute a significant portion of the stars of a globular cluster. I would like to verify that this is true in my $N$-body simulation, but I don't know how to decide whether a star in ...
3
votes
2answers
249 views

How to initialize / bootstrap the Boris algorithm?

The Boris algorithm requires interleaved positions and velocities, i.e. position measured at $t_i$ and velocity at $t_{i-1/2}$ for example. I want to employ the Boris pusher in my application however ...
3
votes
3answers
604 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 $...
3
votes
3answers
2k 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
465 views

Sampling a distribution (from a galaxy model)

I am reading the following article: http://www.kof.zcu.cz/st/dis/schwarzmeier/galaxy_models.html and am currently at section 5.6 (positions of bodies in a galaxy). I am trying to redo the simulations ...
3
votes
2answers
331 views

How can I explicitly express the Ising Hamiltonian in matrix form?

I am reading this book about numerical methods in physics. It has the following question: Consider the Ising Hamiltonian defined as following $$H=-\sum_ {i=1}^{N-1} \sigma_i^x \sigma_ {i+1} ^x + h ...
3
votes
2answers
238 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 $\text{max}\{B_{...
3
votes
2answers
3k 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 ...
3
votes
4answers
2k views

Nanorobots. What stops us from producing them yet? [closed]

If we can already predicts accuratelly motion on molecular levels, what stops us from developing small robots to, for instance, navigate through our blood vessels looking for cancerous cells and ...
3
votes
2answers
389 views

Mixed conductive and convective heat transfer

I want to make a finite element analysis of a cold airflow through warmer pipes. In particular I want to see how the pipes cool down and the air heats up, as it travels through the pipes. Wich are the ...
3
votes
1answer
67 views

Where's the Database in Grover's Algorithm

I am a little confused by parts of Quantum Computation and Quantum Information by Nielsen and Chuang on Grover's algorithm: Grover's algorithm searches a number of quantum states for a "marked ...
3
votes
1answer
455 views

How to numerically solve the Schrodinger equation for the hydrogen atom?

From the analytical solution for the hydrogen atom I remember that the energies were obtained by boundary conditions and then the energy depends on n*n in which n is interpreted as different energy ...
3
votes
2answers
215 views

Updating link variables in lattice $SU(N)$ gauge theory

I'm currently writing a basic program in python to simulate a 1 + 1 dimensional yang mills gauge theory with symmetry group $SU(2)$. On the lattice you work with link variables, which are $SU(N)$ ...
3
votes
1answer
297 views

Tensor network Renormalization group

I am working on Tensor network Renormalization group. While I am trying to implement TRG on Honeycomb lattice I got stuck on pairing two three rank tensors and represent it as D^2*D^2 matrix. Can you ...
3
votes
1answer
617 views

Finding the Riemann tensor for the surface of a sphere with sympy.diffgeom

I have implemented a SymPy program that can calculate the Riemann curvature tensor for a given curve element. However, I am encountering problems solving for the case when the curve element is the ...
3
votes
1answer
366 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 ...
3
votes
2answers
136 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 ...