Abstract: A reactive ReaxFF force field is developed for studying the influence of intrinsic point defects on the chemistry with TiO2 condensed phases. The force field parameters are optimized to ab initio data for the equations of state, relative phase stabilities for titanium and titanium dioxide, potential energy differences for (TiO2)n-clusters (n = 116). Also data for intrinsic point defects in anatase were added. These data contain formation energies for interstitial titanium and oxygen vacancies, diffusion barriers of the oxygen vacancies and molecular oxygen adsorption on a reduced anatase (101) surface. Employing the resulting force field, we study the influence of concentration of oxygen vacancies and expansion or compression of an anatase surface on the diffusion of the oxygen vacancies. Also the barrier for oxygen diffusion in the subsurface region is evaluated using this force field. This diffusion barrier of 27.7 kcal/mol indicates that the lateral redistribution of oxygen vacancies on the surface and in the subsurface will be dominated by their diffusion in the subsurface, since both this barrier as well as the barriers for diffusion from the surface to the subsurface and vice versa (17.07 kcal/mol and 21.91 kcal/mol, respectively, as calculated with DFT), are significantly lower than for diffusion on the surface (61.12 kcal/mol as calculated with DFT).

Abstract: A marker function designed to indicate in which regions of a generic flow field the results from linear eddy-viscosity turbulence models are plausibly inaccurate is introduced. The marker is defined to identify regions that deviate from parallel shear flow. For two different flow fields it is shown that these regions largely coincide with regions where the prediction of the Reynolds stress divergence is inaccurate. The marker therefore offers a guideline for interpreting results obtained from Reynolds-averaged Navier-Stokes simulations and provides a basis for the further development of turbulence model-form uncertainty quantification methods. (C) 2014 AIP Publishing LLC.

Abstract: The conversion of methane to value-added chemicals and fuels is considered to be one of the challenges of the 21st century. In this paper we study, by means of fluid modeling, the conversion of methane to higher hydrocarbons or oxygenates by partial oxidation with CO2 or O2 in a dielectric barrier discharge. Sixty-nine different plasma species (electrons, ions, molecules, radicals) are included in the model, as well as a comprehensive set of chemical reactions. The calculation results presented in this paper include the conversion of the reactants and the yields of the reaction products as a function of residence time in the reactor, for different gas mixing ratios. Syngas (i.e. H2 + CO) and higher hydrocarbons (C2Hx) are typically found to be important reaction products.

Abstract: Coupled shallow impurity states in a freestanding semiconductor nanowire and in a semiconductor nanowire surrounded by a metallic gate are studied within the effective-mass approximation. Bonding and antibonding states are found due to the coupling of the two impurities, and their energy converges with increasing distance di between the two impurities. The dependences of the binding energy on the wire radius R, the distance di between the two impurities, and the impurity radial position in the nanowire are examined.

Abstract: Using the Ginzburg-Landau theory extended to the next-to-leading order, we determine numerically the healing lengths of the two order parameters at the two-gap superconductor/normal metal interface. We demonstrate on several examples that those can be different even in the strict domain of applicability of the Ginzburg-Landau theory. This justifies the use of this theory to describe relevant physics of two-gap superconductors, distinguishing them from their single-gap counterparts. The calculational degree of complexity increases only slightly with respect to the conventional Ginzburg-Landau expansion, thus the extended Ginzburg-Landau model remains numerically far less demanding compared to the full microscopic approaches.

Abstract: The diffusion of charged particles interacting through a repulsive Yukawa potential, exp(-r/lambda)/r, confined by a parabolic potential in the y direction and subjected to a periodic substrate potential in the x direction is investigated. Langevin dynamic simulations are used to investigate the effect of the particle density, the amplitude of the periodic substrate, and the range of the interparticle interaction potential on the diffusive behavior of the particles. We found that in general the diffusion is suppressed with increasing the amplitude of the periodic potential, but for specific values of the strength of the substrate potential a remarkable increase of the diffusion is found with increasing the periodic potential amplitude. In addition, we found a strong dependence of the diffusion on the specific arrangement of the particles, e. g., single-chain versus multichain configuration. For certain particle configurations, a reentrant behavior of the diffusion is found as a function of the substrate strength due to structural transitions in the ordering of the particles.

Abstract: We expand on a recent study of a lattice model of interacting particles [Phys. Rev. Lett. 111, 110601 (2013)]. The adsorption isotherm and equilibrium fluctuations in particle number are discussed as a function of the interaction. Their behavior is similar to that of interacting particles in porous materials. Different expressions for the particle jump rates are derived from transition-state theory. Which expression should be used depends on the strength of the interparticle interactions. Analytical expressions for the self-and transport diffusion are derived when correlations, caused by memory effects in the environment, are neglected. The diffusive behavior is studied numerically with kinetic Monte Carlo (kMC) simulations, which reproduces the diffusion including correlations. The effect of correlations is studied by comparing the analytical expressions with the kMC simulations. It is found that the Maxwell-Stefan diffusion can exceed the self-diffusion. To our knowledge, this is the first time this is observed. The diffusive behavior in one-dimensional and higher-dimensional systems is qualitatively the same, with the effect of correlations decreasing for increasing dimension. The length dependence of both the self-and transport diffusion is studied for one-dimensional systems. For long lengths the self-diffusion shows a 1/L dependence. Finally, we discuss when agreement with experiments and simulations can be expected. The assumption that particles in different cavities do not interact is expected to hold quantitatively at low and medium particle concentrations if the particles are not strongly interacting.

Abstract: The chemical complexity of non-equilibrium plasmas poses a challenge for plasma modeling because of the computational load. This paper presents a dimension reduction method for such chemically complex plasmas based on principal component analysis (PCA). PCA is used to identify a low-dimensional manifold in chemical state space that is described by a small number of parameters: the principal components. Reduction is obtained since continuity equations only need to be solved for these principal components and not for all the species. Application of the presented method to a CO2 plasma model including state-to-state vibrational kinetics of CO2 and CO demonstrates the potential of the PCA method for dimension reduction. A manifold described by only two principal components is able to predict the CO2 to CO conversion at varying ionization degrees very accurately.

Abstract: Germanium is a widely used material for tunnel FETs because of its small band gap and compatibility with silicon. Typically, only the indirect band gap of Ge at 0.66 eV is considered. However, direct band-to-band tunneling (BTBT) in Ge should be included in tunnel FET modeling and simulations since the energy difference between the Ge conduction band edges at the L and G valleys is only 0.14 eV at room temperature. In this paper, we theoretically calculate the parameters A and B of Kane's direct and indirect BTBT models at different tunneling directions ([100], [110], and [111]) for Si, Ge and unstrained Si1-xGex. We highlight how the direct BTBT component becomes more important as the Ge mole fraction increases. The calculation of the band-to-band generation rate in the uniform electric field limit reveals that direct tunneling always dominates over indirect tunneling in Ge. The impact of the direct transition in Ge on the performance of two realistic tunnel field-effect transistor configurations is illustrated with TCAD simulations. The influence of field-induced quantum confinement is included in the analysis based on a back-of-the-envelope calculation.