  Purpose: To develop and apply multiscale Monte Carlo (MC) simulations to
assess variations in nucleus and cytoplasm dose enhancement factors (n,cDEFs)
over tumor-scale volumes. Methods: The intrinsic variation of n,cDEFs (due to
variations in local gold concentration and cell/nucleus size) are estimated via
MC modeling of varied cellular GNP uptake and cell/nucleus sizes. Then, a
Heterogeneous MultiScale (HetMS) model is implemented by combining detailed
models of populations of cells containing GNPs with simplified tissue models to
evaluate n,cDEFs. Simulations of tumors with spatially-uniform and -varying
gold concentrations are performed to determine n,cDEFs as a function of
distance from the source. All simulations are performed for three different GNP
configurations: GNPs distributed around the nucleus or packed into one or four
endosome(s). Results: Intrinsic variations in n,cDEFs can be substantial with
changes in GNP uptake and cell/nucleus radii. In HetMS models of tumors,
sub-unity n,cDEFs (dose decreases) can occur for low energies and high gold
concentrations due to attenuation of primary photons through the gold. In HetMS
simulations of tumors with spatially-uniform gold concentrations, n,cDEFs
decrease with depth into the tumor as photons are attenuated, with relative
differences between GNP models remaining approximately constant with depth in
the tumor. Similar initial n,cDEF decreases with radius are seen in tumors with
spatially-varying gold concentrations, but n,cDEFs for all configurations
converge to a single value for each energy as gold concentration reaches zero.
Conclusions: These results demonstrate that cellular doses are highly sensitive
to cell/nucleus size, GNP intracellular distribution, gold concentration, and
cell position in the tumor. This work demonstrates the importance of proper
choice of computational model when simulating GNPT scenarios.
