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.

14
votes
0answers
818 views

How best to simulate n-body systems in a functional way

I want to simulate a N-body system in the most functional way possible. The system I'm modeling are metallic colloidal particles in a solution of peroxide water, as here http://arxiv.org/pdf/1502....
7
votes
0answers
237 views

How to use periodic boundary conditions on the derivative of u[x,y,z]?

I have a structural mechanical problem, where i try to calculate the strain field in the modelled material. I simplified the problem in the following example: There is a material with two phases, ...
6
votes
0answers
200 views

How to improve precision of finding the eigenfunction for this 2D Schrödinger equation?

I'm trying to solve a 2D Schrödinger equation with the following potential: ...
5
votes
0answers
306 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 ...
4
votes
0answers
270 views

Is it possible to use Gaussian cgs units in Mathematica computations?

I was trying to do computations in Gaussian cgs units, and of course I failed due to dimensional issues. For instance, $q_1 q_2/r$ cannot be converted to $\mathrm{erg}$ since Gaussian system uses $4\...
3
votes
0answers
89 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 ...
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 ...
3
votes
0answers
197 views

Presentation of SNEG creation and annihilation operators in Mathematica 11

I am using the Mathematica package SNEG to make some calculations in solid state physics (see http://nrgljubljana.ijs.si/sneg/). The programs in SNEG library are running, but their presentation are ...
3
votes
0answers
125 views

Finding possible lattice planes of a crystal structure

After generating a crystal structure from a crystallographic data and duplicating it to a larger crystal system I would like to find possible lattice planes of this crystal. It is well-known that the ...
3
votes
0answers
342 views

Solve PDE over time-dependent region

I have some experience with C++ programming, but I am quite new to Mathematica and I think the latter requires a different mindset which I am still not used to. I would like to simulate the effect ...
3
votes
0answers
195 views

Differential equations in partial derivatives

I am trying to get numeric solution for equations: $$ \left\{x^{(2,0)}(u,t) \left(1-\frac{y^{(1,0)}(u,t)^2}{\left(x^{(1,0)}(u,t)^2+y^{(1,0)}(u,t)^2\right)^{3/2}}\right)+\frac{x^{(1,0)}(u,t) y^{(1,0)}(...
3
votes
0answers
237 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 ...
3
votes
0answers
158 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, ...
3
votes
0answers
660 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
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 ...
2
votes
0answers
247 views

Wormhole embedding diagram

I am trying to reproduce the embedding diagrams for the evolution of a Schwarzschild wormhole described in this paper. Following the paper notation, we denote the Kruskal coordinates by $(v,u)$. For ...
2
votes
0answers
535 views

Plotting the Magnetic Field

I have a dataset of $ \textbf{B} = {B_x, B_y, B_z}$. (Quite a big dataset :78*100*150) I want to plot magnetic field lines. If I was working with an analytical function, the magnetic field line ...
2
votes
0answers
319 views

Linear tangent steering law and orbital insertion altitude

I've been trying to figure out how to build a rocket launch simulation that allows the user to specify the desired altitude at which a payload will be inserted into a stable circular orbit, but have ...
2
votes
0answers
180 views

Modelling heat exchange between flowing fluid and surroundings

Starting from this question: Couple a PDE and ODE in NDSolve I would like to modify the fluid from static to a flowing one in the y direction. that's the original code that I modified a little bit: ...
2
votes
0answers
153 views

What is a good way to model fluid stream contours around airfoils?

I am an undergraduate student and am doing a project on conformal mapping for my complex analysis class. I have had one heck of a time trying to plot stream function contours to model fluid flow in ...
2
votes
0answers
102 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 ...
2
votes
0answers
205 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 ...
2
votes
0answers
197 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 ...
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 ...
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 ...
1
vote
0answers
72 views

Problem with Wick Contraction

I'm trying to compute wick contractions knowing that: $\langle a_a a_b \rangle = \delta_{ab}$. So far I've the following ...
1
vote
0answers
64 views

NDEigensystem works for 3-D Schrodinger Equation but not 4-D

NDEigensystem works fine for the following 3-D Schrodinger Equation code -- Clear["Global`*"] Rcube = ImplicitRegion[-1 <= x <= 1 && -1 <= y <= 1 && -1 <= z <= 1, {x,...
1
vote
0answers
35 views
1
vote
0answers
95 views

Conformal field theory using Mathematica

Can someone suggest some open source mathematica packages to perform computations in conformal field theories e.g, OPE expansion, finding conformal blocks etc.
1
vote
0answers
79 views

Evaluate the electrical resistance between any two points of a circuit

I have an electric circuit of which the structure is known but the value of the single components (all resistors) is not. I can measure experimentally the resistance between some points of my circuit. ...
1
vote
0answers
87 views

Fisher Forecasting Using Mathematica

I'm need to do a Fisher Forecast for my experimental physics requirements. I'm trying to constrain parameter for a modified Friedmann equation. Is there any package or code that any of you have handy ...
1
vote
0answers
49 views

How to get a conversion factor in new system of unit

There is some numerical relation $I = \frac{c}{8\pi} E^2$ between the wave intensity $I$ and the electric field strength $E$ written in CGS (centimetre, gram, second) system of units, where $c$ is ...
1
vote
0answers
49 views

How do I get the output of a square-wave input? (singularity problem)

I am trying to get the electrical response of a system to an input square wave. Here is the code: ...
1
vote
0answers
97 views

Simpler way to compute gravity?

I have code for a simple n-body problem: ...
1
vote
0answers
205 views

Solving the 2D heat eqn with a rastering q term

I'm brand new to Mathematica so please forgive my lack of knowledge, I'm trying to learn it for my current internship! I'll try give some background to what I'm doing and where I could do with some ...
1
vote
0answers
255 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 ...
1
vote
0answers
206 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 ...
1
vote
0answers
218 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 ...
1
vote
0answers
95 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
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: ...
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: ...
0
votes
0answers
42 views

Graphene layer rotation

I made single layer graphene by using Code as below: ...
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
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. ...
0
votes
0answers
71 views

Need a method to solve a nonlinear PDE with gradients of the same order as the PDE's

I have a system of first order non-linear PDE's, but some of my boundary conditions are also of first order, so mathematica cannot solve it. When I looked online for similar problems , I found some ...
0
votes
0answers
98 views

FiniteElement and NDSolve won't work

I need to solve a system of DE's in 2 dimensions over time, The problem is that mathematica complains when I have a boundary condition of the same order as the differential order of the DE's, so I ...
0
votes
0answers
60 views

Integration involving vectors, in 3D

Since I'm new to Mathematica I was looking for ways to learn vector integration and how to do them efficiently, so I came across this article where they did calculation for field in a point above a "...
0
votes
0answers
34 views

Implement commutativity and associativity for custom operation

Just to give context: I want to implement quantum mechanical ladder operators and for this I want Mathematica to realize when it can just commute my ladder operators with numbers or variables because ...
0
votes
0answers
109 views

NDSolve error: 2D Convection-diffusion equation

I am trying to numerically solve a 2D convection-diffusion equation using NDSolve. This is version 11.2 Here is how I am modelling it: ...