Skip to main content

Questions tagged [physics]

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

198 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
9 votes
0 answers
791 views

Solving the 2D Schrödinger equation with eigensystem, then verifying orthonormality of eigenfunctions with NIntegrate

I am solving the time-independent 2D Schrödinger equation for an interacting electron and hole in the case of anisotropic electron and hole masses, where the interaction is described a modified form ...
8 votes
0 answers
246 views

How can we improve transonic flow visualization?

With this code we can make 2D FEM simulation of transonic flow around airfoil NACA 0012 at Mach number of 0.925. It takes about 5 minutes on the XENIA-15 laptop of 32 GB memory with processor Intel ...
• 45.8k
8 votes
0 answers
782 views

Numeric Solution Hydrogen Atom

I'm trying to solve the Schrödinger equation for the hydrogen atom without made the variable separation of the polar and radial coordinate. It is my test code to extrapolate to another system with its ...
• 491
7 votes
0 answers
825 views

Mathematica packages for numerically solving the Klein-Gordon and Dirac equations?

The book "Advanced Visual Quantum Mechanics" by Thaller includes Mathematica software packages for the numerical solution of the Klein-Gordon equation and the Dirac equation (subject to user-defined ...
• 1,637
6 votes
0 answers
755 views

• 295
4 votes
0 answers
121 views

• 31
3 votes
0 answers
299 views

How can I improve the solution of my PDE?

I want to improve the last code to find a better solution for umi. For now I'm using only MaxStepFraction. Also, when I put no ...
• 33
3 votes
0 answers
162 views

How to control pedestrians to enter the simulation region with random functions?

I am working on simulating crowds of a cross-typed region, given that pedestrians entering from 4 gates of north, east, west and south, ...
• 141
3 votes
0 answers
1k views

Visualising Special Relativity

Mathematica newbie learning by doing (or not so far, in this case). Problem: Special Relativity - consider two inertial frames $A$ & $A'$ in relative motion. At a time $t0$ as measured in $S$, ...
• 2,560
3 votes
1 answer
499 views

2 body orbits around a barycenter

I'm attempting to plot a barycentric orbit of the sun and Jupiter. Initial Conditions ...
• 847
2 votes
0 answers
91 views

How to effectively substitute variables in xAct?

I am trying to find the gauge field equation by varying the Born-Infeld Action with the help of xact package. I use the following ...
• 2,423
2 votes
0 answers
453 views

• 503
2 votes
0 answers
92 views

Seeming hard-limit on image resolution?

I am trying to generate some very large images of the basins of attraction of a dynamical system. However, the resolution of the images seems to cap out at 10001 x 10001 pixels, and no arrangement of ...
• 21
2 votes
1 answer
375 views

Solving two coupled partial differential equations

I am trying to solve the following system of two coupled partial differential equations (both equations equal 0): Here, $V$ and $Y$ are functions that only depend on $r$, i.e., $V=V(r)$ and $Y=Y(r)$. ...
2 votes
0 answers
533 views

Equation of Motion: Perturb Scalar Field and flat FLRW Metric (xAct)

Is there a way to get the equations of motion (eom) of a scalar field from a Lagrangian in a perturbed FLRW metric? I want to perturb the field as well as the metric to 1st order. What I did so far: ...
• 181
2 votes
0 answers
143 views

Speeding up NIntegrate and DensityPlot

I'm making a density plot of a spectral function Sfull[kx, ky, kz, Ω] along a specified path in {kx,ky,kz}-space. The plot looks ...
• 255
2 votes
0 answers
49 views

OscillatorStrength[Helium] and OscillatorStrength[Rubidium] do not return results

I am new to working with the Resonance Absorption Lines Pacgage << ResonanceAbsorptionLines` and trying to write code using the OscillatorStrength command. ...
• 161
2 votes
0 answers
239 views

Plotting the surface spectral function

I'm not 100% sure if I should be posting this in the physics stack exchange but, I am pretty new to Mathematica, just installed it five weeks ago, and I am having a little trouble implementing the ...
2 votes
0 answers
250 views

2 votes
0 answers
72 views

What kind of performance should I expect out of Eigensystem using FEAST?

I'm numerically solving a time-independent Schrödinger equation using Eigensystem's FEAST method. It takes a lot longer than I ...
• 1,844
2 votes
0 answers
169 views

Symbolic perturbation expansion for quantum mechanics using Hellmann-Feynman derivaties

I am interested in some quantum mechanical perturbation expansion for energies. Actually I want to implement these terms $E_n^{(k)}$. As is stated below one can do that using CAS. I would be ...