Abstract: The interaction of monolayer graphene with specific substrates may break its sublattice symmetry and results in unidirectional chiral states with opposite group velocities in the different Dirac cones (Zarenia et al 2012 Phys. Rev. B 86 085451). Taking advantage of this feature, we propose a valley filter based on a transversal mass kink for low energy electrons in graphene, which is obtained by assuming a defect region in the substrate that provides a change in the sign of the substrate-induced mass and thus creates a non-biased channel, perpendicular to the kink, for electron motion. By solving the time-dependent Schrodinger equation for the tight-binding Hamiltonian, we investigate the time evolution of a Gaussian wave packet propagating through such a system and obtain the transport properties of this graphene-based substrate-induced quantum point contact. Our results demonstrate that efficient valley filtering can be obtained, provided: (i) the electron energy is sufficiently low, i.e. with electrons belonging mostly to the lowest sub-band of the channel, and (ii) the channel length (width) is sufficiently long (narrow). Moreover, even though the transmission probabilities for each valley are significantly affected by impurities and defects in the channel region, the valley polarization in this system is shown to be robust against their presence.

Abstract: Electrically gating of trilayer graphene (TLG) opens a band gap offering the possibility to electrically engineer TLG quantum dots. We study the energy levels of such quantum dots and investigate their dependence on a perpendicular magnetic field B and different types of stacking of the graphene layers. The dots are modeled as circular and confined by a truncated parabolic potential which can be realized by nanostructured gates or position-dependent doping. The energy spectra exhibit the intervalley symmetry E-K(e) (m) = -E (h)(K') (m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number and K and K' label the two valleys. The electron and hole spectra for B = 0 are twofold degenerate due to the intervalley symmetry E-K (m) = E-K' [-(m + 1)]. For both ABC [alpha = 1.5 (1.2) for large (small) R] and ABA (alpha = 1) stackings, the lowest-energy levels show approximately a R-alpha dependence on the dot radius R in contrast with the 1/R-3 one for ABC-stacked dots with infinite-mass boundary. As functions of the field B, the oscillator strengths for dipole-allowed transitions differ drastically for the two types of stackings.

Abstract: Microwave plasmas are one of the most promising techniques for CO2 conversion into value-added chemicals and fuels since they are very energy efficient. Nevertheless, experiments show that this high energy efficiency is only reached at low pressures and significantly drops toward atmospheric pressure, which is a clear limitation for industrial applications. In this paper, we use a zerodimensional reaction kinetics model to simulate a CO2 microwave plasma in a pressure range from 50 mbar to 1 bar, in order to evaluate the reasons for this decrease in energy efficiency at atmospheric pressure. The code includes a detailed description of the vibrational kinetics of CO2, CO, and O2 as well as the energy exchanges between them because the vibrational kinetics is known to be crucial for energy efficient CO2 splitting. First, we use a self-consistent gas temperature calculation in order to assess the key performance indicators for CO2 splitting, i.e., the CO2 conversion and corresponding energy efficiency. Our results indicate that lower pressures and higher power densities lead to more vibrational excitation, which is beneficial for the conversion. We also demonstrate the key role of the gas temperature. The model predicts the highest conversion and energy efficiencies at pressures around 300 mbar, which is in agreement with experiments from the literature. We also show the beneficial aspect of fast gas cooling in the afterglow at high pressure. In a second step, we study in more detail the effects of pressure, gas temperature, and power density on the vibrational distribution function and on the dissociation and recombination mechanisms of CO2, which define the CO2 splitting efficiency. This study allows us to identify the limiting factors of CO2 conversion and to propose potential solutions to improve the process.

Abstract: Knowing the exact number of particles N, and taking this knowledge into account, the quantum canonical ensemble imposes a constraint on the occupation number operators. The constraint particularly hampers the systematic calculation of the partition function and any relevant thermodynamic expectation value for arbitrary but fixed N. On the other hand, fixing only the average number of particles, one may remove the above constraint and simply factorize the traces in Fock space into traces over single-particle states. As is well known, that would be the strategy of the grand-canonical ensemble which, however, comes with an additional Lagrange multiplier to impose the average number of particles. The appearance of this multiplier can be avoided by invoking a projection operator that enables a constraint-free computation of the partition function and its derived quantities in the canonical ensemble, at the price of an angular or contour integration. Introduced in the recent past to handle various issues related to particle-number projected statistics, the projection operator approach proves beneficial to a wide variety of problems in condensed matter physics for which the canonical ensemble offers a natural and appropriate environment. In this light, we present a systematic treatment of the canonical ensemble that embeds the projection operator into the formalism of second quantization while explicitly fixing N, the very number of particles rather than the average. Being applicable to both bosonic and fermionic systems in arbitrary dimensions, transparent integral representations are provided for the partition function Z(N) and the Helmholtz free energy F-N as well as for two- and four-point correlation functions. The chemical potential is not a Lagrange multiplier regulating the average particle number but can be extracted from FN+1 – F-N, as illustrated for a two-dimensional fermion gas. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: We generalize the Gielis Superformula by extending the R. Chacon approach, but avoiding the use of Jacobi elliptic functions. The obtained results are extended to the three-dimensional case. Several new shapes are derived by using the computer algebra system Mathematica(C).

Abstract: The definition of Gielis curves allows for the representation of self intersecting curves. The analysis and the understanding of these representations is of major interest for the analytical representation of sectors bounded by multiple subsets of curves (or surfaces), as this occurs for instance in many natural objects. We present a construction scheme based on R-functions to build signed potential fields with guaranteed differential properties, such that their zero-set corresponds to the outer, the inner envelop, or combined subparts of the curve. Our framework is designed to allow for the definition of composed domains built upon Boolean operations between several distinct objects or some subpart of self-intersecting curves, but also provides a representation for soft blending techniques in which the traditional Boolean union and intersection become special cases of linear combinations between the objects' potential fields. Finally, by establishing a connection between R-functions and Lame curves, we can extend the domain of the p parameter within the R-p-function from the set of the even positive numbers to the real numbers strictly greater than 1, i.e. p is an element of]1, +infinity[.

Abstract: In-situ fabrication of palladium(0) nanoparticles inside zeolitic imidazolate frameworks (ZIF-8) has been established via one-step facile spray-dry technique. Crystal structures and morphologies of the Pd@ZIF-8 samples are investigated by powder XRD, TEM, SAED, STEM, and EDX techniques. High angle annular dark field scanning transmission electron microscopy (HAAD-STEM) and 3D tomographic analysis confirm the presence of palladium nanoparticles inside the ZIF-8 structure. The porosity, surface area and N-2 physisorption properties are evaluated for Pd@ZIF-8 with various palladium contents. Furthermore, Pd@ZIF-8 samples are effectively applied as heterogeneous catalysts in alkenes hydrogenation. This straightforward method is able to speed up the synthesis of encapsulation of metal nanoparticles in metal organic frameworks. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: In this work, a mathematical model was developed to describe the dynamics of fermentation products in sludge alkaline fermentation systems for the first time. In this model, the impacts of alkaline fermentation on sludge disintegration, hydrolysis, acidogenesis, acetogenesis, and methanogenesis processes are specifically considered for describing the high-level formation of fermentation products. The model proposed successfully reproduced the experimental data obtained from five independent sludge alkaline fermentation studies. The modeling results showed that alkaline fermentation largely facilitated the disintegration, acidogenesis, and acetogenesis processes and severely inhibited methanogenesis process. With the pH increase from 7.0 to 10.0, the disintegration, acidogenesis, and acetogenesis processes respectively increased by 53%, 1030%, and 30% while methane production decreased by 3800%.However, no substantial effect on hydrolysis process was found. The model also indicated that the pathway of acetoclastic methanogenesis was more severely inhibited by alkaline condition than that of hydrogentrophic methanogenesis. (C) 2017 Elsevier Ltd. All rights reserved.

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: The electrodeposition of nickel nanostructures on glassy carbon was investigated in 1:2 choline chloride – urea (1:2 ChCl-U) deep eutectic solvent (DES). By combining electrochemical techniques with ex-situ FE-SEM, XPS, HAADF-STEM and EDX, the electrochemical processes occurring during nickel deposition were better understood. Special attention was given to the interaction between the solvent and the growing nickel nanoparticles. The application of a suffciently negative potential results into the electrocatlytic hydrolisis of residual water in the DES, which leads to the formation of a mixed layer of Ni/Ni(OH)2(ads). In addition, hydrogen bonds between hydroxide species and the DES components could be formed, quenching the growth of the nickel clusters favouring their aggregation. Due to these processes, a highly dense distribution of nickel nanostructures can be obtained within a wide potential range. Understanding the role of residual water and the interactions at the interface during metal electrodeposition from DESs is essential to produce supported nanostructures in a controllable way for a broad range of applications and technologies.

Abstract: This work addresses the dissolution corrosion behaviour of 316L austenitic stainless steels. For this purpose, solution-annealed and cold-deformed 316L steels were simultaneously exposed to oxygen-poor (<10-8 mass%) static liquid lead-bismuth eutectic (LBE) for 253e3282 h at 500 °C. Corrosion was consistently more severe for the cold-drawn steels than the solution-annealed steel, indicating the importance of the steel thermomechanical state. The thickness of the dissolution-affected zone was nonuniform, and sites of locally-enhanced dissolution were occasionally observed. The progress of LBE dissolution attack was promoted by the interplay of certain steel microstructural features (grain boundaries, deformation twin laths, precipitates) with the dissolution corrosion process. The identified dissolution mechanisms were selective leaching leading to steel ferritization, and non-selective leaching; the latter was mainly observed in the solution-annealed steel. The maximum corrosion rate decreased with exposure time and was found to be inversely proportional to the depth of dissolution attack.

Abstract: Quantum spin Hall (QSH) insulators with a large topologically nontrivial bulk gap are crucial for future applications of the QSH effect. Among these, group III-V monolayers and their halides, which have a chair structure (regular hexagonal framework), have been widely studied. Using first-principles calculations, we formulate a new structure model for the functionalized group III-V monolayers, which consist of rectangular GaBi-X-2 (X = I, Br, Cl) monolayers with a distorted hexagonal framework (DHF). These structures have a far lower energy than the GaBi-X-2 monolayers with a chair structure. Remarkably, the DHF GaBi-X-2 monolayers are all QSH insulators, which exhibit sizeable nontrivial band gaps ranging from 0.17 to 0.39 eV. The band gaps can be widely tuned by applying different spin-orbit coupling strengths, resulting in a distorted Dirac cone.

Abstract: The formation process of a microdischarge (MD) in both μm- and nm-sized catalyst pores is simulated by a two-dimensional particle-in-cell/Monte Carlo collision model. A parallel-plate dielectric barrier discharge configuration in filamentary mode is considered in ambient air. The discharge is powered by a high voltage pulse. Our calculations reveal that a streamer can penetrate into the surface features of a porous catalyst and MDs can be formed inside both μm- and nm-sized pores, yielding ionization inside the pore. For the μm-sized pores, the ionization mainly occurs inside the pore, while for the nm-sized pores the ionization is strongest near and inside the pore. Thus, enhanced discharges near and inside the mesoporous catalyst are observed. Indeed, the maximum values of the electric field, ionization rate and electron density occur near and inside the pore. The maximum electric field and electron density inside the pore first increase when the pore size rises from 4 nm to 10 nm, and then they decrease for the 100 nm pore, due to
a more pronounced surface discharge for the smaller pores. However, the ionization rate is highest for the 100 nm pore due to the largest effective ionization region.

Abstract: One of the most challenging and recurring problems when modeling plasmas is the lack of data on the key atomic and molecular reactions that drive plasma processes. Even when there are data for some reactions, complete and validated datasets of chemistries are rarely available. This hinders research on plasma processes and curbs development of industrial applications. The QDB project aims to address this problem by providing a platform for provision, exchange, and validation of chemistry datasets. A new data model developed for QDB is presented. QDB collates published data on both electron scattering and heavy-particle reactions. These data are formed into reaction sets, which are then validated against experimental data where possible. This process produces both complete chemistry sets and identifies key reactions that are currently unreported in the literature. Gaps in the datasets can be filled using established theoretical methods. Initial validated chemistry sets for SF6/CF4/O2 and SF6/CF4/N2/H2 are presented as examples.