Theoretical fractional formulation of a three-dimensional radio frequency ion trap (Paul-trap) for optimum mass separation

Abstract

We investigate the dynamics of an ion confined in a Paul–trap supplied by a fractional periodic impulsional potential. The Cantor–type cylindrical coordinate method is a powerful tool to convert differential equations on Cantor sets from cantorian–coordinate systems to Cantor–type cylindrical coordinate systems. By applying this method to the classical Laplace equation, a fractional Laplace equation in the Cantor–type cylindrical coordinate is obtained. The fractional Laplace equation is solved in the Cantor–type cylindrical coordinate, then the ions is modelled and studied for confined ions inside a Paul–trap characterized by a fractional potential. In addition, the effect of the fractional parameter on the stability regions, ion trajectories, phase space, maximum trapping voltage, spacing between two signals and fractional resolution is investigated and discussed.

Keywords

Paul–trap cantor-type cylindrical coordinate fractional parameter fraction ion motion mass separation

Introduction

A brief summary of ion traps and fractional calculus techniques are presented here.

Ion traps

Quadrupole ion traps were invented at the beginning of the 1950s^1,2,3 by Paul et al., demonstrating to be excellent tools to perform mass spectrometry.^4–10Other applications of quadrupole ion traps include quantum computing, ultraprecise atomic clocks, ion crystals, high–precision spectroscopy, fractional ion traps, and etc.^3–12 Moreover, the combined (Paul and Penning) trap^13,14 or the Kingdon trap¹⁵ can be successfully used to achieve mass spectrometry with very good results. Hu et al.¹⁸ proposed the Orbitrap that can be used as a multi-purpose mass spectrometer to examine different types of chemical systems. High resolution, high-mass accuracy and high dynamic range are interesting features of the Orbitrap.^18–20

The cylindrical geometry Paul–trap is easier to design and machine with respect to the hyperbolic geometry trap, and that is why it is increasingly attracting interests.^21,22 Experiments show that the cylindrical ion trap has a good resolution so as to perform mass separation of ions. In addition, its relatively simple geometry and small dimensions make it very suited for ion trapping experiments. Although it is possible to confine particles with distinct charge-to-mass ratios in a Paul trap, this occurs for weakly confined species that are expelled apart from the trap center. Akerman et al. studied the nonlinear mechanical response of a single laser-cooled ion confined in a linear RF–Paul trap,²³ demonstrating that both linear and the nonlinear damping components can be completely and accurately controlled. Mihalcea and Vişan²⁴ investigate the dynamics of an ion confined in a nonlinear Paul trap, which is shown to behave like a damped parametric oscillator that exhibits fractal properties and complex chaotic orbits.

In this paper, we studied about dynamics of a confined ion in a Paul–trap supplied by a fractional periodic potential. In this regard, the upcoming section studies the fractional Laplace equation in cantor-type cylindrical coordinate. In the next section, the fractional Laplace equation in the cantor-type cylindrical coordinate is modelled and studied. In a further section, ion motion inside a Paul-trap with fractional potential in the cantor-type cylindrical coordinates is modelled. The dynamical system consisting of an ion confined in a Paul–trap is investigated in the penultimate section, where numerical simulations are also performed. The effects of the fractional parameter on the stability regions, ion trajectories, phase space, maximum trapping voltage, spacing between two signals and fractional resolution are reviewed and discussed. In the final section, the results are analyzed and discussed.

History of fractional calculus

By looking at articles published in recent decades in the fields of science and engineering, we get acquainted with the topics of fractional calculus, differential equations with fractional derivatives, and concepts of this kind. So far, many books and papers in this field have been written from theoretical and practical points of view.^25–31 The subject of fractional calculus is more than 300 years old. The idea of fractional calculus dates back to the time of basic or classical calculus, and most theories about it were developed before the twentieth century. This was first introduced by Leibniz and L’Hospital’s in 1653.

In the twentieth century, many efforts were made by various scientists in this field. Caputo, by rewriting Riemann–Liouville formula, introduced a new derivative that is now used under the name Caputo derivative. Notable people who have worked on this topic during this period are: Hardy, Samko, Weyl, Riesz and Blair. Since 1970 until now, many people have studied in this field and also left useful articles and books. In this regard, Spanier, Oldham, Miller, Kilbas, Ross and Podlubny can be mentioned. The best resources for studying fractional calculations are books and articles of Miller and Ross, Kilbas and Podlubny. See Ross³² for a more comprehensive study of the history of fractional calculus.

Basic definitions and theorems of the fractional derivatives

Definitions of the fractional derivative of order $α > 0$ are presented in literature.^31,33–40 The Riemann-Liouville and Caputo fractional derivatives are the most used definitions in our paper.

Definition 1.1. For some $α \in ℝ^{+}$ , let n be the nearest integer greater than α. The Caputo fractional derivative of order α of a function $h (ξ)$ is given by,³¹

D_{*}^{α} h (ξ) = J^{n - α} \frac{d^{n}}{d ξ^{n}} h (ξ) = \frac{1}{Γ (n - α)} \int_{0}^{ξ} {(ξ - u)}^{n - α - 1} h^{(n)} (u) d u

(1)

with

n - 1 < α \leq n, n \in ℕ .

Theorem 1.2. The Riemann-Liouville derivative of order $α > 0$ with $n - 1 < α \leq n$ of the power function $f (ξ) = ξ^{β}$ for $β > 0$ is given by,³¹

D^{α} ξ^{β} = \frac{Γ (β + 1)}{Γ (β - α + 1)} ξ^{β - α}

(2)

Proof. Let $h (ξ) = ξ^{β}$ ( $β > 0$ ) then we have,

$\frac{d}{d ξ} (ξ^{β}) = β ξ^{β - 1} \Rightarrow \frac{d^{α}}{d ξ^{α}} ξ^{β} = \frac{β!}{(β - α)!} ξ^{β - α},$ replacing the factorials with the “gamma” function leads to,

D^{α} ξ^{β} = \frac{d^{α}}{d ξ^{α}} ξ^{β} = \frac{Γ (β + 1)}{Γ (β - α + 1)} ξ^{β - α}

(3)

Theorem 1.3. The Caputo derivative of order

α > 0

with

n - 1 < α \leq n

of the power function

f (ξ) = ξ^{β}

for

β > 0

satisfies,³¹

D_{*}^{α} ξ^{β} = {\begin{array}{l} \frac{Γ (β + 1)}{Γ (β - α + 1)} ξ^{β - α}; β > n - 1 \\ 0; β \leq n - 1 \end{array}

(4)

Proof. (see proof of Theorem (1.2)).

Fractional Laplace equation in the Cantor-type cylindrical coordinate

This section presents the fractional Laplace equation in the Cantor–type cylindrical coordinates. The Cantor–type cylindrical–coordinate method is a powerful tool to convert differential equations on cantor sets from cantorian–coordinate systems to Cantor–type cylindrical coordinate systems.

The cantorian–coordinate system was first described by Yang in 2010.^41,42 Both fractional and classical differential equations in the coordinate system to cartesian, cylindrical and spherical coordinates are convertible.^43,44 Newly, the cantorian–coordinate system is set on the fractals problems to obtain acceptable and accurate results. We consider the cantor–type cylindrical coordinates defined in references⁴²,⁴⁵ as, $x^{α} = r^{α} \cos_{α} θ^{α}, y^{α} = r^{α} \sin_{α} θ^{α}$ and $z^{α} = z^{α}$ , where $r \in (0, + \infty), z \in (- \infty, + \infty), θ \in (0, 2 π)$ and $r^{2 α} = x^{2 α} + x^{2 α}$ . The fractional $\sin_{α} θ^{α}$ and $\cos_{α} θ^{α}$ can be defined as follows,

\begin{array}{l} \sin_{α} θ^{α} = \sum_{n = 0}^{\infty} {(- 1)}^{n} \frac{θ^{(2 n + 1) α}}{Γ [1 + (2 n + 1) α]}, \\ \cos_{α} θ^{α} = \sum_{n = 0}^{\infty} {(- 1)}^{n} \frac{θ^{2 α n}}{Γ [1 + 2 α n]} \end{array}

(5)

Now, according to proposed equations and reference,⁴⁵ we can define the fractional gradient and fractional Laplace operators in the Cantor–type cylindrical coordinate system as follows,

\nabla^{α} Ψ (r, θ, z) = e_{r}^{α} \frac{\partial^{α} Ψ}{\partial r^{α}} + e_{θ}^{α} \frac{1}{r^{α}} \frac{\partial^{α} Ψ}{\partial θ^{α}} + e_{z}^{α} \frac{\partial^{α} Ψ}{\partial z^{α}}

(6)

\nabla^{2 α} Ψ (r, θ, z) = \frac{\partial^{2 α} Ψ}{\partial r^{2 α}} + \frac{1}{r^{2 α}} \frac{\partial^{2 α} Ψ}{\partial θ^{2 α}} + \frac{1}{r^{α}} \frac{\partial^{α} Ψ}{\partial r^{α}} + \frac{\partial^{2 α} Ψ}{\partial z^{2 α}}

(7)

where,

e_{r}^{α} = \cos_{α} θ^{α} e_{1}^{α} + \sin_{α} θ^{α} e_{2}^{α}, e_{θ}^{α} = - \sin_{α} θ^{α} e_{1}^{α} + \cos_{α} θ^{α} e_{2}^{α}

and

e_{z}^{α} = e_{3}^{α}

. Suggested fractional vector was given by,

R = r^{α} \cos_{α} θ^{α} e_{1}^{α} + r^{α} \sin_{α} θ^{α} e_{2}^{α} + z^{α} e_{3}^{α} = R_{r} e_{R}^{α} + R_{θ} e_{θ}^{α} + R_{z} e_{z}^{α}

Fractional Laplace equation in the Cantor-type cylindrical coordinate

This section focuses on the fractional Laplace equation in the Cantor–type cylindrical coordinates. The classical 3D Paul trap has a hyperbolic geometry, consisting of a ring and two end cap electrodes that present axial symmetry. In Figure 1, z₀ denotes the distance from the center of the Paul–trap to either of the endcap electrodes, while r₀ denotes the distance from the center of the Paul-trap to the nearest ring surface. Almost any geometry of trap electrodes with ac voltages applied between them, generating a saddle point in the potential, will cater a pseudo–potential minimum in which charged particles can be trapped.²⁰ All commonly used mass analyzers use electric and magnetic fields to apply force on charged particles.^1,13 This force causes the oscillating particle to move around the equilibrium point due to a fractional parabolic potential as, $Φ^{α} (x, y, z) = A (γ_{1} x^{2 α} + γ_{2} y^{2 α} + γ_{3} z^{2 α})$ . Any potential in free space should satisfy the fractional Laplace equation as,

▽^{2 α} Φ^{α} = \frac{\partial^{2 α}}{\partial x^{2 α}} Φ^{α} + \frac{\partial^{2 α}}{\partial y^{2 α}} Φ^{α} + \frac{\partial^{2 α}}{\partial z^{2 α}} Φ^{α} = 0

(8)

where

\frac{\partial^{2 α}}{\partial x^{2 α}} Φ^{α}, \frac{\partial^{2 α}}{\partial y^{2 α}} Φ^{α}

and

\frac{\partial^{2 α}}{\partial z^{2 α}} Φ^{α}

can be computed using the definitions of the fractional derivatives. According to the Theorem (1.2), when

α \to 2 α

we have,

\begin{array}{l} D^{2 α} ξ^{β} = \frac{Γ (β + 1)}{Γ (β - 2 α + 1)} ξ^{β - 2 α}; \\ \frac{n - 1}{2} < α \leq \frac{n}{2}, β > 0 \end{array}

(9)

let

ξ^{β} = x^{2 α}

(

2 α > 0

\begin{array}{l} D^{2 α} x^{2 α} = \frac{Γ (2 α + 1)}{Γ (2 α - 2 α + 1)} x^{2 α - 2 α} = \frac{Γ (2 α + 1)}{Γ (1)}; \\ \frac{n - 1}{2} < α \leq \frac{n}{2}, 2 α > 0 \end{array}

(10)

therefore, in the following we have,

\begin{array}{l} \frac{\partial^{2 α} Φ^{α}}{\partial x^{2 α}} = \frac{A γ_{1} Γ (2 α + 1)}{Γ (1)}, \\ \frac{\partial^{2 α} Φ^{α}}{\partial y^{2 α}} = \frac{A γ_{2} Γ (2 α + 1)}{Γ (1)}, \\ \frac{\partial^{2 α} Φ^{α}}{\partial z^{2 α}} = \frac{A γ_{3} Γ (2 α + 1)}{Γ (1)} \end{array}

(11)

from which we obtain,

▽^{2 α} Φ^{α} = \frac{\partial^{2 α} Φ^{α}}{\partial x^{2 α}} + \frac{\partial^{2 α} Φ^{α}}{\partial y^{2 α}} + \frac{\partial^{2 α} Φ^{α}}{\partial z^{2 α}} = A \frac{Γ (2 α + 1)}{Γ (1)} [γ_{1} + γ_{2} + γ_{3}]

(12)

Figure 1.

Schematic view of a r.f. Paul–trap.

Equation (13) shows that $γ_{1} + γ_{2} + γ_{3} = 0$ when $▽^{2 α} Φ^{α} = 0$ . For an ion trap, $γ_{1} = γ_{2} = 1$ and $γ_{3} = - 2$ and for a quadrupole mass filter $γ_{1} = - γ_{2} = 1$ and $γ_{3} = 0$ . In this paper we focused on the Paul–ion trap, then we assumed, $γ_{1} = γ_{2} = 1$ and $γ_{3} = - 2$ . Therefore, the fractional potential given as, $Φ^{α} (x, y, z) = A (x^{2 α} + y^{2 α} - 2 z^{2 α})$ . Using the standard transformations $x^{α} = r^{α} \cos_{α} θ^{α}, y^{α} = r^{α} \sin_{α} θ^{α}$ and $z^{α} = z^{α}$ , this equation can be transformed into the Cantor-type cylindrical coordinates. Hence, we can derive the fractional potential in the Cantor-type cylindrical coordinates as, $Φ^{α} (r, z) = A (r^{2 α} - 2 z^{2 α})$ , with $r^{2 α} = r^{2 α} \cos_{α}^{2} θ^{α} + r^{2 α} \sin_{α}^{2} θ^{α} = x^{2 α} + y^{2 α}$ .

This potential can be produced by four hyperbolic electrodes. To obtain this form of electrodes, we can consider the surfaces with same potentials $Φ_{0} / 2$ and $- Φ_{0} / 2$ as, $Φ^{α} (r_{0}, 0) = A (r_{0}^{2 α}) = Φ_{0} / 2$ and $Φ^{α} (0, z_{0}) = A (- 2 z_{0}^{2 α}) = - Φ_{0} / 2$ . With this conditions we can find $A = \frac{Φ_{0}}{2 r_{0}^{2 α}}$ and $A = \frac{Φ_{0}}{4 z_{0}^{2 α}}$ ; therefore, $r_{0}^{2 α} = 2 z_{0}^{2 α}$ . Thereby, the fractional electrodes shape for the fractional potential in the presented Cantor–type cylindrical coordinates given by, $Φ^{α} (r, z) = \frac{Φ_{0}}{2 r_{0}^{2 α}} (r^{2 α} - 2 z^{2 α}) = \pm 1$ . The applied electric potential, $Φ_{0}$ (that is applied to the hyperbolic rod’s) is either an r.f. potential $V_{0} \cos Ω t$ or a combination of a d.c. potential, U, of the form,^1,13 $Φ_{0} = U - V_{0} \cos Ω t$ , where $Ω = 2 π f$ is the angular frequency (in rad s⁻¹) of the r.f. field, and f is the frequency in hertz. Using the given definitions and information, the fractional potential $Φ^{α} (r, z)$ can be defined as, $Φ^{α} (r, z) = \frac{1}{2 r_{0}^{2 α}} (r^{2 α} - 2 z^{2 α}) (U - V_{0} \cos Ω t)$ .

The map of the electric field inside the trap and 3D simulations for $α = 0.9, 0.95, 1$ are shown in Figure 2(a) and (b), while 3D simulation for the classical Paul trap (α = 1) is shown in Figure 3.

Figure 2.

Field lines of electric fields; (a): α = 1, (b): $α = 0.95$ and (c): $α = 0.9$ .

Figure 3.

Ion trap simulation in 3D when α = 1, ring: $a x^{2} + b y^{2} - c z^{2} = 1$ and end cap: $a x^{2} + b y^{2} - c z^{2} = - 1$ with a = 1, b = 1, c = 2.

As can be seen in Figure 2(a) and (b), for α = 1, equation $x^{2} + y^{2} = r^{2}$ shows a circle, for $α = 0.95$ equation $x^{1.9} + y^{1.9} = r^{1.9}$ represents a smaller irregular circle compared with α = 1 and for $α = 0.9$ , equation $x^{1.8} + y^{1.8} = r^{1.8}$ indicates a smaller irregular circle compared to $α = 0.95$ and α = 1. In Figure 3, the ring and end cap equations are obtained from $a x^{2 α} + b y^{2 α} - c z^{2 α} = 1$ and $a x^{2 α} + b y^{2 α} - c z^{2 α} = - 1$ , respectively. Figure 4(a) and (b) indicate the contour lines for ring: $a x^{2α} + b y^{2α} - c z^{2α} = 1$ and end cap: $a x^{2α} + b y^{2α} - c z^{2α} = - 1$ , when $α = 1, 0.8, 0.6$ and a = 1, b = 1, c = 2. According to these figures, it can be concluded that by reducing α from 1 to 0.6, the contour lines along the axis $y = - x$ become more elongated.

Figure 4.

Contour lines when $α = 1, 0.8, 0.6$ and a = 1, b = 1, c = 2; (a): ring: $a x^{2 α} + b y^{2 α} - c z^{2 α} = 1$ and (b): end cap: $a x^{2 α} + b y^{2 α} - c z^{2 α} = - 1$ .

Fractional motion of trapped ions in the Paul-trap

In this section, the motion of ion inside a Paul–trap with the fractional potential in the Cantor-type cylindrical coordinates was modelled. The relationship between force, mass, and the applied fields in Newton’s second law and the Lorentz force law is as follows,

\begin{array}{l} F = ma (Newton' s second law) and \\ F = qE (Lorentz force law without magnetic field) \end{array}

(13)

where, F is the force applied to the ion, m is the mass of the ion, a is the acceleration, q is the ionic charge and E is the electric field. Here, F, a and E are vectorial variables. The electric field components in the trap with the fractional potential are as follows,

\begin{array}{l} [\begin{matrix} E_{r}^{α} \\ E_{z}^{α} \end{matrix}] = E^{α} = - \nabla^{α} Φ (r, z) = [\begin{matrix} - \frac{\partial^{α}}{\partial r^{α}} Φ (r, z) \\ - \frac{\partial^{α}}{\partial z^{α}} Φ (r, z) \end{matrix}] \\ = [\begin{matrix} - \frac{2^{2 α} Γ (α + \frac{3}{2})}{\sqrt{π} (2 α + 1) r_{0}^{2 α}} (U - V \cos Ω t) r^{α} \\ \frac{2^{2 α + 1} Γ (α + \frac{3}{2})}{\sqrt{π} (2 α + 1) r_{0}^{2 α}} (U - V \cos Ω t) z^{α} \end{matrix}] \end{array}

(14)

Therefore, the equations of motion for the only positive ion in the Paul–trap with the fractional potential in the Cantor–type cylindrical coordinates without using the magnetic field are given by,

{\begin{array}{l} \frac{d^{2} z}{d ξ} + (a_{z} - 2 q_{z} \cos 2 ξ) \frac{2^{2 α} Γ (α + \frac{3}{2})}{\sqrt{π} (2 α + 1)} z^{α} = 0, \\ \frac{d^{2} r}{d ξ} + (a_{r} - 2 q_{r} \cos 2 ξ) \frac{2^{2 α} Γ (α + \frac{3}{2})}{\sqrt{π} (2 α + 1)} r^{α} = 0 \end{array}

(15)

with

r_{0}^{2 α} = 2 z_{0}^{2 α}

and the assumptions,

ξ = \frac{Ω t}{2}, a_{z} = - 2 a_{r} = - \frac{4 q U}{m r_{0}^{2 α} Ω^{2}}, q_{z} = - 2 q_{r} = \frac{2 q V}{m r_{0}^{2 α} Ω^{2}}

(16)

Assuming α = 1, basic motion equations are as follows,

{\begin{array}{l} \frac{d^{2} z}{d ξ^{2}} + (a_{z} - 2 q_{z} \cos 2 ξ) z = 0, \\ \frac{d^{2} r}{d ξ^{2}} + (a_{r} - 2 q_{r} \cos 2 ξ) r = 0 \end{array}

(17)

Programming and numerical simulations

In this section, programming and numerical simulations of the dynamical system for the trapped ion inside Paul–trap are investigated and discussed. For programming and numerical simulations, the charge state of +1 was considered. We first plot stability regions in the (a, q) and $(V, - U)$ plans, ion trajectories in time, the evolution of phase space ion path, resolution of the ion trap and fractional resolution of the ion trap. Then, we study and discuss the effect of the fractional potential on the mass resolution. The effect of the fractional potential was examined for ions of ¹³¹Xe and ¹³²Xe.

Ion trajectories

Figure 5 shows the first stability region of Paul–trap with the fractional potential when $α = 1, 0.95, 0.9$ . As can be seen, changing the fractional parameter α from 1 to 0.9, first stability region will be smaller along the a axis and bigger along the q axis. The stability diagrams $(V, - U)$ plane for ¹³¹Xe with $Ω = 2 π \times 1.05 \times 10^{6}$ rad/s, $z_{0} = 0.707$ mm and $α = 1, 0.95, 0.9$ have been shown in Figure 6. When the fractional parameter α decreases from 1 to 0.9, the stability diagrams $(V, - U)$ plane become larger. Figure 7 presents the ion trajectories in time when a_z = 0, $q_{z} = 0.9$ and $α = 0.8, 0.9, 1$ , respectively. This figure show that the ion trajectories are comparable for all values of the fractional parameters $α = 0.8, 0.9$ and 1. However, as the value of the parameter α decreases, the ion rotation space increases.

Figure 5.

The first stability region of Paul–trap when $α = 1, 0.95, 0.9$ .

Figure 6.

The stability diagram in $(V, - U)$ plan for ¹³¹Xe with $Ω = 2 π \times 1.05 \times 10^{6}$ rad/s, $z_{0} = 0.707$ mm and $α = 1, 0.95, 0.9$ .

Figure 7.

The ion trajectories in time for $q_{z} = 0.9$ and $α = 0.8, 0.9, 1$ .

The ion trajectories in $z - \dot{z}$ plane for $q_{z} = 0.32$ , $α = 0.8, 0.9, 1$ are shown in Figure 8. The left side and right sides of this figure show the $\max (z)$ and $\max (\dot{z})$ versus fractional parameter α. As can be seen in the right side, the rotation space of the ions increases as the value of α decreases from 1 to 0.8. As the left side shows, by decreasing the value of α from 1 to 0.6, the values of $\max (z)$ and $\max (\dot{z})$ increase, but $\max (z)$ is increasing faster than $\max (\dot{z})$ . As Figure 8 shows, there are two periodic attractors in the system, which are corresponding to forced oscillations confined to the left or right well. The portraits of the phase obviously reflect the existence of one or two attractors and of fractal basin boundaries for the trapped ion, assimilated with a periodically forced double well oscillator. The system can converge rapidly to one of the two attractors, based on the initial conditions and fractional parameter α. Generally the attraction basins have a complicated shape, and the boundary between them is fractal.^24,46

Figure 8.

Left: ion trajectories in $z - \dot{z}$ plane for $q_{z} = 0.32$ and $α = 0.8, 0.9, 1$ ; Right: $\max (z)$ and $\max (\dot{z})$ vs α.

Figure 9 shows the mechanical properties of the confined ions analyzed through the ion displacements in the phase space. Phase space ion trajectory for different values of r.f. fields with initial phases $ξ_{0} = \frac{π}{4}$ and $ξ_{0} = - \frac{π}{4}$ has been proposed for $α = 0.9$ and α = 1. The computational results in this figure show the comparable phase space for different values of fractional parameter $α = 0.9$ and α = 1.

Figure 9.

Evolution of phase space ion trajectory for different values of the phase $ξ_{0} = \frac{π}{4}$ and $ξ_{0} = - \frac{π}{4}$ when $α = 0.9$ and α = 1.

Effect of the fractional factor on the mass resolution

This section presents the effect of fractional parameter α on the mass resolution of trapped ions. As we know, the resolution of a Paul–trap mass spectrometer is a function of the mechanical accuracy of the hyperboloid of the ion trap, $Δ r_{0}$ and the stability performances of the electronics device, such as variations in voltage amplitude $Δ V$ and the r.f. frequency $Δ Ω$ .⁴⁶ The computational resolution will tell us how accurate the form of the voltage signal is. To derive a theoretical formula for the fractional resolution, according to equation (17), there will be,

q_{z} = \frac{2 e V}{m r_{0}^{2 α} Ω^{2}}

(18)

Considering the partial derivatives on the variables of the stability parameters, expression of the resolution $Δ m$ can be computed as follows,

Δ m = 2 α (\frac{2 e V}{r_{0}^{2 α + 1} Ω^{2} q_{z}}) | Δ r_{0} | + (\frac{2 e}{r_{0}^{2 α} Ω^{2} q_{z}}) | Δ V | + 2 (\frac{2 e V}{r_{0}^{2 α} Ω^{3} q_{z}}) | Δ Ω |

(19)

then, there will be,

Δ m = 2 α (\frac{2 e V}{r_{0}^{2 α} Ω^{2} q_{z}}) | \frac{Δ r_{0}}{r_{0}} | + (\frac{2 e V}{r_{0}^{2 α} Ω^{2} q_{z}}) | \frac{Δ V}{V} | + 2 (\frac{2 e V}{r_{0}^{2 α} Ω^{2} q_{z}}) | \frac{Δ Ω}{Ω} |

(20)

therefore, we have,

Δ m = m (2 α | \frac{Δ r_{0}}{r_{0}} | + 2 | \frac{Δ Ω}{Ω} | + | \frac{Δ V}{V} |)

(21)

Thus, the fractional resolution is given by,

\frac{m}{Δ m} = {(2 α | \frac{Δ r_{0}}{r_{0}} | + 2 | \frac{Δ Ω}{Ω} | + | \frac{Δ V}{V} |)}^{- 1}

(22)

Uncertainties $| \frac{Δ Ω}{Ω} | = 10^{- 7}$ , $| Δ V | = 10^{- 4}$ and $| Δ r_{0} | = 10^{- 3}$ have been used for voltage, r.f. and geometry for fractional mass resolution,⁴⁶ respectively. Assuming that $Ω = 2 π \times 1.05 \times 10^{6}$ rad/s, $z_{0} = 0.707$ , a = 0, maximum values of voltage V, V_max, as a function of the fractional parameter, α, and function of ion mass, m, for ¹³¹Xe and 1¹³²Xe when m = 131, 132 and $α = 0.9, 1$ ,⁴⁶ are presented in Figure 10. As can be seen, by increasing the fractional parameter a from 0.55 to 1, the maximum voltage, V_max, decreases rapidly like a negative exponential function. Figure 11 shows the spacing between two signals, $Δ m$ , as a function of the fractional parameter α. As can be seen, by increasing the fractional parameter α from 0.55 to 1, the spacing between two signals, $Δ m$ , decreases rapidly like an exponential function. This means that, by reducing the fraction parameter α, the separation can be performed more accurately.

Figure 10.

Maximum values of V, V_max, as a function of fractional parameter α, when $Ω = 2 π \times 1.05 \times 10^{6}$ rad/s, $z_{0} = 0.707$ mm and a = 0.

Figure 11.

Spacing between two signals, $Δ m$ , as a function of fractional parameter α.

Figure 12, indicates the fractional mass resolution $m / Δ m$ as a function of fractional parameter α. The results of this figure show that by decreasing the fraction parameter α from 1 to 0.55, the fractional mass resolution values rapidly increase from 400 to 1200. The higher fractional mass resolution indicates better and more accurate separation. In Figures 10 to 12, to find the vertical values, we divided the interval $α = [0.55, 1]$ into N = 45 parts using the stepsize h = 0.01, then V_max, $Δ m$ and $m / Δ m$ values were found in all these points. Then, all the curves were plotted using the 45 found points, but to make the curves easily visible, the markers have been used only in ten points.

Figure 12.

The fractional resolution of ion trap, $m / Δ m$ , as a function of fractional parameter α.

Maximum values of V, V_max, as a function of ion mass, m, for hypothetical values $Ω = 2 π \times 1.05 \times 10^{6}$ rad/s, $z_{0} = 0.707 mm$ , a = 0 and $α = 0.6, 0.8, 1$ is shown in Figure 13. This figure also shows the maximum values of voltage V for the ions $^{131} X e$ and $^{132} X e$ when the fractional values are $α = 0.6, α = 0.8$ and α = 1. As can be seen, the maximum voltage for $α = 0.6$ and $α = 0.8$ is less than the maximum voltage for α = 1; and lower voltage indicates better and more accurate separation. Figure 14 presents the spacing between two signals, $Δ m$ , as a function of ion mass, m, when the fractional parameters are $α = 0.6, α = 0.8$ and α = 1. This figure also shows the values of spacing between two signals for the ions $^{131} X e$ and $^{132} X e$ . As can be seen, ${(Δ m)}_{α = 0.6} < {(Δ m)}_{α = 0.8} < {(Δ m)}_{α = 1}$ , and less $Δ m$ indicates better and more accurate separation. In Figures 13 and 14, to find the vertical values, we divided the interval $m = [130.5, 132.5]$ into N = 40 parts using the stepsize h = 0.05, then the values of V_max and $Δ m$ were found in all these points. All the curves were plotted using the 45 found values, but to make the curves easily visible, the markers have been used only in five points.

Figure 13.

Maximum values of V, V_max, as a function of ion mass, m, when $Ω = 2 π \times 1.05 \times 10^{6}$ rad/s, $z_{0} = 0.707$ mm, a = 0 and $α = 0.6, 0.8, 1$ .

Figure 14.

Spacing between two signals, $Δ m$ , as a function of ion mass, m, when $α = 0.6, 0.8, 1$ .

Conclusion

A modified three-dimensional radio frequency Paul–trap with fractional potential was introduced in this study. The first stability region in (q, a) and $(V, - U)$ planes was also shown. Moreover, effect of fractional parameter α on the mass separation was studied. Maximum values of voltage, V_max, as a function of the fractional parameter α for was derived for ions ¹³¹Xe and ¹³²Xe assuming that $Ω = 2 π \times 1.05 \times 10^{6}$ rad/s, $z_{0} = 0.707$ mm and a = 0. Further, the spacing between two signals, $Δ m$ , and mass fractional resolution, $m / Δ m$ , for ions ¹³¹Xe and ¹³²Xe as a function of the fractional parameter α was studied and discussed. The fractional resolution of ion traps $m / Δ m$ increases when the fractional parameter α decreases. As was observed, with decreasing the fractional parameter a from 1 to 0.55, the fractional mass resolution rapidly increased from 400 to 1200. The high fractional resolution in good separation has high mass accuracy. As shown, the maximum voltage for $α = 0.6$ and $α = 0.8$ was less than the maximum voltage for α; and lower voltage indicates better and more accurate separation. The general results of this paper showed that the fractional parameter α can be an important and effective controller to optimize ion mass separation.

Supplemental Material

sj-pdf-1-ems-10.1177_14690667211026790 - Supplemental material for Theoretical fractional formulation of a three-dimensional radio frequency ion trap (Paul-trap) for optimum mass separation

Supplemental material, sj-pdf-1-ems-10.1177_14690667211026790 for Theoretical fractional formulation of a three-dimensional radio frequency ion trap (Paul-trap) for optimum mass separation by Sarkhosh Seddighi Chaharborj, Shahriar Seddighi Chaharborj, Zahra Seddighi Chaharborj and Pei See Phang in European Journal of Mass Spectrometry

Supplemental Material

sj-pdf-2-ems-10.1177_14690667211026790 - Supplemental material for Theoretical fractional formulation of a three-dimensional radio frequency ion trap (Paul-trap) for optimum mass separation

Supplemental material, sj-pdf-2-ems-10.1177_14690667211026790 for Theoretical fractional formulation of a three-dimensional radio frequency ion trap (Paul-trap) for optimum mass separation by Sarkhosh Seddighi Chaharborj, Shahriar Seddighi Chaharborj, Zahra Seddighi Chaharborj and Pei See Phang in European Journal of Mass Spectrometry

Footnotes

Acknowledgement

The authors are pleased to acknowledge the support of the National Sciences and Engineering Research Council of Canada.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by the National Sciences and Engineering Research Council of Canada.

ORCID iD

Sarkhosh Seddighi Chaharborj

References

Paul

(1990). Electromagnetic traps for charged and neutral particles. Reviews of modern physics, 62(3), 531.

Paul

(1990). Electromagnetic traps for charged and neutral particles (Nobel lecture). Angewandte Chemie International Edition in English, 29(7), 739–748.

Major

F. G.

Major

S. P. F. G.

Gheorghe

V. N.

Werth

(2005). Charged particle traps: physics and techniques of charged particle _eld con_nement (Vol. 37). Springer Science & Business Media.

March

R. E.

Todd

J. F.

(Eds.). (1995). Practical aspects of ion trap mass spectrometry: Chemical, environmental, and biomedical applications (Vol. 3). CRC press.

Knoop

Madsen

Thompson

R. C.

(2014). Chapter 1: Physics with Trapped Charged Particles. In Physics with Trapped Charged Particles: Lectures from the Les Houches Winter School (pp. 1–24). https://hal.archives-ouvertes.fr/hal-01455335/

Thompson

R. C.

Madsen

Knoop

(Eds.). (2016). Trapped Charged Particles: A Graduate Textbook with Problems and Solutions (Vol. 1). World Scienti_c.

Knoop

Madsen

Thompson

R. C.

(2014). Physics with Trapped Charged Particles: Lectures from the Les Houches Winter School. World Scienti_c.

Kozlov

M. G.

Safronova

M. S.

Lpez-Urrutia

J. C.

Schmidt

P. O.

(2018). Highly charged ions: optical clocks and applications in fundamental physics. Reviews of Modern Physics, 90(4), 045005.

Snyder

D. T.

Pulliam

C. J.

Wiley

J. S.

Duncan

Cooks

R. G.

(2016). Experimental characterization of secular frequency scanning in ion trap mass spectrometers. Journal of The American Society for Mass Spectrometry, 27(7), 1243–1255.

10.

Cornejo

J. M.

Gutirrez

M. J.

Ruiz

Bautista-Salvador

Ospelkaus

Stahl

Rodrguez

(2016). An optimized geometry for a micro Penning–trap mass spectrometer based on interconnectedions. International Journal of Mass Spectrometry, 410, 22–30.

11.

Noll

R. J.

Makarov

Hardman

Graham Cooks

(2005). The Orbitrap: a new mass spectrometer. Journal of mass spectrometry, 40(4), 430–443.

12.

Rychlik

Kanawati

Roullier-Gall

Hemmler

Liu

Alexandre

, … & Schmitt-Kopplin

(2019). Foodomics assessed by Fourier transform mass spectrometry. In Fundamentals and applications of Fourier transform mass spectrometry (pp. 651–677). Elsevier.

13.

Latorre

Bollen

Gulyuz

Ringle

Bado

Dugan

, … & Translume Collaboration. (2016, September). Micro penning trap for continuous magnetic _eld monitoring in high radiation environments. In APS division of nuclear physics meeting abstracts (Vol. 2016, pp. EA–085).

14.

Shang

J. J.

Cui

K. F.

Cao

Wang

S. M.

Chao

S. J.

Shu

H. L.

Huang

X. R.

(2016). Sympathetic Cooling of 40Ca+-27Al+ Ion Pair Crystal in a Linear Paul Trap. Chinese Physics Letters, 33(10), 103701.

15.

Kingdon

K. H.

(1923). A method for the neutralization of electron space charge by positive ionization at very low gas pressures. Physical Review, 21(4), 408.

16.

Chaharborj

S. S.

Kiai

S. M. S.

Ari_n

N. M.

Gheisari

(2015). Optimum Radius Size between Cylindrical Ion Trap and Quadrupole Ion Trap. Mass Spectrometry Letters, 6(3), 59–64.

17.

Mather

R. E.

Waldren

R. M.

Todd

J. F. J.

March

R. E.

(1980). Some operational characteristics of a quadrupole ion storage mass spectrometer having cylindrical geometry. International Journal of Mass Spectrometry and Ion Physics, 33(3), 201–230.

18.

Akerman

Kotler

Glickman

Dallal

Keselman

Ozeri

(2010). Single–ion nonlinear mechanical oscillator. Physical Review A, 82(6), 061402.

19.

Mihalcea

B. M.

Visan

G. G.

(2010). Nonlinear ion trap stability analysis. Physica Scripta, 2010( T140), 014057.

20.

Chaharborj

S. S.

Moameni

(2017). Applications of the fractional calculus to study the physical theory of ion motion in a quadrupole ion trap. European Journal of Mass Spectrometry, 23(5), 254–271.

21.

Caputo

(1967). Linear models of dissipation whose Q is almost frequency independentII. Geophysical Journal International, 13(5), 529–539.

22.

Miller

K. S.

Ross

(1993). An introduction to the fractional calculus and fractional di_erential equations. Wiley.

23.

Podlubny

(1999). Fractional di_erential equations, mathematics in science and engineering.

24.

Podlubny

(2001). Geometric and physical interpretation of fractional integration and fractional differentiation. arXiv preprint math/0110241.

25.

Diethelm

(2010). The analysis of fractional di_erential equations: An application-oriented exposition using di_erential operators of Caputo type. Springer Science & Business Media.

26.

Luchko

(2020). Fractional derivatives and the fundamental theorem of fractional calculus. Fractional Calculus and Applied Analysis, 23(4), 939–966.

27.

Ross

(1977). The development of fractional calculus 1695–1900. Historia Mathematica, 4(1), 75–89.

28.

Miller

KS.

Ross

An introduction to the fractional calculuse and fractional di_erential equations, New York: John Wiley and Sons Inc., 1993.

29.

Chaharborj

S. S.

Chaharborj

S. S.

Mahmoudi

(2017). Study of fractional order integrodi_erential equations by using Chebyshev Neural Network. J. Math. Stat, 13(1), 1–13.

30.

Z. W.

Trujillo

J. J.

Q. Y.

(2016). Riemann–Liouville fractional calculus in Morrey spaces and applications. Computers & Mathematics with Applications. https://doi.org/10.1016/j.camwa.2016.04.013.

31.

Bhardwaj

Kumar

(2020). Numerical solution of time fractional tricomi-type equation by an rbf based meshless method. Engineering Analysis with Boundary Elements, 118, 96–107.

32.

Radmanesh

Ebadi

M. J.

(2020). A local mesh-less collocation method for solving a class of time-dependent fractional integral equations: 2D fractional evolution equation. Engineering Analysis with Boundary Elements, 113, 372–381.

33.

Bruzewicz

C. D.

Chiaverini

McConnell

Sage

J. M.

(2019). Trapped–ion quantum computing: Progress and challenges. Applied Physics Reviews, 6(2), 021314.

34.

Brown

K. R.

Kim

Monroe

(2016). Co-designing a scalable quantum computer with trapped atomic ions. arXiv preprint arXiv:1602.02840.

35.

Schmidt

P. O.

Leroux

I. D.

(2016). Trapped-Ion Optical Frequency Standards. In Trapped Charged Particles: A Graduate Textbook with Problems and Solutions, 377-425. World Scienti_c.

36.

Ishteva

M. A. R. I. Y. A.

Boyadjiev

L. Y. U. B. O. M. I. R.

Scherer

R. U. D. O. L. F.

(2005). On the Caputo operator of fractional calculus and C-Laguerre functions. Mathematical Sciences Research Journal, 9(6), 161.

37.

Kolwankar

K. M.

Gangal

A. D.

(1996). Fractional di_erentiability of nowhere di_erentiable functions and dimensions. Chaos: An Interdisciplinary Journal of Nonlinear Science, 6(4), 505–513.

38.

Kolwankar

K. M.

(2013). Local fractional calculus: a review. arXiv preprint arXiv:1307.0739.

39.

Yang

X. J.

(2011). Local fractional functional analysis and its applications. Asian Academic Publisher, Hong Kong.

40.

Yang

X. J.

(2012). Advanced local fractional calculus and its applications. World Science, New York, NY, USA.

41.

Boothby

W. M.

(2002). An introduction to di_erentiable manifolds and Riemannian geometry (Vol. 120). Academic press.

42.

Ganji

Z. Z.

Ganji

D. D.

Ganji

A. D.

Rostamian

(2010). Analytical solution of time-fractional NavierStokes equation in polar coordinate by homotopy perturbation method. Numerical Methods for Partial Di_erential Equations, 26(1), 117–124.

43.

Yang

X. J.

Srivastava

H. M.

J. H.

Baleanu

(2013). Cantor-type cylindrical-coordinate method for di_erential equations with local fractional derivatives. Physics Letters A, 377(28), 1696–1700.

44.

Davila–Rodriguez

Ozawa

Hnsch

T. W.

Udem

(2016). Doppler cooling trapped ions with a UV frequency comb. Physical review letters, 116(4), 043002.

45.

Moon

F. C.

(1992). Chaotic and Fractal Dynamicsan introduction for Applied Scientists and Engineers John Wiley & Sons. Inc., New York.

46.

Kiai

S. S.

Baradaran

Adlparvar

Khalaj

M. M. A.

Doroudi

Nouri

, … & Babazadeh

A. R.

(2005). Impulsional mode operation for a Paul ion trap. International Journal of Mass Spectrometry, 247(1), 61–66.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.54 MB

0.25 MB