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
1
vote
0answers
20 views

What is the best numerical way of solving this Schwinger-Dyson equations?

I have to solve these equations numerically \begin{equation} \begin{cases} G(\omega)=\frac{1}{-i\omega-\Sigma(\omega)} \\ \Sigma(\tau)=J^2G(\tau)^2G(-\tau) \end{cases} \end{equation} These are the ...
0
votes
0answers
29 views

Comparison of room acoustics simulation techniques

What is the state of the art technique for the numerical simulation of the propagation of a sound wave inside a room (room acoustics) regarding the time domain? I am specifically interested in the ...
0
votes
0answers
22 views

Verifying completeness of constants of motion

I can find constants of motion by looking at the null space of the Poisson Bracket operator $ \{H, \cdot\} $ over a polynomial space by brute force with symbolic algebra (code). This scales terribly ...
2
votes
0answers
22 views

Modal filtering with Vandermonde matrix - Artificial viscosity

I am trying to implement a 2D shock detector for an artificial viscosity model to control strong nonlinearities in compressible fluids. The method I am relying on is originally from: Yu, M. L., ...
0
votes
0answers
52 views

Inexpensive nuclear reaction simulation tools

I am interested in seeing if there are any simulation software programs for less than $50 that may be used to simulate how neutrons (energies between thermal and 15 MeV) interact with certain nuclei ...
0
votes
0answers
55 views

Clarification on DMRG computational complexity

I was reading a paper on the density matrix renormalization group (https://arxiv.org/abs/1008.3477). In DMRG, we gradually grow a chain by inserting a unit cell at the center of the chain (for ...
-2
votes
1answer
46 views

Finding free orbit formula by cartesian coordinate

As a simple idea, I want to find a simple orbital function using calculus rather than using the law of Kepler. However, using the recurrence relation is too difficult to solve. Help me to solve...
0
votes
1answer
50 views

Sampling Maxwell-Jüttner distribution for non-unity mass and speed of light

I am trying to sample Maxwell-Jüttner distribution using the Sobol method as described in Zenitani Loading relativistic Maxwell distributions in particle simulations (2015). Equation (2) in the paper ...
2
votes
0answers
20 views

How do I know that the time-step I'm using in my Verlet integrator is small enough?

I am performing a molecular dynamics simulation of a many-body system using the Verlet algorithm in an implicit solvent, i.e. adding noise to the equation of motion via a thermostat term. I want to ...
1
vote
0answers
53 views

Calculating the local energy in neural network quantum state

given a Hamiltonian of Heisenberg 1D model as following: $$H = -J\sum_{I}\sigma_{i}^{z}\sigma_{i+1}^{z}$$ I am trying to solve it with a neural network function given by Restricted Boltzmann machine ...
2
votes
1answer
60 views

Is there a program or a website able to perform all Wick contractions for a given expression?

Imagine I have an expression of the type: $$\langle \phi_{x_1} \phi_{x_1} \phi_{x_2} \phi_{x_2} \phi_{z_1} \phi_{z_1} A_{z_1} \phi_{z_2} \phi_{z_2} A_{z_2} \phi_{z_3} \phi_{z_3} A_{z_3} \phi_{z_4} \...
-1
votes
1answer
25 views

Tracing Path of a Body and Uncertainities Involved [closed]

If for instance suppose we take a situation like this- A satellite just exploded due to some reasons and it's sharpnel are falling in Earth. We have to find out the exact location where they will ...
0
votes
0answers
45 views

Is QCD computable?

Are arbitrarily good approximations of the time evolution of any QCD system, given initial and boundary conditions, Turing computable? Can lattice QCD simulations be used to do so in theory? How ...
1
vote
0answers
23 views

Correct method for error propagation in pseudo-data approach using model

I want to compute the error on the integral over experimental data $(y_i,\delta y_i)$ in a pseudo-data approach using model ($x_i$ are assumed to be known error-less numbers). If I were using direct ...
0
votes
0answers
16 views

How to calculate optical conductivity from numerical eigenstates of tight-binding model?

Let's say we have a 1D spatially inhomogeneous tight-binding model that does not have momentum as a good quantum number. We can numerically diagonalize it to get the spectrum and eigenstates. But how ...
-1
votes
0answers
35 views

Software to simulate molten salt flow and thermodynamic operations

I was curious if there was any software (preferably in C++, Java, and/or python) that could be used to simulate the following: Heat capacity of a fluid Heat transfer through a liquid and a solid ...
0
votes
1answer
30 views

Optical Logic Gates using Polarisers [closed]

I am trying to build an optical logic circuit using the horizontal or vertical orientations of light as the two binary states and laser as the source. Before I ask you my question I wanted to first ...
0
votes
0answers
43 views

(Finite Element Method) Zero thermal conductivity and Non-zero convection boundary condition

I asked this same question in Computational Science Stack but no one proved it mathematically (see the link below). Link of original question: https://scicomp.stackexchange.com/questions/32959/fem-1d-...
1
vote
0answers
20 views

Fluid rendering technique (C++/3D) [closed]

I'd like to render fluid based on density value. Let's say I have a cube with 8 points (each of thoses points has its proper density value). So, If i want represent the fluid form on this cubic ...
1
vote
0answers
38 views

Numerically solving unbounded Fokker-Planck equation

I am wanting to solve a 3D Fokker-Planck equation of the form: $$\partial_{t}p(\mathbf{x}, t) = -\nabla \cdot \mathbf{J}$$ where $\mathbf{J} = \mathbf{v}(\mathbf{x})p(\mathbf{x}, t) - D\nabla p(\...
0
votes
0answers
29 views

What methods can be used to interpolate values of a gradient?

Let's say I have a cubic mesh, the points of which I have obtained a vector value for, thus forming this discrete vector space of sort, I'm not really sure of the proper terminology for this which ...
1
vote
0answers
11 views

Finit element modelling of rigid bodies and machining

I would like to computer model of material machining. Is finite element metod suitable for this? I am imagining 3D grid of voxels, which are tracking matter properties across space. Is this correct ...
0
votes
0answers
40 views

Quantum Monte Carlo Loop Algorithm for quantum spin: why is the freezing graph present in ferromagnetic Ising model?

I study the loop algorithm (Evertz et al). I cannot understand, why the freezing graph type where we have to flip all 4 spins together is not present for the quantum-XY model and the anti-/...
2
votes
0answers
54 views

Finite differences as a variational method

I think I should be able to derive finite differences for the Schroedinger equation by starting from a variational method. More specifically: in finite differences we approximate $$\{\psi(x),x\in\...
-1
votes
1answer
70 views

Modelling lunar motion in Python [closed]

I'm trying to solve the equations of motion of the Moon in the Moon-Sun-Earth system (equations 1139-1141 here:http://farside.ph.utexas.edu/teaching/336k/Newtonhtml/node130.html) My initial ...
0
votes
1answer
80 views

How can I use Runge-Kutta4 to solve this orbit in polar coordinates?

I want to simulate a simple orbit of a planet moving around a star which is fixed in position. I have formulated the ODEs for this problem using Lagrangian Mechanics and have found the equations of ...
2
votes
1answer
26 views

Scattered fields in surface plasmon excitation with grating

We need to use prism-coupling (e.g. Kretschmann or Otto configuration) or use grating to excite surface plasmon between metal and dielectric. What I understand is that when we use a grating or ...
2
votes
0answers
18 views

Combined effects of strain and doping on the electronic structure of semiconductors

I've read papers in DFT studies on band gap tuning in semiconductors and the usually studied methods are either through doping or application of external strain. But its always just one or the other. ...
1
vote
0answers
32 views

Lippmann-Schwinger Equation (Numerical Method)

I recently tried to code up the numerical solution for solving the LS equation, which is based on the "scatt.py" code from R.H. Landau's "Computational Physics". The formulation is kind of like: $$ R_{...
0
votes
0answers
13 views

What are possible causes that can lead a self consistent calculation to diverge in DFT?

I am extending a code written to do self consistent Density Functional Theory calculations to the case of spin polarized systems. Due to the modifications, the calculations are leading to diverging ...
1
vote
1answer
64 views

Discontinuities in a Poincare map for a double pendulum

I'm generating poincare sections of a double pendulum, and they mostly look okay, but some of them have weird discontinuities that seem wrong. The condition for these sections is the standard $\...
3
votes
2answers
102 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 ...
0
votes
1answer
33 views

Metropolis Algorithm Transition-Proposal Probability

I'm working my way through a short section on the Metropolis algorithm in the lecture notes on Computational Quantum Physics by Prof. Troyer. However, I am not sure what probability distribution was ...
0
votes
0answers
49 views

Symplectic time-integrator for Fokker-Planck equation

Is there a way to use a symplectic time-integrator for the numerical solution Fokker-Planck equation of the form: $$\partial_{t}p = -\nabla \cdot (\underline{v}p - D\nabla p)~?$$ where $p = p(\...
1
vote
1answer
50 views

Physics Graphic Simulations

I am interested in modeling physics using graphics simulations. I took a python course with CS majors and one of the assignments was to simulate a planet orbiting a star. That being said, I was ...
1
vote
0answers
23 views

Why do the Binder Cumulants of different system sizes intersect at the critical point?

When Monte Carlo simulations are performed for spin models (Ising model etc.) the critical temperature can be found by simulating for different lattice sizes and plotting the Binder Cumulant for them. ...
1
vote
0answers
38 views

How do I calculate the pressure on the walls of my MD simulation?

So I've built a simulation in vpython where I have multiple hard balls that collide with each other and with the walls of the cube I've placed them in. I then replaced the values for their radius and ...
1
vote
0answers
22 views

Density of states LL in graphene [closed]

I am using the Kernel Polynomial Method to determine the spectral density of a 2DEG system that has been sujected to a perpendicular magnetic field B. I wish to determine (a) What the amplitudes of ...
1
vote
1answer
42 views

Numerical calculation of skyrmion number

I'm having trouble finding a numerical method to calculate the skyrmion number for a certain system. I want to calculate $$ n = \frac{1}{4\pi}\int \textbf{M} \cdot \left( \frac{\partial \textbf{M}}{\...
0
votes
0answers
43 views

$<P_T>$ distribution as function of PL

I have a conceptual problem when working with some Monte Carlo generated data. I have to simulate the $P_T$ and $P_L$ distribution for $N$ particles through Monte Carlo methods as follows: $$\theta \...
0
votes
1answer
18 views

Obtain element boundary forces in Finite Element Method

Suppose that we have a rectangular domain discretized using 16 voxels (4x4 voxels) and I apply a Finite Element Analysis with arbitrary boundary conditions in the domain. As a result, I obtain the ...
0
votes
0answers
25 views

Winner-takes-it-all algorithm in particle decays

How does the winner-takes-it-all algorithm work for two competing processes, specifically in particle decay processes? What are its advantages against veto algorithm, for example?
-3
votes
1answer
221 views

How do we find the equation for the gyrating motion of a particle in a uniform magnetic field and a non-uniform Electric field? [closed]

Considering the gyrating motion is not negligible and also retaining the guiding center drift, how do we get the trajectories x(t),y(t),z(t) of the particle? In this case is the variation in the ...
1
vote
0answers
22 views

FDTD simulations for metamaterials and complex refractive index

I have designed a working FDTD code to simulate electromagnetic wave propagation using Yee's algorithm. I have been able to place materials such as glass slabs, etc which have positive and real ...
1
vote
0answers
24 views

About likelihood ratio computation

My head is stuck into a rather simple problem I don't know why. I'll explain: 1- bosons have an angular distribution in $\theta$ and $\phi$: $$F(\theta, \phi) = \frac{3}{4\pi}[0.5(1-\alpha) + (0.5)(3\...
0
votes
1answer
61 views

Double pendulum Poincaré section issues

I'm generating high-resolution Poincaré sections of a double pendulum, and I'm running into some issues in creating the initial conditions for a given section. In general, I describe my pendulum with ...
1
vote
1answer
44 views

Analyticity of the initial data

Consider an initial data $(q_{i j},K_{i j})$ on a Cauchy hypersurface $\Sigma$ in a globally hyperbolic spacetime. How reasonable is it to impose the condition that these functions are analytic? I am ...
0
votes
1answer
37 views

Program for simplifying expressions of many-body Pauli operators

I have a system of $N$ spins, with the usual $\sigma_i^{x,y,z}$ operators defined on each site (operators at different sites commute). I would like to find a program, perhaps in Mathematica, with the ...
2
votes
1answer
75 views

Choosing initial condition for Hamilton-Jacobi PDE from initial $x$ and $p$

For separable solutions to Hamilton-Jacobi PDE (say in 2D), we treat the Hamilton's principal function $S$ as $$S= W(x) + W(y) - E*t$$ and treat the separate parts as constants and find $W(x)$, $W(y)$...
0
votes
1answer
105 views

Master equation for the mechanical modes

Consider the standard model of optomechanical systems with a single optical cavity mode coupled to a mechanical oscillator, which is canonically modeled as a FP cavity with one fixed mirror and one ...