One of the most challenging and frequently arising problems in many areas of science is to find solutions of a system of multivariate nonlinear equations. There are several numerical methods that can find many (or all if the system is small enough) s
olutions but they each exhibit characteristic problems. Moreover, traditional methods can break down if the system contains singular solutions. Here, we propose an efficient implementation of Newton homotopies, which can sample a large number of the stationary points of complicated many-body potentials. We demonstrate how the procedure works by applying it to the nearest-neighbor $phi^4$ model and atomic clusters.
We study the dependence of the quantum yield of photoluminescence of a dense, periodic array of semiconductor nanocrystals (NCs) on the level of doping and NC size. Electrons introduced to NCs via doping quench photoluminescence by the Auger process,
so that practically only NCs without electrons contribute to the photoluminescence. Computer simulation and analytical theory are used to find a fraction of such empty NCs as a function of the average number of donors per NC and NC size. For an array of small spherical NCs, the quantization gap between 1S and 1P levels leads to transfer of electrons from NCs with large number of donors to those without donors. As a result, empty NCs become extinct, and photoluminescence is quenched abruptly at an average number of donors per NC close to 1.8. The relative intensity of photoluminescence is shown to correlate with the type of hopping conductivity of an array of NCs.
