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

8
votes
0answers
433 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
52 views

Numerical Solution of the Convection 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
123 views

The use of Artificial Intelligence in physics resarch

I have been reading about a machine that observed a double pendulum and created equations that both described its motion and associated conservation laws. The authors claims: We have developed a ...
5
votes
0answers
599 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 ...
3
votes
0answers
121 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
64 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
844 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
140 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
711 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
125 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
43 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
143 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
34 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
319 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
60 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
34 views

Boundary conditions for solving Poisson's Equation with Experimental Data

I want to numerically (with Matlab) solve Poisson's equation : $ \frac{\partial^2u}{\partial x^2} + \frac{\partial^2u}{\partial y^2} = f(x,y)$ On a rectangular domain using experimental data. From ...
1
vote
0answers
22 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
59 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
77 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 ...
1
vote
0answers
45 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
33 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
62 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
203 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
56 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
35 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
368 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
74 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
105 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
110 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
530 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
59 views

Find out ground sates for large 2D classical spin model

Reaching the ground state of a large 2D classical spin model (e.g. classical Heisenberg model) might be a relatively difficult task while using conventional "flip/reject" Monte Carlo method. The ...
1
vote
0answers
186 views

RGEs of the MSSM - problems with Mathematica

I'm having some troubles with the trilinear soft couplings of the MSSM RGEs. I've used the ones written in Martin's supersymmetry primer and I run them using mathematica, if I do so without taking ...
1
vote
0answers
3k 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
246 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
15 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 ...
0
votes
0answers
47 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 ...
0
votes
0answers
15 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 ...
0
votes
0answers
56 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
0answers
36 views

Using a mask with a (2D) Fourier Transform

For a project I'm working on, I (1) take a (gaussian) random field of 256x256 pixels, (2) apply numpy's fast fourier transform (numpy.fft package), (3) filter the k-space image (using a Tukey ...
0
votes
0answers
19 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
0answers
40 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
29 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
0answers
32 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 ...
0
votes
0answers
57 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 ...
0
votes
0answers
45 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, ...
0
votes
0answers
19 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 ...
0
votes
0answers
27 views

Optimize locations of vertices and edge length using physics without simulation

I have a directed graph with set of vertices and edges. Each vertices don't have a position coordinate to be displayed on 2D computer screen. I want to add unique position coordinate to each vertices ...
0
votes
0answers
35 views

Adams-Moulton and BDF methods

1.What are the differences between Adams-Moulton and BDF methods. Which one is better and which one computes the solution faster? I think Adams-Moulton is a better method as it can get to the ...
0
votes
0answers
96 views

How to calculate a Mooney-Rivlin material

Well, I'm working on a fsi (fluid-solid/structure interaction) model. I've got ready the fluid part, now I need to get ready the solid part. Since I'm working with tissues, I need to use an ...
0
votes
0answers
41 views

need some references in molecular dynamics simulation

I Just start we a new model on my project and I will need to work with Molecular dynamic simulation. I have to simulate dynamics of many filament inside a network (actomyosin network).The problem is ...