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
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
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:
...
