---
title: "Photonic crystal surface modes for trapping and waveguiding of ultracold atoms in quantum sensors"
authors: ["Valery Konopsky 1"]
affiliation: "% Institute of Spectroscopy, Fizicheskaya, 5, Troitsk, Moscow, 108840, Russia; konopsky@isan.troitsk.ru"
journal: "Sensors"
year: 2023
volume: ""
issue: ""
article_number: ""
pages: ""
doi: "10.3390/s23218812"
type: journal-article
site_group: ""
url_abstract: ""
url_pdf: "https://valery.konopsky.com/kvnlocal/Konopsky_Sensors-23-08812.pdf"
language: en
source_tex: "Z:\\ValeryData\\Valery_New\\my_articles\\Sensors2023\\Sensors\\manuscript\\Konopsky2Sensors2023.tex"
source_pdf: "Z:\\ValeryData\\Valery_New\\my_articles\\Sensors2023\\Sensors\\published\\Konopsky_Sensors-23-08812.pdf"
---
## Abstract

The design of a photonic system for trapping and waveguiding of ultracold atoms far above a dielectric surface is proposed and analysed. The system consists of an optical rib waveguide deposited on a planar one-dimensional photonic crystal, which sustains two wavelengths of photonic crystal surface modes tuned in the red and blue sides relative to the atomic transition of the neutral atom. The addition of a third blue-tuned wavelength to the system allows the neutral atoms to be stabilised in the lateral dimension above the rib waveguide. Trapping atoms at relatively large distances above the dielectric surface allows to reduce the influence of Casimir-Polder forces in this system. The detailed design methodology and specifications of the photonic system are provided.

# Introduction

The development of quantum sensors [1] based on coherent atomic circuits is important for the realisation of many potential technological applications [2], such as gravimeters [3, 4, 5, 6, 7], gradiometers [8] and gyroscopes [9, 10]. Atom chips [11, 12, 13] provide a versatile technique for the generation and coherent manipulation of ultracold atoms for such sensors on the micrometer scale.

For the successful development of compact and robust atom chips integrated with guided atomic interferometers, the optimal choice of the appropriate interface for the light-atom interaction is of crucial importance. In this study, we examine the trapping and waveguiding of ultracold atoms employing two-colour optical surface waves propagating along the external boundary of a rib waveguide located on top of a specially designed planar one-dimensional photonic crystal (1D PC).

The use of optical waveguides to guide atoms was proposed in 1993 by Letokhov *et al.* [14]. They proposed to create a potential well inside a hollow core fiber by a laser with a red-detuned light directed into the fiber. Since then, the use of optical waveguides has been identified as a key element of atomic chip interferometry, and many types of optical waveguide have been considered as candidates for the best interfaces for light-atom interactions in them. The employment of integrated optical waveguides to trap and guide ultracold atoms above planar surfaces has been analysed in [15, 16]. They exploit an earlier proposal by Ovchinnikov *et al.* [17] to use two colors of light, with red and blue detuning, and with different evanescent decay lengths, to obtain a planar (one-dimensional) trap with a potential minimum above a dielectric prism. That is, the insertion of optical waveguides into planar surface should add lateral confinement and restrict the ensemble of cold atoms in two dimensions.

The main problem here is that the evanescent decay length into an external medium in such waveguides cannot be very large, since the effective refractive index of an optical wave $\mathrm{n_{eff}}$ in a standard waveguide cannot be smaller than the refractive index of the waveguide substrate $n_0$. The penetration length of the evanescent wave intensity (i.e., $I = |E|^2$) in the external medium, is
``` math
\begin{equation}
   L_e=\frac{\lambda}{4\pi \sqrt{\rho^2-n_{e}^2}} \; .

\end{equation}
```
Here $\rho$ is either the numerical aperture of the light beam in the prism, $\rho=n_0\sin(\theta_0)$, or the effective refractive index of the optical mode in the corresponding planar 1D slab waveguide, $\rho=\mathrm{n_{eff}(w_{\infty})}$. The refractive index of the external medium (i.e., vacuum, for atom optics applications) is $n_{e}=1$. The value of $L_e$ can be very large only in the case when the difference $(\rho^2-1)$ is small.

It will be shown below that this condition $\rho=\mathrm{n_{eff}(w_{\infty})}\simeq 1$ can be satisfied for a photonic crystal (PC) waveguide, in contrast to a standard optical waveguide. It is this unique property of PC waveguides that makes it possible to design the photonic system that sustains the ultracold atomic ensembles at relatively large distances above the dielectric surface.

# Materials and Methods

## Planar photonic crystal waveguides and their differences from standard waveguides

The effective index of a two-dimensional optical waveguide $\mathrm{n_{eff}(w)}$ with finite lateral width $\mathrm{w}$ (i.e., rib or ion-implanted waveguide) is somewhat smaller than the effective index of the corresponding planar 1D slab waveguide $\mathrm{n_{eff}(w_{\infty})}$, where the width of $\mathrm{w}=\infty$ (i.e., extends to infinity in the lateral direction). This diminution of $\mathrm{n_{eff}(w)}$ is due to the zigzag propagation (from left to right) of the wave in the plane of the waveguide (in the ray model approximation), resulting in a reduction in the projection of the wavevector on the z-direction. But the penetration length of the waveguide mode intensity into the external medium for a two-dimensional waveguide is still determined by the formula (1), where the effective index of the corresponding planar 1D slab waveguide $\mathrm{n_{eff}(w_{\infty})}$ is inserted as $\rho$.

According to equation (1) only modes with $\mathrm{n_{eff}(w_{\infty})}\simeq 1$ can have large penetration length in the medium with $n_{e}=1$. In standard waveguides, the optical field confinement in the vertical direction is the result of total internal reflection (TIR) on both sides of the waveguide layer. As a result, $\mathrm{n_{eff}(w_{\infty})}>n_0 \geqslant n_e$ (i.e., $\mathrm{n_{eff}(w_{\infty})} > 1$) and the maximum possible penetration length cannot be greater than $L_e(n_0)=\lambda/(4\pi [n_0^2-n_e^2]^{1/2})$. For a waveguide on a quartz substrate with $n_0=1.46$ and $\lambda=850$ nm, this gives $L_e(1.46) \leqslant 64$ nm. In [15], the authors consider sodium fluoride (the lowest-index common optical mineral, $n_0=1.32$) as a possible substrate, but it yields only $L_e(1.32)\leqslant 79$ nm.

<figure id="fig1" data-latex-placement="h!">
<span class="image placeholder" data-original-image-src="Figures/Fig1.eps" data-original-image-title="" width="75%"></span>
<figcaption><span id="fig1" data-label="fig1"></span> Concept of the photonic system under consideration. The left inset shows the cross-section of the 1D PC structure, and the right inset shows the outline of the trap-well potential for ultracold atoms.</figcaption>
</figure>

One possible way to mitigate this limitation is to use a freely suspended thin dielectric film [18, 19] to achieve $n_0=n_e=1$. With this approach, the following penetration length values for the *intensity* of evanescent waves were obtained [18]: $L_e(1.0) = 87.5$ nm and 69.75 nm for modes excited in a 300 nm-thick suspended silica film with waveguide rib height $\mathrm{h=15}$ nm and width $\mathrm{w=2}\,\upmu$m (at wavelengths $\lambda=850$ nm and 720 nm, respectively).

This article presents a different approach to overcome this limitation, in which the uniform waveguide substrate is replaced by a planar 1D photonic crystal, which is designed to create a photonic bandgap in the spectral region of interest. Thus, in this case, the optical field confinement in the vertical direction is due to reflection from the photonic bandgap on one side (bottom) and total internal reflection – as usual – on the another side (top) of the waveguide as shown in Fig. 1.

The planar 1D PC is a simple dielectric stack. It is durable and compatible with existing technologies for producing dielectric mirrors, so its fabrication is much easier than suspended thin dielectric films. In the vicinity of the interface between the planar 1D PC and the external medium, optical surface waves can be excited subject to an appropriate choice of the thicknesses of the double layers and the thickness of the final truncated layer [20].

These optical surface waves are excitations of photonic crystal surface modes (PC SMs), also known as ‘(asymmetric) planar Bragg waveguide modes’, ‘surface waves in periodic layered media’, ‘photonic bandgap surface modes’, ‘optical Bloch surface waves’, ‘photonic crystal surface waves’ and ‘surface waves in multilayer coatings’. They were first studied both theoretically [21, 22] and experimentally [23] in the 1970s. Twenty years later, the excitation of optical surface modes in the Kretschmann-like configuration was demonstrated [24, 25]. In recent years, PC SMs have found an increasing number of applications in the fields of optical sensors [26, 27, 28, 29, 30], optical biosensors [31, 32, 33, 34, 35, 36]. and in other fields [37, 38, 39, 40, 41, 42].

The penetration length in the external medium can be very large for these modes as their effective refractive index can be infinitesimally close to the refractive index (RI) of the external medium. This is a unique property of PC SMs that makes their excitation with the $\rho=\mathrm{n_{eff}(w_{\infty})}\simeq 1$ quite possible. Already in the first optical biosensor based on these modes [31], this property was used to distinguish between a surface adsorption and bulk RI (i.e., RI far from the surface).

<figure id="fig2" data-latex-placement="h!">
<span class="image placeholder" data-original-image-src="Figures/Fig2.eps" data-original-image-title="" width="100%"></span>
<figcaption>Dispersion of the planar 1D PC structure (2), for <span class="math inline"><em>N</em> = 5</span>, with <span class="math inline"><em>L</em><sup>′</sup></span> thickness <span class="math inline"><em>d</em><sub>3</sub> = 207</span> nm (A), and with <span class="math inline"><em>d</em><sub>3</sub> = 182</span> nm (B). Intensity profiles of modes at <span class="math inline"><em>d</em><sub>3</sub> = 207</span> nm for <span class="math inline"><em>λ</em><sub>1</sub> = 850</span> nm, <span class="math inline"><em>ρ</em> = n<sub>eff</sub>(w<sub>∞</sub>) = 1.0056</span> (C), and <span class="math inline"><em>λ</em><sub>2</sub> = 640</span> nm, <span class="math inline"><em>ρ</em> = n<sub>eff</sub>(w<sub>∞</sub>) = 1.0105</span> (D). </figcaption>
</figure>

## Design methodology for a 2D waveguide placed on a 1D PC

### Design of a 1D PC 

The design of a 1D PC, sustaining two PC SMs at wavelengths $\lambda_1=850$ nm and $\lambda_2=640$ nm, was performed using a free Windows program available at [43]. Results are presented in Fig. 2. The wavelengths are chosen to be detuned far enough away from the $D_1$, $D_2$ transitions of the $^{87}$Rb atom ($\lambda_{D_1}=795$ nm, $\lambda_{D_2}=780$ nm) on the blue and on the red sides. The light potential (which causes the attraction and repulsion forces) scales as $I/\Delta$, while the light scattering rate (which causes the heating rate) scales as $I/\Delta^2$ [44], where $\Delta$ is a detuning from the $D$ transitions (see Eqs. (3), (4) below). Therefore, to increase the trap lifetime, it is better to be away from the resonance if the intensities still allow the desired trap depth to be reached. The larger detuning in the blue side is due to the expectation that a larger value of intensity of the repulsive (blue-detuned) light beam will be needed to form a well trap potential of the desired shape.

The following 1D PC structure is designed to sustain both $\lambda_1=850$ nm and $\lambda_2=640$ nm PC SMs:
``` math
\begin{equation}
\mathrm{prism / (LH)^{N}L''HL' / vacuum} \,,

\end{equation}
```
where $L$ is the $SiO_2$ layer (thickness $d_1=170. 73$ nm), $H$ is the $TiO_2$ layer (thickness $d_2=86. 46$ nm), $L''$ is the extended $SiO_2$ layer (thickness $d_{12}=304$ nm), and $L'$ is the final $SiO_2$ layer (thickness $d_3=207$ nm). The RIs of $SiO_2$, $TiO_2$ layers used in the calculations are given in Table 1. The obtained parameters of PC SMs in this planar 1D PC structure are also presented in Table 1 and in Fig.2.

\|C\|CC\|CC\| **wavelength** & $n_1 (SiO_2) \textsuperscript{\textborn}$ & $n_2 (TiO_2) \textsuperscript{\textborn}$ &$\rho=\mathrm{n_{eff}(w_{\infty})}$ & $L_e$, nm\
**850 nm** & 1.4666& 2.3137 &1.0056 & 638\
**640 nm** & 1.4679& 2.3674 &1.0105 & 351\

<sup> </sup>+ $i\cdot10^{-5}$

It can be seen that this 1D PC is designed so that both modes at $d_3=207$ nm have an excitation angle very close to the TIR angle, and hence the effective RI is close to unity ($\rho=\mathrm{n_{eff}(w_{\infty})}\simeq1$). As a consequence, the penetration lengths of their intensities into the vacuum are very long (see $L_e$ in Fig. 2(c,d) and in Table 1).

### From 1D PC waveguide to 2D waveguide

Note that if we reduce the thickness of the final $SiO_2$ layer by 25 nm (or more), from $d_3=207$ nm to 182 nm, then no PC SMs will be excited on such planar 1D PC near $\rho\simeq 1$ (up to $\rho=1.17$ for $\lambda_2=640$ nm), see Fig.2(b). Therefore, if we will use the planar 1D PC with $d_3=182$ nm and a waveguide of height $\mathrm{h=25}$ nm above it, the excitation of the PC SMs will take place inside the two-dimensional waveguide only, subject that its width $\mathrm{w}$ will be sufficient. Thus, if excitation is performed by beams with excitation angles close to $\theta_0$, then only PC SMs with $\rho=\mathrm{n_{eff}(w_{\infty})}=n_0\sin(\theta_0)$ will be excited, i.e., only modes in two-dimensional waveguides of width $\mathrm{w}$, but not modes in the surrounding planar slab waveguide. Moreover, the penetration of such waveguide wave into the external medium (vacuum, in our case) will be the same as for the wave in the planar slab waveguide in Fig.2(a, c, d), where the waveguide width $\mathrm{w_\infty}$ is extended to infinity in the lateral direction. The equation (1) still gives the correct value of $L_e$ when $\rho=\mathrm{n_{eff}(w_{\infty})}$ is substituted. This is because the same wave propagates in the 2D waveguide as in the planar slab 1D waveguide, but in a zigzag fashion from left to right, if considered in the ray model approximation.

### Design of a 2D waveguide placed on a 1D PC 

The mode analysis of the 2D waveguide on the surface of the 1D PC was performed by using both a semi-analytical mode solver based on the Variational Effective Index Method (VEIM) [45, 46] and a numerical Finite Difference Eigenmode (FDE) solver [47] with the Lumerical Inc. MODE software. Both approaches give very similar results for effective refractive indices and mode profiles.

## Calculating the optical dipole potential at large wavelength detuning

For very large detunings that significantly exceed the splitting of the D-line doublet, the following approximate expression for the dipole potential can be used [44]:
``` math
\begin{equation}
   U_\mathrm{dip}(x,y,\omega)=\frac{3\pi c^2}{2 \omega_D^3} \frac{\gamma_D}{(\omega-\omega_D)}I(x,y) \; ,

\end{equation}
```
where $\omega_D= 2\pi \times 380.77$ THz is the frequency of the center of the D-line doublet, and $\gamma_D= 37$ MHz is the mean value of natural linewidths of the transitions $D_1$, $D_2$. Using the two-dimensional spatial distributions of light intensities $I(x,y)$ shown in Fig.3, we obtain the trap well potential shown in Fig.4 (in temperature units $U/k_B$). Atomic light scattering rate also may be calculated in this approximation:
``` math
\begin{equation}
   \hslash\, \Gamma_\mathrm{sc}(x,y,\omega)= \frac{\gamma_D}{(\omega-\omega_D)} U_\mathrm{dip}(x,y,\omega)\; .

\end{equation}
```
Each scattered photon imparts a recoil energy $E_\mathrm{rec}=( \hslash\omega/c)^2/(2m)$ to the atom of mass $m$. As a result, taking into account the approximated expressions (3) and (4), a trap lifetime $\tau$ can be estimated as:
``` math
\begin{equation}
   \frac{1}{\tau}\simeq \frac{ 2 E_\mathrm{rec} }{U_\mathrm{dip}} \Gamma_\mathrm{sc}  = \frac{2 E_\mathrm{rec} }{\hslash |\omega-\omega_D|}\gamma_D\; .

\end{equation}
```
In addition to atom losses due to photon scattering, the trap lifetime can be reduced by background gas collision losses, which outweigh photon scattering losses at large detuning values [48].

The total potential of the atom $U_\mathrm{total}(x,y)$ near the surface is the sum of the optical potentials $U_\mathrm{dip}(x,y,\omega)$ of all participating light beams and the Casimir-Polder potential together with the gravitational potential:
``` math
\begin{equation}
   U_\mathrm{total}(x,y)= \sum^{\omega_3}_{\omega=\omega_1}{U_\mathrm{dip}(x,y,\omega)} - \frac{C_3 \lambda_\mathrm{eff}}{x^3 (x + \lambda_\mathrm{eff})}   -  mgx \; ,

\end{equation}
```
where $C_3 = 5.7\times 10^{-49}$ J m$^3$ is the van der Waals coefficient, $\lambda_\mathrm{eff} = 710/(2\pi)$ nm is the reduced wavelength [49, 18] describing the attractive interaction between atoms and the glass surface, and $g=9.8$ m/s$^2$ is the acceleration of free fall.

<figure id="fig3" data-latex-placement="h!">
<span class="image placeholder" data-original-image-src="Figures/Fig3.eps" data-original-image-title="" width="75%"></span>
<figcaption> Electric field intensity distribution for red-detuned attracting light beam, <span class="math inline"><em>λ</em><sub>1</sub> = 850</span> nm, <span class="math inline">n<sub>eff<sub>1</sub></sub> = 1.0012</span> (A), blue-detuned repulsive light beam, <span class="math inline"><em>λ</em><sub>2</sub> = 640</span> nm, <span class="math inline">n<sub>eff<sub>2</sub></sub> = 1.0074</span> (B), and auxiliary blue-detuned repulsive light beam for lateral support, <span class="math inline"><em>λ</em><sub>3</sub> = 633</span> nm, <span class="math inline">n<sub>eff<sub>3</sub></sub> = 1.0076</span> (C). </figcaption>
</figure>

# Results

## Example of specification of a 2D waveguide placed on a 1D PC

After selecting a waveguide height in 25 nm (above the last layer $L'$ with a thickness of $d_3=182$ nm, see Fig.2), the width of the waveguide has to be selected. Both the semi-analytical VEIM solver and the numerical FDE solver show that the waveguide with $\mathrm{w=3\,\upmu{m}}$ sustains single TE00 mode propagation for $\lambda_1=850$ nm and $\lambda_2=640$ nm with effective RI $\mathrm{n_{eff_1}}=1.0012$ and $\mathrm{n_{eff_2}}=1.0074$, respectively. In addition, at wavelength $\lambda_3=633$ nm, the TE01 mode with $\mathrm{n_{eff_3}}=1.0076$ is also supported in the waveguide with $\mathrm{3\,\upmu{m}}$ width. The next [TE00] mode at $\lambda_3=633$ nm has an effective index of $\mathrm{n_{eff_{30}}}=1.0163$. Therefore, only three modes with the electric field intensity distribution shown in Fig 3 are excited at excitation angles of $\rho\leqslant 1.016$. The third mode with $\lambda_3=633$ nm is needed to provide "lateral support" for the ultracold atoms in the optical trap above the waveguide with $\mathrm{3\,\upmu{m}}$ width.

\|C\|CCCC\|CC\| **wavelength** &$\mathrm{n_{eff}(w,h)}$ & vacuum fraction, $S^{e}_z / S^\mathrm{total}_z$ & mode effective area, $\upmu\mathrm{m}^2$ & loss, dB/cm & optical power in Fig.4, mW & intensity in Fig.4, GW/m$^2$\
**850 nm** &1.0012<sup>†</sup> & 76% & 6.0 & 2 & 8.4 & 1.4\
**640 nm** &1.0074<sup>†</sup> & 46% & 3.5 & 7 & 29.1 & 8.3\
**633 nm** &1.0076 & 39% & 3.7 & 8 & 7.4 & 2.0\

<sup> †</sup>TE00; <sup> </sup>TE01; <sup> </sup>for $Im(n_j)=10^{-5}, \forall j$

The calculated mode parameters for the rib waveguide with $\mathrm{w=3\,\upmu{m}, h=25\,nm}$ are presented in Table 2. Their effective RI, the fraction of Poynting energy transmitted through the vacuum, the effective area of the modes and the losses are given. The optical power in these modes required to form a 100 $\upmu\mathrm{K}$ trap, as in Fig.4(a,b,c), is also indicated. Note that in these waveguides a significant fraction of the Poynting energy $S_z$ propagates through the external medium (vacuum). This is a unique property of the presented 1D PC waveguides that is unattainable for standard waveguides, where most of the Poynting energy is transmitted through the waveguide core and only a small fraction of the evanescent wave is available to manipulate atoms in the vacuum.

The mode losses were calculated assuming that the imaginary part of the RI of $SiO_2$ and $TiO_2$ is $10^{-5}$. According to [50], values as low as $Im(n_j)<10^{-5}$ in this spectral range can be obtained for these dielectric materials using ion-assisted electron beam deposition. We set $Im(n_j)=10^{-5}$ to account for scattering losses, which can be quite prominent for PC SMs with effective refractive indices close to unity.

## Optical trap above the waveguide located on the 1D PC

The total trap well potentials of the atom above the rib glass waveguide on 1D PC, which are calculated using the equation (6), are shown in Fig.4. Figure 4(a) shows a potential with two wavelengths: attracting (red-detuned, $\lambda_1=850$ nm) and repulsive (blue-detuned, $\lambda_2=640$ nm). One can see in this figure a common problem for this type of potential: weak localisation in the lateral dimensions and a decrease of the height of the trap well near the edge of the waveguide (i.e., near $y=\pm 1.5\,\upmu$m in our case). In the vicinity of these points, atoms will escape from the trap. To solve this problem, the third beam with the wavelength of $\lambda_3=633$ nm [TE01] is added and the resulting total potential is shown in Fig.4(b). That is, all the modes shown in Fig.3 are used in Fig.4(b). It is clear that such "lateral support" dramatically improves lateral localisation and eliminates atom escape paths to the surface through points $y=\pm 1.5\,\upmu$m.

<figure id="fig4" data-latex-placement="h!">
<span class="image placeholder" data-original-image-src="Figures/Fig4.eps" data-original-image-title="" width="100%"></span>
<figcaption> Optical potential well when two wavelengths of 850 nm and 640 nm are involved (A), and when the third auxiliary wavelength of 633 nm is added (B). Vertical cross-section of the potential well at y=0 (C), and horizontal cross-sections of the potential well at x=665 nm for different intensities of the 633 nm-wave (D).</figcaption>
</figure>

In Fig.4(c) a vertical cross section of the potential well at $y=0$ is shown. The minimum of the well is at $h=665$ nm from the surface. This should provide a very stable confinement of atoms at large distances from the dielectric surface, where the influence of the van der Waals force is negligible. The light intensities (see Table 2) are chosen so that the trap depth is 100 $\upmu$K. The sum of all light intensities gives the total atomic light scattering rate at the centre of the trap, $\Gamma_\mathrm{sc}=8$ Hz, by the equation (4).

The extent of the lateral localisation can be easily controlled by varying the intensity of the $\lambda_3=633$ nm [TE01] lightwave, as shown in Fig.4(d), where cross sections of the potential well at $x=665$ nm is shown. The dotted cyan line corresponds to the absence of the auxiliary lateral localisation (i.e., this is the horizontal cross-section of Fig.4(a)). The solid magenta line is the horizontal cross-section of Fig.4(b), while the other two lines show how the lateral confinement increases as the intensity of the TE01 auxiliary lightwave increases.

## Excitation of waveguide modes in a Kretschmann-like scheme

For excitation of only two/three necessary modes and for avoiding excitation of undesired additional modes (e.g. $\rho=1.17$ at $\lambda_2=640$ nm on the flat surface, see Fig.2(b), or $\rho=1.016$ at $\lambda_3=633$ nm [TE00] in the rib waveguide), the scheme of Kretschmann type is quite convenient. In this scheme, mode excitation takes place via a coupling prism on which the optimal number of double layers $N$ (see (2)) is deposited. For the presented photonic structure, this number for the optimal excitation is $N\simeq 5$. With optimal excitation of a surface wave, there is no reflection of an incident beam due to destructive interference of the reflected wave with the optical surface wave re-radiated back to the prism. Only modes with $\mathrm{n_{eff}}$ close to unity can be excited if the beam excitation angle $\theta_0$ (see Fig. 1) is chosen near the total internal angle (i.e. near $\rho\sim 1$). The prism coating can be arranged in such a way that half of the prism surface is covered by a PC structure (2) with $N\simeq 5$ (for excitation), while the other half has the same structure with $N>10$ (for propagation without re-radiation into the prism).

# Discussion

Engineering the optimal interface for the light-atom interaction is critically important for developing compact and robust quantum sensors based on guided atom interferometers. In this article, a novel waveguide structure based on planar 1D PC has been presented, which is designed to trap and waveguide neutral atoms by optical dipole force. This photonic structure makes it possible to increase penetration depths of the evanescent light waves, which provides very stable trapping of atoms at relatively large distances from the dielectric surface of the waveguide. Additional stabilisation of the atom in the lateral dimension by the auxiliary TE01 mode has been considered. By modulating the intensity of this mode, the lateral width of the trap well can be easily modulated, which can be used, e.g., for adiabatic cooling and other manipulations of ultracold atomic ensembles.

-0cm

999

Bongs, K.; Bennett, S.; Lohmann, A. Quantum sensors will start a revolution — if we deploy them right. **2023**, *617*, 672–675.

Bongs, K.; Holynski, M.; Vovrosh, J.; Bouyer, P.; Condon, G.; Rasel, E.; Schubert, C.; Schleich, W.P.; Roura, A. Taking atom interferometric quantum sensors from the laboratory to real-world applications. **2019**, *1*, 731–739.

Cladé, P.; Guellati-Khélifa, S.; Schwob, C.; Nez, F.; Julien, L.; Biraben, F. A promising method for the measurement of the local acceleration of gravity using Bloch oscillations of ultracold atoms in a vertical standing wave. **2005**, *71*, 730.

Ferrari, G.; Poli, N.; Sorrentino, F.; Tino, G. Long-lived Bloch oscillations with bosonic Sr atoms and application to gravity measurement at the micrometer scale. **2006**, *97*, 060402.

Tino, G.M. Testing gravity with cold atom interferometry: results and prospects. **2021**, *6*, 024014.

Bolotin, Y.; Golovan, A.; Vyazmin, V.; Evstifeev, M.; Vershovskii, A. Advanced gravity field survey methods. In *Methods and Technologies for Measuring the Earth’s Gravity Field Parameters*; Springer, 2022; pp. 237–301.

Ge, G.; Chen, X.; Li, J.; Zhang, D.; He, M.; Wang, W.; Zhou, Y.; Zhong, J.; Tang, B.; Fang, J.; et al. Accuracy improvement of a compact $^{85}$Rb atom gravimeter by suppressing laser crosstalk and light shift. **2023**, *23*, 6115.

Mao, D.K.; Deng, X.B.; Luo, H.Q.; Xu, Y.Y.; Zhou, M.K.; Duan, X.C.; Hu, Z.K. A dual-magneto-optical-trap atom gravity gradiometer for determining the Newtonian gravitational constant. **2021**, *92*.

Dutta, I.; Savoie, D.; Fang, B.; Venon, B.; Alzar, C.G.; Geiger, R.; Landragin, A. Continuous cold-atom inertial sensor with 1 nrad/sec rotation stability. **2016**, *116*, 183003.

Garrido Alzar, C.L. Compact chip-scale guided cold atom gyrometers for inertial navigation: Enabling technologies and design study. **2019**, *1*.

Keil, M.; Amit, O.; Zhou, S.; Groswasser, D.; Japha, Y.; Folman, R. Fifteen years of cold matter on the atom chip: promise, realizations, and prospects. **2016**, *63*, 1840–1885.

McGilligan, J.P.; Moore, K.; Dellis, A.; Martinez, G.; de Clercq, E.; Griffin, P.; Arnold, A.; Riis, E.; Boudot, R.; Kitching, J. Laser cooling in a chip-scale platform. **2020**, *117*.

Li, D.; He, W.; Shi, S.; Wu, B.; Xiao, Y.; Lin, Q.; Li, L. Review of atom chips for absolute gravity sensors. **2023**, *23*, 5089.

Ol’shanii, M.; Ovchinnikov, Y.B.; Letokhov, V. Laser guiding of atoms in a hollow optical fiber. **1993**, *98*, 77–79.

Barnett, A.; Smith, S.; Olshanii, M.; Johnson, K.; Adams, A.; Prentiss, M. Substrate-based atom waveguide using guided two-color evanescent light fields. **2000**, *61*, 023608.

Burke Jr, J.P.; Chu, S.T.; Bryant, G.W.; Williams, C.J.; Julienne, P.S. Designing neutral-atom nanotraps with integrated optical waveguides. **2002**, *65*, 043411.

Ovchinnikov, Y.B.; Shul’ga, S.; Balykin, V. An atomic trap based on evanescent light waves. **1991**, *24*, 3173.

Ovchinnikov, Y.B.; Ayi-Yovo, F.E. Towards all-optical atom chips based on optical waveguides. **2020**, *22*, 053003.

Ovchinnikov, Y.B. A perspective on integrated atomo-photonic waveguide circuits. **2022**, *120*, 010502.

Konopsky, V. Design of 1D photonic crystals sustaining optical surface modes. **2022**, *12*, 1489.

Arnaud, J.A.; Saleh, A.A.M. Guidance of surface waves by multilayer coatings. **1974**, *13*, 2343–2345.

Yeh, P.; Yariv, A.; Hong, C.S. lectromagnetic propagation in periodic stratified media. I. General theory. **1977**, *67*, 423–438.

Yeh, P.; Yariv, A.; Cho, A.Y. Optical surface waves in periodic layered media. **1978**, *32*, 104–105.

Robertson, W.M.; May, M.S. Surface electromagnetic waves on one-dimensional photonic band gap arrays. **1999**, *74*, 1800–1802.

Shinn, A.; Robertson, W. Surface plasmon-like sensor based on surface electromagnetic waves in a photonic band-gap material. **2005**, *105*, 360–364.

Konopsky, V.N.; Alieva, E.V. Long-range propagation of plasmon polaritons in a thin metal film on a one-dimensional photonic crystal surface. **2006**, *97*, 253 904.

Hamidi, S.; Ramezani, R.; Bananej, A. Hydrogen gas sensor based on long-range surface plasmons in lossy palladium film placed on photonic crystal stack. **2016**, *53*, 201–208.

Ignatyeva, D.O.; Knyazev, G.A.; Kapralov, P.O.; Dietler, G.; Sekatskii, S.K.; Belotelov, V.I. Magneto-optical plasmonic heterostructure with ultranarrow resonance for sensing applications. **2016**, *6*, 28077.

Li, J.; Tang, T.; Zhang, Y.; Luo, L.; Sun, P. Magneto-plasmonic sensor with one dimensional photonic crystal for methane detection. **2018**, *155*, 74–80.

Ignatyeva, D.; Kapralov, P.; Golovko, P.; Shilina, P.; Khramova, A.; Sekatskii, S.; Nur-E-Alam, M.; Alameh, K.; Vasiliev, M.; Kalish, A.; et al. Sensing of surface and bulk refractive index using magnetophotonic crystal with hybrid magneto-optical response. **2021**, *21*, 1984.

Konopsky, V.N.; Alieva, E.V. Photonic crystal surface waves for optical biosensors. **2007**, *79*, 4729–4735.

Guo, Y.; Ye, J.Y.; Divin, C.; Huang, B.; Thomas, T.P.; Baker, Jr., J.R.; Norris, T.B. Real-time biomolecular binding detection using a sensitive photonic crystal biosensor. **2010**, *82*, 5211–5218.

Konopsky, V.N.; Karakouz, T.; Alieva, E.V.; Vicario, C.; Sekatskii, S.K.; Dietler, G. Photonic crystal biosensor based on optical surface waves. **2013**, *13*, 2566–2578.

Rivolo, P.; Michelotti, F.; Frascella, F.; Digregorio, G.; Mandracci, P.; Dominici, L.; Giorgis, F.; Descrovi, E. Real time secondary antibody detection by means of silicon-based multilayers sustaining Bloch surface waves. **2012**, *161*, 1046–1052.

Khodami, M.; Hirbodvash, Z.; Krupin, O.; Wong, W.R.; Lisicka-Skrzek, E.; Northfield, H.; Hahn, C.; Berini, P. Fabrication of Bloch long range surface plasmon waveguides integrating counter electrodes and microfluidic channels for multimodal biosensing. **2021**, *30*, 686–695.

Kalas, B.; Ferencz, K.; Saftics, A.; Czigany, Z.; Fried, M.; Petrik, P. Bloch surface waves biosensing in the ultraviolet wavelength range – Bragg structure design for investigating protein adsorption by in situ Kretschmann-Raether ellipsometry. **2021**, *536*, 147869.

Delfan, A.; Liscidini, M.; Sipe, J.E. Surface enhanced Raman scattering in the presence of multilayer dielectric structures. **2012**, *29*, 1863–1874.

Konopsky, V.N.; Alieva, E.V.; Alyatkin, S.Y.; Melnikov, A.A.; Chekalin, S.V.; Agranovich, V.M. Phase-matched third-harmonic generation via doubly resonant optical surface modes in 1D photonic crystals. **2016**, *5*, e16168.

Fong, N.R.; Menotti, M.; Lisicka-Skrzek, E.; Northfield, H.; Olivieri, A.; Tait, N.; Liscidini, M.; Berini, P. Bloch long-range surface plasmon polaritons on metal stripe waveguides on a multilayer substrate. **2017**, *4*, 593–599.

Degli-Eredi, I.; Sipe, J.; Vermeulen, N. -polarized graphene modes sustained by photonic crystal structures. **2015**, *40*, 2076–2079.

Konopsky, V.; Prokhorov, V.; Lypenko, D.; Dmitriev, A.; Alieva, E.; Dietler, G.; Sekatskii, S. Electrical excitation of long-range surface plasmons in PC/OLED structure with two metal nanolayers. **2020**, *12*, 1–8.

Kovalevich, T.; Belharet, D.; Robert, L.; Ulliac, G.; Kim, M.S.; Herzig, H.P.; Grosjean, T.; Bernal, M.P. Bloch surface waves at the telecommunication wavelength with lithium niobate as the top layer for integrated optics. **2019**, *58*, 1757–1762.

Konopsky, V., 2022. <https://www.pcbiosensors.com/1DPC4all.htm>, Version 1.1.8401.20377.

Grimm, R.; Weidemüller, M.; Ovchinnikov, Y.B. Optical dipole traps for neutral atoms. In *Advances in atomic, molecular, and optical physics*; Elsevier, 2000; Vol. 42, pp. 95–170.

Ivanova, A.; Stoffer, R.; Hammer, M. A variational mode solver for optical waveguides based on quasi-analytical vectorial slab mode expansion. **2013**.

Ivanova, O.; Hammer, M.; Stoffer, R.; Van Groesen, E. A variational mode expansion mode solver. **2007**, *39*, 849–864.

Zhu, Z.; Brown, T.G. Full-vectorial finite-difference analysis of microstructured optical fibers. **2002**, *10*, 853–864.

Miller, J.; Cline, R.; Heinzen, D. Far-off-resonance optical trapping of atoms. **1993**, *47*, R4567.

Thompson, J.D.; Tiecke, T.; de Leon, N.P.; Feist, J.; Akimov, A.; Gullans, M.; Zibrov, A.S.; Vuletić, V.; Lukin, M.D. Coupling a single trapped atom to a nanoscale optical cavity. **2013**, *340*, 1202–1205.

Sidqi, N.; Clark, C.; Buller, G.S.; Thalluri, G.K.V.; Mitrofanov, J.; Noblet, Y. Comparative study of dielectric coating materials for micro-cavity applications. **2019**, *9*, 3452–3468.

## References

1. Bongs, K.; Bennett, S.; Lohmann, A. newblock Quantum sensors will start a revolution --- if we deploy them right. em Nature bf 2023, em 617, 672--675.

2. Bongs, K.; Holynski, M.; Vovrosh, J.; Bouyer, P.; Condon, G.; Rasel, E.; Schubert, C.; Schleich, W.P.; Roura, A. newblock Taking atom interferometric quantum sensors from the laboratory to real-world applications. em Nature Reviews Physics bf 2019, em 1, 731--739.

3. Clad'e, P.; Guellati-Kh'elifa, S.; Schwob, C.; Nez, F.; Julien, L.; Biraben, F. newblock A promising method for the measurement of the local acceleration of gravity using Bloch oscillations of ultracold atoms in a vertical standing wave. em Europhysics Letters bf 2005, em 71, 730.

4. Ferrari, G.; Poli, N.; Sorrentino, F.; Tino, G. newblock Long-lived Bloch oscillations with bosonic Sr atoms and application to gravity measurement at the micrometer scale. em Physical Review Letters bf 2006, em 97, 060402.

5. Tino, G.M. newblock Testing gravity with cold atom interferometry: results and prospects. em Quantum Science and Technology bf 2021, em 6, 024014.

6. Bolotin, Y.; Golovan, A.; Vyazmin, V.; Evstifeev, M.; Vershovskii, A. newblock Advanced gravity field survey methods. In em Methods and Technologies for Measuring the Earth’s Gravity Field Parameters; Springer, 2022; pp. 237--301.

7. Ge, G.; Chen, X.; Li, J.; Zhang, D.; He, M.; Wang, W.; Zhou, Y.; Zhong, J.; Tang, B.; Fang, J.; et al. newblock Accuracy improvement of a compact $^85$Rb atom gravimeter by suppressing laser crosstalk and light shift. em Sensors bf 2023, em 23, 6115.

8. Mao, D.K.; Deng, X.B.; Luo, H.Q.; Xu, Y.Y.; Zhou, M.K.; Duan, X.C.; Hu, Z.K. newblock A dual-magneto-optical-trap atom gravity gradiometer for determining the Newtonian gravitational constant. em Review of Scientific Instruments bf 2021, em 92.

9. Dutta, I.; Savoie, D.; Fang, B.; Venon, B.; Alzar, C.G.; Geiger, R.; Landragin, A. newblock Continuous cold-atom inertial sensor with 1 nrad/sec rotation stability. em Physical Review Letters bf 2016, em 116, 183003.

10. Garrido Alzar, C.L. newblock Compact chip-scale guided cold atom gyrometers for inertial navigation: Enabling technologies and design study. em AVS Quantum Science bf 2019, em 1.

11. Keil, M.; Amit, O.; Zhou, S.; Groswasser, D.; Japha, Y.; Folman, R. newblock Fifteen years of cold matter on the atom chip: promise, realizations, and prospects. em Journal of Modern Optics bf 2016, em 63, 1840--1885.

12. McGilligan, J.P.; Moore, K.; Dellis, A.; Martinez, G.; de Clercq, E.; Griffin, P.; Arnold, A.; Riis, E.; Boudot, R.; Kitching, J. newblock Laser cooling in a chip-scale platform. em Applied Physics Letters bf 2020, em 117.

13. Li, D.; He, W.; Shi, S.; Wu, B.; Xiao, Y.; Lin, Q.; Li, L. newblock Review of atom chips for absolute gravity sensors. em Sensors bf 2023, em 23, 5089.

14. Ol'shanii, M.; Ovchinnikov, Y.B.; Letokhov, V. newblock Laser guiding of atoms in a hollow optical fiber. em Optics Communications bf 1993, em 98, 77--79.

15. Barnett, A.; Smith, S.; Olshanii, M.; Johnson, K.; Adams, A.; Prentiss, M. newblock Substrate-based atom waveguide using guided two-color evanescent light fields. em Physical Review A bf 2000, em 61, 023608.

16. Burke Jr, J.P.; Chu, S.T.; Bryant, G.W.; Williams, C.J.; Julienne, P.S. newblock Designing neutral-atom nanotraps with integrated optical waveguides. em Physical Review A bf 2002, em 65, 043411.

17. Ovchinnikov, Y.B.; Shul'ga, S.; Balykin, V. newblock An atomic trap based on evanescent light waves. em Journal of Physics B: Atomic, Molecular and Optical Physics bf 1991, em 24, 3173.

18. Ovchinnikov, Y.B.; Ayi-Yovo, F.E. newblock Towards all-optical atom chips based on optical waveguides. em New Journal of Physics bf 2020, em 22, 053003.

19. Ovchinnikov, Y.B. newblock A perspective on integrated atomo-photonic waveguide circuits. em Applied Physics Letters bf 2022, em 120, 010502.

20. Konopsky, V. newblock Design of 1D photonic crystals sustaining optical surface modes. em Coatings bf 2022, em 12, 1489.

21. Arnaud, J.A.; Saleh, A.A.M. newblock Guidance of surface waves by multilayer coatings. em Appl. Opt. bf 1974, em 13, 2343--2345.

22. Yeh, P.; Yariv, A.; Hong, C.S. Electromagnetic propagation in periodic stratified media. I. General theory. em J. Opt. Soc. Am. bf 1977, em 67, 423--438.

23. Yeh, P.; Yariv, A.; Cho, A.Y. newblock Optical surface waves in periodic layered media. em Appl. Phys. Lett. bf 1978, em 32, 104--105.

24. Robertson, W.M.; May, M.S. newblock Surface electromagnetic waves on one-dimensional photonic band gap arrays. em Appl. Phys. Lett. bf 1999, em 74, 1800--1802.

25. Shinn, A.; Robertson, W. newblock Surface plasmon-like sensor based on surface electromagnetic waves in a photonic band-gap material. em Sens. Actuator B-Chem. bf 2005, em 105, 360--364.

26. Konopsky, V.N.; Alieva, E.V. newblock Long-range propagation of plasmon polaritons in a thin metal film on a one-dimensional photonic crystal surface. em Phys. Rev. Lett. bf 2006, em 97, 253,904.

27. Hamidi, S.; Ramezani, R.; Bananej, A. newblock Hydrogen gas sensor based on long-range surface plasmons in lossy palladium film placed on photonic crystal stack. em Optical Materials bf 2016, em 53, 201--208.

28. Ignatyeva, D.O.; Knyazev, G.A.; Kapralov, P.O.; Dietler, G.; Sekatskii, S.K.; Belotelov, V.I. newblock Magneto-optical plasmonic heterostructure with ultranarrow resonance for sensing applications. em Scientific Reports bf 2016, em 6, 28077.

29. Li, J.; Tang, T.; Zhang, Y.; Luo, L.; Sun, P. newblock Magneto-plasmonic sensor with one dimensional photonic crystal for methane detection. em Optik bf 2018, em 155, 74--80.

30. Ignatyeva, D.; Kapralov, P.; Golovko, P.; Shilina, P.; Khramova, A.; Sekatskii, S.; Nur-E-Alam, M.; Alameh, K.; Vasiliev, M.; Kalish, A.; et al. newblock Sensing of surface and bulk refractive index using magnetophotonic crystal with hybrid magneto-optical response. em Sensors bf 2021, em 21, 1984.

31. Konopsky, V.N.; Alieva, E.V. newblock Photonic crystal surface waves for optical biosensors. em Anal. Chem. bf 2007, em 79, 4729--4735.

32. Guo, Y.; Ye, J.Y.; Divin, C.; Huang, B.; Thomas, T.P.; Baker, Jr., J.R.; Norris, T.B. newblock Real-time biomolecular binding detection using a sensitive photonic crystal biosensor. em Anal. Chem. bf 2010, em 82, 5211--5218.

33. Konopsky, V.N.; Karakouz, T.; Alieva, E.V.; Vicario, C.; Sekatskii, S.K.; Dietler, G. newblock Photonic crystal biosensor based on optical surface waves. em Sensors bf 2013, em 13, 2566--2578.

34. Rivolo, P.; Michelotti, F.; Frascella, F.; Digregorio, G.; Mandracci, P.; Dominici, L.; Giorgis, F.; Descrovi, E. newblock Real time secondary antibody detection by means of silicon-based multilayers sustaining Bloch surface waves. em Sensors and Actuators B: Chemical bf 2012, em 161, 1046--1052.

35. Khodami, M.; Hirbodvash, Z.; Krupin, O.; Wong, W.R.; Lisicka-Skrzek, E.; Northfield, H.; Hahn, C.; Berini, P. newblock Fabrication of Bloch long range surface plasmon waveguides integrating counter electrodes and microfluidic channels for multimodal biosensing. em Journal of Microelectromechanical Systems bf 2021, em 30, 686--695.

36. Kalas, B.; Ferencz, K.; Saftics, A.; Czigany, Z.; Fried, M.; Petrik, P. newblock Bloch surface waves biosensing in the ultraviolet wavelength range -- Bragg structure design for investigating protein adsorption by in situ Kretschmann-Raether ellipsometry. em Applied Surface Science bf 2021, em 536, 147869.

37. Delfan, A.; Liscidini, M.; Sipe, J.E. newblock Surface enhanced Raman scattering in the presence of multilayer dielectric structures. em JOSA B bf 2012, em 29, 1863--1874.

38. Konopsky, V.N.; Alieva, E.V.; Alyatkin, S.Y.; Melnikov, A.A.; Chekalin, S.V.; Agranovich, V.M. newblock Phase-matched third-harmonic generation via doubly resonant optical surface modes in 1D photonic crystals. em Light: Science & Applications bf 2016, em 5, e16168.

39. Fong, N.R.; Menotti, M.; Lisicka-Skrzek, E.; Northfield, H.; Olivieri, A.; Tait, N.; Liscidini, M.; Berini, P. newblock Bloch long-range surface plasmon polaritons on metal stripe waveguides on a multilayer substrate. em ACS Photonics bf 2017, em 4, 593--599.

40. Degli-Eredi, I.; Sipe, J.; Vermeulen, N. TE-polarized graphene modes sustained by photonic crystal structures. em Optics Letters bf 2015, em 40, 2076--2079.

41. Konopsky, V.; Prokhorov, V.; Lypenko, D.; Dmitriev, A.; Alieva, E.; Dietler, G.; Sekatskii, S. newblock Electrical excitation of long-range surface plasmons in PC/OLED structure with two metal nanolayers. em Nano-Micro Letters bf 2020, em 12, 1--8.

42. Kovalevich, T.; Belharet, D.; Robert, L.; Ulliac, G.; Kim, M.S.; Herzig, H.P.; Grosjean, T.; Bernal, M.P. newblock Bloch surface waves at the telecommunication wavelength with lithium niobate as the top layer for integrated optics. em Applied Optics bf 2019, em 58, 1757--1762.

43. Konopsky, V., 2022. newblock https://www.pcbiosensors.com/1DPC4all.htm, Version 1.1.8401.20377.

44. Grimm, R.; Weidem"uller, M.; Ovchinnikov, Y.B. newblock Optical dipole traps for neutral atoms. In em Advances in atomic, molecular, and optical physics; Elsevier, 2000; Vol. 42, pp. 95--170.

45. Ivanova, A.; Stoffer, R.; Hammer, M. newblock A variational mode solver for optical waveguides based on quasi-analytical vectorial slab mode expansion. em arXiv preprint arXiv:1307.1315 bf 2013.

46. Ivanova, O.; Hammer, M.; Stoffer, R.; Van Groesen, E. newblock A variational mode expansion mode solver. em Optical and Quantum Electronics bf 2007, em 39, 849--864.

47. Zhu, Z.; Brown, T.G. newblock Full-vectorial finite-difference analysis of microstructured optical fibers. em Optics Express bf 2002, em 10, 853--864.

48. Miller, J.; Cline, R.; Heinzen, D. newblock Far-off-resonance optical trapping of atoms. em Physical Review A bf 1993, em 47, R4567.

49. Thompson, J.D.; Tiecke, T.; de Leon, N.P.; Feist, J.; Akimov, A.; Gullans, M.; Zibrov, A.S.; Vuleti'c, V.; Lukin, M.D. newblock Coupling a single trapped atom to a nanoscale optical cavity. em Science bf 2013, em 340, 1202--1205.

50. Sidqi, N.; Clark, C.; Buller, G.S.; Thalluri, G.K.V.; Mitrofanov, J.; Noblet, Y. newblock Comparative study of dielectric coating materials for micro-cavity applications. em Optical Materials Express bf 2019, em 9, 3452--3468. thebibliography.

---
Generated: 2026-08-30 from Z:\ValeryData\Valery_New\my_articles\Sensors2023\Sensors\manuscript\Konopsky2Sensors2023.tex, validated against Z:\ValeryData\Valery_New\my_articles\Sensors2023\Sensors\published\Konopsky_Sensors-23-08812.pdf.

