Questions tagged [quantum-computing]
Questions about simulation of quantum circuits and quantum algorithms in Mathematica.
53 questions
4
votes
1
answer
143
views
How to implement a 2x2 operator in Mathematica?
I am trying to write a Mathematica program to compute the following:
For a given Hermitian matrix $\rho$, the operator $L_\theta$ with respect to a parameter $\theta$ is defined as:
\begin{equation}
...
- 765
2
votes
1
answer
166
views
Why does NDSolve fail or slow down when simulating a driven JC model with time-dependent terms, even when the drive is off?
I'm performing a numerical check of an effective Hamiltonian transformation applied to a driven Jaynes–Cummings (JC) model. After applying a rotating frame transformation and rotating-wave ...
- 165
1
vote
0
answers
146
views
Overload Issue with Solving Schrödinger Equation in Mathematica [closed]
When attempting to calculate the Franck Condon Factor for a molecule, I am encountering an issue calculating the higher vibrational levels of the molecule (which is integral to my analysis of the ...
1
vote
1
answer
106
views
Issue with Interpolated Potential for Time-Independent Schrodinger Equation
I am attempting to use an interpolated potential for a molecule, initially in the form of (x, V(x)), as input for solving the Time Independent Schrodinger Equation.
After defining the data points, ...
2
votes
1
answer
153
views
Seeking Recommendation on Time-Independent Schrödinger Equation Solving Method in Mathematica
I am seeking to solve a Time-Independent Schrödinger Equation using the Mathematica Platform. I will detail the problem below:
Scaling Constant is "K." This is calculated by the formula (2 m ...
0
votes
0
answers
38
views
How to Expand and Simplify a Complex Hamiltonian in Mathematica?
I am working on a quantum system and need to expand and simplify a Hamiltonian in Mathematica. The Hamiltonian is given as:
$$
H_{I} = \hbar \frac{\chi}{2} a_{1}^{\dagger 2} a_{1}^{2} + \left(i\hbar \...
6
votes
1
answer
328
views
How to Numerically Solve the Steady State of a Lindblad Equation with Parity Symmetry?
I am trying to reproduce some figures from this paper. like below:
The effective Hamiltonian of the system is
\begin{align}
H_{\text{eff}} = \frac{3g_x^2}{\omega} |e\rangle \langle e| + \frac{2g_x^...
- 165
1
vote
0
answers
73
views
Entropy Production Rates
I am trying to reproduce the results shown in Fig. 2(a) of the paper https://journals.aps.org/prresearch/pdf/10.1103/PhysRevResearch.6.033037 to calculate entropy production rates in a cavity ...
- 11
6
votes
4
answers
863
views
Is it possible for Mathematica to output the name of a matrix as opposed to its matrix form?
I am interested in working with the names of matrices as opposed to the actual matrices themselves. I would like Mathematica to recognise the content of a matrix and output its name rather than its ...
- 997
0
votes
1
answer
131
views
How to confirm that these tensor-product matrices commute [closed]
So I want to show that these "stabilizer operators" commute.
These operators are the tensor products of $8$ $2 \times 2$ matrices. As such the operators themselves are $256 \times 256$ ...
- 997
0
votes
0
answers
121
views
Using CrossEntropyLoss function for basic neural network
Hi I am interested in understanding how to use the neural network functions in Mathematica. I have just started learning hence this is a basic query. The simple example I am looking at is as follows: ...
- 291
1
vote
0
answers
343
views
Modern way to do quantum optics calculations in 2024?
Quantum Optics can often require messy algebraic simplifications of noncommuting operators (usually ladder operators and commutators). I've asked here about how to do it before.
What are some modern ...
- 1,066
2
votes
0
answers
84
views
Deriving Expressions for R1, R2, and R12 , Normalized Coincident Detection Probability, and LG Mode Formalism
I am currently reviewing a research paper on orbital angular momentum entanglement between two photons (signal and idler) generated from a pump photon via SPDC. The paper introduces the following ...
- 45
3
votes
1
answer
530
views
Solving a Lindblad differential equation in matrix form
I am trying to solve the following Lindblad equation analytically with Mathematica:
$$\frac{d\rho}{dt}=-i[H,\rho]+\sum\nolimits_{i=0}^1 K_i\rho K_i^+-\frac12\cdot\{K_i^+K_i,\rho\},$$
where
$$H=\begin{...
- 53
1
vote
0
answers
101
views
Why must the "&" symbol be placed inside function instead of outside? [closed]
The following code is taken from Wolfram Community "Quantum Error Correction: bit flip".
...
- 997
0
votes
0
answers
102
views
Solving a system of complex matrix equation using Reduce function
I am trying to solve a system of four matrix equations using Reduce function. The system arises in Quantum Information. The system is solvable with a solution i.e. phi, but the following code gives ...
1
vote
0
answers
141
views
Gaussian Matrix Model - Large N Expansion Physics with Mathematica
Update: Attempted 2nd time by replica method
I'm trying to implement basic large N expansion in Mathematica.
I'm trying to understand what this means on the fly as well.
The model I am trying to ...
- 165
4
votes
2
answers
582
views
How to generate all permutation matrices for 4 qubits?
In quantum mechanics, a qubit can be understood as a 2 by 1 vector denoted by Dirac notation as $$|0\rangle \equiv \left( \begin{array}{c}
1\\
0\\
\end{array} \right) ,|1\rangle \equiv \left( \...
- 521
6
votes
1
answer
236
views
Minimizing computational time for a quantum walk problem
Any suggestions as to how to speed up the computational time for this quantum walk problem which is coded using a normalized SparseArray coin operator as follows:
...
- 3,303
1
vote
2
answers
148
views
Defining and executing large matrix operators on vectors efficiently
For an $N$-qubit system, I am interested in coding, in an efficient way, the following set of Kraus operators $\{A_c\}_c$ defined in the computational basis as
\begin{align}
\hat{A}_c := \sum_{m=...
- 291
2
votes
1
answer
127
views
Confusion in Eigenvalue and eigenfunction distribution using NDEigensystem [duplicate]
I have been using NDEigensystem to solve the Schrodinger equation for different potentials. However, whenever the potential has its minimum in the negative y-axis, ...
- 2,836
2
votes
1
answer
528
views
Decomposing Hamiltonian in terms of Pauli strings
Does anyone know of a native Mathematica function or a package that allows one to decompose an arbitrary Hamiltonian into Pauli strings? There are a few packages for Python that do this (e.g., one ...
1
vote
0
answers
341
views
How to define momentum operator for use with creation and annihilation operators in quantum mechanics
I am having trouble writing the momentum operator in the differential form to complete a piece of code shown here so that I can show a commutator [P,Q] = I .
Here is what I have done so far, and I ...
- 165
-4
votes
1
answer
167
views
Constructing NxN tridiagonal block diagonal matrix with N is a variable [closed]
Hi please can someone help me writing a program of generating a tridiagonal block matrix with dimension of NxN where N is a variable can be 50, 100, 200, 300...etc , Well I don't know how I'm a ...
- 117
0
votes
1
answer
88
views
How to export integration values to a data file? [closed]
I'm calculating quantum mechanical expectation value for a very high order truncation value. But due to low ram, it gets hung each time. So I am trying to write those data to a .dat file into the HDD ...
- 2,836