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

0
votes
0answers
3 views

Atomic flux on the surface?

I have a surface which is made up of two triangles of unequal areas. I impinge a single atom on the surface at a certain angle at a certain time interval. As flux is atoms/unit area/unit time. In ...
0
votes
0answers
14 views

Graphene plasmonics, Integro-differential eigenvalue equation

In my research around plasmonics in 2-Dimensional nanostructures like graphene, I encounter an eigenvalue integro-differential equation of the form: ...
2
votes
1answer
173 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 ...
1
vote
2answers
70 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, ...
0
votes
1answer
198 views

Position based dynamics constraint scaling factor

Reading through Müller et al., Position Based Dynamics, 2007 I got lost when passing from equation (5) $$\Delta p = \frac{C(p)}{|\nabla_pC(p)|^2}\nabla_pC(p)$$ to equation (6) (and applying the ...
0
votes
1answer
82 views

Leapfrog method in Particle-in-cell

Recently I wanted to write a 3D electromagnetic Particle-in-cell code with c++.I know that I should use leapfrog method.For example,when I calculates the position and velocity of particles,i should ...
0
votes
2answers
35 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
57 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
54 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
47 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 ...
1
vote
1answer
62 views

Determing Velocity of Moons

I have a question that I believe is relatively easy to answer, I am working on an $N$-body simulation of a fictional star system and am having trouble finding the velocity of moons so that they will ...
0
votes
1answer
216 views

My python trajectory calculator for a small mass in a gravitational field - why is energy conservation violated?

I wrote this thing to calculate the next position of a small mass, given its current position and velocity, and a vector field determines acceleration. It sort of works - except I can never get the ...
1
vote
0answers
15 views

Modeling simple laser induced population transfer via adiabatic passage in python [migrated]

I'm trying to model adiabatic passage between two levels in a three-level atom interacting with two laser fields using Scipy and Numpy.. I'm not sure if my model is wrong due to my incorrectly ...
2
votes
1answer
55 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 ...
1
vote
0answers
31 views

Mixed spin Ising Model

As we know ferrimagnets can be modeled by the Ising model. I came across this equation in "Compensation Temperature of the Mixed-Spin Ising Model on the Hexagonal Lattice" by W. Figueiredo, M. Godoy, ...
7
votes
5answers
769 views

Church-Turing hypothesis as a fundamental law of physics

The Church-Turing hypothesis says one can not build a computing device which has more computing power (in terms of computability) than the abstract model of Turing machine. So, there is something in ...
0
votes
0answers
25 views

Mechanics - Rotational Inertia [migrated]

I am facing a problem solving a mechanical system which is shown in the scheme below. We assume that coupling of the two shafts is done with a clutch Equations: ...
1
vote
0answers
49 views

Tight binding model in a magnetic field

The standard way to treat a tight binding method in a magnetic is to replace the hopping matrix element: $t_{i,j}\rightarrow e^{i\int_i^j\mathbf{A(x)}.d\mathbf{x}}$ the so called "Peierls ...
0
votes
1answer
72 views

How to solve highly oscillating differential equation [closed]

The equation looks like: $$x''(t)+bx'(t)+c x(t)+dx^3(t)=0.$$ This is the motion of a particle in a potential $cx^2/2+dx^4/4$ with friction force $bx'$. In my case, the friction term is very small and ...
0
votes
0answers
30 views

What is the largest atom with a reliable configuration-interaction (CI) calculation?

The simplest approximation for calculating the ground state of an atom is the Hartree-Fock approximation. To get accurate result for the ground state energy, one has to do configuration-interaction ...
0
votes
0answers
26 views

Q: Any standard terminology for “gold standard” Monte Carlo solutions

The "standard" way in our field to compute scattering statistics in wave propagation problems is via a "deterministic" wave prop code operating in a Monte Carlo context where we iterate over many ...
0
votes
1answer
54 views

Connections between Density Matrix Renormalization Group and Conformal Field Theory

Can we use the density matrix renormalization group (DMRG) method to understand problems in conformal field theory? I have been trying to find some connections, but nothing is coming up when I search. ...
0
votes
0answers
21 views

What is the difference between finite displacement and linear response for calculating vibrational properties?

I see these concepts appearing in the context of calculating phonons and other vibrational properties, but I can't find a concrete explanation of the differences between linear response (DFPT) and the ...
1
vote
2answers
63 views

How to verify/falsify the existence of localised edge states numerically?

I have to consider a Hamiltonian given in second quantized form in real space $$H = \sum c_k^\dagger h_{kl} c_l \, ,$$ describing fermions on a 2d hypercubic lattice. The concrete form of the matrix ...
1
vote
0answers
60 views

How do I model this Quantum Cascade Laser (QCL) using COMSOL? [closed]

I'm trying to model this: GaAs on-top of a heat-sink made of Copper. It's trying to simulate a type of Quantum Cascade Laser, I'm looking at the effect the temperature has by running a 10V voltage ...
0
votes
1answer
31 views

Amount of unknown parameters in compressible Euler equations

I'm looking at this page for the compressible Euler equations. To me it seems, in the 1-dimensional case, there should be 3 unknowns: density, velocity, and pressure. This is because the energy $E$ ...
7
votes
2answers
4k views

Solving one dimensional Schrodinger equation with finite difference method

Consider the one-dimensional Schrodinger equation $$-\frac{1}{2}D^2 \psi(x)+V(x)\psi(x)=E\psi(x)$$ where $D^2=\dfrac{d^2}{dx^2},V(x)=-\dfrac{1}{|x|}$. I want to calculate the ground state ...
-4
votes
1answer
68 views

Computational physics using mathematica [closed]

So I was confused about this question on how to exactly begin to answer it. I am a novice in mathematica and I am teaching myself thus I require help in this question. From what I think I should do, ...
1
vote
1answer
77 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 ...
3
votes
2answers
466 views

Numerical solution to Schrödinger equation - eigenvalues

This is my first question on here. I'm trying to numerically solve the Schrödinger equation for the Woods-Saxon Potential and find the energy eigenvalues and eigenfunctions but I am confused about how ...
1
vote
1answer
393 views

Does the universe follow a general physical and computational principle of evolution?

Evolution is a principle in biology, whereby organisms evolve their ability to replicate and multiply in numbers over successive generations. From a computational point of view, the organisms employ a ...
0
votes
1answer
51 views

Using the Fourier transform to find the natural frequencies of coupled oscillators

How can I find the natural frequencies of a system consisting of a pair of coupled oscillators using Fourier transforms? The System consists of two masses and three springs. One of the springs ...
1
vote
0answers
21 views

Introduction material for lattice field simulation

I am looking for introductory material of lattice field theory simulation. It is better start from simple example (e.g. \phi^4) and include some source code. Is there any on-line or book which is ...
12
votes
2answers
397 views

How do you simulate chiral gauge theories on a computer?

David Tong and Lubos Motl have argued that our universe can't possibly be a digital computer simulation because chiral gauge theories can't be discretized, and the Standard Model is a chiral gauge ...
0
votes
0answers
40 views

Ising Model 2D Correlation Length

I'm using Metropolis and Wolff Clustering algorithms to estimate the spin-spin correlation function $$<s_os_r>$$ I know that this is related to the correlation length but how do we determine ...
1
vote
1answer
116 views

How to measure the spin-spin correlation in a Monte Carlo simulation of the Ising model?

I'm simulating the Ising Model in 2D up to 5D and I want to calculate the spin-spin correlation, correlation length, and critical exponent of the system. What is a good way to go about doing this? ...
1
vote
0answers
50 views

Ising model Monte Carlo simulations in 4D and 5D

I'm going to be simulating the Ising Model in 4D and above to calculate spin-spin correlations and critical exponents and am wondering how to tackle this algorithmically. For example, in 1D, use an ...
0
votes
1answer
58 views

Probability density function of a particle for computation [closed]

I'm writing a program, part of which relies on a particle being able to change location similar to a how a real particle would behave (pardon my physics). For example, on a grid of 100x100, a ...
0
votes
0answers
39 views

Why are functional representations of systems important in physics or computational physics?

This was an addendum to a previous question I asked, but I figured I should make it it's own discussion. Assuming I am able derive a functional representation for any dynamical system (dissipative, ...
1
vote
0answers
176 views

Problems while numerically computing band structure using k.p theory

I want to use k.p theory to numerically compute the band structure of a bulk semiconductor. The band I like to include are the lowest conduction band (cb), the heavy-hole (hh), the light-hole (lh) and ...
3
votes
1answer
34k views

How does force relate to velocity

I had originally asked this question on math overflow and it was suggested that I ask it here. So I know that a force will change the magnitude of velocity if it is at an angle other that 90 degrees. ...
2
votes
3answers
83 views
2
votes
1answer
35 views

Hard processes in particle collisions

This is in reference to a MC generator called Pythia, but has wider applicability because the phrases are used widely in particle phyiscs, even though I have not managed to find a clear, ...
1
vote
0answers
25 views

Cellular automata and time T

Reading around about various CA models I saw that there are even those who are following the track to provide a single mechanicist model of the universe. For example Gerard van 't Hooft claims that ...
1
vote
3answers
111 views

Textbook recommendation for computational physics

Can anyone recommend some good textbooks for undergraduate level computational physics course? I am using numerical recipe but find it not a very good textbook.
3
votes
1answer
70 views

Calculating a two-dimensional orbital path with infinite granularity (non-Euler integration)

For a game I am making, I am trying to calculate the position of an orbiting object around one or more bodies. I have successfully implemented this gravity simulation by calculating the force, then ...
0
votes
1answer
23 views

How do we define capillary number in 2D (two dimensions)

In 3D the capillary number $Ca$ is defined as: $$Ca=\frac{\nu \rho U}{\gamma}$$ where $\nu$ is the kinematic viscosity ($m^2/s$), $U$ is the velocity ($m/s$)and $\gamma$ is the interfacial tension ...
3
votes
1answer
89 views

Why is the objective function used in Nudged Elastic Band method reasonable?

In Nudged Elastic Band (NEB) method, which is used to find reaction pathways when both initial and final states are know, an objective function is first constructed and then minimized to find reaction ...
0
votes
0answers
18 views

Why breakdown voltage decreases with increasing FWHM of an asperity?

I am running a simulation for a dielectric breakdown. I am getting breakdown voltage decreases with FWHM, which does not make sense as field enhancement is higher for a small width. Simulations solve ...
1
vote
1answer
27 views

Effects of terrain on outdoor sound propagation?

I'm looking for a computationally fast model for capturing the effects of varying terrain topology in far-field outdoor acoustic propagation, and was wondering if anyone had any insight or could point ...