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 $\chi^{(3)}$0 and radial envelope functions $\chi^{(3)}$1. The lowest interband transition (confined exciton) energy, generalizing the Brus formula, includes finite barrier and dielectric mismatch effects: $\chi^{(3)}$2 where $\chi^{(3)}$3, $\chi^{(3)}$4 is the core/shell average permittivity, $\chi^{(3)}$5 is the MOF host permittivity, and $\chi^{(3)}$6–$\chi^{(3)}$7. The dipole matrix element for the exciton transition is

$\chi^{(3)}$8

with $\chi^{(3)}$9 the dimensionless envelope overlap. The two-level saturation intensity is

$R$0

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 $R$1 (ground, exciton, biexciton), with dipole coupling and coherent light-matter interaction. The system Hamiltonian under a time-dependent electric field is: $R$2 The Liouville–von Neumann equation (with dephasing $R$3 and relaxation $R$4 via superoperator $R$5) is expanded order-by-order in the driving field: $R$6 At steady-state, the third-order polarization at the fundamental frequency is

$R$7

A closed-form expression for $R$8 in the degenerate Kerr configuration is

$R$9

where $t$0 is the QD density, the denominators $t$1 encode detuning and dephasing, and the sum over permutations enforces causality. Near resonance, the leading behaviors are Lorentzian: $t$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 $t$3, intrinsic $t$4) in a MOF host (permittivity $t$5, susceptibility $t$6), at volume fraction $t$7. The Maxwell–Garnett (MG) and Bruggeman (Br) mixing formulas are invoked. For the MG case,

$t$8

with cubic nonlinear susceptibility

$t$9

The Bruggeman formula gives a self-consistent relation: $R_{\rm tot}=R+t$0 with local-field factors

$R_{\rm tot}=R+t$1

yielding

$R_{\rm tot}=R+t$2

The choice of homogenization model determines sensitivity to $R_{\rm tot}=R+t$3 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 $R_{\rm tot}=R+t$4 emerge from the interplay between QD quantum structure, local-field factors, and composite geometry: $R_{\rm tot}=R+t$5 Critical qualitative trends are:

• Shell thickness ($R_{\rm tot}=R+t$6): Increasing $R_{\rm tot}=R+t$7 enhances envelope-state overlap $R_{\rm tot}=R+t$8 and lowers dephasing rate $R_{\rm tot}=R+t$9, boosting $$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$0.
• Core radius ($$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$1): Decreasing $$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$2 blue-shifts the exciton resonance ($$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$3) and modifies the absorption/scattering features.
• Host permittivity ($$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$4) and fill fraction ($$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$5): Higher $$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$6 and $$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$7 increase $$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$8 and linearly scale the overall $$-\frac{\hbar^2}{2\,m^*(r)}\nabla^2\psi(\mathbf r)+V(r)\psi(\mathbf r) =E\,\psi(\mathbf r),$$9, with limits set by percolation.

Optimization for large third-order nonlinearity thus involves:

• Choosing $$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}$$0–$$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}$$1 nm to position the two-photon resonance in the NIR regime.
• Applying a ZnS shell with $$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}$$2 nm to maximize $$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}$$3 and minimize $$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}$$4.
• Selecting a MOF with $$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}$$5 and moderate $$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}$$6 ($$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}$$7–$$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}$$8).
• Minimizing polydispersity to narrow the resonance (limiting inhomogeneous Voigt broadening).

5. Validation of Analyticity and Causality

The physicality of the computed $$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}$$9 is validated via a Kramers–Kronig (KK) consistency check. The real part of $$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}$$0 is reconstructed from its imaginary part by

$$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}$$1

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

$$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}$$2

remains $$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}$$3 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 ($$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}$$4 nm, $$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}$$5 nm, $$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}$$6, $$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}$$7, $$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}$$8, $$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}$$9 meV, $\psi$0 meV) yields the following effective $\psi$1:

$\psi$2 (nm)	$\psi$3 (m$\psi$4/V$\psi$5)	$\psi$6 (m$\psi$7/V$\psi$8)	$\psi$9 (m$r = R$0/V$r = R$1)
900	$r = R$20.0831	0.0034	0.0832
1000	$r = R$30.144	0.0179	0.1441
1100	$r = R$40.203	0.0667	0.2137
1200	$r = R$50.236	0.160	0.283
1300	$r = R$60.235	0.301	0.389
1400	$r = R$70.208	0.483	0.525

The nonlinear response shows a pronounced $r = R$8 peak around $r = R$9 nm, with the peak shifting to shorter wavelengths as $\chi^{(3)}$00 decreases. Increased shell thickness $\chi^{(3)}$01 narrows the resonance and enhances amplitude; varying $\chi^{(3)}$02 and $\chi^{(3)}$03 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 $\chi^{(3)}$04 to target the operational wavelength (e.g., NIR, telecom).
• Use a ZnS shell $\chi^{(3)}$05 nm to enhance $\chi^{(3)}$06 and suppress $\chi^{(3)}$07.
• Employ high-$\chi^{(3)}$08 MOFs with moderate $\chi^{(3)}$09 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)}$10) 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.