Optimised surface-electrode ion-trap junctions for experiments with cold molecular ions

TL;DR

This paper presents the design and optimization of surface-electrode ion-trap junctions for experiments with cold molecular ions, enabling advanced chip-based trapping techniques for various scientific applications.

Contribution

It introduces a multi-objective optimization procedure for planar ion-trap junctions and proposes a multi-functional chip for cold molecular ion experiments.

Findings

Optimized junction designs for quadrupole and octupole configurations.

Proposed a multi-functional ion-trap chip for millikelvin temperature molecular ions.

Opened pathways for complex chip-based trapping in quantum technology and chemistry.

Abstract

We discuss the design and optimisation of two types of junctions between surface-electrode radiofrequency ion-trap arrays that enable the integration of experiments with sympathetically cooled molecular ions on a monolithic chip device. A detailed description of a multi-objective optimisation procedure applicable to an arbitrary planar junction is presented, and the results for a cross junction between four quadrupoles as well as a quadrupole-to-octupole junction are discussed. Based on these optimised functional elements, we propose a multi-functional ion-trap chip for experiments with translationally cold molecular ions at temperatures in the millikelvin range. This study opens the door to extending complex chip-based trapping techniques to Coulomb-crystallised molecular ions with potential applications in mass spectrometry, spectroscopy, controlled chemistry and quantum technology.

Tables

Table 1. Table 1: Optimal electrode configurations for a typical five-wire surface-electrode trap as the asymptotic quadrupolar channel for the cross junction. a 𝑎 a and b 𝑏 b are the RF electrode separation and width, respectively, h ℎ h is the trapping height above the surface. See text for discussion.

Reference	a	$b/a$	Relative trap depth at fixed $h$	Relative channel curvature at fixed $h$	Method
[40]	1.09 $h$	1.19	68%	83%	maximum depth at fixed $a$
[36, 41]	0.69 $h$	3.68	100%	94%	maximum depth at fixed $h$
[4, 36], this work	0.83 $h$	2.41	95%	100%	maximum curvature at fixed $h$

Table 2. Table 2: Specification of the optimised cross junction (figure 2 (a)), see the labelled cubic splines in figure  9 (a). a 𝑎 a and b 𝑏 b are the RF electrode separation and width for the optimal SE quadrupole, respectively (table 1 ). Note that { x a , y a } subscript 𝑥 𝑎 subscript 𝑦 𝑎 \{x_{a},y_{a}\} and { x b , y b } subscript 𝑥 𝑏 subscript 𝑦 𝑏 \{x_{b},y_{b}\} were constrained to remain on the line of equality (i.e., x = y 𝑥 𝑦 x=y ), and thus the geometry is specified by 12 degrees of freedom. The length scales are in units of the trapping hight h ℎ h , and the unit of angles is radian.

Cubic spline 1	$\{x_a,y_a,v_a,\alpha \}$	$\{x_c,y_c,v_c,\gamma \}$
	$\{0.004,0.004,2.242,0.353\}$	$\{0.643,$ fixed at $a/2,0.633,3.091\}$
Cubic spline 2	$\{x_b,y_b,v_b,\beta \}$	$\{x_d,y_d,v_d,\delta \}$
	$\{1.562,1.562,0.064,-0.172\}$	$\{3.690,$ fixed at $a/2+b,1.25098,1.637\}$

Table 3. Table 3: Specification of the quadrupole-to-octupole junction, see figure  9 (b). y i subscript 𝑦 𝑖 y_{i} are the RF-electrode edges positions obtained from the optimisation of the SE octupole discussed in appendix B, and a 𝑎 a and b 𝑏 b are the RF electrode separation and width for the quadrupole guide, respectively. Note that due to the constraints used, e.g., x a = x d = x e subscript 𝑥 𝑎 subscript 𝑥 𝑑 subscript 𝑥 𝑒 x_{a}=x_{d}=x_{e} , the geometry was specified by 14 degrees of freedom. The length scales are in units of the trapping hight h ℎ h , and the unit of angles is radian.

Cubic spline 1	$\{x_a,y_a,v_a,\alpha \}$	$\{x_c,y_c,v_c,\gamma \}$
	$\{3.969,$ fixed at $y_4,18.420,2.008\}$	$\{$fixed at $0,$ fixed at $a/2+b,2.539,3.141\}$
Cubic spline 2	$\{x_b,y_b,v_b,\beta \}$	$\{x_d,y_d,v_d,\delta \}$
	$\{$fixed at $0,$ fixed at $a/2,5.204,0.794\}$	$\{3.969,$ fixed at $y_3,4.107,3.024\}$
Cubic spline 3	$\{x_e,y_e,v_e,\eta \}$	$\{x_f,y_f,v_f,\varphi \}$
	$\{3.969,$ fixed at $y_2,0.0,1.452\}$	$\{0.215,$ fixed at $y_1,1.634,-1.127\}$

Equations

$$F_1$$

$$F_2$$

$$F_3$$

$$U$$

$${\xi }_i=\frac{1}{F_i^0},$$

$$Q_{ii}(x)=\frac{2q{V}_{\mathrm{RF}}}{m{\Omega }_{\mathrm{RF}}^2}\frac{{\partial }^2{\Theta }_{\mathrm{RF}}}{{\partial x_i}^2}({\mathbf{r}}_{\min }(x)),$$

$$\begin{gathered} \text{cases}\text{eqalign}\underset{{\mathrm{y}}_1,\dots ,{\mathrm{y}}_{\mathrm{p}}}{\mathrm{maximise}}\frac{{\partial }^n{\Theta }_{\mathrm{RF}}(y_1,...,y_p)}{\partial z^n} \\ \mathrm{subject}\mathrm{to}\frac{{\partial }^{i+j+k}{\Theta }_{\mathrm{RF}}(y_1,...,y_p)}{\partial x^i\partial y^j\partial z^k}=0, \end{gathered}$$

$$F(d_1,d_2,d_3,d_4)={\lambda }_1{f}^{(1)}+{\lambda }_2{f}^{(2)}+{\lambda }_3{f}^{(3)}-{\lambda }_4{f}^{(4)},$$

Full text

1 Introduction

Two-dimensional ion-trap arrays have become important tools for the development of large-scale quantum information processing [1, 2, 3] and quantum simulation [4, 5, 6] based on atomic ions. In particular, surface-electrode (SE) radiofrequency (RF) ion traps in which all electrodes lie in a plane [7] offer the possibility of integrating multiple trapping zones and of precisely shaping electrodes, and therefore, trapping potentials. The scalability of SE structures offered by microfabrication methods paves the way for novel large-scale experiments [8].

Translationally cold and quantum-state controlled molecular ions, on the other hand, are promising platforms for fundamental studies such as the test of the time invariance of physical constants [9], as well as for applications ranging from precision spectroscopy [10, 11, 12, 13] and metrology [14] to collision studies [15] and controlled chemistry [16, 17]. The workhorse technique for preparing such systems is sympathetic cooling with laser-cooled co-trapped atomic ions [18] in RF traps [19]. This method enables the formation of Coulomb-crystallised molecular ions at temperatures in the millikelvin range[18]. Recently, the sympathetic cooling of molecular ions in a SE RF trap has been demonstrated [20]. The extension of SE trapping technology to cold molecular ions opens up exciting possibilities, e.g., the implementation of miniaturised guided-ion-beam experiments which could combine mass spectrometry, spectroscopy, reaction studies, and cold chemistry experiments on a chip device.

In order to develop an ion-trapping network, the key ingredient is junctions. Usually, a junction is thought of as a connection between identical trap arrays, e.g., a T-, Y-, or X- (cross) junction connecting quadrupolar channels [21, 22, 23, 24]. In addition, a junction can be seen as an element that enables the modification of trapping fields between two different trapping configurations, e.g., a quadrupole-to-octupole junction which can be of interest for, e.g., connecting a standard quadrupole ion transport channel with a large-volume octupolar ion storage region. Such field-modifying junctions allow the manipulation of the ion-crystal structures as well as the separation of a single harmonic well into a double-well trapping configuration with potential applications in quantum information and quantum matter-wave experiments [25].

The main challenge to the use of SE junctions arises from the fact that the lowest-order multipole component of the RF field at the centre of a given intersection is a hexapole determined by the tangents of the intersecting channels [26]. In RF traps, the trapping potential is usually calculated within the adiabatic approximation [27] in which a time-independent pseudopotential derived from the RF field governs the confinement of ions [19]. Straight RF electrodes intersecting at right angles do not provide pseudopotential confinement perpendicular to the paths due to field cancellation. The pseudopotential at the centre becomes stiffer as the angle between straight arms is reduced and also as the number of intersecting channels decreases. For this reason, SE Y-junctions [23, 28] were initially preferred to SE X-junctions [24].

X-junctions are particularly important for applications which need to economize space, e.g., highly integrated experiments. Recently, efforts have been invested in developing such junctions resulting in their successful demonstration and in new techniques for better controlling ion shuttling operations [24, 29]. One of the key elements for these achievements is the optimisation of the RF-electrode geometry in the vicinity of the intersections such that ions are confined along the direction perpendicular to the surface [24].

The ion-trap junctions that have been developed thus far were designed to enable the shuttling of single-species atomic ion crystals while minimising ion heating during transport [23, 28, 24, 30]. Besides the applications in molecular physics mentioned above, the extension of this technology to bicomponent Coulomb crystals is also important for quantum technology and quantum computation with mixed-species crystals [31]. The challenge for the transportation of bicomponent crystals through junctions arises from the mass dependence of the pseudopotential.

In this paper, we describe a general method for the optimisation of an arbitrary planar junction in accordance with the requirements for the transportation of mixed-species ion crystals (section 2). This approach is applied to design a cross between four quadrupoles and to a quadrupole-to-octupole junction discussed in detail in section 3. Based on this study, a design of a multi-functional ion-trap array which enables the integration of several experiments on a single-layer chip is presented (section 4).

2 Multi-objective optimisation of electrode structures

2.1 Objective functions and weighted sum method

Transporting ions through SE junctions requires the minimisation of several objectives. A variety of objective functions has been considered [23, 28, 24, 32], and thus, the question arises which physical quantities are relevant when optimising an electrode structure for a given application. First and foremost, the suppression of pseudopotential barriers that impede ion transport is crucial [26]. Note that the position of pseudopotential minima, where the RF field vanishes, solely depends on the trap geometry, and thus is identical for all ion masses. However, the mass dependence of the pseudopotential may cause a segregation of laser-cooled and sympathetically cooled ions during the transportation over barriers which could lead to insufficient cooling and high loss rates of molecular ions. Thus, for shuttling mixed-species ion crystals the height of pseudopotential barriers is the most important quantity to be minimised (objective 1).

For reliable shuttling of ion qubits, axial pseudopotential gradients that cause ion heating must be suppressed [23, 28, 24, 30, 29], since at intersections, where noise on the RF potential is present, axial motional modes of ions might be excited [22]. Thus, for applications in quantum computation the minimisation of the pseudopotential gradient is crucial (objective 2). An alternative approach to near-motional-ground-state shuttling of ions is the minimisation of the variation of the axial potential curvature [29]. We note that this optimisation objective is less important in the present context as the molecular ions are continuously sympathetically cooled by laser-cooled ions, also during shuttling.

Furthermore, the variation of the trapping height needs to be controlled for two reasons. First, weakly confining potentials at SE intersections might cause significant ion losses. This issue can be mitigated by using static fields that forces ions to circumnavigate the junction centre at the expense of moving the trapping minima away from the pseudopotential null [24]. This technique is unsuited for bicomponent crystals due to the risk of the drastic segregation of ionic species under the influence of such a static field [33]. Second, for the transportation of molecular ions together with laser cooled atomic ions in the presence of cooling beams large variations of the trapping height should be avoided (objective 3).

These three objectives evaluated at the ion channel ${\mathbf{r}_{\rm min}(x)=\mathop{\mathrm{arg\,min}}\limits_{y,z}({\Phi}_{\rm ps}(x,y,z))}$, where $\mathop{\mathrm{arg\,min}}$ represents the values of y and z that minimise ${\Phi}_{\rm ps}(x,y,z)$ at a certain x, can be cast into the following mathematical form:

[TABLE]

$F_{i}$ ($i$=1, 2, 3) denotes the $i$th objective function, $\mathit{\Phi}_{\rm ps}$ is the pseudopotential, $L$ is the coordinate of the initial trapping zone, $z_{\rm min}$ is the component of ${\mathbf{r}_{\rm min}(x)}$ perpendicular to the chip surface, and $h$ is the intended trapping height. Note that the minimisation of pseudopotential barriers might occur at the cost of large trapping height variations and vice versa.

An efficient method of solving multi-objective optimisation problems is the weighted sum method in which a function $U$ is minimised [34]:

[TABLE]

Here, $w_{i}>0$ ($i$=1, 2, 3) is the weighting factor for the $i$th objective function $F_{i}$ which determines a comparative contribution specifing the user’s preference. The lower-bound transformation method [34] was used to ensure Fi are of comparable magnitudes in numerical calculations. Thereby, the transformation coefficient $\xi_{i}$ ($i$=1, 2, 3) was given by:

[TABLE]

where $F_{i}^{0}$ = minimum $\{F_{i}(d_{1},...,d_{p})\}$, and $d_{1-p}$ denote the $p$ degrees of freedom, see section 2.2. As highlighted by equation (2.1), it is important to notice a subtle distinction between using weights to determine the user’s preferences and using weights to transform objective functions. The $w_{i}$ are arbitrary values systematically varied in the interval (0,1) to set the relative importance of the objectives; see, e.g., the specified values of $w_{i}$ and $\xi_{i}$ presented in appendix A.

In our calculations, we used normalized (non-dimensional) objective functions $F_{i}$ in equation (1). A non-dimensional pseudopotential can be written as $\widetilde{\mathit{\Phi}}_{\rm ps}={{\mathit{\Phi}}_{\rm ps}}({\frac{\ q^{2}V_{\rm RF}^{2}}{4m\Omega_{\rm RF}^{2}h^{2}}})^{-1}={\|\nabla\mathit{\Theta}_{\rm RF}\|}^{2}$, where $\mathit{\Theta}_{\rm RF}$ is the RF-electrodes basis function [35] which depends solely on the geometry. Here, $m$ and $q$ are the mass and charge of the trapped ions, $h$ is the trapping height, and $V_{\rm RF}$ and $\Omega_{\rm RF}$ denote the RF drive amplitude and angular frequency, respectively. Building on this formalism, the objective functions and, hence, solutions of the optimisation are independent of the mass of the ions, the trapping parameters ($V_{\rm RF}$ and $\Omega_{\rm RF}$) as well as the trapping height. Note that owing to the normalisation of $F_{3}$ in equation (1c) the trapping height, which is a characteristic length of the system, does not contribute to our optimisation calculations, and can be substituted in the final result. As a result, the geometries obtained (see appendix A) present universal SE structures applicable to any mass or set of trapping parameters.

2.2 Design space

In this study, the design space refers to a set of parameters defining RF electrode geometries to be used in optimisation processes. The number and locations of control points specifying RF electrode shapes were varied according to the three methods shown with the example of a cross junction in figure 1. The first approach utilized four control points, where two points on each perimeter are connected simply via a line (figure 1(a)). The choice of only four parameters accelerated optimisation processes and led to satisfactory solutions despite the low number of degrees of freedom.

In the second approach, the inner and outer RF electrode edges were broken up into several line elements (figure 1(b)). We began with four control points and increased them up to 20 points while aiming for improved results. Careful attention must be paid when increasing the number of degrees of freedom because it does not necessarily result in better solutions due to the higher probability of getting trapped in local minima of the optimisation function. This method of parametrisation has been widely used [23, 28, 24, 29, 25, 32], however, electrode structures obtained often exhibited sharp features arising from underdetermined geometrical constraints. Fourier components of the surface potential of wavelength $\lambda$ exponentially decrease as exp($-\frac{2\pi z}{\lambda}$) [36], and consequently, the field is not significantly influenced by geometric elements of size $\lambda\ll z$. Such sharp features [28, 24] might also hinder precise manufacturing and would reduce breakdown voltages through the bulk or surface flashover.

Thus, we employed a third approach to mitigate these issues by using cubic splines to define the design. Control points on each electrode perimeter were bound to a cubic spline specified by eight parameters, see figure 1(c)-(d). Because of symmetry, the number of parameters for each spline was eventually reduced to six, and thus the geometry of the cross junction was defined by 12 parameters. This method of parametrisation was found to be flexible as well as efficient and was mainly used in our calculations.

2.3 Implementation

In our calculations, a two-step optimisation process was applied. In the first optimisation, the potential generated by SE electrodes was analytically calculated in the “gapless plane approximation” [37] using the SurfacePattern software package [38]. These fast calculations enabled many iterations, and hence the exploration of some of the generic aspects of junction geometries. Due to weakly confining potentials at ion channel intersections, the effect of gaps cannot be neglected, and therefore, re-optimisation of the structure in the presence of the gaps is essential. The geometries obtained from the first optimisation served as starting points for the second, in which electric potentials generated by three-dimensional structures including gaps were computed using the finite element method (FEM) [39]. In a multi-objective optimisation process, a single set of parameters that simultaneously minimises all objective functions usually does not exist, hence, a Pareto optimal solution is searched among many possible solutions [34].

3 Optimisation results and analyses

3.1 Optimisation of the cross junction

First, it is worth discussing the geometries of intersecting RF rails for infinite quadrupolar channels. In reference [40], the optimal ratio of the width of the RF electrodes $b$ to their separation $a$ in a five-electrode design (DC-RF-DC-RF-DC) with RF electrodes of equal width (see figure 1) was calculated for a given $a$ while maximising the trapping depth in the absence of control static fields (referred to as the intrinsic depth, that is, the value of the pseudopotential at escape points). Because of the importance of the trapping height $h$ as a constraint in many experiments, the intrinsic depth as well as the electric potential curvature was optimised for a fixed $h$ in references [4, 36, 41]. Table 1 presents a detailed comparison between these three calculations, highlighting the fact that the optimal value of $\frac{b}{a}$ depends on the optimisation goal.

We fixed the ratio $\frac{b}{a}$ to an optimised value for a specific pseudopotential curvature in the transverse directions, and used the geometric degrees of freedom in the vicinity of the junction to suppress barriers. We chose the value $\frac{b}{a}=2.414$, which maximises the curvature of the RF field in the transverse plane and additionally leads to the trapping depth of 95% of the maximum achievable one at a given trapping height (table 1).

We chose the trapping height to be $\approx$ 700 $\mu$m, and hence, the width of the RF electrodes $b=$ 1400 $\mu$m and their separation $a=$ 580 $\mu$m. Thus, a sufficiently large trapping volume for a few hundred ions as well as stable trapping conditions for a variety of interesting ion species can be achieved (see below and figure 2). Moreover, a suitable choice of the trap size is crucial to eliminate effects of anharmonic terms in the trapping potential which increase with decreasing trap dimensions. These effects are pronounced in SE traps owing to their asymmetric geometries, and their consequences on bicomponent Coulomb crystals, e.g., the segregation of atomic from molecular ions, have been discussed in [33].

The minimisation of pseudopotential barriers is assessed based on a suppression factor which is defined as the ratio between the height of the suppressed barriers and either those of a straight junction or the intrinsic depth of a straight channel far away from the junction. The result of the first optimisation of the cross junction is presented in figure 2. The normalised pseudopotential is evaluated for the RF amplitude VRF=300 V, and the RF angular frequency $\Omega_{\rm RF}$=$2\pi\times 10.0$ MHz. In this case, the highest barrier for 40Ca*+* is 0.36 meV (figure 2(a)), 614 times less than the depth of a straight channel, and 125 times less than the barriers present in a design with completely straight electrodes.

The non-confining characteristic of the trapping potential at the intersection region manifests itself in a drastic change of the stability parameters as illustrated in figure 2(b). The $Q_{ii}$ ($i\in\{x,y,z\}$) in figure 2(b) are the diagonal elements of the multi-dimensional form of the stability parameter $q_{\rm M}$ of the Mathieu equation [42], which are defined by

[TABLE]

evaluated along the ion channel ${\mathbf{r}_{\rm min}(x)=\mathop{\mathrm{arg\,min}}\limits_{y,z}({\Phi}_{\rm ps}(x,y,z))}$, where $\mathop{\mathrm{arg\,min}}$ represents the values of y and z that minimise ${\Phi}_{\rm ps}(x,y,z)$ at a certain x, $\mathit{\Theta}_{\rm RF}$ is the RF-electrodes basis function [35], and $m$ and $q$ are the mass and charge of the trapped ions, respectively.

3.2 Second optimisation of the cross junction and characterisation of a central bridge

The three-dimensional structure of the cross junction with the inclusion of 120 $\mu$m wide gaps was modelled using FEM (figure 3(a) (i-ii)). Due to the 330 $\mu$m diameter hole at the centre, the pseudopotential exhibited a spatially broad, loosely confining trapping region along the vertical direction $z$ (figure 3(a) (ii)).

To address this issue, we investigated the incorporation of a bridge between the central DC electrodes with a carefully refined geometry. To explore the consequences of this modification, over 100 geometries were simulated using FEM. Bridges formed by merging the tips of RF electrodes at the centre cause extremely high pseudopotential barries. This situation is opposed to junctions in three-dimensional electrode geomeries, e.g., the one demonstrated in [22], where RF bridges were located symmetrically above and below the ion channel. Our calculations show that as the bridge widens, pseudopotential barriers drastically increase. Further calculations were performed using a 70-${\mu}$m wide DC bridge in order to enhance mechanical stability of the structure [43]. For a fixed DC bridge width, the height of pseudopotential barriers reduces as the size of the gaps in the vicinity of the junction centre decreases, and thus a further decrease of the barriers depends on how narrow the gaps can be fabricated. These calculations enabled the suppression of pseudopotential barriers by a factor of three compared to an initial straight bridge of the same width while the variation of the trapping height was controlled in a range of $\pm$10%. The resulting geometry depicted in figure 3(b) (i) features a confining trapping potential perpendicular to the ion channel at the intersection (figure 3(b) (ii)).

Figure 4 illustrates the influence of the incorporation of the optimised DC bridge on the characteristics of the cross junction. A significant increase of the trap stiffness at the centre has been achieved (figure 4(a)), which also manifests itself in the larger $Q_{zz}$ parameter improved by a factor of two (figure 4(c)). This junction exhibits pseudopotential barriers $<$ 5.5 meV for 40Ca*+* (figure 4(b)), which is less than 2.5% of the trap depth far from the junction in the quadrupolar channel, whilst the ion channel is elevated at a well-controlled height over the ion shuttling path (figure 4(d)).

3.3 Quadrupole-to-octupole junction

To design a quadrupole-to-octupole junction that smoothly converts the trapping fields, similar optimisation processes were implemented. The geometry was parametrised using three cubic splines (section 2.2) corresponding to 14 degrees of freedom, see appendix A. Two different solutions were obtained: patterns with gradually enlarging RF electrodes and non-trivial electrode geometries with rather complex patterns, which are labelled with J1 and J2 , respectively, in figure 5(a). The corresponding pseudopotential profiles of the two junctions are presented in figure 5(b), calculated for 40Ca*+* at the RF amplitude VRF=300 V and the RF angular frequency $\Omega_{\rm RF}$=$2\pi\times 10.0$ MHz, where the asymptotic quadrupolar channel minimum is at 700 $\mu$m above the electrode plane.

J1 was obtained from optimisation processes in which pseudopotential barriers were treated with large weighting factors in equation (2.1). Such an electrode pattern can also be intuitively envisioned for a quadrupole-to-octupole junction as it slowly interpolates the two asymptotic channel electrode patterns. In contrast, this transformation occurs on a shorter distance in J2, which was obtained from an optimisation in which the length of the junction was constrained. This can be seen in figure 6 showing one-dimensional cuts through the pseudopotential along the vertical direction $z$ at six different locations along the ion transport path for J1 and J2.

The pseudopotential barriers calculated along the transport channel for the two structures are compared for N${}_{2}^{+}$, Ca*+, and CaO+* in figure 7. Although the pseudopotential barriers in J2 are larger, the shuttling distance to the octupolar region is shorter making it the preferable geometry.

The electric fields generated by J2 with the inclusion of 120 $\mu$m-wide gaps were modelled using FEM where only slight refinements were applied in order to avoid sharp electrode shapes. The result ensured a smooth transition between the two trapping regions, and thus no further optimisation was undertaken. Note that in the field-modifying linear junctions the pseudopotential barriers are typically much lower than those of multichannel intersections.

4 Design of a multi-functional ion-trap chip for integrated experiments with cold molecular ions

Building on this study, we have developed a design for a multi-functional SE RF ion-trap chip to enable the performance of various tasks such as loading and preparation of ions, mass spectrometry, spectroscopy, and manipulation of molecular Coulomb crystals on a chip device. Figure 8 illustrates the layout of this chip. The chip consists of four quadrupolar ion channels and one octupolar ion channel connected via the two optimised SE junctions discussed in the previous sections. The optimisation of the SE octupole trap integrated on the chip is discussed in appendix B.

The proposed trapping architecture features the functionality of a large guided ion-beam apparatus, e.g., see references [44, 45], at a fraction of the production and operation costs. The quadrupolar zone Q1 can be used to prepare cold molecular ions and to characterise them using, e.g., a resonant-excitation mass spectrometry technique [46]. The experimental sites labelled with Q2 and Q3 are designed for spectroscopy and ion-neutral reaction studies with suitable accesses for laser and molecular beams. The octupolar region along the fourth arm can be exploited for reaction studies in a large RF-field-free trapping volume as well as for storing ions. The monolithic, miniaturised design also facilitates the operation of cryogenic systems. The concrete experimental implementation of this device will be reported elsewhere.

5 Conclusions

We have discussed the design and optimisation of surface-electrode ion-trap junctions of two types, an X-junction between four quadrupoles and a quadrupole-to-octupole junction, in view of developing integrated, monolithic trapping architectures for cold molecular ions. A detailed description of the multi-objective optimisation exploited to design arbitrary SE junctions has been presented. The results are modular, scalable components that can be used as design libraries for developing a complex ion-trap network. These optimal geometries are independent of the ion mass and trapping parameters. We have shown that for a cross junction between four quadrupolar channels, the incorporation of an optimised bridge results in a significant increase of the trap stiffness as well as a better control over the trapping height at the intersection. Our optimisation yields a non-trivial geometry for the quadrupole-to-octupole junction which features an effectively reduced distance between the two trapping regions.

Based on these results, a multi-functional surface-electrode ion-trap chip has been proposed that combines various tasks such as loading and preparation of ions, mass spectrometry, spectroscopy, and manipulation of the structure of crystals in a miniaturised device using precisely shaped trapping potentials. The carefully designed ion channel intersections are intended for the transport of sympathetically cooled molecular ions in the form of bicomponent crystals.

A. M. thanks Dr. Bjoern Lekitsch from the University of Sussex for helpful discussions about the microfabrication techniques. This work has been supported by the University of Basel, the COST Action MP1001 “Ion Traps for Tomorrow’s Applications”, the Swiss National Science Foundation through the National Centre of Competence in Research “Quantum Science and Technology” and the Swiss Nanoscience Institute.

Appendix A: Specifications of the optimised junctions

The electrode geometries of the optimised cross junction (figure 2(a)) and quadrupole-to-octupole junction J2 (figure 5(a)) obtained are fully defined by the cubic splines specified in tables 2 and 3, respectively. For clarity, each spline is labelled in figure 9, see also figure 1(d) for the definition of the parameters used. Note that the dimensions were normalised, e.g., lengths are in units of the trapping height, and therefore, the values can be accordingly scaled to suit any required size of junction architectures.

For the optimisation of the cross junction, the transformation coefficients in equation (2.1) were calculated to be $\{\xi_{1}$, $\xi_{2}$, $\xi_{3}\}$ = $\{$1.6$\times 10^{-6}$, 6.1$\times 10^{-12}$, 1.0$\times 10^{-4}\}$ associated with pseudopotential barriers and gradients, and the trapping height variation, respectively. The corresponding weighting factors were altered in the interval $(0,1)$ such that the three objective functions were significantly represented without deterministic priority, e.g., in this case $\{w_{1}$, $w_{2}$, $w_{3}\}$ = $\{0.2$, $0.6$, $0.2\}$. For the quadrupole-to-octupole junction, these values were given by $\{\xi_{1}$, $\xi_{2}$, $\xi_{3}\}$ = $\{$2.2$\times 10^{-5}$, 9.9$\times 10^{-11}$, 8.9$\times 10^{-4}\}$ and $\{w_{1}$, $w_{2}$, $w_{3}\}$ = $\{0.375$, $0.5$, $0.125\}$.

Appendix B: Optimisation of a SE octupole channel

Three approaches can be considered to find the optimal geometry of a SE multipole trap within the ”gapless plane approximation” [37]. In the first approach, the OptimalFinitePattern routine from the SurfacePattern package enables the computation of an optimal surface geometry for a set of given local constraints, i.e., fixed values of the potential and derivatives at a certain point based on the relaxed Linear Programming method [47]. In the second approach, the optimum position of the electrode edges is analytically calculated. We maximised the strength of the RF-electrodes basis function $\mathit{\Theta}_{\rm RF}$ at a given height $h$ while the lower order multipole coefficients were set to zero:

[TABLE]

where $\{i,j,k\}\in\{\mathbb{N}~{}|~{}\exists n:(i+j+k)<n\}$ and $n$ denotes the order of the multipole and $p$ is the number of degrees of freedom. For an optimal octupole ($n=4$) with symmetric design of four RF electrodes ($p=4$), the solution of equation (1e) is given by $\{y_{1},y_{2},y_{3},y_{4}\}$ = $h\times\{-1-\sqrt{2}+\sqrt{2\left(2+\sqrt{2}\right)},1-\sqrt{2}+\sqrt{4-2\sqrt{2}},-1+\sqrt{2}+\sqrt{4-2\sqrt{2}},1+\sqrt{2}+\sqrt{2\left(2+\sqrt{2}\right)}\}$. The third approach is based on conformal mapping of conventional cylindrical trap geometries to strip electrodes on a plane [36]. Figure 10(a) depicts the configuration of an optimal octupole trap featuring maximum absolute value of $\partial^{4}\mathit{\Theta}_{\rm RF}/\partial z^{4}$. All three approaches generate the same result [43].

The geometry obtained was analysed using FEM with the inclusion of the realistic inter-electrode spacings. The result showed a non-vanishing RF field at the trapping centre which yield three pseudopotential minima (figure 10(c)). The ramification would be disordered structures of Coulomb crystals with ions driven to unwanted regions, those in which atomic ions are not efficiently cooled by lasers. Attempts to compensate this effect by the application of static fields would be inadequate because it is solely caused by incorporating gaps to the geometry.

To address this issue, electric potentials obtained from FEM calculations were used to generate numerical objective functions for a multi-objective optimisation using the Nelder-Mead algorithm. The merit function was given by:

[TABLE]

where the objective functions $f^{(n)}$ were formulated as $f^{(n)}=\frac{\partial^{n}\mathit{\Theta}_{\rm RF}(d_{1},d_{2},d_{3},d_{4})}{{\partial z}^{n}}$. The parameters $d_{1},d_{2},d_{3}$, and $d_{4}$ are schematically shown in figure 10 (b). $\lambda_{i}$ stands for the weighting factor including a transformation coefficient, see section 2.1. The result of the second optimisation is presented in figure 10 (d). As the gap aspect ratio $\frac{t}{g}$ decreases particularly for $\frac{t}{g}<2$, the residual RF field is enhanced [43], where $t$ and $g$ are the inter-electrode spacing and the trench depth, respectively, see figure 10(b). The explanation is provided by understanding the role of the gap polarization and interpolation potentials that alter the surface potential [37]. The effect is pronounced in the octupolar channel because of the sensitivity of the field-free region to the variation of the surface potential.