Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [physics]

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

2
votes
1answer
50 views

center of mass, total momentum

Has someone an idea how to solve the 1. question for ANY number of mass points and a clue to the 2. question? center of mass system, momentum
0
votes
1answer
76 views

center of mass system, momentum

when I have systems with 3 mass points ($m_i, i=1,2,3$) to a time $t=t_0$ on positions $\vec r_i(t_0)$ and velocity $\vec v_i(t_0)$. How could a syntax look like to calculate 1) center of mass 2) ...
0
votes
0answers
57 views

Plotting the first Brillouin zone

I'm trying to plot the first Brillouin zone from an equation I solved for this case. This is the code I used: ...
1
vote
0answers
47 views

Decay chain browser

This https://demonstrations.wolfram.com/IsotopeBrowser/ used to be an excellent tool to demonstrate decay chains to my students. Unfortunately it does not work under newer versions of mma. Can anybody ...
3
votes
2answers
100 views

Initial values of position (x) and speed (v) of a particle visualizing using Mathematica

$$\vec{F}(\vec{r})=-m\omega^2\begin{pmatrix}x\\4y\end{pmatrix}$$ I have the force $F$ shown above. How could I specify the initial values of position ($x$) and speed ($v$) in Mathematica using the <...
0
votes
1answer
65 views

trajectory of a thrown object [closed]

I have a equation of motion v(t), that is integrated (t=0) to x(t). This symbolizes a trajectory of a thrown object. I want to create compared Plots for x(t) with different values for beta. When I try ...
0
votes
1answer
75 views

Solving a difficult differential equation

I've the following DE, describing a physical phenomenon. And the prupose is to solve that DE using Mathematica $$x(t)\cdot r+x'(t)\cdot l+a\cdot\ln\left(1+\frac{x(t)}{b}\right)=0\space\...
4
votes
3answers
107 views

Plotting an Equation Using ParametricNDSolve

I'm trying to plot the solution of a set of differential equations and see how the solution changes when the value of a certain variable de is changed. The equations have solutions to the values of de ...
0
votes
0answers
72 views

Potential with two variables [duplicate]

How would a procedure look like to calculate a potential$$ V(r)=-\int_0^r F(s) \,\mathrm{d} s $$ with example data: ...
1
vote
2answers
173 views

Time-dependent Schrödinger equation and density plot

I'm interested in the density plot of the solution of a 1D time-dependent Schrödinger equation with a given potential. So I have: ...
3
votes
3answers
197 views

Work of a force field [closed]

I have a parameter curve in 3d: r1[t_] := {Cos[2 Pi t], Sin[2 Pi t], t} r2[t_] := {1, 0, t} r3[t_] := {1 - 1/2 Sin[Pi t], 0, 1/2 - Cos[Pi t]} with a vector field:...
0
votes
1answer
53 views

Finding the minimum points of a discrete trigonometric equation

Basically I'm trying to solve the equation $$A\Big(\sin(\theta_a-\theta_{a+1}) + \sin(\theta_a-\theta_{a-1})\Big) + B\Big(n_a^z \sin(\theta_a) - n_a^x \cos(\theta_a)\Big) = 0$$ subjected to the ...
3
votes
1answer
86 views

Using ImageTransformation to implement relativistic aberration effect

****The desired effect should be images that look like something out of Star Wars like when you are boosting to high speeds in the Millennium Falcon.**** This problem was tackled by Daniel Weiskopf ...
4
votes
1answer
144 views

Solving the eigenvalue problem for a double well potential using a 1D particle in a box as a basis set

My first question is how would I go about getting the 1D particle in a box eigenfunctions using matrix techniques and how would I use the particle in a box eigenfunctions as a basis set for the ...
4
votes
1answer
226 views

Dynamic, nonlinear, damped Euler–Bernoulli beam equation

I would like to solve the 3 coupled PDEs describing a damped, nonlinear (i.e displacements in the $x$ direction along the beam need to be considered along with the $y$ displacements normally ...
0
votes
1answer
32 views

NSolve cannot work with MinValue?

I am trying to solve the following problem: find the smallest value of $b$ such that the inverse of the solution u[p] given by ...
2
votes
0answers
38 views

HypExp and HPL packages for hypergeometric functions: Evaluating a function HPL[{minus,plus},x]?

I am currently using the HypExp and HPL packages, which are useful for expanding hypergeometric functions in series around integer or half-integer values, as in common in dimensional regularization ...
1
vote
1answer
70 views

Is NDSolve being precise here? Light deflection angle

I am currently trying to calculate deflection angles for the path of photons near Schwarzschild black holes, using the formula $$\Delta\phi=2\int_{r_0}^\infty \frac{dr}{r^2\sqrt{1/b^2-1/r^2(1-1/r)}}$...
1
vote
1answer
118 views

Plotting and Animating Spherical Spring Pendulum

first post and somewhat new to Mathematica. Extremely new to making animations in Mathematica. My goal is to make a 3D animation of a Spherical Spring Pendulum system for a personal project. I ...
0
votes
1answer
44 views

Solve nonlinear differential equation system (Plot solution)

When I solve the system of differential equations I get two solutions for each variable (See the picture). My question is how I can neglect one of the solutions at the plot (The yellow and the blue). ...
0
votes
0answers
36 views

Solve nonlinear differential equation system - Plot

When I solve the system of differential equations I get two solutions for each variable (See the picture). My question is how I can neglect one of the solutions at the plot. Thanks. Code: ...
1
vote
1answer
49 views

Taking the limit of a parametric function produced by ParametricNDSolveValue

I'm working on a project for hard-sphere scattering, and to solve the for the phase function $\overline{\delta}_l(k,r)$ I use the variable phase equation $$ \frac{d\overline{\delta}_l(r,k)}{dr} = -\...
0
votes
0answers
42 views

Graphene layer rotation

I made single layer graphene by using Code as below: ...
2
votes
2answers
79 views

Spinning top: how to plot position on a sphere?

I am modelling a spinning top on a sphere and I want to obtain plots such as these I have solved the equations of motion ...
5
votes
2answers
193 views

Black hole orbit in Mathematica

I am currently investigating the motion of a particle around a black hole. The Lagrangian for the system is $$ \mathcal{L}=F(r)\dot{t}^2-\frac{\dot{r}^2}{F(r)}-(r\dot{\phi})^2 $$ Where $ F(r) = 1 - ...
0
votes
0answers
59 views

Jupiter as a Circular Mass Ring

I am trying to model the planet Jupiter as a uniform circular ring of matter and calculate the central force inside of it. This ring has the Sun as its center and defined by its linear mass density <...
0
votes
1answer
47 views

Trying to evaluate NDSolve inside of FindRoot

I want to find the equation for the range of a projectile as a function of its elevation angle at launch. Currently, I am stuck in my effort to calculate the time. What I have doesn't work right ...
0
votes
1answer
60 views

Equation of motion of a system

I am studying the problem in Vibration effects. I would like to obtain the simplified equation of motion in (2.7). ...
0
votes
0answers
74 views

How can I get a numerical value from a function?

I am a Mathematica beginner working on a problem in the area of quantum mechanics. I adapted this source code to work with the problem I'm trying to solve. At present, the graphics produced by the ...
0
votes
0answers
32 views

Assumptions in integration

I am trying to evaluate the integral $$ \int_0^{2\pi}\frac{D_{11}^2+2D_{12}^2+D_{22}^2+D_{33}^2}{1+e\cos x(t)}dx $$ Where $D_{ij}$ are third derivatives wrt $t$ as defined below. ...
1
vote
1answer
105 views

Having trouble working with two mass, three spring dynamic system

I am following a dynamic analysis example https://www.math24.net/mass-spring-system/ and trying to implement it in Mathematica. However, I am having trouble even getting simple properties like the ...
0
votes
2answers
109 views

How can I know how many points are used to by the interpolating function returned by NDSolve?

In a problem in which the solution has oscillations that I solved, the amplitude of the oscillations is much lower than I expected. It looks as if the interpolating function doesn't use enough points ...
0
votes
1answer
70 views

Gravity inside the earth with non-uniform density [closed]

I want to compute gravity inside the earth in different layers I have their thickness and density, can you please help me which formula I can use for this?
0
votes
1answer
101 views

Why does Monte Carlo integration underestimate my integral as integration range increases?

I'm having a huge problem with Monte Carlo integration. This simple example adequately shows the problem. ...
6
votes
2answers
207 views

Conductance from Landauer-Buttiker approach

I have a single-rectangular potential barrier in graphene. The transmission $T(E,\theta)$ is given by $$T(E,\theta) = \left[1+\left(\frac{V}{\hbar v_Fk_x}\right)^2\tan^2\theta\sin^2(k_xD)\right]^{-1}$$...
3
votes
1answer
73 views

WhenEvent - energy conservation not valid [closed]

I would like to simulate the following conservative system, with purely elastic impacts (coefficient of restitution is 1). Applying the WhenEvents the energy of the system is not conserved. How can I ...
0
votes
1answer
33 views

How to combine a ContourPlot with a Manipulate[Parametric Plot]

I have this code of an orbit at a Lagrange Point of the Sun-Jupiter system ...
1
vote
2answers
94 views

Stability of collinear Lagrange points

I am currently simulating the restricted 3-body problem in Mathematica. I have identified the Lagrange points and their coordinates. Now I want to see how particles move in their vicinity, to check ...
1
vote
1answer
93 views

Cutting off some terms from a matrix

I'm working with perturbation in differential geometry. I got a huge matrix for Christoffel symbol, and I need to simplify it using the consideration that $\phi^2$, $\psi^2$, $\phi\partial_{\mu}\phi$, ...
3
votes
0answers
90 views

Stationary 1-D Schrodinger equation for calculating energy of a -CH3 group quantum rotor in a metal organic framework

I am trying to calculate the energy of a quantum mechanical rotor by plugging a potential function predicted from DFT calculations in this paper: DOI:10.1039/C8CC02650E The potential can be ...
8
votes
2answers
155 views

Proper use of NDSolve in the context of the reduced 3 body problem

I am currently trying to solve the reduced three body problem in Mathematica. I have the equations of motion $\ddot{x}-2\dot{y}=-\frac{\partial\Omega}{\partial x } $ $\ddot{y}+2\dot{x}=-\frac{\...
1
vote
1answer
75 views

Contour Plot of planetary motion in Mathematica

I was trying to make the contour plot on $(v,r)$ axis for planetary motion for different orbits to check for bounded and unbounded motion. This is what I tried ...
0
votes
1answer
43 views

Syntax Error in NDSolve Code [closed]

I am trying to solve the motion equations for the Trojan Asteroids in Sun-Jupiter system. However, i am getting this error Syntax: '(' cannot be followed by 'x[t],y[t]' in this code ...
1
vote
0answers
29 views

How can I put the gradient of a scalar field in xAct? [closed]

[ for example the gradient of scalar field to square in the Gauss-Bonnet action?
1
vote
1answer
54 views

Integrating a incredibly long expression, whenever I try to use Monte Carlo I got complex answers

I am trying to integrate a recursion relation that is formulated via conformal invariance. The expression if copy and pasted from Mathematica to Word is over 20,000 pages long. Below is an ...
1
vote
0answers
70 views

Animating a 3D plot of electromagnetic waves

Could an expert on this site guide me in making an animated 3D plot of electromagnetic waves? I would like the user to be able to change the polarization interactively. Also, I would like to run ...
0
votes
1answer
112 views

What is wrong with this MonteCarlo? [closed]

I am making a trial for a MonteCarlo. In this code I simulate 10 protons hitting on a slab of gold: ...
3
votes
0answers
47 views

Sneg Library - Number operator simplify

I am trying to use the SNEG Library for evaluating some commutators. For instance I am trying to evaluate $$ \left[n_1^2,a^\dagger_1\right]$$ with the code in Mathematica ...
0
votes
1answer
74 views

Finding the limit of a function

I have the following code: ...
10
votes
3answers
346 views

DFT for Helium atom

I'm trying to solve this problem I'm not a Mathematica expert, but my program works perfect when vhartree={0,0,..,0} This is the program I wrote: ...