ASU Scholarship Showcase

Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm³ to 62 mm³, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.

Background:
Data assimilation refers to methods for updating the state vector (initial condition) of a complex spatiotemporal model (such as a numerical weather model) by combining new observations with one or more prior forecasts. We consider the potential feasibility of this approach for making short-term (60-day) forecasts of the growth and spread of a malignant brain cancer (glioblastoma multiforme) in individual patient cases, where the observations are synthetic magnetic resonance images of a hypothetical tumor.

Results:
We apply a modern state estimation algorithm (the Local Ensemble Transform Kalman Filter), previously developed for numerical weather prediction, to two different mathematical models of glioblastoma, taking into account likely errors in model parameters and measurement uncertainties in magnetic resonance imaging. The filter can accurately shadow the growth of a representative synthetic tumor for 360 days (six 60-day forecast/update cycles) in the presence of a moderate degree of systematic model error and measurement noise.

Conclusions:
The mathematical methodology described here may prove useful for other modeling efforts in biology and oncology. An accurate forecast system for glioblastoma may prove useful in clinical settings for treatment planning and patient counseling.

The development of non-volatile logic through direct coupling of spontaneous ferroelectric polarization with semiconductor charge carriers is nontrivial, with many issues, including epitaxial ferroelectric growth, demonstration of ferroelectric switching and measurable semiconductor modulation. Here we report a true ferroelectric field effect—carrier density modulation in an underlying Ge(001) substrate by switching of the ferroelectric polarization in epitaxial c-axis-oriented BaTiO₃ grown by molecular beam epitaxy. Using the density functional theory, we demonstrate that switching of BaTiO₃ polarization results in a large electric potential change in Ge. Aberration-corrected electron microscopy confirms BaTiO₃ tetragonality and the absence of any low-permittivity interlayer at the interface with Ge. The non-volatile, switchable nature of the single-domain out-of-plane ferroelectric polarization of BaTiO₃ is confirmed using piezoelectric force microscopy. The effect of the polarization switching on the conductivity of the underlying Ge is measured using microwave impedance microscopy, clearly demonstrating a ferroelectric field effect.

The emission properties of GeSn heterostructure pin diodes have been investigated. The devices contain thick (400–600 nm) Ge_1-y Sn_y i-layers spanning a broad compositional range below and above the crossover Sn concentration y_c where the Ge_1-y Sn_y alloy becomes a direct-gap material. These results are made possible by an optimized device architecture containing a single defected interface thereby mitigating the deleterious effects of mismatch-induced defects. The observed emission intensities as a function of composition show the contributions from two separate trends: an increase in direct gap emission as the Sn concentration is increased, as expected from the reduction and eventual reversal of the separation between the direct and indirect edges, and a parallel increase in non-radiative recombination when the mismatch strains between the structure components is partially relaxed by the generation of misfit dislocations. An estimation of recombination times based on the observed electroluminescence intensities is found to be strongly correlated with the reverse-bias dark current measured in the same devices.

Predicting the timing of a castrate resistant prostate cancer is critical to lowering medical costs and improving the quality of life of advanced prostate cancer patients. We formulate, compare and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). We accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). While these models are simplifications of a previously published model, they fit data with similar accuracy and improve forecasting results. Both models describe the progression of androgen resistance. Although Model 1 is simpler than the more realistic Model 2, it can fit clinical data to a greater precision. However, we found that Model 2 can forecast future PSA levels more accurately. These findings suggest that including more realistic mechanisms of androgen dynamics in a two population model may help androgen resistance timing prediction.

The compositional dependence of the lowest direct and indirect band gaps in Ge_1-ySn_y alloys has been determined from room-temperature photoluminescence measurements. This technique is particularly attractive for a comparison of the two transitions because distinct features in the spectra can be associated with the direct and indirect gaps. However, detailed modeling of these room temperature spectra is required to extract the band gap values with the high accuracy required to determine the Sn concentration y_c at which the alloy becomes a direct gap semiconductor. For the direct gap, this is accomplished using a microscopic model that allows the determination of direct gap energies with meV accuracy. For the indirect gap, it is shown that current theoretical models are inadequate to describe the emission properties of systems with close indirect and direct transitions. Accordingly, an ad hoc procedure is used to extract the indirect gap energies from the data. For y < 0.1 the resulting direct gap compositional dependence is given by ΔE₀ = −(3.57 ± 0.06)y (in eV). For the indirect gap, the corresponding expression is ΔE_ind = −(1.64 ± 0.10)y (in eV). If a quadratic function of composition is used to express the two transition energies over the entire compositional range 0 ≤ y ≤ 1, the quadratic (bowing) coefficients are found to be b₀ = 2.46 ± 0.06 eV (for E0) and b_ind = 1.03 ± 0.11 eV (for E_ind). These results imply a crossover concentration y_c = $0.073 [+0.007 over -0.006], much lower than early theoretical predictions based on the virtual crystal approximation, but in better agreement with predictions based on large atomic supercells.

Earth-abundant sustainable inorganic thin-film solar cells, independent of precious elements, pivot on a marginal material phase space targeting specific compounds. Advanced materials characterization efforts are necessary to expose the roles of microstructure, chemistry, and interfaces. Herein, the earth-abundant solar cell device, Cu₂ZnSnS_(4-x)Se_x, is reported, which shows a high abundance of secondary phases compared to similarly grown Cu₂ZnSnSe₄.

The ongoing Zika virus (ZIKV) epidemic in the Americas poses a major global public health emergency. While ZIKV is transmitted from human to human by bites of Aedes mosquitoes, recent evidence indicates that ZIKV can also be transmitted via sexual contact with cases of sexually transmitted ZIKV reported in Argentina, Canada, Chile, France, Italy, New Zealand, Peru, Portugal, and the USA. Yet, the role of sexual transmission on the spread and control of ZIKV infection is not well-understood. We introduce a mathematical model to investigate the impact of mosquito-borne and sexual transmission on the spread and control of ZIKV and calibrate the model to ZIKV epidemic data from Brazil, Colombia, and El Salvador. Parameter estimates yielded a basic reproduction number R_0 = 2.055 (95% CI: 0.523–6.300), in which the percentage contribution of sexual transmission is 3.044% (95% CI: 0.123–45.73). Our sensitivity analyses indicate that R_0 is most sensitive to the biting rate and mortality rate of mosquitoes while sexual transmission increases the risk of infection and epidemic size and prolongs the outbreak. Prevention and control efforts against ZIKV should target both the mosquito-borne and sexual transmission routes.

A data set of observed daily precipitation, maximum and minimum temperature, gridded to a 1/16° (~6 km) resolution, is described that spans the entire country of Mexico, the conterminous U.S. (CONUS), and regions of Canada south of 53° N for the period 1950-2013. The dataset improves previous products in spatial extent, orographic precipitation adjustment over Mexico and parts of Canada, and reduction of transboundary discontinuities. The impacts of adjusting gridded precipitation for orographic effects are quantified by scaling precipitation to an elevation-aware 1981-2010 precipitation climatology in Mexico and Canada. Differences are evaluated in terms of total precipitation as well as by hydrologic quantities simulated with a land surface model. Overall, orographic correction impacts total precipitation by up to 50% in mountainous regions outside CONUS. Hydrologic fluxes show sensitivities of similar magnitude, with discharge more sensitive than evapotranspiration and soil moisture. Because of the consistent gridding methodology, the current product reduces transboundary discontinuities as compared with a commonly used reanalysis product, making it suitable for estimating large-scale hydrometeorologic phenomena.

We present two-dimensional Mg(OH)₂ sheets and their vertical heterojunctions with CVD-MoS₂ for the first time as flexible 2D insulators with anomalous lattice vibration and chemical and physical properties. New hydrothermal crystal growth technique enabled isolation of environmentally stable monolayer Mg(OH)₂ sheets. Raman spectroscopy and vibrational calculations reveal that the lattice vibrations of Mg(OH)₂ have fundamentally different signature peaks and dimensionality effects compared to other 2D material systems known to date. Sub-wavelength electron energy-loss spectroscopy measurements and theoretical calculations show that Mg(OH)₂ is a 6 eV direct-gap insulator in 2D, and its optical band gap displays strong band renormalization effects from monolayer to bulk, marking the first experimental confirmation of confinement effects in 2D insulators. Interestingly, 2D-Mg(OH)₂ sheets possess rather strong surface polarization (charge) effects which is in contrast to electrically neutral h-BN materials. Using 2D-Mg(OH)₂ sheets together with CVD-MoS₂ in the vertical stacking shows that a strong change transfer occurs from n-doped CVD-MoS₂ sheets to Mg(OH)₂, naturally depleting the semiconductor, pushing towards intrinsic doping limit and enhancing overall optical performance of 2D semiconductors. Results not only establish unusual confinement effects in 2D-Mg(OH)₂, but also offer novel 2D-insulating material with unique physical, vibrational, and chemical properties for potential applications in flexible optoelectronics.