Self-Consistent Theoretical Framework for Third-Order Nonlinear Susceptibility in CdSe/ZnS--MOF Quantum Dot Composites

TL;DR

This work addresses predicting the third-order nonlinear optical response, $\chi^{(3)}(\omega)$, of CdSe/ZnS–MOF composites by linking microscopic confinement to macroscopic nonlinear optics. It proposes a self-consistent chain: finite-potential quantum confinement via EMA, a density-matrix expansion to third order, Maxwell–Garnett/Bruggeman homogenization with local-field factors, and Kramers–Kronig consistency to map to observable Z-scan quantities. Key findings include a confinement-induced bandgap blue shift of $\Delta E \approx 0.2$–$0.3$ eV for $R=2.5$–$4.0$ nm, a dominant two-photon resonance near $\lambda \approx 1.2\,\mu$m, and a macroscopic response $χ^{(3)}$ of order $10^{-22}$ m$^2$/V$^2$ with $n_2<0$ and rising $\beta$ with wavelength, consistent with experimental trends. The framework yields design maps showing how shell thickness, dielectric environment, and loading fraction tune $χ^{(3)}$, $n_2$, and $β$, offering a path to engineering strong, tunable nonlinearities in hybrid QD–MOF materials.

Abstract

This work presents a fully theoretical and self consistent framework for calculating the third-order nonlinear susceptibility of CdSe/ZnS--MOF composite quantum dots. The approach unifies finite-potential quantum confinement,the Liouville von Neumann density matrix expansion to third order, and effective-medium electrodynamics (Maxwell--Garnett and Bruggeman) within a single Hamiltonian-based model, requiring no empirical fitting. Electron hole quantized states and dipole matrix elements are obtained under the effective-mass approximation with BenDaniel--Duke boundary conditions; closed analytic forms for(including Lorentzian/Voigt broadening) follow from the response expansion. Homogenization yields macroscopic scaling laws that link microscopic descriptors (core radius, shell thickness, dielectric mismatch) to bulk coefficients and. A Kramers--Kronig consistency check confirms causality and analyticity of the computed spectra with small residuals. The formalism provides a predictive, parameter-transparent route to engineer third-order nonlinearity in hybrid quantum materials,clarifying how size and environment govern the magnitude and dispersion of.

Figures (6)

• Figure 1: Finite spherical potential well and quantized levels in a CdSe/ZnS–MOF quantum dot. The 3D rendering shows the CdSe core, ZnS shell, and MOF matrix; internal guides mark the confined electron and hole states $E_{1e}$ and $E_{1h}$.
• Figure 2: Contour map of the confinement-enhanced bandgap $E_g(R,t)$ for a CdSe/ZnS–MOF core–shell dot. Parameters: $E_g^{\mathrm{bulk}}{=}1.74~\mathrm{eV}$, $m_e^\ast{=}0.13m_0$, $m_h^\ast{=}0.45m_0$, $\varepsilon_i{=}6.0$ (CdSe), $\varepsilon_{\mathrm{shell}}{=}8.0$ (ZnS), $\varepsilon_h{=}2.1$ (MOF). Screening model $\varepsilon_{\mathrm{eff}}(t)=\varepsilon_i+(\varepsilon_{\mathrm{shell}}-\varepsilon_i)\,[1-\exp(-t/\lambda)]$, $\lambda=0.4~\mathrm{nm}$.
• Figure 3: Spectral dependence of $\chi^{(3)}(\omega)$ for radii $R=2.5,\,3.0,\,4.0~\mathrm{nm}$; solid: $\Re[\chi^{(3)}]$, dashed: $\Im[\chi^{(3)}]$. Local-field factor $L^4=(3\varepsilon_h/(\varepsilon_i+2\varepsilon_h))^4$ with $\varepsilon_h=2.1$, $\varepsilon_i=6.0$.
• Figure 4: Mapping of $\chi^{(3)}(\lambda,t)$: surface shows normalized $\Re[\chi^{(3)}]$ vs. $\lambda$ and $t$; contours denote $\Im[\chi^{(3)}]$. Parameters: $\varepsilon_i=6.0$, $\varepsilon_{\mathrm{shell}}=8.0$, $\varepsilon_h=2.1$, $\lambda_0=0.4~\mathrm{nm}$, $\mu=7\times10^{-29}~\mathrm{C\,m}$, $\gamma=3\times10^{13}~\mathrm{s^{-1}}$, $N=2.5\times10^{25}~\mathrm{m^{-3}}$.
• Figure 5: Local-field enhancement factor $L^4$ (log$_{10}$ scale) vs. loading fraction $\phi$ and host permittivity $\varepsilon_h$ for CdSe/ZnS–MOF composites. $L^4=(3\varepsilon_h/(\varepsilon_i^{\mathrm{eff}}+2\varepsilon_h))^4$ with $\varepsilon_i^{\mathrm{eff}}=\varepsilon_{\mathrm{CdSe}}+(\varepsilon_{\mathrm{ZnS}}-\varepsilon_{\mathrm{CdSe}})(1-e^{-t/\lambda_0})$, $t=0.8\,\mathrm{nm}$, $\lambda_0=0.4\,\mathrm{nm}$.