Abstract: The adsorption and diffusion of Li, Na, K and Ca atoms on a Mo2C monolayer are systematically investigated by using first principles methods. We found that the considered metal atoms are strongly bound to the Mo2C monolayer. However, the adsorption energies of these alkali and earth alkali elements decrease as the coverage increases due to the enhanced repulsion between the metal ions. We predict a significant charge transfer from the ad-atoms to the Mo2C monolayer, which indicates clearly the cationic state of the metal atoms. The metallic character of both pristine and doped Mo2C ensures a good electronic conduction that is essential for an optimal anode material. Low migration energy barriers are predicted as small as 43 meV for Li, 19 meV for Na and 15 meV for K, which result in the very fast diffusion of these atoms on Mo2C. For Mo2C, we found a storage capacity larger than 400 mA h g(-1) by the inclusion of multilayer adsorption. Mo2C expands slightly upon deposition of Li and Na even at high concentrations, which ensures the good cyclic stability of the atomic layer. The calculated average voltage of 0.68 V for Li and 0.30 V for Na ions makes Mo2C attractive for low charging voltage applications.

Abstract: In this work, we study the abatement of dilute trichloroethylene (TCE) in air with a negative direct current corona discharge. A numerical model is used to theoretically investigate the underlying plasma chemistry for the removal of TCE, and a reaction pathway for the abatement of TCE is proposed. The Cl atom, mainly produced by dissociation of COCl, is one of the controlling species in the TCE destruction chemistry and contributes to the production of chlorine containing by-products. The effect of humidity on the removal efficiency is studied and a good agreement is found between experiments and the model for both dry (5 % relative humidity (RH)) and humid air (50 % RH). An increase of the relative humidity from 5 % to 50 % has a negative effect on the removal efficiency, decreasing by ±15 % in humid air. The main loss reactions for TCE are with ClO·, O· and CHCl2. Finally, the by-products and energy cost of TCE abatement are discussed.

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: In this paper, an overview is given of modeling activities going on in our research group, for describing the plasma chemistry and plasmasurface interactions in reactive plasmas. The plasma chemistry is calculated by a fluid approach or by hybrid Monte Carlo (MC)fluid modeling. An example of both is illustrated in the first part of the paper. The example of fluid modeling is given for a dielectric barrier discharge (DBD) in CH4/O2, to describe the partial oxidation of CH4 into value-added chemicals. The example of hybrid MCfluid modeling concerns an inductively coupled plasma (ICP) etch reactor in Ar/Cl2/O2, including also the description of the etch process. The second part of the paper deals with the treatment of plasmasurface interactions on the atomic level, with molecular dynamics (MD) simulations or a combination of MD and MC simulations.

Abstract: Because of the unique properties of plasma technology, its use in gas conversion applications is gaining significant interest around the globe. Plasma-based CO2 and CH4 conversion has become a major research area. Many investigations have already been performed regarding the single-component gases, that is, CO2 splitting and CH4 reforming, as well as for two-component mixtures, that is, dry reforming of methane
(CO2/CH4), partial oxidation of methane (CH4/O2), artificial photosynthesis (CO2/H2O), CO2 hydrogenation (CO2/H2), and even first steps toward the influence of N2 impurities have been taken, that is, CO2/N2 and CH4/N2. In this Feature Article we briefly discuss the advances made in literature for these different steps from a plasma chemistry modeling point of view. Subsequently, we present a comprehensive plasma chemistry set, combining the knowledge gathered in this field so far and supported with extensive experimental data. This set can be used for chemical kinetics plasma modeling for all possible combinations of CO2, CH4, N2, O2, and H2O to investigate the bigger picture of the underlying plasmachemical pathways for these mixtures in a dielectric barrier discharge plasma. This is extremely valuable
for the optimization of existing plasma-based CO2 conversion and CH4 reforming processes as well as for investigating the influence of N2, O2, and H2O on these processes and even to support plasma-based multireforming processes.

Abstract: Natural deep eutectic solvents (NADES) are mixtures of naturally derived compounds with a significantly decreased melting point due to the specific interactions among the constituents. NADES have benign properties (low volatility, flammability, toxicity, cost) and tailorable physicochemical properties (by altering the type and molar ratio of constituents), hence they are often considered as a green alternative to common organic solvents. Modeling the relation between their composition and properties is crucial though, both for understanding and predicting their behavior. Several efforts were done to this end, yet this review aims at structuring the present knowledge as an outline for future research. First, we reviewed the key properties of NADES and relate them to their structure based on the available experimental data. Second, we reviewed available modeling methods applicable to NADES. At the molecular level, density functional theory and molecular dynamics allow interpreting density differences and vibrational spectra, and computation of interaction energies. Additionally, properties at the level of the bulk media can be explained and predicted by semi-empirical methods based on ab initio methods (COSMO-RS) and equation of state models (PC-SAFT). Finally, methods based on large datasets are discussed; models based on group contribution methods and machine learning. A combination of bulk media and dataset modeling allows qualitative prediction and interpretation of phase equilibria properties on the one hand, and quantitative prediction of melting point, density, viscosity, surface tension and refractive indices on the other hand. In our view, multiscale modeling, combining the molecular and macroscale methods, will strongly enhance the predictability of NADES properties and their interaction with solutes, yielding truly tailorable solvents to accommodate (bio)chemical reactions.

Abstract: Antibiotics such as sulfamethoxazole (SFX) are environmentally hazardous after being released into the aquatic environment and challenges remain in the development of engineered prevention strategies. In this work, a mathematical model was developed to describe and evaluate cometabolic biotransformation of SFX and its transformation products (TPs) in an enriched ammonia oxidizing bacteria (AOB) culture. The growth-linked cometabolic biodegradation by AOB, non-growth transformation by AOB and nongrowth transformation by heterotrophs were considered in the model framework. The production of major TPs comprising 4-Nitro-SFX, Desamino-SFX and N-4-Acetyl-SFX was also specifically modelled. The validity of the model was demonstrated through testing against literature reported data from extensive batch tests, as well as from long-term experiments in a partial nitritation sequencing batch reactor (SBR) and in a combined SBR + membrane aerated biofilm reactor performing nitrification/denitrification. Modelling results revealed that the removal efficiency of SFX increased with the increase of influent ammonium concentration, whereas the influent organic matter, hydraulic retention time and solid retention time exerted a limited effect on SFX biodegradation with the removal efficiencies varying in a narrow range. The variation of influent SFX concentration had no impact on SFX removal efficiency. The established model framework enables interpretation of a range of experimental observations on SFX biodegradation and helps to identify the optimal conditions for efficient removal. (C) 2017 Elsevier Ltd. All rights reserved.

Abstract: Synchrotron small angle X-ray scattering (SAXS) was used to characterize silicalite-1 zeolite crystallization from TEOS/TPAOH/water clear sol. SAXS patterns were recorded over a broad range of length scales, enabling the simultaneous monitoring of nanoparticles and crystals occurring at various stages of the synthesis. A simple two-population model accurately described the patterns. Nanoparticles were modeled by polydisperse coreshell spheres and crystals by monodisperse oblate ellipsoids. These models were consistent with TEM images. The SAXS results, in conjunction with in situ light scattering, showed that nucleation of crystals occurred in a short period of time. Crystals were uniform in size and shape and became increasingly anisotropic during growth. In the presence of nanoparticles, crystal growth was fast. During crystal growth, the number of nanoparticles decreased gradually but their size was constant. These observations suggested that the nanoparticles were growth units in an aggregative crystal growth mechanism. Crystals grown in the presence of nanoparticles developed a faceted habit and intergrowths. In the final stages of growth, nanoparticles were depleted. Concurrently, the crystal growth rate decreased significantly.

Abstract: Given the ecological problems associated to the CO2 emissions of fossil fuels, plasma technology has gained
interest for conversion of CO2 into value-added products. Microwave plasmas operating at atmospheric pressure
have proven to be especially interesting, due to the high gas temperatures inside the reactor (i.e. up to 6000 K)
allowing for efficient thermal dissociation of CO2 into CO and O2. However, the performance of these high
temperature plasmas is limited by recombination of CO back into CO2 once the gas cools down in the afterglow.
In this work, we computationally investigated several quenching nozzles, developed and experimentally tested
by Hecimovic et al., [1] for their ability to quickly cool the gas after the plasma, thereby quenching the CO
recombination reactions. Using a 3D computational fluid dynamics model and a quasi-1D chemical kinetics
model, we reveal that a reactor without nozzle lacks gas mixing between hot gas in the center and cold gas near
the reactor walls. Especially at low flow rates, where there is an inherent lack of convective cooling due to the
low gas flow velocity, the temperature in the afterglow remains high (between 2000 and 3000 K) for a relatively
long time (in the 0.1 s range). As shown by our quasi-1D chemical kinetics model, this results in a important loss
of CO due to recombination reactions. Attaching a nozzle in the effluent of the reactor induces fast gas quenching
right after the plasma. Indeed, it introduces (i) more convective cooling by forcing cool gas near the walls to mix
with hot gas in the center of the reactor, as well as (ii) more conductive cooling through the water-cooled walls of
the nozzle. Our model shows that gas quenching and the suppression of recombination reactions have more
impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: Mesoporous strontium hydroxyapatite (SrHAp) nanorods (NRs) have been successfully synthesized using a simple and efficient chemical route, i.e. the hydrothermal method. Structural and morphological characterization of the as-synthesized SrHAp NRs have been performed by transmission electron microscopy (TEM) and high angle annular dark field scanning transmission electron microscopy (HAADF-STEM). TEM and HAADF-STEM measurements of the NRs reveal the coexistence of longer and shorter particles with the length ranging from 50 nm to 400 nm and a diameter of about 20-40 nm. Electron tomography measurements of the NRs allow us to better visualize the mesopores and their facets. Two model drugs, hydrophobic risperidone and hydrophilic pramipexole, were loaded into the SrHAp NRs. These nanorods were coated using a modified chitosan (CS) with poly(2-hydroxyethyl methacrylate) (PHEMA), in order to encapsulate the drug-loaded SrHAp nanoparticles and reduce the cytotoxicity of the loaded materials. The drug release from neat and encapsulated SrHAp NRs mainly depends on the drug hydrophilicity. Importantly, although neat SrHAp nanorods exhibit some cytotoxicity against Caco-2 cells, the Cs-g-PHEMA-SrHAp drug-loaded nanorods show an acceptable cytocompatibility.

Abstract: MOF-5 nanoparticles were prepared by mixing a solution of [Zn4O(C6H5COO)(6)] with a solution of benzene-1,4-dicarboxylic acid in DMF at ambient conditions. The former species mimics as a secondary building unit (SBU), and the latter acts as linker. Mixing of the two solutions induced the formation of MOF-5 nanoparticles in dilute suspension. The applied conditions were identified as suitable for a closer investigation of the particle formation process by combined light and small angle neutron scattering (SANS). Scattering analysis revealed a significant impact of the molar ratio of the two components in the reaction mixture. Excessive use of the building unit slowed down the process. A similar effect was observed upon addition of 4n-decylbenzoic acid, which is supposed to act as a modulator. The formation mechanism leads to initial intermediates, which turn into cubelike nanoparticles with a diameter of about 60-80 nm. This initial stage is followed by an extended formation period, where nucleation proceeds over hours, leading to an increasing number of nanoparticles with the same final size of 60-80 nm.

Abstract: A molecular dynamics (MD) study of MgxAlyOz thin films grown by magnetron sputtering is presented using an ionic model and comparing two potential sets with formal and partial charges. The applicability of the model and the reliability of the potential sets for the simulation of thin film growth are discussed. The formal charge potential set was found to reproduce the thin film structure in close agreement with the structure of the experimentally grown thin films. Graphical abstract A molecular dynamics study of growth of MgxAlyOz thin films is presented using an ionic model and comparing two potential sets with formal and partial charges. The simulation results with the formal charge potential set showed a transition in the film from a crystalline to an amorphous structure, when the Mg metal content decreases below 50% in very close agreement with the structure of the experimentally deposited films.

Abstract: Biomembranes are under constant attack of free radicals that may lead to lipid oxidation in conditions of oxidative stress. The products generated during lipid oxidation are responsible for structural and dynamical changes which may jeopardize the membrane function. For instance, the local rearrangements of oxidized lipid molecules may induce membrane rupture. In this study, we investigated the effects of mechanical stress on oxidized phospholipid bilayers (PLBs). Model bilayers were stretched until pore formation (or poration) using nonequilibrium molecular dynamics simulations. We studied single-component homogeneous membranes composed of lipid oxidation products, as well as two-component heterogeneous membranes with coexisting native and oxidized domains. In homogeneous membranes, the oxidation products with —OH and —OOH groups reduced the areal strain required for pore formation, whereas the oxidation product with ]O group behaved similarly to the native membrane. In heterogeneous membranes composed of oxidized and non-oxidized domains, we tested the hypothesis according to which poration may be facilitated at the domain interface region. However, results were inconclusive due to their large statistical variance and sensitivity to simulation setup parameters. We pointed out important technical issues that need to be considered in future simulations of mechanically-induced poration of heterogeneous membranes. This research is of interest for photodynamic therapy and plasma medicine, because ruptured and intact plasma membranes are experimentally considered hallmarks of necrotic and apoptotic cell death.

Abstract: The reaction behavior of species that may affect the growth of ultrananocrystal line and nanocrystalline diamond ((U)NCD) films is investigated by means of molecular dynamics simulations. Impacts of CHx (x = 0 – 4), C2Hx (x=0-6), C3Hx (x=0-2), C4Hx (x = 0 – 2), H, and H-2 on clean and hydrogenated diamond (100)2 x 1 and (111) 1 x 1 surfaces at two different substrate temperatures are simulated. We find that the different bonding structures of the two surfaces cause different temperature effects on the sticking efficiency. These results predict a temperature-dependent ratio of diamond (100) and (111) growth. Furthermore, predictions of which are the most important hydrocarbon species for (U)NCD growth are made.

Abstract: We have studied the dependence of the binding energy of a cubane dimer on the mutual orientation of and the distance between the composing monomers employing the second-order Møller-Plesset perturbation scheme (MP2) with the cc-pVDZ molecular basis set. We have found that the MP2 contribution from the molecular correlations is responsible for the bound state of the cubane dimer, whereas the Hartree-Fock contribution remains anti-bonding at all intermolecular distances. Starting with two molecules in the standard orientation and centers of mass at (0,0,0) and (0,0,d), respectively, the maximal binding energy is found at d = 5.125 Å and one of the monomers rotated by 45° about the z-axis. This configuration implies that the hydrogen atoms belonging to different monomers tend to repel each other. The results are in agreement with experimental data on the optimal packing of cubane molecules in the solid state.