Modeling high-order harmonic generation in quantum dots using a real-space tight-binding approach

This paper introduces a computationally efficient three-dimensional real-space tight-binding model, parameterized via DFT and Wannierization, to accurately simulate the size-dependent suppression of high-order harmonic generation in quantum dots, thereby filling a theoretical gap for medium-sized nanostructures that cannot be described by existing methods for periodic solids or small molecules.

The Gist

Here is an explanation of the paper, translated from complex physics jargon into a story about tiny, glowing marbles and a new way to predict how they dance.

The Big Idea: The "Goldilocks" Problem of Tiny Lights

Imagine you have a giant, solid block of crystal (like a diamond). If you hit it with a super-powerful laser, it doesn't just glow; it starts singing. It takes the low note of the laser and spits out a whole choir of higher, sharper notes. Scientists call this High-Order Harmonic Generation (HHG). It's like a guitar string that, when plucked hard, suddenly produces a perfect symphony of overtones.

Recently, scientists tried this with Quantum Dots (QDs). Think of these as tiny, isolated islands of crystal, so small they are only a few nanometers wide (thousands of times smaller than a human hair).

The Mystery:
When they shined the laser on these tiny islands, something weird happened.

• Big islands (3nm+): They sang beautifully, just like the big crystal block.
• Tiny islands (<3nm): They went completely silent. The "singing" stopped.

The problem? Scientists didn't have a computer model that could explain why the tiny dots went silent.

• The "Molecule" tools were too slow and clunky for these medium-sized dots (too many atoms to count individually).
• The "Solid Crystal" tools assumed the material was infinite and endless, which doesn't work for a tiny, finite dot.

It was like trying to measure the wind in a hurricane using a thermometer meant for a teacup, or trying to predict the weather of a single room using a model for the whole planet. Neither worked.

The Solution: A New "Digital Lego" Model

The authors of this paper built a new, super-fast computer model to solve this. They call it a Real-Space Tight-Binding Model.

Here is how it works, using a creative analogy:

1. The "Digital Lego" Approach (Tight-Binding)
Imagine the quantum dot is built out of Lego bricks. Instead of trying to simulate every single atom's movement from scratch (which is like trying to simulate the physics of every single molecule of air in a room), this model treats the Lego bricks as connected blocks.

• They used a "map" (derived from complex supercomputer calculations called DFT) to know exactly how these Lego bricks connect.
• They built the model in 3D, allowing them to cut out a perfect sphere of these bricks to represent the tiny dot.

2. The "Wannier" Translator
To make the math work, they used a special translator called Wannier functions.

• Analogy: Imagine you have a song written in a complex musical score (Bloch functions, used for infinite crystals). It's hard to read for a single instrument. The Wannier translator rewrites that song into simple, localized notes that only belong to specific instruments (atoms). This lets them see how the "music" (the electrons) behaves when the instrument is isolated in a small room (the dot).

3. The "Ghost" vs. The "Real Thing"
The authors built two versions of their model to prove it works:

• The "Finite" Model: This simulates the actual tiny dot with hard edges. It's like simulating a single, isolated drum.
• The "Periodic" Model: This simulates an endless wall of drums. They used this to check if their math was right by comparing it to known results for big crystals.

What They Discovered

Once they turned on their new model, the mystery was solved.

Why did the tiny dots go silent?
The model showed that the "singing" (HHG) relies on a three-step dance:

1. Kick: The laser kicks an electron loose.
2. Run: The electron runs away, gaining speed.
3. Crash: The electron crashes back into its hole, releasing a burst of light (the harmonic).

In a big crystal, the electron has a long runway to build up speed before it crashes.
In a tiny dot (<3nm), the walls are too close!

• The electron tries to run, but it hits the "wall" of the dot too early.
• It gets scattered and confused (dephasing).
• It never builds up enough speed to make a loud "crash" sound.
• Result: The light is suppressed. The dot goes silent.

The Ellipticity Test
They also tested what happens if the laser light spins (elliptical polarization) instead of just shaking back and forth.

• In big crystals, spinning the light kills the singing very quickly.
• Their model showed that tiny dots behave similarly, but the "killing" effect is even more sensitive to the size of the dot. This confirmed that the model captures the real physics of these nano-worlds.

Why This Matters

This paper is a bridge.

• Before, we could model atoms (too small) or big crystals (too big).
• Now, we have a tool that perfectly handles the "Goldilocks" zone: medium-sized nanostructures (like quantum dots).

The Impact:

1. Speed: Their model is incredibly fast. It can simulate a complex quantum dot in minutes on a graphics card, whereas previous methods would take days or were impossible.
2. Design: Engineers can now use this model to design better quantum dots for future technologies (like ultra-fast computers or new types of lasers) without having to build them in a lab first. They can simulate how changing the size or shape will affect the light output.

In a nutshell: The authors built a fast, accurate "digital twin" for tiny quantum dots. They used it to prove that when these dots get too small, the electrons hit the walls too hard to sing, explaining a mystery that experimentalists had been scratching their heads over for years.

Technical Summary

Here is a detailed technical summary of the paper "Modelling HHG in QDs using a real-space tight-binding approach" by Thümmler et al.

1. Problem Statement

High-order harmonic generation (HHG) in semiconductor quantum dots (QDs) has recently been investigated experimentally, revealing a strong suppression of HHG yield in dots smaller than 3 nm, particularly for longer driving wavelengths. However, a comprehensive theoretical framework to explain these size-dependent effects has been missing due to computational limitations:

• Small Molecules/Atoms: Methods like direct Schrödinger equation solutions or real-time time-dependent density functional theory (rt-TDDFT) are too computationally expensive for QDs containing hundreds of atoms.
• Periodic Solids: Standard models like Semiconductor Bloch Equations (SBEs) or periodic tight-binding models assume infinite periodicity and cannot account for the finite size and boundary effects of QDs.

The authors aim to bridge this gap by developing a computationally efficient model capable of simulating HHG in medium-sized nanostructures (up to ~500 atoms) while capturing finite-size effects.

2. Methodology

The authors propose a three-dimensional real-space tight-binding (TB) model tailored for confined systems. The methodology consists of the following key components:

A. Parameterization from First Principles

• DFT Calculations: The model parameters are derived from Density Functional Theory (DFT) calculations for the semiconductor bulk (Silicon and CdSe wurtzite) using Quantum Espresso.
• Wannierization: Maximally localized Wannier functions (MLWFs) are generated using Wannier90. The authors separate the valence and conduction bands to create orthogonal sets of Wannier functions.
• Hamiltonian Construction: The field-free Hamiltonian ($\hat{H}_0$) and dipole matrix elements are constructed in the Wannier basis. This allows the model to naturally include the periodic limit while being adaptable to finite geometries.

B. The Finite-Size Model (Core Contribution)

• Geometry: The QD is modeled as a spherical cut-out of the periodic crystal structure, defined by unit cell boundaries.
• Operators: The system is described using electron ($\hat{e}^\dagger$) and hole ($\hat{h}^\dagger$) creation operators acting on Wannier states.
• Density Matrix Propagation: The time evolution is governed by the von-Neumann equation for the density matrix $\hat{\rho}$. The authors neglect Coulomb interactions (justified by the dominance of field-induced energy shifts over exciton binding energies in strong fields) and ionization effects.
• Numerical Implementation:
  • The Liouvillian operator is split into field-free and field-dependent parts.
  • The density matrix is vectorized, and the equation of motion is solved using a Runge-Kutta 4/5 algorithm on GPUs (using vexcl and boost.odeint).
  • Optimization: To handle the large matrix sizes, the authors exploit the sparsity of the system and the hermiticity of the density matrix, avoiding explicit calculation of redundant matrix elements. This reduces computational cost significantly compared to naive implementations.

C. Real-Space Periodic Model (Benchmark)

To validate the finite-size approach and compare it with established methods, the authors developed a real-space periodic model.

• This model tracks coherences between Wannier sites that shift simultaneously (preserving translational symmetry).
• It is shown to be mathematically equivalent to the Semiconductor Bloch Equations (SBEs) in the moving-frame Wannier basis in the limit of infinite coherence range.

3. Key Results

A. Validation and Benchmarking

• Linear Response: The periodic model's optical conductivity (absorption spectra) for Si and CdSe matches well with results from the Kubo formula and rt-TDDFT calculations, validating the Wannier parameterization.
• Nonlinear Response (Bulk): For bulk Silicon, the periodic model reproduces HHG spectra qualitatively similar to rt-TDDFT (Octopus code) but with a significant speedup (20 minutes vs. several days).
• Nonlinear Response (Small Clusters): For a small Cd$_{11}$Se$_{11}$ cluster, the finite-size model reproduces the rt-TDDFT result: a complete suppression of high-order harmonics, consistent with experimental observations.

B. Quantum Dot Size Dependence

The model successfully reproduces the experimentally observed trend where HHG yield drops drastically for QDs smaller than 3 nm:

• Mechanism: The suppression is attributed to the finite size effect on the three-step model of HHG. While electron-hole pair excitation is unaffected, the acceleration and recombination steps are hindered by the dot boundaries.
• Wavelength Dependence: The suppression is more pronounced for longer driving wavelengths (e.g., 5 $\mu$m vs. 3 $\mu$m) because longer wavelengths induce larger electron quiver amplitudes, leading to more frequent scattering and dephasing at the QD boundaries.

C. Ellipticity Dependence

• The authors simulated the HHG yield for elliptically polarized pulses.
• Result: The yield decreases as ellipticity increases (approaching circular polarization), following a Gaussian-like dependence.
• Comparison: This behavior mirrors that of bulk solids, but the specific width of the ellipticity dependence varies with driving wavelength and energy range, indicating a combination of bulk-like behavior and finite-size boundary effects.

4. Significance and Impact

• Bridging the Gap: This work provides the first computationally feasible model to simulate HHG in medium-sized nanostructures (up to ~3.3 nm diameter, ~1000 Wannier functions), a regime inaccessible to rt-TDDFT and unsuitable for purely periodic models.
• Physical Insight: It offers a microscopic explanation for the experimental suppression of HHG in small QDs, linking it directly to boundary-induced dephasing and scattering.
• Efficiency: The GPU-accelerated implementation allows for the simulation of ensemble averages over random orientations, which is crucial for comparing with experimental data where QDs are randomly oriented.
• Transferability: Since the model is parameterized from ab initio DFT calculations, it is transferable to other semiconductor materials and nanostructure geometries, facilitating the design of nanostructures for specific nonlinear optical applications.

In conclusion, the paper presents a robust, efficient, and physically grounded theoretical tool that successfully explains the size-dependent nonlinear optical response of quantum dots, filling a critical void in the theoretical description of strong-field physics in confined nanostructures.