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
42 views

What is wrong with my formulas for a mathematical model of a double pendulum? [on hold]

I tried mathematics but I've had no response in a month. I wanted to create a computer simulation on Matlab, using a model for a pendulum from this study (A double pendulum model of tennis strokes. ...
1
vote
1answer
38 views

Single slit diffraction simulation (mathematica) [on hold]

I'm trying to simulate the Fraunhoffer diffraction at slits(single,double,triple) with Mathematica. In the picture, the red one is analytical result and the green one is numerical result. The ...
2
votes
0answers
19 views

Boundary conditions and density when using smoothed particle hydrodynamics/corrective smoothed particle method for two-temperature models

I am trying to solve the dual-hyperbolic two temperature model (TTM) using a modified smoothed particle hydrodynamics (SPH) scheme known as CSPM (corrective smoothed particle method), as demonstrated ...
0
votes
0answers
22 views

Efficient weak form for numerically generating a static magnetic field from arbitrary current distribution (FEM)

This may stray into computational physics, but I'm curious to know if there is some more efficient way for what I am trying to do. Perhaps a better gauge formulation is needed. Consider two ...
6
votes
1answer
180 views

Potential function - numerical simulation

Using MATLAB, I fixed the potential in a region inside a rectangular plate (100 V) and in the border (50 V). I got the following result of the potential along the plate: I can't find an intuitive ...
-1
votes
0answers
26 views

SU(2) Coherent State and Probabilities

The SU(2) Corherent State for a for more than 1 two level system is given by: $$ | \eta, J \rangle = (1 + |\eta|^2)^{-J} \sum_{m=-J}^{J} \sqrt{\binom{2J}{J + m}} \times \eta^{J + m} |J, m \rangle, ...
-2
votes
0answers
26 views

Imaginary part in calculation [closed]

I want to simulate (I use mathematica) a matrix that has imaginary number. Given the matrix, $$A = \begin{pmatrix} e^{2I+3}Sin{(2\Pi t)} & e^{I (\Pi)}Cos{(2\Pi t)} \\ Cos{(2\Pi t)} & -e^{...
0
votes
2answers
289 views

Modelling multiple heat sources?

I am working on a project to model heat transfer in a data center. I have basic models for a one dimensional heat transfer using conduction and convection. Such as, taking a server and modelling the ...
15
votes
4answers
328 views

Why don't hexapods gallop?

I haven't observed any terrestrial hexapods that gallop like agile tetrapods(cheetahs/horses). Instead, I've observed that all of them scuttle which constrains their maximum velocity. And after ...
-3
votes
0answers
26 views

Book on galaxy formation and theoretical astrophysics for physics graduate [closed]

I am looking for a graduate level book on galaxy formation and theoretical astrophysics and also numerical relativity for a graduate student who has completed quantum field theory and general ...
2
votes
2answers
65 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
17 views

Capturing superfluid condensation with exact diagonalization

Doing exact diagonalization on bosonic systems is tricky, because the possibility of multiple occupancy means that even the single-site Hilbert space is infinite-dimensional. It's my understanding ...
-2
votes
0answers
20 views

What should I study if I am interested in creating innovative gadgets? [closed]

I’m interested in gadgets and I know it requires programming skills. But on the other hand I really like experimental physics. So, I don’t know how to combine these science areas together.
-2
votes
0answers
12 views

Newton's law - satellite position (in 2D) for specific moment of time [closed]

I’m trying designate function that calculate satellite position (2 dimensional) for specific moment of time for better understanding Newton's law of gravity. To verify my calculations I created second ...
5
votes
0answers
112 views

Computational scaling of quantum and classical Monte Carlo algorithms

How does the computational complexity of finding an equilibrium thermal state for a given Hamiltonian at a given temperature scale with system size under classical and quantum Monte Carlo? I know ...
3
votes
0answers
42 views

Calculating 2 particle Partial Trace for Density Matrix in Zeeman basis for a large number of Spins

I want to trace out all spins but 2 from a density matrix in the zeeman basis for N spins. For N=3 for example I would have the basisvectors: $ |S=1.5, m=-1.5\rangle =|000\rangle, |S=1.5, m=1.5\rangle ...
-6
votes
0answers
50 views

Please explain this project to me [closed]

I've been struggling with this project for a few months now and I frankly can't break out of my way of thinking. The questions I'm struggling with understanding are written down below. The questions ...
2
votes
1answer
43 views

Meaning of $\nabla_{\mathbf{p}_k} W(\mathbf p, h)$ in PBF

I'm reading this paper on Position Based Fluids and I couldn't understand the meaning of $\nabla_{\mathbf{p}_k} W(\mathbf{ p_i - p_j}, h)$ in the equation 7 (see below). …the gradient of the ...
-3
votes
0answers
27 views
0
votes
0answers
20 views

Young's modulus graphene nanoribbons [closed]

How to calculate volume and young's modulus in graphene nanoribbons? How to calculate second derivative Total energy on the strain? Band structure of graphene nano-ribbons is the code below: ...
-2
votes
2answers
38 views

Data transmission in optical fibers

Im reading though my Network Architecture lecture notes and it states optical fibers need a repeater every >10km. Repeaters are used to amplify digital signals but the notes do not explain or state ...
0
votes
1answer
43 views

Heat transfer in solidification of cylindrical billet

I am trying to plot IN MATLAB the temperature profile of solidification of metal in manufacture of cylindrical billets. The metal is solidifying and the billet moves downwards as a constant water ...
2
votes
1answer
89 views

Gravitational field of thin 2D ring - numerical simulation

I'm aware of Newton's Shell Theorem, which states that inside of a thin ring of uniform density, the gravitational force exerted on a point mass should be zero. I wrote a quick field simulation to ...
0
votes
1answer
127 views

Is there any way to find difference between two same sounds of different people?

I am trying to understand and find a way to distinguish two same sounds of different people by some physics formula, so could you guys help me? OK I'll try to explain my question in this way that, ...
7
votes
2answers
2k views

Should acceleration be included in state vector of a Kalman filter?

I'm developing (actually adopting existing solution) a Kalman filter to model motion of a vehicle (UAV or automobile). The state vector will include position, velocity, and, possibly, acceleration. ...
3
votes
2answers
566 views

System that is more efficient to simulate than to compute

I know some physics calculations require alot of computing to solve. Are there currently systems where it is more logical to do an experiment to "let the universe simulate the experiment for us" and ...
0
votes
1answer
52 views

Is Schmidt Decomposition well defined for periodic matrix product representation(MPS)?

We are able to perform Schmidt decomposition for open-boundary MPSs with dimension of boundary bond $m_0 = 1$, $|\psi\rangle=\sum w_{a_l}|a_l\rangle_L|a_l\rangle_R$. Because we can make $|a_l\rangle_{...
0
votes
1answer
120 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
1answer
122 views

Molecular dynamics and detailed balance

In developing methods to perform Monte Carlo simulations one sufficient condition to preserve the stationarity of the target probability distribution is to impose detailed balance i.e. [Gardiner page ...
1
vote
1answer
318 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 ...
1
vote
1answer
145 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
73 views

How to plot numerically the wave functions according to the Hamiltonian?

It is often difficult to analytically solve the Schrodinger equation, and so we need to obtain a solution numerically. An example plot is shown below. Here, the wave functions for a three junction ...
3
votes
2answers
132 views

Numerical Tools to find Braiding Statistics of Quasiparticles

While certain classes of systems that exhibit topological order can be solved exactly (such as the Toric Code, Abelian FQH Edges, etc.) there also exist systems (think of perturbed versions of the ...
0
votes
1answer
247 views

Tight binding Hamiltonian for 2D finite dimensional lattice and nanowire

The Hamiltonian of a 1D lattice having finite $N$ atoms, (if we consider one basis per atom) is given by the following $N\times N$ matrix- Here $E$ is the onsite energy and $t$ is the hopping ...
9
votes
2answers
282 views

In a Monte Carlo NVT simulation How do I determine equilibration

I'm running an NVT (constant number of particles, volume and temperature) Monte Carlo simulation (Metropolis algorithm) of particles in two dimensions interacting via Lennard-Jonse potential ($U = 4(\...
1
vote
0answers
41 views

What is the best order of learning mathematics before start Monte Carlo Methods for Physics simulation? [closed]

I doing medical physics masters and hope to start my research in Monte Carlo methods to simulate some particle interaction and hence to calculate some unknown physical quantities ( Stopping power and ...
1
vote
1answer
254 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
0answers
36 views

Mathematica code for self-consistent field loop in Hartree-Fock method

I want to write a Mathematica code for Hartree-Fock treatment of helium (or hydrogen molecule-> namely two electrons systems) , but I don't know how to write a loop for self consistent field procedure....
5
votes
1answer
315 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 ...
2
votes
0answers
52 views

Efficient Method for Multiplying Angular Momentum Operators

I'm doing a calculation that involves canonical symmetrization of angular momentum. For example: $H_{\text{classical}} = J_x J_y \rightarrow \hat{H}_{\text{quantum}} = \frac{1}{2}(\hat{J_x}\cdot \...
4
votes
1answer
61 views

Finding $\pm 2 \pi$ defects in 2-D lattice nematic simulation

I'm working on a Monte Carlo simulation of a two-dimensional nematic system (XY-like model with even-order Legendre polynomial interactions, such that the director angle $\theta$ obeys $\theta \equiv \...
1
vote
0answers
14 views

Modelling a flow due to an oscillating hemisphere?

Say I have a long cylinder filled with air, closed at one end and open at the other. Now say I place a compressible hemisphere at the closed end and make it oscillate (compression and expansion of the ...
0
votes
1answer
46 views

Lattice QFT: Non-homogenous lattice spacing

I am interested why we fix the lattice spacing, $a$ to be homogenous in all dimensions. After a Wick rotation, $a=i \epsilon$ where $\epsilon=t_{i+1} - t_{i}$ and with euclidean time given by $\tau = ...
0
votes
0answers
39 views

C libraries for Solid State Physics

During last year I've attended courses of solid state physics, and I studied methods for bands structure calculations in the Hartree Fock framework. However, I know only theoretical aspects of these ...
3
votes
0answers
58 views

Simulation of plasma in tokamak

I am reading some papers on numerical algorithms for simulation of plasma in the context of nuclear fusion in a tokamak. I am getting a little lost as there is a huge number of references, and it is ...
0
votes
0answers
38 views

“Instantaneous” time stepping with time dependent Hamiltonian Schrodinger equation

The Schrodinger equation for time-dependent Hamiltonian is $$i\hbar\frac{d}{dt}\psi(t) = H(t)\psi(t) \, .$$ I know that the "instantaneous" solution of this equation is $$\psi(t+dt) = e^{-\frac{i}{...
0
votes
0answers
39 views

How to make an equal flow rate in 4 glass tubes with different inner diameters?

I would like to ask a question regarding a flow rate project I've been working on that required pneumatic software to solve. Unfortunately, I don't have the software, and I don't know how to control ...
3
votes
0answers
33 views

The definition of fidelity for fermion

The definition of fidelity for two mixed ensembles is $F=Tr\sqrt{\sqrt{\rho_1}\rho_2\sqrt{\rho_1}}$. Now I came across a problem in numerical calculation, Systems A,B are identical, but attached to ...