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
1answer
83 views

What is the time imaginary method? [closed]

I have to submit homework about the scheme which solves the time-independent Schrödinger equation and finds the ground state by the imaginary time method. I know the substitution $-\mathrm{i}\tau\...
0
votes
0answers
31 views

Numerically evaluating momentum sums

For many-body systems, calculations often involve momentum sums, which are ultimately evaluated numerically for a finite system size $N$. The larger $N$ is, the more accurate the results. Is there ...
0
votes
1answer
28 views

Numerical computation of electric multipoles

I'm trying to do a plot of the first multipole terms in Mathematica. My plot isn't what I expected, so maybe my problem is in the mathematics. Following the Griffiths' book on Electrodynamics, the ...
1
vote
0answers
76 views

Experiments to resolve dillema between continuity and discrete

Which experiments/experimental methods are suggested to resolve an alternative about the structure of our universe space and time - is it continuous or is it discrete in a very small scale, especially ...
1
vote
1answer
77 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 ...
0
votes
0answers
40 views

canonical ensemble density matrix numerical integration of the von Neumann equation

I am working with a numerically calculated Eigensystem of a given non-linear Hamiltonian. As a test I integrated Schrödinger equation for each eigenstate in time with Dormand Prince 54 and also Runge ...
30
votes
6answers
6k views

Is there any aspect of an explosion resulting from a nuclear weapon test that cannot be simulated using super computers?

This Washington Post news article states that with the advent of computer simulation of nuclear tests, live tests are no longer needed. Generally speaking there are 3 aspects of an explosion ...
1
vote
1answer
48 views

Single fluid simulation algorithm

Every fluid simulation I've seen implicitly involves a fluid that isn't mentioned, usually air. This is clear because the simulation usually has an ambient temperature, and always an ambient velocity. ...
1
vote
1answer
28 views

Why will a Fourier transform of a pickup placed at the middle of a string only result in odd harmonics?

My physics teacher explained a lot of information about strings but he never went in depth with my question. What happens to the even harmonics and why are the odd ones the only ones picked up?
0
votes
0answers
49 views

recovering numerical wave function using path integral after Wick rotation

I have written two different path-integral codes, PATHINT and PATHTREE, to numerically solve some classical-physics problems in nonlinear systems, finance and neuroscience. They work just fine. My ...
1
vote
2answers
188 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 set,...
0
votes
0answers
22 views

Heat transfer in liquid helium channel connected to “infinite” reservoir

I am trying to solve the temperature distribution in liquid helium using numerical analysis (OpenFOAM for those who are familiar with it). One part of the fluid is contained in a narrow channel, one ...
0
votes
0answers
80 views

Graduate level orbital mechanics book

I recently finished an undergraduate course in classical mechanics and really enjoyed the subject, particularly the sections regarding the mechanics of orbits. I am considering pursuing a graduate ...
2
votes
1answer
59 views

Difference between finite volume, characteristic method and plug flow models of a pipe

I have to model pipes (of a district heating network) with ODE's. My background is computer-science, so it is not that easy for me to understand different approaches. Finite volume approach. Method ...
2
votes
0answers
45 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 ...
0
votes
0answers
166 views

How do you fourier transform a tight binding hamiltonian numerically?

The task is to do a fourier transformation of a tight binding hamiltonian of a 1D-chain with unit cell size 2, but even after many tries and googling I still don't have a idea how to do it correctly. ...
0
votes
2answers
93 views

Obtaining Velocity from Acceleration

I'm following a research paper (PDF via RG), part of which is about using the accelerometer of a smartphone to assist the user's positioning via WLAN. Accelerometer is used to determine if the user ...
2
votes
2answers
129 views

Exact Diagonalization of a BdG Hamiltonian on a Finite Lattice

I would like to numerically find the edge modes of a $p_x$ + $i p_y$ BdG Hamiltonian. The lattice version is given by H = $\sum\left[-t \left(c_{m+1,n}^{\dagger} c_{m,n} + \text{h.c} \right) - t\left(...
0
votes
0answers
12 views

Numeric fermiomic creation operators and unit cell

I have to do some numerics (e.g. FFT in Maple/Octave) on a 1D fermionic chain without forces between the particles. The description says that two sites build an unit cell. What does this mean? And ...
0
votes
2answers
86 views

Molecular Dynamics simulation question

In molecular dynamics simulations you apply and solve Newton’s equations of motion for a set of atoms to obtain real dynamical information on how individual atomic positions change as a function of ...
1
vote
1answer
93 views

Modeling transient heat transfer between two 1-D materials

So I have been working on creating a 1-D model of a material with a layer of heat resistant material stuck to it. For this I have been trying to use the implicit finite-differential method. The ...
0
votes
0answers
41 views

RF cavity Field Maps using Poisson Superfish

This may be a slightly obscure question but I'm trying to use Poisson Superfish to generate field maps for single cell of an ILC SRF Cavity. There are examples in the documentation for radio ...
2
votes
1answer
94 views

Why is modelling water so difficult?

There are several models of water that have come into existence so far. Nevertheless, there is no one particular model which effectively gives all the properties of water correctly. A model like TIP3P ...
0
votes
1answer
84 views

Is it possible to understand physics and make new discoveries using computer simulation? [closed]

I'm a computer science major and I want to learn Physics. I can create computer simulations of any type. I'm not good at math that is required to learn QFT or GR,but I'm thinking is it possible to ...
2
votes
1answer
174 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 ...
1
vote
1answer
111 views

How do I calculate the effect of a prograde, retrograde, radial and anti-radial burn on the orbital elements of a two-dimensional orbit?

My relevant knowledge: This question works in two dimensions, with a two-body problem (see bottom for context). To my understanding, a two-dimensional orbit has the following orbital elements (3D ...
-1
votes
2answers
84 views

What is the output of a simulation of quantum particles? [closed]

When physicists simulate a quantum system, what is the output of such computation? Is it the probability wave function of the system? What if the system consists of enough particles to be considered ...
0
votes
1answer
138 views

Calculating Reynolds number of ellipse in fluid

I know Reynolds number is calculated from the formula by including length or diameter of an object in fluid. I have a case where there is ellipse with 2a major and 2b axis. I hv learnt when the ...
0
votes
2answers
184 views

Must-know topics for Computational Fluid Dynamics (CFD)

My undergrad major was Electrical Engineering. My MS thesis is related to simulation of hydrogen movement in a complex geometry. I did take a course on Fluid Dynamics but I have to say it wasn't one ...
1
vote
0answers
91 views

Obtaining propagating solutions for Schrodinger equation from known bound states (in 2 and 3 dimensions)?

If I found all the bound states for a certaing potential in 2 or 3 dimensions (numerically), can I immediately obtain some information about the propagating solutions for the same potential (such as ...
0
votes
3answers
177 views

How to propagate $E_x (x,t) = \exp(-t^2/\tau^2-i\omega_0 t) \exp(-x^2/w_0^2)$ in finite difference time domain (FDTD) analysis

Finite difference time domain (FDTD) allows to solve differential equations for time evolution. For example, we can analyze ultra-short pulses in free space by solving the Maxwell's equations. The ...
0
votes
0answers
51 views

Numerical approach to Mean Field Hubbard Model

I am trying to produce the phase diagram of a 2D Hubbard model on a Lattice for spineless bosons. The hamiltonian of the problem is: $$ H = -t\sum_{<i,j>}(b_i^+b_j + c.c.) +U\sum_i \frac{n_i(n_i-...
2
votes
1answer
195 views

Problems implementing Lattice Boltzmann method

I want to make a basic 2D implementation of the Lattice Boltzmann method (in javascript) to simulate gas/fluid flow, but I have run into some problems and I can't seem to find the cause of it. ...
0
votes
2answers
142 views

How Much Computing Power would be Required to Fully Simulate a Cubic Meter? [closed]

Imagine you want to simulate a cubic meter down to the particle. By following the Standard Model and other basic physical equations, how much computing power would be required to do this, in say, a ...
0
votes
2answers
210 views

Diffusion of gas into vaccum

I'm interested in solving the diffusion equation for gas in vacuum. I have a general question and a more specific questions. What I know: The Diffusion Equation: For density function $\phi(\vec{\...
1
vote
1answer
45 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 ...
3
votes
1answer
202 views

Calculation of long-range forces in Molecular Dynamics - Ewald summation

I am trying to write a code to calculate the potential and forces, for the same using ewald summation.For this purpose, the formula for potential and force I have used is : $$ U = U^{(r)} + U^{(k)} +...
1
vote
1answer
125 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 ...
0
votes
1answer
48 views

How does Parallel tempering work for optimisation procedures?

Markov Chain Monte Carlo procedures can be used in global optimisation procedures. Parallel tempering is a particular type of stochastic sampling, but how does it work?
3
votes
1answer
47 views

Symmetry in program for ewald summation

The formula for Ewald summation as given in Allen and Tildesley - $$ U = U^{(r)} + U^{(k)} + U^{(bc)} + U^{self} $$ where the k-space contribution of potential is given by $$ U^{(k)} = \frac{1}{\pi ...
1
vote
3answers
147 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
31 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
48 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
42 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
63 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
19 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
116 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
74 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
41 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 ...