CdSe/ZnS–MOF Composite Quantum Dots

• The paper demonstrates a rigorous theoretical framework linking microscopic quantum descriptors to macroscopic third-order nonlinear susceptibility without empirical fitting.
• It applies quantum confinement modeling and density-matrix expansion to quantify size- and structure-dependent exciton transitions and saturation intensities.
• The effective-medium theory and design optimization guidelines provide practical parameters for engineering optical switching, modulation, and photonic devices.

CdSe/ZnS–MOF composite quantum dots (QDs) are nanostructured hybrid materials consisting of a spherical cadmium selenide (CdSe) core encapsulated by a concentric zinc sulfide (ZnS) shell, embedded within a metal–organic framework (MOF) host matrix. These composites exploit the strong quantum confinement of semiconductor nanocrystals and the tunable dielectric/electrostatic environment provided by the MOF scaffold to engineer highly nonlinear optical responses, notably a large third-order nonlinear susceptibility, $\chi^{(3)}$. A rigorous, self-consistent theoretical framework enables quantitative, parameter-transparent prediction and control of $\chi^{(3)}$ in these hybrid systems, linking microscopic quantum descriptors to macroscopic optical observables without empirical fitting (Wu et al., 4 Nov 2025).

1. Quantum Confinement Model and Electronic Structure

The electronic states of a CdSe core of radius $R$, surrounded by a ZnS shell of thickness $t$ (total radius $R_{\rm tot}=R+t$), are modeled using the envelope-function effective-mass approximation (EMA) with BenDaniel–Duke boundary conditions. The single-particle Schrödinger equation for electrons or holes in spherical coordinates is: $$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$ with position-dependent effective mass

$$m^*(r)= \begin{cases} m_1^*, & 0\le r<R \quad (\text{CdSe core}) \ m_2^*, & R\le r<R_{\rm tot} \quad (\text{ZnS shell}) \end{cases}$$

and radial potential

$$V(r)= \begin{cases} 0, & r<R \ V_b^{(e/h)}, & R\le r<R_{\rm tot} \ \infty, & r\ge R_{\rm tot} \end{cases}$$

The transcendental eigenvalue equation derived from continuity of $\psi$ and mass-flux at $r = R$ determines quantized levels $E_{e,0},\,E_{h,0}$ and radial envelope functions $R_{e,h}(r)$. The lowest interband transition (confined exciton) energy, generalizing the Brus formula, includes finite barrier and dielectric mismatch effects: $$E_X(R,t) = E_g^{\rm bulk} + E_{e,0} + E_{h,0} - \frac{1.786\,e^2}{4\pi\varepsilon_0\,\varepsilon_{\rm in}\,R_{\rm eff}} +\frac{e^2}{8\pi\varepsilon_0\,R_{\rm eff}} \left(\frac1{\varepsilon_{\rm in}}-\frac1{\varepsilon_{\rm out}}\right)^\xi,$$ where $R_{\rm eff}=R+\alpha t$, $\varepsilon_{\rm in}$ is the core/shell average permittivity, $\varepsilon_{\rm out}$ is the MOF host permittivity, and $\xi\simeq0.2$–$0.4$. The dipole matrix element for the exciton transition is

$$\mu_{01} = e\int_0^{R_{\rm tot}} R_e(r)\,r\,R_h(r)\,4\pi r^2\,dr \equiv\mu_{cv}^{\rm bulk}S(R,t),$$

with $S$ the dimensionless envelope overlap. The two-level saturation intensity is

$$I_{\rm sat} = \frac{\varepsilon_0\,c\,n_0\,\hbar^2 \left[\gamma^2+(\omega-\omega_{01})^2\right]} {2|\mu_{01}|^2}.$$

Numerical solution of the transcendental equation yields size- and structure-dependent quantum levels and transition strengths.

2. Nonlinear Susceptibility via Density-Matrix Expansion

The nonlinear optical response is modeled by promoting the QD to a three-level system $|0\rangle, |1\rangle, |2\rangle$ (ground, exciton, biexciton), with dipole coupling and coherent light-matter interaction. The system Hamiltonian under a time-dependent electric field is: $$H(t) = H_0 - \hat{\mu}E(t),\qquad H_0 = \sum_{j=0}^2 \hbar\omega_j |j\rangle\langle j|.$$ The Liouville–von Neumann equation (with dephasing $\gamma_{ij}$ and relaxation $\Gamma_i$ via superoperator $\mathcal D[\rho]$) is expanded order-by-order in the driving field: $$\dot\rho = -\frac{i}{\hbar}[H(t),\rho] + \mathcal D[\rho],\quad \rho = \rho^{(0)} + \rho^{(1)} + \rho^{(2)} + \rho^{(3)} + \cdots.$$ At steady-state, the third-order polarization at the fundamental frequency is

$$P^{(3)}(\omega) = \varepsilon_0 \chi^{(3)}(\omega) |E(\omega)|^2 E(\omega).$$

A closed-form expression for $\chi^{(3)}$ in the degenerate Kerr configuration is

$$\chi^{(3)}(\omega) = \frac{N}{\varepsilon_0\,\hbar^3} \sum_{\mathcal P\in S_3} \frac{\mu_{01}^{2}\mu_{12}\mu_{20}} {D_{\mathcal P(10)}(\omega)\,D_{\mathcal P(20)}(2\omega)\,D_{\mathcal P(10)}(\omega)},$$

where $N$ is the QD density, the denominators $D$ encode detuning and dephasing, and the sum over permutations enforces causality. Near resonance, the leading behaviors are Lorentzian: $$\Re[\chi^{(3)}]\propto\frac{\Delta}{\Delta^2+\gamma^2},\qquad \Im[\chi^{(3)}]\propto\frac{\gamma}{\Delta^2+\gamma^2}.$$ Inhomogeneous broadening (e.g., size polydispersity) is incorporated as Gaussian convolution (Voigt profile), preserving analytic structure.

3. Homogenization and Effective-Medium Theory

For the bulk composite, each QD is treated as an inclusion (permittivity $\varepsilon_i(\omega)$, intrinsic $\chi_i^{(3)}$) in a MOF host (permittivity $\varepsilon_h$, susceptibility $\chi_h^{(3)}$), at volume fraction $\phi$. The Maxwell–Garnett (MG) and Bruggeman (Br) mixing formulas are invoked. For the MG case,

$$\varepsilon_{\rm eff} = \varepsilon_h\, \frac{\varepsilon_i+2\varepsilon_h+2\phi(\varepsilon_i-\varepsilon_h)} {\varepsilon_i+2\varepsilon_h-\phi(\varepsilon_i-\varepsilon_h)},\qquad L(\omega) = \frac{3\varepsilon_h}{\varepsilon_i+2\varepsilon_h}$$

with cubic nonlinear susceptibility

$$\chi^{(3)}_{\rm eff}(\omega) \simeq \phi|L(\omega)|^4 \chi^{(3)}_i(\omega) + (1-\phi)\chi^{(3)}_h(\omega).$$

The Bruggeman formula gives a self-consistent relation: $$\phi\,\frac{\varepsilon_i-\varepsilon_{\rm eff}}{\varepsilon_i+2\varepsilon_{\rm eff}} + (1-\phi)\frac{\varepsilon_h-\varepsilon_{\rm eff}}{\varepsilon_h+2\varepsilon_{\rm eff}} = 0,$$ with local-field factors

$$L_i^{(B)} = 3\varepsilon_{\rm eff}/(\varepsilon_i+2\varepsilon_{\rm eff}),$$

yielding

$$\chi^{(3)}_{\rm eff} \simeq \phi|L_i^{(B)}|^4 \chi^{(3)}_i + (1-\phi)|L_h^{(B)}|^4 \chi_h^{(3)}.$$

The choice of homogenization model determines sensitivity to $\phi$ and dielectric contrast; Maxwell–Garnett is accurate for dilute systems, while Bruggeman applies near percolation.

4. Scaling Laws and Design Optimization

The parameter dependencies of $\chi^{(3)}_{\rm eff}$ emerge from the interplay between QD quantum structure, local-field factors, and composite geometry: $$\chi^{(3)}_{\rm eff}(\omega)\propto \phi|L(\omega)|^4\,N\,\frac{|\mu_{01}(R,t)|^4}{\hbar^3\gamma(R,t)^3} \sim \phi|L|^4\frac{S(R,t)^4\,R^{-3}}{\gamma(R,t)^3}.$$ Critical qualitative trends are:

• Shell thickness ($t$): Increasing $t$ enhances envelope-state overlap $S$ and lowers dephasing rate $\gamma$, boosting $\chi^{(3)}$.
• Core radius ($R$): Decreasing $R$ blue-shifts the exciton resonance ($E_X\propto R^{-2}$) and modifies the absorption/scattering features.
• Host permittivity ($\varepsilon_h$) and fill fraction ($\phi$): Higher $\varepsilon_h$ and $\phi$ increase $|L|^4$ and linearly scale the overall $\chi^{(3)}_{\rm eff}$, with limits set by percolation.

Optimization for large third-order nonlinearity thus involves:

• Choosing $R \sim 2.5$–$3.5$ nm to position the two-photon resonance in the NIR regime.
• Applying a ZnS shell with $t \gtrsim 0.8$ nm to maximize $S$ and minimize $\gamma$.
• Selecting a MOF with $\varepsilon_h \in [2,4]$ and moderate $\phi$ ($0.1$–$0.3$).
• Minimizing polydispersity to narrow the resonance (limiting inhomogeneous Voigt broadening).

5. Validation of Analyticity and Causality

The physicality of the computed $\chi^{(3)}_{\rm eff}(\omega)$ is validated via a Kramers–Kronig (KK) consistency check. The real part of $\chi^{(3)}$ is reconstructed from its imaginary part by

$$\Re[\chi^{(3)}_{\rm KK}(\omega)] = \frac{1}{\pi}\mathcal P\int_{-\infty}^{\infty} \frac{\Im[\chi^{(3)}(\omega')]}{\omega'-\omega}\,d\omega',$$

with numerical zero-padding and tapered windows. The normalized error

$$K(\omega)=\frac{|\Re[\chi^{(3)}(\omega)]-\Re[\chi^{(3)}_{\rm KK}(\omega)]|} {\sqrt{\Re^2+\Im^2}}$$

remains $\lesssim5\%$ in the central spectral region, confirming analyticity and causality in the resulting nonlinear spectra.

6. Quantitative Spectral Predictions

Numerically applying the full model to typical parameters ($R=3.0$ nm, $t=0.8$ nm, $\varepsilon_h=2.1$, $\varepsilon_i=6.0$, $\phi=0.15$, $\gamma_{10}=20$ meV, $\gamma_{20}=30$ meV) yields the following effective $\chi^{(3)}_{\rm eff}(\lambda)$:

$\lambda$ (nm)	$\Re[\chi^{(3)}]\times10^{22}$ (m$^2$/V$^2$)	$\Im[\chi^{(3)}]\times10^{22}$ (m$^2$/V$^2$)	$|\chi^{(3)}|\times10^{22}$ (m$^2$/V$^2$)
900	$-$0.0831	0.0034	0.0832
1000	$-$0.144	0.0179	0.1441
1100	$-$0.203	0.0667	0.2137
1200	$-$0.236	0.160	0.283
1300	$-$0.235	0.301	0.389
1400	$-$0.208	0.483	0.525

The nonlinear response shows a pronounced $|\chi^{(3)}|$ peak around $1200$ nm, with the peak shifting to shorter wavelengths as $R$ decreases. Increased shell thickness $t$ narrows the resonance and enhances amplitude; varying $\varepsilon_h$ and $\phi$ modulates the response up to an order of magnitude.

7. Photonic Application Guidelines and Engineering Implications

Designing CdSe/ZnS–MOF nano-composites for maximal third-order nonlinearity is guided by parameter-transparent scaling relations:

• Select $R$ to target the operational wavelength (e.g., NIR, telecom).
• Use a ZnS shell $t\gtrsim 0.8$ nm to enhance $S$ and suppress $\gamma$.
• Employ high-$\varepsilon_h$ MOFs with moderate $\phi$ for large local-field enhancement and connectivity without film percolation.
• Minimize size polydispersity to preserve narrow resonances.
• Confirm Kramers–Kronig self-consistency to rule out nonphysical artifacts in the computed spectra.

This formalism enables direct, quantitative prediction of macroscopic observables ($\chi^{(3)}, n_2, \beta$) for engineering optical switching, modulation, and wavelength conversion in hybrid quantum-plasmonic and photonic platforms. The methodology provides analytical access and design rules for the synthesis and deployment of QD–MOF nanocomposites with tailored nonlinear optical properties, connecting nanoscale structure to device-relevant figures of merit.