# Tagged Questions

Questions related to the application of Mathematica to problems in physics.

1k views

### Efficient implementation of cubic equation of state

...In this case the Peng Robinson Equation of State. Equations of state are empirical equations with parameters derived from experimental data. They are used to predict pure component and mixture ...
41 views

### How to combine two matrices.? [duplicate]

let A, B are any 4 by 2 Matrices , i want to know whether it is possible to obtain a 4 by 4 matrix by using equation A=B Thanks in advance
45 views

### Polarization sum in FeynCalc

I think I know how to sum over final state fermion spins (FermionSpinSum[...]?) but what is the analog expression (in FeynCalc) if you have two or more polarization ...
57 views

### How to implement time dependence in the DoFun package?

DoFun is a very usefull mathematica package that helps to write Dyson-Schwinger and FRG flow equations when these are to complex to write down by hand. You can download it here and read it's manual ...
68 views

### Calculating the normalization constant for a wavefunction

How can I compute the normalization constant for a quantum mechanics wave-function, like $\Psi(x) = N \exp(-\lambda x^2/2)$ by using Mathematica? (The normalization constant is $N$).
34 views

### How to combine similar terms in a density matrix expression? [Ref: Quantum Computing]

I tried the following : Simplify, FullSimplify, and Expand. Due to this problem, I am not ...
7k views

### Find eigen energies of time-independent Schrödinger equation

I'm trying to get the eigenvalues of a one dimensional time-independent Schrödinger equation, $-\frac{h^2}{2m_0}\frac{d^2\psi}{dx^2}+U(x)~\psi=Ei~\psi$ where U(x) is some potential and Ei is the ...
866 views

### NDSolve::ndcf: Repeated convergence test failure. How to solve?

I am trying to simulate a system of $n$ pendulums with some friction in Mathematica 9. This is the code I am using: ...
84 views

### Divergence of $1/r^2$: Mathematica yields 0, while the correct answer is a DiracDelta

In spherical coordinates, the divergence of a point charge is $4 \pi \delta (r)$, the three-dimensional Dirac delta function. In ...
163 views

### Attempting to use NDSolve to plot harmonic oscillator solutions

NDSolve[{ (-y''[r]/1880) + (470 (0.04077)^2 r^2 - 48 + 1/(1880 r^2)) y[r] == 0, y[0] = 0, y'[0] = 0}, y, {r, -4, 4}] I use this but get errors and am ...
173 views

### Help with NonLinearModelFit

I have a question about fitting data with a model. I have a quite complicated model, based on Arctans, modified Bessel functions, etc, that I want to fit to my data, a table of couples (minimal ...
81 views

### Singularity while solving PDE

I have to solve a PDE (in the context of the functional renormalization group in physics). I have a function of two variables, $U(l,p)$. I know $U(0,p)=U_0(p)$ and then I have an equation (eq925) that ...
113 views

### A 1D numerical integral Mathematica cannot compute, from physics

A well know result in theoretical physics is that a sum over Matsubara fermionic frequencies, i.e.: $$S = \sum_{n=-\infty}^{\infty} h(\omega_n) \hspace{32pt} \omega_n=(2n+1)\frac{\pi}{\beta}$$ can ...
364 views

### Plotting 2-D Electron Density Carbon Atoms

I need to create a 2-D electron energy density plot in Mathematica to compare with my STM experimental results in my lab class. This would be done by plotting the sum of the symmetric and ...
188 views

### Boson commutation relations

I have tried to adapt this answer to my problem of calculating some bosonic commutation relations, but there are still some issues. The way I'm implementing the commutator is straightforward: ...
389 views

### Time evolution of a bispinor [duplicate]

The effective Hamiltonian for polariton graphene is given by: ...
42 views

### Plotting a parameterized vector in the x-y plane

I am trying to plot the trajectory of a parachutist descending at a constant vertical rate of 8 m/s while a south wind also blows at a rate of 8 m/s. He initially jumps from a height of 300 m. I am ...
291 views

### Time evolution/dynamics of circular plate with hole (Biharmonic equation and stiffness)

On this Mathematica.SE form, there exists information on how to use Mathematica to demonstrate the vibration of a circular membrane and the deflection of an orifice plate (the latter I had raised ...
275 views

### FEM Solution desired for “Plate with orifice” deflection: Application of Boundary Conditions and use of Regions

I found the deflection of an orifice plate (circular plate with a hole) subject to uniform pressure using Mathematica's NDSolve functionality. The plate is fixed at its outer circumference and is ...
99 views

### How to reduce an expression into user-defined variables?

In the following code, I define a set of matrices. Then, I define a function that calculates the commutation relation relation between these matrices ...
3k views

### Is there something like DensityPlot3D to visualize atomic orbitals?

I'm visualizing some hydrogen like atomic orbitals. For looking at plane slices of the probability density, the DensityPlot function works well, and with something ...
143 views

### Synchrotron Radiation and ListDensityPlot

I'm trying to reproduce the following plots: As described by this blog entry: Synchrotron radiation As you can see in the linked article, the task is essentially to make a ...
142 views

### Interfacing Red Pitaya board with Mathematica [closed]

I just discovered the Red Pitaya board and borrowed a unit to play during this weekend. Surprisingly I did not find any reference to Mathematica on the Red Pitaya forum, nor did I find any reference ...
128 views

### Help nedeed for understanding NDSolve architecture -Vacuum Heating Simulation

I am a beginner in mathematica and I am facing a problem that is way bigger than my knowledge.. I am tryng to simulate a semplificated version of Brunel Effect, wich is a Laser-Plasma interaction ...
125 views

### Creating an orrery

I am interested in creating an orrery, ideally with SystemModeler. I have been experimenting with the MultiBody library since that allows you to do animations. The most relevant example seems to be ...
121 views

### Lens profile in a grid [closed]

So I have a FDTD code which works well. The propagator of the code itself is: ...
117 views

### Equal flow boundary conditions

I'm having the following issue that can not find a solution , by using this boundary condition v[0, t] == u[1, t] that Mathematica does not understand, someone ...
148 views

### FeynRules issue when trying to export UFO file

I am having problems with the UFO file generated by a very simple model using FeynRules. The following error occurs when MadGraph tries to import the "Model_UFO" file: ...
768 views

### How to plot an emission spectrum?

If I have a list of data with various wavelengths in nanometers, how would I plot them on a graph so it looks like this: So far I have managed to plot a spectrum in ...
2k views

### ColorData[“VisibleSpectrum”] is wrong?

I am trying to visualize the visible spectrum using the built-in ColorData["VisibleSpectrum"] function which "colors based on light wavelength in nanometers". But I ...
62 views

### Slater Koster parameters generated using Mathematica

Here is a link where people describe how one can generate the entries of the Slater Koster parameter table seen here as Table of interatomic matrix elements. The ...
258 views

### Sampling a phase space uniformly for a given energy?

The energy for a 1D harmonic oscillator can be written: $$E = \frac{1}{2 m} p^2 + \frac{m \omega^2}{2} x^2$$ where x is position and p is the momentum. I would like to sample phase space $(x,p)$ ...
312 views

### Crank-Nicolson scheme for the Schrödinger equation

I am trying to solve the 1-D Schrödinger equation with the Crank-Nicholson method. I use the atomic units: time < 10^-14 => 413 a.u.t. length 30 * 10^-9 m => 567 Bohr radii ...
195 views

### Potential flow over cylinder with FEM: Pressure field is flipped

Much in the fashion of a few problems already answered on Mathematica.SE (1, 2), I am trying to solve the Laplace equation in velocity potential $\phi$ to simulate irrotational, inviscid flow over a ...
703 views

### How to solve fluid flow problem based on Navier-Stokes equations?

Does anyone know or can provide any examples how fluid flow problem can be formulated and solved in Wolfram Language? Simplest cases of 1D or 2D flows based on Navier-Stokes equations or even their ...
1k views

### Trying to graph trajectory function in Mathematica? Won't Graph?

The exercise is: Define the trajectory function y(x) in Mathematica. Pick the initial speed (in meters per second) equal to ...
90 views

### How to obtain the magnetic field of the following set up [closed]

This is the experimental set up: A ferromagnetic cylindrical disk (of irrelevant dimensions) has either sides set at different electric potentials (left side at 0 volts and right side at 1 volt). A ...
118 views

### Alternatives to Mathematica that are not computer algebra systems [closed]

Are there alternatives to Mathematica that are not computer algebra systems? I know there have been questions on this, but I am focusing on oprn source projects and instead of focusing on computer ...
77 views

### Solving 4 coupled differential equation for modelling boson star

I'm trying to study about modelling star from a journal. To get the characteristics of the star, I have to solve 4 differential equations which are coupled each other. I want to use Mathematica 9 to ...
77 views

### Find all degenerate eigenvalues of a cubic equation

I have an equation that is cubic in w. The three solutions correspond to bands in a bandstructure, and are a function of wavevector ...
208 views

### Finding the conditions for degenerate solutions to a cubic equation

I have a cubic equation (corresponding to the band structure of a physical system) given by w in the code below. ...
154 views

### Specifying NeumannValue for Axial load instead of transverse load (Plane Stress situation)

I am trying to simulate a bar/link element under axial tensile load. The link is of unit thickness and I may use a plane stress model. A plane stress model is an approximation that may be used for ...
788 views

### A better “VisibleSpectrum” function?

Since ColorData["VisibleSpectrum"] is wrong, I would like to have a more accurate function to use. Can this information be extracted from Mathematica itself?
83 views

### Eigenvalues function not producing expected results?

I'm running into a strange issue when trying to compute eigenvalues of a 2x2 matrix for graphene. If I "manually diagonalise" by multiplying the two off-diagonal terms together and taking the square ...
92 views

### FullSimplify not giving me a physically useful answer for Mott scattering

I am using Mathematica to find the cross-section for Mott scattering over an extended target. I evaluated ...
2k views

### Convert spectral distribution to RGB color

Presume I have a spectrum as a function of wavelength (an example being the blackbody spectrum): I want to convert that to a single RGB color to display on-screen, i.e. the "color" of that object ...
180 views

### Modeling a stellar atmosphere

I'm trying to model a stellar atmosphere using Mathematica. Basically I'm working in a physical and chemical model which combined should give me a profile of the chemical abundance in the atmosphere. ...
353 views

### Pull out scalars from NonCommutativeMultiply in commutator of SU2 spin algebra

I am trying to work on spin chains in Mathematica. That is, N spin-1/2 on a one-dimensional lattice with periodic boundary condition (N+1=1). I need therefore the SU2 spin algebra $S_z, S_+, S_-$ with ...