Abstract: The strain fields of periodically buckled graphene induce a periodic pseudomagnetic field (PMF) that modifies the electronic band structure. From the geometry, amplitude, and period of the periodic pseudomagnetic field, we determine the necessary conditions to access the regime of correlated phases by examining the band flattening. As compared to twisted bilayer graphene the proposed system has the advantages that (1) only a single layer of graphene is needed, (2) one is not limited to hexagonal superlattices, and (3) narrower flat bandwidth and larger separation between flat bands can be induced. We, therefore, propose that periodically strained graphene single layers can become a platform for the exploration of exotic many-body phases.

Abstract: External fields are a powerful tool to probe optical excitations in a material. The linear energy shift of an excitation in a magnetic field is quantified by its effective g factor. Here we show how exciton g factors and their sign can be determined by converged first-principles calculations. We apply the method to monolayer excitons in semiconducting transition metal dichalcogenides and to interlayer excitons in MoSe2/WSe2 heterobilayers and obtain good agreement with recent experimental data. The precision of our method allows us to assign measured g factors of optical peaks to specific transitions in the band structure and also to specific regions of the samples. This revealed the nature of various, previously measured interlayer exciton peaks. We further show that, due to specific optical selection rules, g factors in van der Waals heterostructures are strongly spin and stacking-dependent. The calculation of orbital angular momenta requires the summation over hundreds of bands, indicating that for the considered two-dimensional materials the basis set size is a critical numerical issue. The presented approach can potentially be applied to a wide variety of semiconductors.

Abstract: The dielectric properties of confined water is of fundamental interest and is still controversial. For water confined in channels with height smaller than h = 8 angstrom, we found a commensurability effect and an extraordinary decrease in the out-of-plane dielectric constant down to the limit of the dielectric constant of optical water. Spatial resolved polarization density data obtained from molecular dynamics simulations are found to be antisymmetric across the channel and are used as input in a mean-field model for the dielectric constant as a function of the height of the channel for h > 15 angstrom. Our results are in excellent agreement with a recent experiment [L. Fumagalli et al., Science 360, 1339 (2018)].

Abstract: Significant structural evolution occurs during the deposition of CuInSe2 solar materials when the Cu content increases. We use in situ heating in a scanning transmission electron microscope to directly observe how grain boundaries migrate during heating, causing nondefected grains to consume highly defected grains. Cu substitutes for In in the near grain boundary regions, turning them into a Cu-Se phase topotactic with the CuInSe2 grain interiors. Together with density functional theory and molecular dynamics calculations, we reveal how this Cu-Se phase makes the grain boundaries highly mobile.

Abstract: We report charge transfer up to a single electron per interfacial unit cell across nonpolar heterointerfaces from the Mott insulator LaTiO3 to the charge transfer insulator LaCoO3. In high-quality bi- and trilayer systems grown using pulsed laser deposition, soft x-ray absorption, dichroism, and scanning transmission electron microscopy-electron energy loss spectroscopy are used to probe the cobalt-3d electron count and provide an element-specific investigation of the magnetic properties. The experiments show the cobalt valence conversion is active within 3 unit cells of the heterointerface, and able to generate full conversion to 3d7 divalent Co, which displays a paramagnetic ground state. The number of LaTiO3/LaCoO3 interfaces, the thickness of an additional, electronically insulating “break” layer between the LaTiO3 and LaCoO3, and the LaCoO3 film thickness itself in trilayers provide a trio of control knobs for average charge of the cobalt ions in LaCoO3, illustrating the efficacy of O−2p band alignment as a guiding principle for property design in complex oxide heterointerfaces.

Abstract: The wide measured range of thermal conductivities (k) for monolayer MoS2 and the corresponding incongruent calculated values in the literature all suggest that extrinsic defect thermal resistance is significant and varied in synthesized samples of this material. Here we present defect-mediated thermal transport calculations of MoS2 using interatomic forces derived from density functional theory combined with Green's function methods to describe phonon-point-defect interactions and a Peierls-Boltzmann formalism for transport. Conductivity calculations for bulk and monolayer MoS2 using different density functional formalisms are compared. Nonperturbative first-principles methods are used to describe defect-mediated spectral functions, scattering rates, and phonon k, particularly from sulfur vacancies (VS), and in the context of the plethora of measured and calculated literature values. We find that k of monolayer MoS2 is sensitive to phonon-VS scattering in the range of experimentally observed densities, and that first-principles k calculations using these densities can explain the range of measured values found in the literature. Furthermore, measured k values for bulk MoS2 are more consistent because VS defects are not as prevalent.

Abstract: We used atom probe tomography to complement electron microscopy for the investigation of spinodal decomposition in alkali feldspar. To this end, gem-quality alkali feldspar of intermediate composition with a mole fraction of a(K) = 0.43 of the K end-member was prepared from Madagascar orthoclase by ion-exchange with (NaK)Cl molten salt. During subsequent annealing at 550 degrees C and close to ambient pressure the ion-exchanged orthoclase unmixed producing a coherent lamellar intergrowth of Na-rich and K-rich lamellae. The chemical separation was completed, and equilibrium Na-K partitioning between the different lamellae was attained within four days, which was followed by microstructural coarsening. After annealing for 4 days, the wavelength of the lamellar microstructure was approximate to 17 nm and it increased to approximate to 30 nm after annealing for 16 days. The observed equilibrium compositions of the Na-rich and K-rich lamellae are in reasonable agreement with an earlier experimental determination of the coherent solvus. The excess energy associated with compositional gradients at the lamellar interfaces was quantified from the initial wavelength of the lamellar microstructure and the lamellar compositions as obtained from atom probe tomography using the Cahn-Hilliard theory. The capability of atom probe tomography to deliver quantitative chemical compositions at nm resolution opens new perspectives for studying the early stages of exsolution. In particular, it helps to shed light on the phase relations in nm scaled coherent intergrowth.

Abstract: The penetration of reactive oxygen and nitrogen species (RONS) across the stratum corneum (SC) is a necessary and crucial process in many skin‐related plasma medical applications. To gain more insights into this penetration behavior, we combined experimental measurements of the permeability of dry and moist SC layers with computer simulations of model lipid membranes. We measured the permeation of relatively stable molecules, which are typically generated by plasma, namely H2O2, NO3−, and NO2−. Furthermore, we calculated the permeation free energy profiles of the major plasma‐generated RONS and their derivatives (i.e., H2O2, OH, HO2, O2, O3, NO, NO2, N2O4, HNO2, HNO3, NO2−, and NO3−) across native and oxidized SC lipid bilayers, to understand the mechanisms of RONS permeation across the SC. Our results indicate that hydrophobic RONS (i.e., NO, NO2, O2, O3, and N2O4) can translocate more easily across the SC lipid bilayer than hydrophilic RONS (i.e., H2O2, OH, HO2, HNO2, and HNO3) and ions (i.e., NO2− and NO3−) that experience much higher permeation barriers. The permeability of RONS through the SC skin lipids is enhanced when the skin is moist and the lipids are oxidized. These findings may help to understand the underlying mechanisms of plasma interaction with a biomaterial and to optimize the environmental parameters in practice in plasma medical applications.

Abstract: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

Abstract: Multicellularity is a key evolutionary innovation, leading to coordinated activity and resource sharing among cells, which generally occurs via the physical exchange of chemical compounds. However, filamentous cable bacteria display a unique metabolism in which redox transformations in distant cells are coupled via long-distance electron transport rather than an exchange of chemicals. This challenges our understanding of organismal functioning, as the link among electron transfer, metabolism, energy conservation, and filament growth in cable bacteria remains enigmatic. Here, we show that cells within individual filaments of cable bacteria display a remarkable dichotomy in biosynthesis that coincides with redox zonation. Nanoscale secondary ion mass spectrometry combined with 13 C (bicarbonate and propionate) and 15 N-ammonia isotope labeling reveals that cells performing sulfide oxidation in deeper anoxic horizons have a high assimilation rate, whereas cells performing oxygen reduction in the oxic zone show very little or no label uptake. Accordingly, oxygen reduction appears to merely function as a mechanism to quickly dispense of electrons with little to no energy conservation, while biosynthesis and growth are restricted to sulfide-respiring cells. Still, cells can immediately switch roles when redox conditions change, and show no differentiation, which suggests that the “community service” performed by the cells in the oxic zone is only temporary. Overall, our data reveal a division of labor and electrical cooperation among cells that has not been seen previously in multicellular organisms.

Abstract: High rate activated sludge (HRAS) processes have a high potential for carbon and energy recovery from sewage, yet they suffer frequently from poor settleability due to flocculation issues. The process of flocculation is generally optimized using jar tests. However, detailed jar hydrodynamics are often unknown, and average quantities are used, which can significantly differ from the local conditions. The presented work combined experimental and numerical data to investigate the impact of local hydrodynamics on HRAS flocculation for two different jar test configurations (i.e., radial vs. axial impellers at different impeller velocities) and compared the hydrodynamics in these jar tests to those in a representative section of a full scale reactor using computational fluid dynamics (CFD). The analysis showed that the flocculation performance was highly influenced by the impeller type and its speed. The axial impeller appeared to be more appropriate for floc formation over a range of impeller speeds as it produced a more homogeneous distribution of local velocity gradients compared to the radial impeller. In contrast, the radial impeller generated larger volumes (%) of high velocity gradients in which floc breakage may occur. Comparison to local velocity gradients in a full scale system showed that also here, high velocity gradients occurred in the region around the impeller, which might significantly hamper the HRAS flocculation process. As such, this study showed that a model based approach was necessary to translate lab scale results to full scale. These new insights can help improve future experimental setups and reactor design for improved HRAS flocculation.

Abstract: Waste treatment taxation is a popular policy instrument in many European countries and regions. Its impact on household waste has extensively been researched. However, only little research exists which looks into the impact of waste treatment taxation on industrial waste generation. Nevertheless, industrial waste constitutes more than ninety percent of waste generated in the European Union. This study assesses the impact of an incineration tax on the generation of industrial plastic waste in Flanders, Belgium. We conduct different types of econometrical panel analyzes and provide statistical evidence that firms show lagged behavior, which means that the previous year’s waste generation partly determines the current year’s. The dynamic panel estimations show robust results, indicating that a growth of incineration taxes exert significant negative effects on the growth of industrial plastic waste generation. This result offers no argument to iteratively raise incineration taxes. We conclude that incineration taxation is meaningful if tax rates are set according to the prevailing market conditions, i.e. taking into account the marginal costs of alternatives for incineration. In the short run, the effectiveness of taxation will quickly diminish due to the rapidly rising marginal costs of waste reduction. In the long run, extra recycling capacity is needed to recycle the minimized waste fraction. The role of taxation in the long run is to maintain an equilibrium in which recycling is preferred by the market.

Abstract: In a circular economy (CE), materials, components and products should be kept at the highest level of functionality, while phenomena like dilution, mixing and contamination, often referred to as the loss of resources, should be avoided. One method that can assess the performance of systems to concentrate or avoid dilution of resources is Statistical Entropy Analysis (SEA). Up till now, the method has been applied on the substance level (elements and compounds) only, but showed its applicability to various scales and a variety of systems. Further development of the method allowed to consider information on the product, component and material levels, which makes the method applicable to different combinations of CE strategies, both destructive (e.g. recycling) and non-destructive (e.g. reuse). The method is demonstrated on a simplified vehicle life-cycle, which is modeled through four component groups and six materials. It shows that the method allows to evaluate different CE strategies and identify critical stages which lead to the most severe resource and functionality losses. Based on the methods results, it is possible to determine a perfect circularity reference level, representing a system state that preserves functionality and avoids resource losses. The introduction of a circularity reference level enables the establishment of a framework for resource effectiveness in which diluting and concentrating effects of activities (e.g. sorting) are quantified. The distance of a system to an ideal circular state determines the deviation from a resource-effective system that maintains the original product functionality over a maximum period of time, with minimal efforts.

Abstract: We present KITE, a general purpose open-source tight-binding software for accurate real-space simulations of electronic structure and quantum transport properties of large-scale molecular and condensed systems with tens of billions of atomic orbitals (N similar to 10(10)). KITE's core is written in C++, with a versatile Python-based interface, and is fully optimized for shared memory multi-node CPU architectures, thus scalable, efficient and fast. At the core of KITE is a seamless spectral expansion of lattice Green's functions, which enables large-scale calculations of generic target functions with uniform convergence and fine control over energy resolution. Several functionalities are demonstrated, ranging from simulations of local density of states and photo-emission spectroscopy of disordered materials to large-scale computations of optical conductivity tensors and real-space wave-packet propagation in the presence of magneto-static fields and spin-orbit coupling. On-the-fly calculations of real-space Green's functions are carried out with an efficient domain decomposition technique, allowing KITE to achieve nearly ideal linear scaling in its multi-threading performance. Crystalline defects and disorder, including vacancies, adsorbates and charged impurity centres, can be easily set up with KITE's intuitive interface, paving the way to user-friendly large-scale quantum simulations of equilibrium and non-equilibrium properties of molecules, disordered crystals and heterostructures subject to a variety of perturbations and external conditions.