Giacomo Miceli

The minimum energy paths for the migration of interstitial Mg in wurtzite GaN are studied through density functional calculations. The study also comprises Li, Na, and Be dopants to examine the dependence on size and charge of the dopant species. In all cases considered, the impurities diffuse like ions without any tendency of localizing charge. Li, Mg, and to some extent Na, diffuse almost isotropically in GaN, with average diffusion barriers of 1.1, 2.1, and 2.5 eV, respectively. Instead Be shows a marked anisotropy with energy barriers of 0.76 and 1.88 eV for diffusion paths perpendicular and parallel to the c-axis. The diffusion barrier generally increases with ionic charge and ionic radius, but their interplay is not trivial. The calculated migration barrier for Mg is consistent with the values estimated in a recent beta- emission channeling experiment.

We study the formation energies of native point defects in GaN through density-functional theory. In our first-principles scheme, the band edges are positioned in accord with hybrid density functional calculations, thus yielding a band-gap in agreement with experiment. With respect to previous semilocal calculations, the calculated formation energies and charge transition levels are found to be significantly different in quantitative terms, while the overall qualitative trend remains similar. In Ga-rich conditions, the nitrogen vacancy corresponds to the most stable defect for all Fermi energies in the band gap, but its formation energy is too high to account for autodoping. Our calculations also indicate that the gallium vacancy does not play any compensating role in n-type GaN.

We investigate the structural properties of liquid water at near ambient conditions using first-principles molecular dynamics simulations based on a semilocal density functional augmented with nonlocal van der Waals interactions. The adopted scheme offers the advantage of simulating liquid water at essentially the same computational cost of standard semilocal functionals. Applied to the water dimer and to ice Ih, we find that the hydrogen-bond energy is only slightly enhanced compared to a standard semilocal functional. We simulate liquid water through molecular dynamics in the NpH statistical ensemble allowing for fluctuations of the system density. The structure of the liquid departs from that found with a semilocal functional leading to more compact structural arrangements. This indicates that the directionality of the hydrogen-bond interaction has a diminished role as compared to the overall attractions, as expected when dispersion interactions are accounted for. This is substantiated through a detailed analysis comprising the study of the partial radial distribution functions, various local order indices, the hydrogen-bond network, and the selfdiffusion coefficient. The explicit treatment of the van der Waals interactions leads to an overall improved description of liquid water.

We determine and compare structural, dynamical, and electronic properties of liquid water at near ambient conditions through density-functional molecular dynamics simulations, when using either plane-wave or atomic-orbital basis sets. In both frameworks, the electronic structure and the atomic forces are self-consistently determined within the same theoretical scheme based on a nonlocal density functional accounting for van der Waals interactions. The overall properties of liquid water achieved within the two frameworks are in excellent agreement with each other. Thus, our study supports that implementations with plane-wave or atomic-orbital basis sets yield equivalent results and can be used indiscriminately in study of liquid water or aqueous solutions.