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)

8
votes
0answers
476 views

Intuition for when the replica trick should work and why it works

I am a graduate student in mathematics working in probability (without a very good background in physics honestly) and I've started to see arguments based on computations derived from the replica ...
5
votes
0answers
144 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 ...
5
votes
0answers
659 views

Interpretation of Stiffness Matrix and Mass Matrix in Finite Element Method

I would like to have a general interpretation of the coefficients of the stiffness matrix that appears in FEM. For instance if we are solving a linear elasticity problem and we modelize the relation ...
4
votes
0answers
39 views

How to calculate dispersion relation from a Finite Difference (FD) wave simulation

I have a python code that calculates the solution of the inhomogeneous acoustic wave equation for a 2D medium with any velocity and source configuration. It was implemented using Finite Differences ...
4
votes
0answers
1k 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} ...
3
votes
0answers
70 views

Are there resources for simulating and/or theoretically describing solitons?

Recently there are striking new ideas emerging on "lower level" dynamics with respect to quantum mechanics involving fluid mechanics principles, including hints of soliton-like aspects to particle ...
3
votes
0answers
133 views

black body simulation

black body radiation is typically understood from Planck's argument of light resonance in a box, from which the density of states is computed. Now, suppose I want to simulate a black body ...
3
votes
0answers
69 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 ...
3
votes
0answers
86 views

Non-linear Wave Equation - Numerical Methods

Motivation: I'm working with a highly non-linear spherical wave-like equation (second order PDE). The equation can be written on the form $$\ddot{u} = f(t, \dot{u},\dot{u}',u',u'')$$ where ...
3
votes
0answers
147 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 ...
3
votes
0answers
867 views

Comparison of different ab-initio codes

One may find on the web a lot of different computational packages to perform "ab-initio" calculations of electron structure of the solids. Usually, the documentation is not quite transparent about the ...
2
votes
0answers
26 views

Open-source code for computing response functions

Summing Feynman diagrams to compute the response functions of a microscopic model is common in many areas of physics. While conceptually straightforward, the task can be computationally intensive. ...
2
votes
0answers
45 views

Numerical problem with Hartree-Fock equations for dilute Bose gas

I have to solve the following set of equations self-consistently: $$\begin{align} n_c(\mathbf{r}) & = \frac{1}{g}\left[\mu - V_{\rm ext}(\mathbf{r}) - 2 g n_{T}(\mathbf{r}) \right] \\[3mm] ...
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 ...
2
votes
0answers
45 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 ...
2
votes
0answers
243 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: ...
2
votes
0answers
59 views

How can the fictitious mass in the Car-Parrinello method reproduce the “real” dynamics?

In the Car-Parrinello method, to solve simultaneously the classical equations of motion for the atoms and the Kohn-Sham equations for the electrons, the following effective Lagrangian is used: $$ ...
2
votes
0answers
83 views

Modelling gravitational potential of a galaxy

I am interested in modelling the gravitational potential of a disc-shaped galaxy with radius $R$, i.e. solving the 2D Poisson equation numerically by Gauss-Seidel relaxation: $$\nabla^2 \phi = 4\pi G ...
2
votes
0answers
210 views

How to compute matsubara frequency summation over computer?

Matsubara frequency sum usually takes the following form: $S_\eta = \frac{1}{\beta}\sum_{i\omega_n} g(i\omega_n).$ But in my problem, $g(i\omega_n)$ is a lengthy expression, which can not be ...
2
votes
0answers
37 views

Coefficients and Parameters for contracted Gaussian basis sets

This is a repost from Chemistry.stackexchange in the hopes that someone here will be able to help me. Any help at all would be greatly appreciated. As far as I understand, an STO-NG contracted ...
2
votes
0answers
400 views

What is the drag coefficient of an open wedge?

To check my two dimensional CFD calculation I am looking for reference data on the drag coefficient of an open wedge. The geometry is shown below, together with the flow direction. I have found ...
2
votes
0answers
63 views

Trotter splitting and entanglement entropy

I have heard that a numerical solution to the Schrodinger equation using the Trotter splitting formula for a many-body Hamiltonian can cause an artificial increase in the entanglement entropy. I was ...
1
vote
0answers
20 views

Vortex in a Pipe / Solutions or approximations?

I want to study a vortex flow in a pipe; in other words helical flow. A couple introductory texts on fluid mechanics (e.g. Shames, Mechanics of Fluids) describe the solution to vortex flow in 2D: ...
1
vote
0answers
16 views

Efficient ways to solve complicated kinematics problems

I'm curious about the range of approaches one can take to perform motion/kinematic simulations on many 3D objects. Currently I am familiar with the iterative process of doing such a simulation, that ...
1
vote
0answers
20 views

Monte Carlo simulation of Breit-Wigner

I am using a Monte Carlo approach to integrate a scattering amplitude which has a propagator-squared term like $f(s)=\frac{1}{(s - 4)^2 + 1}$, between $s=0$ and $s=20$. Using a uniform random ...
1
vote
0answers
67 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
0answers
71 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 ...
1
vote
0answers
34 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
0answers
99 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
0answers
94 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 ...
1
vote
0answers
43 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
0answers
85 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, ...
1
vote
0answers
42 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 ...
1
vote
0answers
97 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 ...
1
vote
0answers
124 views

Numerical Solution of 1D Boltzmann Transport Equation

I need to solve the one-dimensional Boltzmann transport equation in a semiconductor numerically, and I want to take a deterministic approach toward the problem (i.e. not use Monte-Carlo or similar ...
1
vote
0answers
274 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 ...
1
vote
0answers
58 views

Does Planck's constant imply limits to computing *results*

... I don't mean quantum effects limiting hardware fabrication sizing. Such small scales have for some time exhibited issues. Rather, along the lines of imagining the smallest possible divisions of ...
1
vote
0answers
41 views

Computer coding of Perdew Burke and Ernzerho (PBE) method

I have been assigned the task to write a computer code to implement the Perdew Burke and Ernzerhof (PBE) method. Does anyone know a good reference which can make the coding of this method easier?
1
vote
0answers
137 views

Lattice Gas Cellular Automata - HPP model square lattice

In the HPP model of LGCA, a square lattice is used and there is only one collision configuration as mentioned in figure (taken from the book Lattice Gas Cellular Automata and Lattice Boltzmann models ...
1
vote
0answers
580 views

Hysteresis Curve and how to implement it using Preisach model or other models

As a homework I need to draw a hysteresis curve (preferably an interactive one) using Matlab or any other programming language. The problem is I have trouble finding a good algorithm to do so. I need ...
1
vote
0answers
79 views

Application of wavelets in computational electromagnetics

I am a first-year graduate student (electrical eng.); wavelets is one of the courses that I am taking presently. As part of the course, I plan to explore application of wavelets in computational ...
1
vote
0answers
125 views

How does the Lennard Jones Potential changes for interaction between particles of different sizes?

I am interested in incorporating a Lennard-Jones potential in a simulation. When the interaction only involves the same type of particle, with same characteristics, we can use reduced units, scaling ...
1
vote
0answers
120 views

Gerstner Wave Formula to Vertex information $(x,y,z)$

I have created a program, that draws a plane mesh, and allows to me edit or animate all the separate vertices along it. The idea behind it, is for me to be able to create a gridMesh based wave ...
1
vote
0answers
605 views

Is Bremermann's limit redundant with Landauer's limit for all practical purposes?

Trying to understand the physical limits to computation, I notice that among these we have two types of limits that constrain the minimum allowable energy for a computation. Limits that constrain ...
1
vote
0answers
5k views

How to calculate the velocity of fluid at the outlet when density and the pressure drop are known?

I have a U- like pipe. Its inlet has atmospheric pressure $p_o=10^{5} \, Pa$. Vacuum is applied to the other end with a pressure gradient $\nabla p_v=-30 \cdot 10^{3} \, kPa/s$. The total time of the ...
1
vote
0answers
281 views

2D quantum well energy spectrum (analytical vs numerical)

I am trying to understand the energy spectrum difference between the analytical and the approximated solution for a quantum well. The particle is inside a box with domain $\Omega=(0,0)$X$(1,1)$. For ...
0
votes
0answers
14 views

Numerical simulation of strong field gravitational vacuum solutions colliding

I am interested in the current state of knowledge of strong field General Relativity learned from numerical investigations of gravitational wave packets colliding with each other or black holes. If ...
0
votes
0answers
19 views

photometric redshifts - the origin of 1 / (1+ zspec) in photometric accuracy

i have a training sample of sources with known spectroscopic redshifts ($z_{\rm spec}$) and i run a photometric fitting code that uses a template to determine the photometric redshifts ($z_{\rm ...
0
votes
0answers
11 views

Explicit solvent simulations in NVE ensemble: Molecular Dynamics

I am simulating a SDS micelle system in NVE ensemble at room temperature. I am not sure if the system is equilibrated, since the pressure variations are around 100's of atm, where the average pressure ...
0
votes
0answers
21 views

Books on Spectral Jacobian

Can anyone point me a book to better understand "spectral Jacobian". I am studying X-parameters, I have the book X-Parameters: Characterization, Modeling, and Design of Nonlinear RF and Microwave ...