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

learn more… | top users | synonyms

1
vote
0answers
92 views

Tracking solutions over a Brillouin zone

Thanks to @bbgodfrey I have a function that finds the two solutions to a set of equations bounded by a hexagonal region (corresponding to the Brillouin zone). It takes the angle ...
0
votes
0answers
84 views

Difficulty solving an equation involving a couple of different powers

I have this semi-complicated equation to which I want to find an analytical solution for $\epsilon$. The equation is a s following: $$a^2 (2 b k \epsilon \sqrt{1 - \epsilon^2/b^2} - 2 c k \...
8
votes
1answer
373 views

Diverging solution to coupled second order ODEs from NDSolve

For a physics application I am considering the radial Bogoliubov-de Gennes equations in two dimensions. In dimensionless form they read as follows $$ \begin{cases} \bigg[\frac{1/4-m^2}{r^2}+2(E-1)\...
7
votes
2answers
563 views

Integrating a function over a surface integral

From a first principles bandstructure calculation I get an energy scalar field in three dimensions $E(x,y,z)$. It's now easy to plot a constant energy (contour)-surface for dedicated values $\epsilon=...
1
vote
0answers
63 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 ...
0
votes
1answer
193 views

Eigen energies of 2D time independent Schrödinger equation

Firstly, I am aware that a similar question of this type was asked, but it was not helpful as it was only one dimensional (so please don't mark my question as a duplicate) The trouble I am having is ...
3
votes
1answer
154 views

Controlling scaling of multiple vector fields in VectorPlot3D

I'm trying to plot the electric and magnetic fields of a relativistically moving electric charge using VectorPlot3D, and apply "Manipulate" to change the velocity. The code I'm using is as follows: <...
0
votes
1answer
85 views

Problem simulating a wave on a string by solving the wave equation as a difference equation.

This is a question I posted earlier however the purpose of what I'm doing is going to be expounded much more clearly. ...
1
vote
1answer
68 views

Implicit region is not “valid” in ParamatricNDSolveValue function?

I am trying to solve a geometric problem with relation to my Schrodinger equation and its boundaries. Here is my code: ...
1
vote
0answers
51 views

How to evaluate Poisson Bracket? [closed]

I'am new in Mathematica so please don't yell at me. I want to define Poisson Bracket in Mathematica 9 and evaluate it for two functions,but it didn't worked. Could you help me please? ...
0
votes
1answer
257 views

How to plot and visualize a single linear vector in 3D? [closed]

I have a vector (in physic) designated as F1=250cos(60)i+250cos(60)j+250cos(45)k, and i would like to see it in a 3D graphic with the axis centered at the origin, after what i would include other ...
1
vote
1answer
217 views

Particle motion with discrete values for the potential

Sorry if the question seems trivial, but I am new to physical simulations in Mathematica. I have a numerical estimation of an electric potential $V(x)$ in 3D, to be more precise I have a matrix ...
1
vote
1answer
243 views

NSolve fails to solve a nonlinear equation

Mathematica can not solve this: g = 9.82; ω = 0.5; h = 5; y0 = 1; v = 0; τ = 0; NSolve[h + v t - (g t^2)/2 == y0 Sin[ω t], t]; The error code is: NSolve::...
31
votes
3answers
977 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?
1
vote
1answer
58 views

Variables within functions within DSolve

Hi I am just beginning to learn Mathematica and this is my first time I have been exposed to any type of coding. I am encountering a problem in a basic physics problem. For example we are always ...
3
votes
0answers
170 views

Numerical integration: complicated 2D integral seems to be poorly estimated

In the course of some physics research I've been working on, a very annoying integral has appeared that I'm having difficulty evaluating numerically. Any help you could offer would be greatly ...
3
votes
0answers
150 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, ...
0
votes
1answer
86 views

How can I use manipulate on multiple lists within a function?

So I have my N-Body Problem set up like so, and I can plot it in ListPointPlot3D, but I cannot figure out how to stick this into a manipulate function. ...
4
votes
2answers
322 views

Embedded Graphics: General Relativity

Is there a way to make the following graph with two dimples? Like the central indent and another indent to the side to simulate a central mass like the sun and another smaller mass like the Earth in ...
2
votes
0answers
141 views

Position from 3 dimensional acceleration data? [closed]

Per a previous post on Mathematica StackExchange here I'm trying to utilize 3-axis accelerometer data to determine approximate 3D position. I wanted to start with a simple case. To this end I've ...
12
votes
2answers
302 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)$ ...
3
votes
1answer
377 views

Introduction to Vectors and NDSolve

I want to use NDSolve for a complicated equation that contains vectors. I am still new to Mathematica, and have really enjoyed what I have learned thus far. Concerning this question, I have already ...
0
votes
0answers
192 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: ...
5
votes
3answers
593 views

Numerically Simulating a Simple Differential Equation: RC Circuit

I need to learn to simulate a differential equation in Mathematica. The differential equation that I need to simulate is complicated, without an analytical solution. However, I am brand new to ...
5
votes
2answers
203 views

Multidimensional integral to compute the end-to-end distribution of a FENE ideal chain

How do I evaluate the mean end-to-end squared distance of a FENE ideal chain at fixed inverse temperature $\beta$ in the canonical ensemble? This quantity is defined as the mean value $$\left<\...
2
votes
0answers
143 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 ...
2
votes
3answers
792 views

draw part of parabola using graphics

today I was doing my physics homework. one of the questions was draw velocity-time and position-time graphs from the acceleration-time graph. So for the first step I draw acceleration-time using ...
13
votes
2answers
513 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 ...
13
votes
2answers
782 views

How to calculate position from 3 dimensional acceleration data?

I discussed this problem yesterday with Pickett here. We concluded that the best solution may be just Interpolate and then ...
2
votes
1answer
214 views

2D Fourier transform of a few (4) disjoint discs on a plane

I'd really appreciate some advice. Short Version I'm trying to calculate the following $$ \psi(X,Y,z=d)=\underset{aperature}{\iint}\psi(x,y,z=0)e^{-i\frac{k}{d}\left(xX+yY\right)}dxdy\\\psi(x,y,z=0)=...
9
votes
3answers
678 views

Calculate inertia tensors

Moments of inertia and inertia tensors are well-known characteristics of rigid bodies in physics and applied mathematics. How to calculate them in Mathematica? Is possible to use regions which was ...
5
votes
1answer
1k views

NDSolve and Piecewise functions

I want to clarify upfront, This is not a please help me get the answer to my homework question. I know what I need to do I just do not know how to get mathematica to do it. I am working on an ...
8
votes
3answers
719 views

Biot-Savart Law/Magnetostatics Solution

I am working on solving the Biot-Savart Law equation for the magnetic field around a charged ring of uniform current density. The expression that Mathematica gives is rather nasty, as expected. Notice ...
8
votes
1answer
528 views

How to draw quantum circuit with mathematica?

How to draw quantum circuit as the following one with mathematica? I only find the quantum package which is designed for MMA 7.0 as the author said. But the version of my Mathematica is 9.0 and when ...
7
votes
1answer
1k 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 ...
21
votes
4answers
3k views

Convert spectral distribution to RGB color

Presume that 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 ...
8
votes
3answers
192 views

Substituting values in 3j symbols

I believe these three expressions should give the same answer of -1/Sqrt[3], but the first gives +1/Sqrt[3] whereas the second ...
3
votes
1answer
179 views

Interpolation function rejected by ParametricPlot3D

I was playing around with the bouncing ball example, and decided to make a 3d version. I decided to use ParametricPlot3D for a simple independent 3D presentation....
1
vote
0answers
207 views

Solving PDEs with complicated boundary conditions [duplicate]

The system I'm trying to solve is $$\nabla^2 C_{(r,\theta)} =0$$ $$C_{(\infty,\theta)}=C_0$$ $$ [ \frac{\partial C_{(r,\theta)}}{\partial r} \cos(\theta)+\frac{1}{r}\frac{\partial C_{(r,\theta)}}{\...
1
vote
1answer
79 views

How to make a discretized NMinimize more precise

I am using Mathematica for physics research and I want to minimize a Hamiltonian equation with respect to two variables (I have also discretized the problem). I have a single constraint. When I plot ...
9
votes
1answer
1k views

Schrödinger eigenvalue problem in two dimensions (Harmonic Oscillator)

I read here the discussion about how to solve a one-dimensional eigenvalue problem. I am wondering, how can one generalize these methods to two dimensions? For example, how to solve this equation: $$...
0
votes
1answer
139 views

Help with Computation-heavy triple integration

I'm Trying to use Mathematica to symbolically integrate a function with respect to x, y, and z. For some reason, this computation either takes about, 60s and returns the input, or takes many hours ...
1
vote
0answers
91 views

Help wanted in evaluating an integral

Trying to evaluate an integral. To be specific, I'm trying to find the expectation value of $x^2$ for the ground state of hydrogen. Here is the code I used: ...
0
votes
1answer
372 views

Solving for the time-evolution operator in a periodically driven system

I am looking at the Hamiltonian $$H(t)=\begin{pmatrix} 0 & e^{i\Omega t}\\ e^{-i\Omega t}& 0\end{pmatrix}$$ I am trying to solve for the unitary operator $U(t,0)=\mathcal{T}\exp(-i\int_0^t ...
22
votes
2answers
3k views

Animation of double pendulum

Sadly, I am a completely newbie. I am studying Physics and in our theoretical physics class we got the task to solve the double pendulum using Mathematica. We just got the program, but no introduction ...
9
votes
3answers
2k views

How to create an inverse pole figure color map

To interpret electron backscatter diffraction (EBSD) results, inverse pole figures are used. A three dimensional space of directions ([100], [110] and [111] for those who know about crystallography) ...
3
votes
1answer
154 views

How to simulate the 1-D motion of a particle under constant acceleration? [duplicate]

I am trying to represent the position of an object under freefall, but am completely lost. Here is my attempt so far: ...
3
votes
1answer
610 views

Traces of products of Pauli matrices

Given a matrix $M$ of shape $2^L*2^L$, I would like to compute all the traces $\text{Tr}( M.(\sigma^{n_1}\otimes\sigma^{n_2}\otimes\ldots\otimes\sigma^{n_L})) $ for $n_1=1...4$, $n_2=1...4$, ..., $n_{...
1
vote
2answers
167 views

Making mathematica do regulated integrals

Consider the integration, $\int _0 ^\infty dx\ x \tanh( \pi x) \sqrt {x^2 + a^2 } $ where $a$ is a real number. This integral is divergent. We note that an asymptotic expansion of the integrand ...
2
votes
1answer
161 views

Operator operation

There are two basic operators r={x,y,z} and p=-I*h*{D[#,x] , D[#,y], D[#,z]} &. The angular momentum operator is defined as $...