We derive a nonsymmetrized 8-band effective-mass Hamiltonian for quantumdot heterostructures (QDHs) in Burt's envelope-function representation. The 8×8 radial Hamiltonian and the boundary conditions for the Schrödinger equation are obtained for spherical QDHs. Boundary conditions for symmetrized and nonsymmetrized radial Hamiltonians are compared with each other and with connection rules that are commonly used to match the wave functions found from the bulk k ·p Hamiltonians of two adjacent materials. Electron and hole energy spectra in three spherical QDHs: HgS/CdS, InAs/GaAs, and GaAs/AlAs are calculated as a function of the quantum dot radius within the approximate symmetrized and exact nonsymmetrized 8×8 models. The parameters of dissymmetry are shown to influence the energy levels and the wave functions of an electron and a hole and, consequently, the energies of both intraband and interband transitions.PACS numbers: 73.20. Dx, 73.40.Kp, 73.40.Lq This model explicitly includes eight bands around the Γ point of the Brillouin zone, namely, electron, heavy-, light-, and spin-orbit split-off hole bands (each of them is twice-degenerate due to the spin), and treats all other bands as remote. Along with more simple models, the 8×8 k ·p Hamiltonian has been used to investigate different QDs (see, e. g. .Recently, one has begun to apply multiband effective-mass Hamiltonians to investigate elastic, electronic, and optical properties of multilayer nanostructures such as quantumdot heterostructures (QDHs): CdS/HgS 11 , InAs/GaAs 12,13 , GaAs/Al x Ga 1−x As 14,15 , and CdS/HgS/CdS/H 2 O 16,17 . However, it should be emphasized, that multiband k ·p Hamiltonians are derived for homogeneous bulk materials, i.e. under the assumption that all effective-mass parameters are constant. This is important, because at a certain step of the derivation, wavenumbers k are declared as operatorsp/h that do not commute with the functions of coordinates. But, at the heterointerfaces of the multilayer nanostructures, there occurs an abrupt change of effective-mass parameters from their values in one material to those in the adjacent material. Inside a thin transitional layer that contains the heterointerface, the ordering of the differential operators and coordinate-dependent effective-mass parameters in the multiband Hamiltonian becomes crucial. In QDs with an infinitely high confining potential for electrons and holes, all components of the wave function vanish at the heterointerface, and there remains a possibility of applying the bulk multiband k ·p Hamiltonian straightforwardly. [3][4][5][6][8][9][10][11] There are two ways to proceed from QDs to QDHs.(i) The first way is to use an appropriate bulk multiband Hamiltonian for each constituent material separately, and then to match the obtained homogeneous solutions at the abrupt heterojunctions applying the connection rules (CRs) that are usually obtained by imposing the continuity of the wave function envelopes and of the normal to the heterointerface component of the velocity. 11,16...