{"id":170,"date":"2025-01-13T12:09:10","date_gmt":"2025-01-13T12:09:10","guid":{"rendered":"https:\/\/qplai.liacs.nl\/?p=170"},"modified":"2025-01-13T12:35:06","modified_gmt":"2025-01-13T12:35:06","slug":"autonomous-extraction-of-semiconductor-quantum-dot-parameters-using-a-basic-fourier-transform-approach","status":"publish","type":"post","link":"https:\/\/qplai.liacs.nl\/index.php\/2025\/01\/13\/autonomous-extraction-of-semiconductor-quantum-dot-parameters-using-a-basic-fourier-transform-approach\/","title":{"rendered":"Autonomous Extraction of Semiconductor Quantum Dot Parameters Using a Basic Fourier Transform Approach"},"content":{"rendered":"\n<p>The pursuit of practical quantum computing relies on scalable architectures, and transistor-based quantum dots offer a promising path forward. By trapping single electrons in these semiconductor devices, we hope to process information using their internal degree of freedom &#8211; spin. However, achieving this requires precise control and tuning<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Tuning Quantum Dots: A Complex Prelude to Quantum Operations<\/h3>\n\n\n\n<p>Before engaging in quantum operations, quantum dots must be carefully tuned to achieve the desired number of electrons. This involves adjusting voltages on metallic gates, a process fraught with challenges like parameter drift and the difficulty of sensing charge configurations. Tuning even a small number of quantum dots manually is labor-intensive, making the task infeasible as the number of dots grows. <strong>This underlines the necessity for autonomous tuning methods to advance the field<\/strong>. (See: <a href=\"https:\/\/arxiv.org\/abs\/2407.20061\">Autonomous Bootstrapping of Quantum Dot Devices<\/a> or <a href=\"https:\/\/arxiv.org\/abs\/2402.03931\">Fully autonomous tuning of a spin qubit<\/a>)<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Charge Stability Diagrams: The Map to Electron Configuration<\/h3>\n\n\n\n<p>A critical tool in quantum dot tuning is the charge stability diagram (CSD), a two-dimensional voltage map depicting charge configurations. <\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"442\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-30-1024x442.png\" alt=\"\" class=\"wp-image-204\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-30-1024x442.png 1024w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-30-300x129.png 300w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-30-768x331.png 768w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-30.png 1126w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption class=\"wp-element-caption\">Left: typical device scheme, Right: model capacitance diagram. From: <a href=\"https:\/\/ieeexplore.ieee.org\/stamp\/stamp.jsp?arnumber=7751431\">https:\/\/ieeexplore.ieee.org\/stamp\/stamp.jsp?arnumber=7751431<\/a><\/figcaption><\/figure>\n\n\n\n<figure class=\"wp-block-image aligncenter size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" width=\"820\" height=\"432\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-28.png\" alt=\"\" class=\"wp-image-202\" style=\"width:480px;height:auto\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-28.png 820w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-28-300x158.png 300w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-28-768x405.png 768w\" sizes=\"auto, (max-width: 820px) 100vw, 820px\" \/><figcaption class=\"wp-element-caption\">Example of charge stability diagram, with mapped occupation of the electrons [Q0, Q1]. From: <a href=\"https:\/\/scipost.org\/SciPostPhysCodeb.43\">https:\/\/scipost.org\/SciPostPhysCodeb.43<\/a><\/figcaption><\/figure>\n\n\n\n<p>These diagrams are governed by the capacitance model, where electrostatic energy determines charge states. Clean and well-defined CSDs enable the extraction of model parameters, such as the size and slopes of Coulomb diamonds, which are essential for understanding and optimizing the system.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Overcoming Challenges with Scalable Strategies<\/h3>\n\n\n\n<p>As the number of quantum dots increases, the number of CSD cuts grows exponentially, making manual analysis impractical. Current strategies include:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>Autonomous identification of features like diamond edges and slopes.<\/li>\n\n\n\n<li>Using convolutional neural networks (CNNs) in a supervised fashion to extract parameters.<\/li>\n<\/ul>\n\n\n\n<p>We propose a simpler yet effective alternative: leveraging Fourier transforms to analyze CSDs.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Capacitance Model: The Physics Behind CSDs<\/h3>\n\n\n\n<p>The system can be described using a capacitance matrix that links charges and voltages on the dots and gates, analogous to the relationship \\( Q = CV \\):<\/p>\n\n\n\n<p>$$<br>\\begin{pmatrix}<br>Q_D \\\\<br>Q_G<br>\\end{pmatrix} =<br>\\begin{pmatrix}<br>C_{DD} &amp; C_{DG} \\\\<br>C_{DG}^T &amp; C_{GG}<br>\\end{pmatrix}<br>\\begin{pmatrix}<br>U_D \\\\<br>V_G<br>\\end{pmatrix}<br>$$<\/p>\n\n\n\n<p>Here:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li> \\(C_{DD}\\) : Dot-dot capacitance matrix.<\/li>\n\n\n\n<li>\\( C_{DG} \\) : Dot-gate cross-capacitance matrix.<\/li>\n\n\n\n<li>\\(C_{GG}\\) : Gate-gate capacitance matrix.<\/li>\n<\/ul>\n\n\n\n<p>Minimising the electrostatic energy, \\( E \\), determines the charge configuration. This is visualised in a CSD, with gate voltages as axes and charge states represented by color. Parameters such as the size and slope of Coulomb diamonds are directly derived from this model.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Key Parameters for Simplified Analysis<\/h3>\n\n\n\n<p>Four critical parameters derived from the capacitance model:<\/p>\n\n\n\n<ol class=\"wp-block-list\">\n<li><strong>Diamond Sizes<\/strong>: \\( \\Delta V_1, \\Delta V_2 \\)<br>\\[<br>\\Delta V_i = \\frac{|e| C_{DD}^{-1}[i,i]}{(C_{DD}^{-1}C_{DG})[i,i]}<br>\\]<\/li>\n\n\n\n<li><strong>Diamond Slopes<\/strong>: \\( \\theta_1, \\theta_2 \\)<br>\\[<br>\\theta_1 = \\text{arctan}\\left(\\frac{\\alpha[1,1]}{\\alpha[1,2]}\\right),\\quad \\theta_2 = \\text{arctan}\\left(\\frac{\\alpha[2,2]}{\\alpha[2,1]}\\right) \\quad \\alpha = C_{DD}^{-1}C_{DG}<br>\\]<\/li>\n<\/ol>\n\n\n\n<figure class=\"wp-block-image aligncenter size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" width=\"814\" height=\"812\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-25.png\" alt=\"\" class=\"wp-image-199\" style=\"width:532px;height:auto\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-25.png 814w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-25-300x300.png 300w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-25-150x150.png 150w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-25-768x766.png 768w\" sizes=\"auto, (max-width: 814px) 100vw, 814px\" \/><\/figure>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Our Procedure: Fourier Transform Analysis of CSDs<\/h3>\n\n\n\n<p>Our method focuses on two-dot systems but can be generalized to larger setups. Here\u2019s the step-by-step process:<\/p>\n\n\n\n<h4 class=\"wp-block-heading\">1. Data Acquisition<\/h4>\n\n\n\n<p>Charge stability diagrams are obtained experimentally by sweeping gate voltages and measuring responses, such as currents through sensing dots. Alternatively, synthetic data can be generated using the QDarts simulator to explore different noise levels. Below we present 10 CSDs obtained for a different level of noise (x-axis)<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"546\" height=\"404\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-12.png\" alt=\"\" class=\"wp-image-184\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-12.png 546w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-12-300x222.png 300w\" sizes=\"auto, (max-width: 546px) 100vw, 546px\" \/><\/figure>\n\n\n\n<h4 class=\"wp-block-heading\">2. Edge Extraction<\/h4>\n\n\n\n<p>We use the Sobel filter, a simple edge-detection tool with 3&#215;3 kernels, to identify horizontal and vertical edges. The resulting gradient magnitudes highlight the edges of Coulomb diamonds.<\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"560\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-14-1024x560.png\" alt=\"\" class=\"wp-image-186\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-14-1024x560.png 1024w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-14-300x164.png 300w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-14-768x420.png 768w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-14.png 1396w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/figure>\n\n\n\n<h4 class=\"wp-block-heading\">3. Fourier Transform<\/h4>\n\n\n\n<p>Applying a two-dimensional Fourier transform (fft2) to the filtered image reveals periodic structures in the CSD. The transform captures prominent \u201cjets&#8221; corresponding to the directions perpendicular to the diamond edges.<\/p>\n\n\n\n<figure class=\"wp-block-image aligncenter size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" width=\"754\" height=\"714\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-17.png\" alt=\"\" class=\"wp-image-189\" style=\"width:300px;height:auto\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-17.png 754w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-17-300x284.png 300w\" sizes=\"auto, (max-width: 754px) 100vw, 754px\" \/><\/figure>\n\n\n\n<h4 class=\"wp-block-heading\">4. Jet Identification<\/h4>\n\n\n\n<p>To identify jets, we:<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>Smooth the Fourier transform image using a Gaussian filter.<\/li>\n\n\n\n<li>Locate peaks along a specific cut of the Fourier-transformed space. These peaks define the jets, allowing us to extract angles \\( \\theta_1 \\) and \\( \\theta_2 \\).<\/li>\n<\/ul>\n\n\n\n<figure class=\"wp-block-image size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"345\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-15-1024x345.png\" alt=\"\" class=\"wp-image-187\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-15-1024x345.png 1024w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-15-300x101.png 300w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-15-768x259.png 768w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-15.png 1181w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><\/figure>\n\n\n\n<h4 class=\"wp-block-heading\">5. Period Extraction<\/h4>\n\n\n\n<p>The period of oscillations along the jets is determined by the distance from the origin to the first peak along each jet. This period is inversely proportional to the size of the Coulomb diamonds.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"691\" height=\"298\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-16.png\" alt=\"\" class=\"wp-image-188\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-16.png 691w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-16-300x129.png 300w\" sizes=\"auto, (max-width: 691px) 100vw, 691px\" \/><\/figure>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Results: Analyzing Synthetic Data<\/h3>\n\n\n\n<p>Using synthetic data generated by the QDarts simulator, we validated our method\u2019s accuracy under varying noise levels. The Fourier-based approach proved robust, efficiently extracting key parameters even in noisy conditions.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" width=\"702\" height=\"439\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-24.png\" alt=\"\" class=\"wp-image-196\" style=\"width:620px;height:auto\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-24.png 702w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-24-300x188.png 300w\" sizes=\"auto, (max-width: 702px) 100vw, 702px\" \/><\/figure>\n\n\n\n<p>Clearly below certain noise level, the relative prediction error in any of the quantites is less than 10%. In general it is likely to be dominated by the &#8220;unlucky&#8221; CSD.<\/p>\n\n\n\n<figure class=\"wp-block-image size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" width=\"540\" height=\"416\" src=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-23.png\" alt=\"\" class=\"wp-image-195\" style=\"width:608px;height:auto\" srcset=\"https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-23.png 540w, https:\/\/qplai.liacs.nl\/wp-content\/uploads\/image-23-300x231.png 300w\" sizes=\"auto, (max-width: 540px) 100vw, 540px\" \/><\/figure>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Conclusion<\/h3>\n\n\n\n<p>The Fourier transform-based method offers a simplified, yet effective way to extract parameters underlying charge stability diagrams in the autonomous way. <\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Code availability<\/h3>\n\n\n\n<p>The code is available at: <a href=\"https:\/\/github.com\/jan-a-krzywda\/ftt-charge-stability-diagram\">https:\/\/github.com\/jan-a-krzywda\/ftt-charge-stability-diagram<\/a><br>Data: <a href=\"https:\/\/leidenuniv1-my.sharepoint.com\/:u:\/g\/personal\/krzywdaja_vuw_leidenuniv_nl\/EccqnFnrldFDreY5xaaeRVEBOGKyDFh6RZMv1QUsemxrHQ?e=RX60mY\">data_csd_fft.zip<\/a><\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n\n\n\n<h3 class=\"wp-block-heading\">Outlook and limitations<\/h3>\n\n\n\n<p>TODO<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\"\/>\n","protected":false},"excerpt":{"rendered":"<p>The pursuit of practical quantum computing relies on scalable architectures, and transistor-based quantum dots offer a promising path forward. By trapping single electrons in these semiconductor devices, we hope to process information using their internal degree of freedom &#8211; spin. However, achieving this requires precise control and tuning Tuning Quantum Dots: A Complex Prelude to [&hellip;]<\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"class_list":["post-170","post","type-post","status-publish","format-standard","hentry","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/posts\/170","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/comments?post=170"}],"version-history":[{"count":18,"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/posts\/170\/revisions"}],"predecessor-version":[{"id":207,"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/posts\/170\/revisions\/207"}],"wp:attachment":[{"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/media?parent=170"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/categories?post=170"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/qplai.liacs.nl\/index.php\/wp-json\/wp\/v2\/tags?post=170"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}