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)

0
votes
0answers
24 views

Mixed Neumann/Dirichlet BCs in Poisson-Schrodinger self-consistent solver

I'm currently in the process of writing a self-consistent Schrodinger-Poisson solver for a device heterstructure (High Mobility Electron Transistor). I am having a few difficulties however. It seems ...
0
votes
0answers
26 views

Partial Derivative in FeynCalc [migrated]

I try to compute some loop diagram in mathematica with FeynCalc package. In the nominator, I have a huge trace and it shows the right solution but in the next step I need to find the solution of a ...
0
votes
0answers
17 views

How can use genetic algorithm to regression curve? [on hold]

i have a equation $f(z,a,b,c)=\sqrt{a+b*z+c*z^2}$ and have 20 experiment points $f(z),z$. I tried least square method to fitting $a,b,c$. Its good, but i want to try other method. It's call genetic ...
0
votes
0answers
11 views

Computing the dynamics with a bilateral constraint in position

Let's says I have a dynamical system of generalized coordinates $x$, subject to a bilateral constraint $g(x(t),t)=x_1(t)-d(t)=0$ where $d$ is a prescribed displacement. I am looking for the ...
0
votes
1answer
63 views

Finite Difference Method (FDM) solution to heat equation in material having two different conductivity

I am not from a mechanical engineering background and I have not taken any courses in PDE so this may seem trivial for many. and I am writing a Matlab code with the objective to solve for the ...
0
votes
1answer
39 views

Laser resonator modes

For a resonator, plot the intensity as a function of time for (i) when the modes have identical phases, and (ii) when the modes have random phases. Attempt Here is the equation for intensity (p. ...
1
vote
0answers
27 views

Numerically integrate equations of motion with noise

Consider the equations of motion $$\begin{cases} \dot{x}(t) & = v(t) \\ \dot{v}(t) & = -\frac{\lambda}{m} v(t) + \frac{1}{m}F^{c}(x(t)) = a(x(t), v(t)) \end{cases},$$ where $x$ is the ...
0
votes
0answers
26 views

Keplerian orbits and radial velocities: problem with true anomalies

I am writing a pure Python code to calculate radial velocities of stars with massive objects orbiting them. I am following the formalism from Murray & Correia (2011) (hereafter MC11). In a given ...
0
votes
0answers
30 views

Generalized Eigenvalue Problem from linear stability analysis

I am trying to solve a generalized eigenvalue problem raised by linear stability analysis $$AV=\lambda BV.$$ $A$ and $B$ are non-symmetric complex valued matrices. The set of equations I am trying to ...
0
votes
0answers
18 views

2-D Ising model temperature of equilibrium [duplicate]

I have simulated 2-D Ising model in C#, but my scatter plot for magnetization versus temperature is an decreasing line (you know that the correct plot decreases quickly in a temperature near 2.3). ...
0
votes
0answers
24 views

Why is more difficult for a DFT to converge with non-periodic conditions?

I am simulating a cubic amorphous silica system. When all boundaries are periodic the SCF iteration converges rapidly. However, when the boundary is not periodic in the $x$ direction, the SCF ...
0
votes
2answers
74 views

Use Runge Kutta method to solve schrodinger equation

The schrodinger equation in spherical coordinates after seperation of variables as a solution of hydrogen atom is given by $$ \frac{-\hbar^{2}}{2 \ m} \left[ \frac{1}{r^{2}} \frac{d}{dr} \left(r^{2} ...
2
votes
0answers
19 views

Derivation of generation of time between two subsequent particle enters into computational domain

I have problem with understanding derivation of one equation in following problem. You have 1D computational domain (it is not 1D but because it is symmetrical and we are watching only radial ...
1
vote
1answer
83 views

Simulation of relativistic probe passing through an external solar system

I recently read about the Breakthrough Initiative to launch "StarShot", a nano-probe that is designed to travel to Alpha Centauri at $0.2c$. One of the challenges to be solved involves the precise ...
0
votes
0answers
25 views

Matrix representation of the radial Laplace operator isn't symmetric (or hermition as a result)

I'm working with the cylindrical coordinates. I'm using the central difference to convert the radial part of Laplace operator into a matrix. $\nabla^2 u = \frac{\partial ^2u}{\partial ...
4
votes
1answer
72 views

Can I re-wind gravity?

Suppose I perform an arbitrary simulation where I integrate the motions of a collection of particles which interact only gravitationally. Suppose I use a time reversible integrator (to be specific, ...
4
votes
1answer
67 views

How to conserve momentum when numerically-integrating the path of a charged particle through a magnetic field

I have the python script attached below, which is intended to track the trajectory of a charged particle in a static, uniform magnetic field. It is very simple to calculate the instantaneous force ...
6
votes
1answer
79 views

FCC-to-BCC phase transition in NaCl, Buckingham or Lennard-Jones potentials?

Background The transformation from B1 (face centered cubic (FCC) type) to B2 (body centered cubic (body centered cubic (BCC) type) structures is one of the best documented high pressure phase ...
1
vote
0answers
24 views

How is time evolution done in numerical GR?

Suppose we're simulating what happens when a fairly massive object falls into a black hole. Say the object starts far away, so that the initial condition is that the metric is the Schwarzschild metric ...
0
votes
0answers
14 views

Fermi-level of lattice model

I have a chain of $N$ lattice points with periodic boundary conditions. Every lattice point only has one orbital for an electron to occupy and spin is not included for simplicity. The lattice points ...
5
votes
1answer
66 views

Methods for handling close approaches in $N$-body simulations

In direct gravitational $N$-body simulations, what are the preferred methods for handling close approaches between bodies in order to preserve the accuracy of the evolution of the system?
1
vote
0answers
50 views

Simulating TE011 mode in cylindrical cavity, Comsol [closed]

I have been trying to simulate TE011 mode in a cylindrical cavity by modifying the example given in Comsol ...
1
vote
0answers
24 views

Choosing the boundary conditions knowing the potential

So I have to apply Numerov's algorithm to solve the Schrödinger equation in spherical coordinates using the potential $V(r)=-A e^{-r/a}$, where $A=32.7MeV$ and $a=2.18fm$. I have already ...
1
vote
1answer
72 views

Schrödinger equation in spherical coordinates using Numerov's method?

So I need to solve the Schrödinger equation using Numerov's method. However, the problem I have assigned to me is in spherical coordinates. As it seems, I only need to take into account the dependance ...
0
votes
0answers
16 views

Spherical Harmonic fitting excess polar magnitudes

I am trying to fit an expansion of spherical harmonic functions to a dataset distributed over the surface of a sphere using the least squares method. Each data point is in terms of ...
2
votes
1answer
44 views

Current Status of the Monte Carlo Sign Problem

I've been reading about the Monte Carlo sign problem, and I am a little confused about its current status. Specifically, after reading this post When is the "minus sign problem" in quantum ...
1
vote
0answers
30 views

Calculate 2D Effective mass from bulk effective mass

I am trying to create a self consistent Shrodinger Poisson Solver for various semiconductors. There is already one done by Professor Hu from UC Berkeley - QM CV Simulator. Looking at the code, they ...
0
votes
1answer
50 views

Numerically finding the energy diagram of the hamiltonian

I'm looking at a collection of three two level systems (qubits) coupled to each other (with known bare state energies and couplings). The hamiltonian is given by ...
0
votes
0answers
19 views

Standing waves and acoustic levitation

i am doing a new personnal project about acoustic levitation with standing waves into ambiant air using a piezo ultrasonic tranceiver and a reflective pannel. I modelised both of the emission and the ...
0
votes
0answers
21 views

Plotting EM fields

I am currently studying transformation optics and learning about electromagnetic cloaking. To plot my results(Electric and magnetic fields) i have to solve maxwell's equations given the permeability ...
1
vote
1answer
34 views

Simulating the Electric Field inside a Hollow Sphere of Uniform Charge

I've recently been studying Gauss's Law and have come across some results that I want to verify. For one, I am trying to verify that the electric field inside a conductor (in this case it could be ...
0
votes
1answer
84 views

Quadratic drag projectile motion

I have calculated formulas with 1 dimensional trajectory motion (free-fall) including quadratic drag, and have created the following equations. These equations of motion are not of much use on its ...
0
votes
1answer
14 views

Which are the best books I can refer for turbulent heat transfer and its computation?

I am a graduate student. I wanted to do computation on heat transfer problem mainly heat convection. I am working on turbulent regime. Can you kindly suggest some good books so that I can get ...
0
votes
0answers
64 views

Which formula was used to get Gravitational Waves animation?

I would like to make animation like this by myself It is from wikipedia article "Gravitational wave". But if I use formula from wikipedia article "Gravitational wave" subsection "7.2.2 Wave ...
1
vote
1answer
35 views

N-Body Simulation for Photons

For my undergrad optics class I am taking we have to select final projects to work on for second half of the semester. I have spent the last few weeks working on my own project which was a GPU ...
0
votes
1answer
45 views

Review article recommendation in the field of 2d materials

I am new to 2d materials. I tried to search related review articles in review of modern physics, but did not find anyone covering the whole of the 2d material area. Anybody can recommend some latest ...
2
votes
1answer
165 views

Is Machine Learning a useful tool for analyzing experimental data? [closed]

Apart from being a student of Physics I'm also a developer and recently in the context of cloud development I've been hearing a lot about cloud based Machine Learning. After looking at some examples ...
4
votes
1answer
77 views

Do black hole merger simulations include regions inside event horizons?

Inspired by this question, I would like to ask the following specific point. In numerical simulations of general relativity that involve black holes, like the ones used to understand the black-hole ...
0
votes
2answers
72 views

About unit definition in numerical computation [closed]

I attend an undergraduate level course of quantum mechanics. I am thinking of computer computation on a quantum mechanical problem. I do not have much experience in numerical computation, so I read ...
1
vote
1answer
48 views

Calculating orbits with varying laws of physics? [closed]

The general way planets orbit stars in our universe is determined by the inverse square law and the three laws of motion. Both the inverse square law and the three laws of motion are the reason that ...
0
votes
0answers
23 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
23 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 ...
1
vote
0answers
35 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: ...
2
votes
2answers
70 views

time step in kinetic montecarlo simulation

I have a doubt about how the time step is calculated in kinetic montecarlo simulations. One state with index $i$ is connected to other $N$ states, indexed by $j=1...N$, by transitions that happen ...
0
votes
0answers
17 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 ...
1
vote
0answers
22 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 ...
0
votes
0answers
24 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 ...
0
votes
0answers
26 views

Initial guess for calculating the Universal Variable Formulation of the two-body Kepler problem

I am trying to simulate orbits using the Universal Variable Formulation of Keplers two-body problem, using the algorithm provided by Howard curtis in Orbital Mecanichs for Engineering Students as per ...
2
votes
0answers
30 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. ...
1
vote
1answer
45 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 ...