Abstract: Piezo- and flexoelectricity are manifestations of electromechanical coupling in solids with potential applications in nanoscale materials. Naumov etal. [Phys. Rev. Lett. 102, 217601 (2009)] have shown by first principles calculations that a monolayer BN sheet becomes macroscopically polarized in-plane when in a corrugated state. Here, we investigate the interplay of layer corrugation and in-plane polarization by atomistic lattice dynamics. We treat the coupling between static flexural modes and in-plane atomic ion displacements as an anharmonic effect, similar to the membrane effect that is at the origin of negative thermal expansion in layered crystals. We have derived analytical expressions for the corrugation-induced static in-plane strains and the optical displacements with the resulting polarization response functions. Beyond h-BN, the theory applies to transition metal dichalcogenides and dioxides. Numerical calculations show that the effects are considerably stronger for 2D h-BN than for 2H-MoS2.

Abstract: The continuous miniaturization of nanodevices, such as transistors, solar cells, and optical fibers, requires the controlled synthesis of (ultra)thin gate oxides (<10 nm), including Si gate-oxide (SiO2) with high quality at the atomic scale. Traditional thermal growth of SiO2 on planar Si surfaces, however, does not allow one to obtain such ultrathin oxide due to either the high oxygen diffusivity at high temperature or the very low sticking ability of incident oxygen at low temperature. Two recent techniques, both operative at low (room) temperature, have been put forward to overcome these obstacles: (i) hyperthermal oxidation of planar Si surfaces and (ii) thermal or plasma-assisted oxidation of nonplanar Si surfaces, including Si nanowires (SiNWs). These nanooxidation processes are, however, often difficult to study experimentally, due to the key intermediate processes taking place on the nanosecond time scale.
In this Account, these Si nano-oxidation techniques are discussed from a computational point of view and compared to both hyperthermal and thermal oxidation experiments, as well as to well-known models of thermal oxidation, including the Deal−Grove, Cabrera−Mott, and Kao models and several alternative mechanisms. In our studies, we use reactive molecular dynamics (MD) and hybrid MD/Monte Carlo simulation techniques, applying the Reax force field. The incident energy of oxygen species is chosen in the range of 1−5 eV in hyperthermal oxidation of planar Si surfaces in order to prevent energy-induced damage. It turns out that hyperthermal growth allows for two growth modes, where the ultrathin oxide thickness depends on either (1) only the kinetic energy of the incident oxygen species at a growth temperature below Ttrans = 600 K, or (2) both the incident energy and the growth temperature at a growth temperature above Ttrans. These modes are specific to such ultrathin oxides, and are not observed in traditional thermal oxidation, nor theoretically considered by already existing models. In the case of thermal or plasma-assisted oxidation of small Si nanowires, on the other hand, the thickness of the ultrathin oxide is a function of the growth temperature and the nanowire diameter. Below Ttrans, which varies with the nanowire diameter, partially oxidized SiNW are formed, whereas complete oxidation to a SiO2 nanowire occurs only above Ttrans. In both nano-oxidation processes at lower temperature (T < Ttrans), final sandwich c-Si|SiOx|a-SiO2 structures are obtained due to a competition between overcoming the energy barrier to penetrate into Si subsurface layers and the compressive stress (∼2−3 GPa) at the Si crystal/oxide interface. The overall atomic-simulation results strongly indicate that the thickness of the intermediate SiOx (x < 2) region is very limited (∼0.5 nm) and constant irrespective of oxidation parameters. Thus, control over the ultrathin SiO2 thickness with good quality is indeed possible by accurately tuning the oxidant energy, oxidation temperature and surface curvature.
In general, we discuss and put in perspective these two oxidation mechanisms for obtaining controllable ultrathin gate-oxide films, offering a new route toward the fabrication of nanodevices via selective nano-oxidation.

Abstract: Utilization of Au and nanocrystalline diamond ( NCD) as interlayers noticeably modifies the microstructure and field electron emission ( FEE) properties of hexagonal boron nitride nanowalls ( hBNNWs) grown on Si substrates. The FEE properties of hBNNWs on Au could be turned on at a low turn-on field of 14.3V mu m(-1), attaining FEE current density of 2.58mAcm(-2) and life-time stability of 105 min. Transmission electron microscopy reveals that the Au-interlayer nucleates the hBN directly, preventing the formation of amorphous boron nitride ( aBN) in the interface, resulting in enhanced FEE properties. But Au forms as droplets on the Si substrate forming again aBN at the interface. Conversely, hBNNWs on NCD shows superior in life-time stability of 287 min although it possesses inferior FEE properties in terms of larger turn-on field and lower FEE current density as compared to that of hBNNWs-Au. The uniform and continuous NCD film on Si also circumvents the formation of aBN phases and allows hBN to grow directly on NCD. Incorporation of carbon in hBNNWs from the NCD-interlayer improves the conductivity of hBNNWs, which assists in transporting the electrons efficiently from NCD to hBNNWs that results in better field emission of electrons with high life-time stability. (C) 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Abstract: To unravel the effects of phonon confinement, the influence of size and morphology on the atomic vibrations is investigated in Sn nanoislands and cluster-assembled films. Nuclear resonant inelastic x-ray scattering is used to probe the phonon densities of states of the Sn nanostructures which show significant broadening of the features compared to bulk phonon behavior. Supported by ab initio calculations, the broadening is attributed to phonon scattering and can be described within the damped harmonic oscillator model. Contrary to the expectations based on previous research, the appearance of high-energy modes above the cutoff energy is not observed. From the thermodynamic properties extracted from the phonon densities of states, it was found that grain boundary Sn atoms are bound by weaker forces than bulk Sn atoms.

Abstract: Systems which rapidly evolve through symmetry-breaking transitions on timescales comparable to the fluctuation timescale of the single-particle excitations may behave very differently than under controlled near-ergodic conditions. A real-time investigation with high temporal resolution may reveal insights into the ordering through the transition that are not available in static experiments. We present an investigation of the system trajectory through a normal-to-superconductor transition in a prototype high-temperature superconducting cuprate in which such a situation occurs. Using a multiple pulse femtosecond spectroscopy technique we measure the system trajectory and time evolution of the single-particle excitations through the transition in La1.9Sr0.1CuO4 and compare the data to a simulation based on the time-dependent Ginzburg-Landau theory, using the laser excitation fluence as an adjustable parameter controlling the quench conditions in both experiment and theory. The comparison reveals the presence of significant superconducting fluctuations which precede the transition on short timescales. By including superconducting fluctuations as a seed for the growth of the superconducting order we can obtain a satisfactory agreement of the theory with the experiment. Remarkably, the pseudogap excitations apparently play no role in this process.

Abstract: To gain insight into the nature of the adhesion mechanism between hydroxyapatite (HA) and rutile (rTiO(2)), the mutual affinity between their surfaces was systematically studied using density functional theory (DFT). We calculated both bulk and surface properties of HA and rTiO(2), and explored the interfacial bonding mechanism of amorphous HA (aHA) surface onto amorphous as well as stoichiometric and nonstoichiometric crystalline rTiO(2). Formation energies of bridging and subbridging oxygen vacancies considered in the rTiO(2)(110) surface were evaluated and compared with other theoretical and experimental results. The interfacial interaction was evaluated through the work of adhesion. For the aHA/rTiO(2)(110) interfaces, the work of adhesion is found to depend strongly on the chemical environment of the rTiO(2)(110) surface. Electronic analysis indicates that the charge transfer is very small in the case of interface formation between aHA and crystalline rTiO(2)(110). In contrast, significant charge transfer occurs between aHA and amorphous rTiO(2) (aTiO(2)) slabs during the formation of the interface. Charge density difference (CDD) analysis indicates that the dominant interactions in the interface have significant covalent character, and in particular the Ti-O and Ca-O bonds. Thus, the obtained results reveal that the aHA/aTiO(2) interface shows a more preferable interaction and is thermodynamically more stable than other interfaces. These results are particularly important for improving the long-term stability of HA-based implants.

Abstract: The Wigner-Liouville equation is reformulated using a spectral decomposition of the classical force field instead of the potential energy. The latter is shown to simplify the Wigner-Liouville kernel both conceptually and numerically as the spectral force Wigner-Liouville equation avoids the numerical evaluation of the highly oscillatory Wigner kernel which is nonlocal in both position and momentum. The quantum mechanical evolution is instead governed by a term local in space and non-local in momentum, where the non locality in momentum has only a limited range. An interpretation of the time evolution in terms of two processes is presented; a classical evolution under the influence of the averaged driving field, and a probability-preserving quantum-mechanical generation and annihilation term. Using the inherent stability and reduced complexity, a direct deterministic numerical implementation using Chebyshev and Fourier pseudo-spectral methods is detailed. For the purpose of illustration, we present results for the time evolution of a one-dimensional resonant tunneling diode driven out of equilibrium. (C) 2017 Elsevier Inc. All rights reserved.

Abstract: We study the oxygen thermal stability of two binary
systems. The larger particles are magnetic amorphous Co (7.2 nm) or
Fe3O4 (7.5 nm) nanocrystals, whereas the smaller ones (3.7 nm) are
Au nanocrystals. The nanocrystal ordering as well as the choice of the
magnetic nanoparticles very much influence the stability of the binary
system. A perfect crystalline structure is obtained with the Fe3O4/Au
binary supracrystals. For the Co/Au binary system, oxidation of Co
results in the chemical transformation from Co to CoO, where the size
of the amorphous Co nanoparticles increases from 7.2 to 9.8 nm in
diameter. During the volume expansion of the Co nanoparticles, Au
nanoparticles within the binary assemblies coalesce and are at the
origin of the instability of the binary nanoparticle supracrystals. On the
other hand, for the Fe3O4/Au binary system, the oxidation of Fe3O4 to
γ-Fe2O3 does not lead to a size change of the nanoparticles, which
maintains the stability of the binary nanoparticle supracrystals. A similar behavior is observed for an AlB2-type Co−Ag binary
system: The crystalline structure is maintained, whereas in disordered assemblies, coalescence of Ag nanocrystals is observed.

Abstract: Epitaxially grown ultra-thin Si layers are often used to passivate Ge surfaces in the high-k gate module of (strained) Ge FinFET and Gate All Around devices. We use Si4H10 as Si precursor as it enables epitaxial Si growth at temperatures down to 330 degrees. C-V characteristics of blanket capacitors made on Ge virtual substrates point to the presence of an optimal Si thickness. In case of compressively strained Ge fin structures, the Si growth results in non-uniform and high strain levels in the strained Ge fin. These strain levels have been calculated for different shapes of the Ge fin and in function of the grown Si thickness. The high strain is the driving force for potential (unwanted) Ge surface reflow during Si deposition. The Ge surface reflow is strongly affected by the strength of the H-passivation during Si-capping and can be avoided by carefully selected process conditions. (C) The Author(s) 2018. Published by ECS.

Abstract: Discrepancies exist between the theoretically predicted and experimentally measured performance of band-to-band tunneling devices, such as Esaki diodes and tunnel field-effect transistors (TFETs). We resolve this discrepancy for highly-doped, direct-bandgap Esaki diodes by successfully calibrating a semi-classical model for high-doping-induced ballistic band-tails tunneling currents at multiple temperatures with two In0.53Ga0.47As Esaki diodes using their SIMS doping profiles, C-V characteristics and their forward-bias current density in the negative differential resistance (NDR) regime. The current swing in the NDR regime is shown not to be linked to the band-tails Urbach energy. We further demonstrate theoretically that the calibrated band-tails contribution is also the dominant band-tails contribution to the subthreshold swing of the corresponding TFETs. Lastly, we verify that the presented procedure is applicable to all direct-bandgap semiconductors by successfully applying it to InAs Esaki diodes in literature.

Abstract: Supersonic flows provide a high thermodynamic
nonequilibrium, which is crucial for energy-efficient conversion of
CO 2 in microwave plasmas and are therefore of great interest.
However, the effect of the flow on the chemical reactions is poorly
understood. In this work, we present a combined flow and plasma
chemical kinetics model of a microwave CO 2 plasma in a Laval
nozzle setup. The effects of the flow field on the different dissociation
and recombination mechanisms, the vibrational distribution, and the
vibrational transfer mechanism are discussed. In addition, the effect
of experimental parameters, like position of power deposition, outlet
pressure, and specific energy input, on the CO 2 conversion and
energy efficiency is examined. The short residence time of the gas in
the plasma region, the shockwave, and the maximum critical heat,
and thus power, that can be added to the flow to avoid thermal
choking are the main obstacles to reaching high energy efficiencies.

Abstract: When slow ions incident on a surface are neutralized, the excess potential energy is passed on to an electron inside the surface, leading to emission of secondary electrons. The microscopic description of this process, as
well as the calculation of the secondary electron yield, is a challenging problem due to its complexity as well
as its sensitivity to surface properties. One of the first quantitative descriptions was articulated in the 1950s by
Hagstrum, who based his calculation on a parametrization of the density of states of the material. In this paper, we
present a model for calculating the secondary electron yield, derived from Hagstrum’s initial approach. We use
first-principles density functional theory calculations to acquire the necessary input and introduce the concept of
electron cascades to Hagstrum’s model in order to improve the calculated spectra, as well as remove its reliance
on fitting parameters. We apply our model to He+ and Ne+ ions incident on Ge(111) and Si(111) and obtain
yield spectra that match closely to the experimental results of Hagstrum.

Abstract: The deployment of Gen-IV lead-cooled fast reactors requires a good compatibility between the selected structural/cladding steels and the inherently corrosive heavy liquid metal coolant. An effective liquid metal corrosion mitigation strategy involves the in-situ steel passivation in contact with the oxygen-containing Pb-alloy coolant. Transmission electron microscopy was used in this work to study the multi-layered oxide scales forming on an austenitic stainless steel fuel cladding exposed to oxygen-containing (CO ≈ 10−6 mass%) static liquid leadbismuth eutectic (LBE) for 1000 h between 400 and 500 °C. The oxide scale constituents were analyzed, including the intertwined phases comprising the innermost biphasic layer.

Abstract: Aquaporins (AQPs) are transmembrane proteins that conduct not only water molecules across the cell membrane but also other solutes, such as reactive oxygen and nitrogen species (RONS), produced (among others) by cold atmospheric plasma (CAP). These RONS may induce oxidative stress in the cell interior, which plays a role in cancer treatment. The underlying mechanisms of the transport of RONS across AQPs, however, still remain obscure. We apply molecular dynamics simulations to investigate the permeation of both hydrophilic (H₂O₂and OH) and hydrophobic (NO₂and NO) RONS through AQP1. Our simulations show that these RONS can all penetrate across the pores of AQP1. The permeation free energy barrier of OH and NO is lower than that of H₂O₂and NO₂, indicating that these radicals may have easier access to the pore interior and interact with the amino acid residues of AQP1. We also study the effect of RONS-induced oxidation of both the phospholipids and AQP1 (i.e., sulfenylation of Cys₁₉₁) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H₂O₂and NO₂permeation, while for OH and NO, we do not observe a strong effect of oxidation. The simulation results help to gain insight in the underlying mechanisms of the noticeable rise of CAP-induced RONS in cancer cells, thereby improving our understanding on the role of AQPs in the selective anticancer capacity of CAP.

Abstract: In van der Waals heterostructures consisting of graphene and a substrate, lattice mismatch often leads to a moire pattern with a huge supercell, preventing its treatment within first- principles calculations. Previous theoretical works considered mostly simple stacking models such as AB, AA with straining the lattice of graphene to match that of the substrate. Here, we propose a moire superlattice build from graphene and porous graphene or graphyne like monolayers, having a lower interlayer binding energy, needing little strain in order to match the lattices. In contrast to the results from the simple stacking models, the present ab initio calculations for the moire superlattices show different properties in lattice structure, energy, and band structures. For example, the Dirac cone at the K point is preserved and a linear energy dispersion near the Fermi level is obtained.