Abstract: The electron states in axially symmetric quantum wires are computed by means of the effective-mass Schrodinger equation, which is written in cylindrical coordinates phi, rho, and z. We show that a direct discretization of the Schrodinger equation by central finite differences leads to a non-symmetric Hamiltonian matrix. Because diagonalization of such matrices is more complex it is advantageous to transform it in a symmetric form. This can be done by the Liouville-like transformation proposed by Rizea et al. (2008), which replaces the wave function psi(rho) with the function F(rho) = psi(rho)root rho and transforms the Hamiltonian accordingly. Even though a symmetric Hamiltonian matrix is produced by this procedure, the computed wave functions are found to be inaccurate near the origin, and the accuracy of the energy levels is not very high. In order to improve on this, we devised a finite-difference scheme which discretizes the Schrodinger equation in the first step, and then applies the Liouville-like transformation to the difference equation. Such a procedure gives a symmetric Hamiltonian matrix, resulting in an accuracy comparable to the one obtained with the finite element method. The superior efficiency of the new finite-difference scheme (FDM) is demonstrated for a few p-dependent one-dimensional potentials which are usually employed to model the electron states in free-standing and core shell quantum wires. The new scheme is compared with the other FDM schemes for solving the effective-mass Schrodinger equation, and is found to deliver energy levels with much smaller numerical error for all the analyzed potentials. It also gives more accurate results than the scheme of Rizea et al., except for the ground state of an infinite rectangular potential in freestanding quantum wires. Moreover, the PT symmetry is invoked to explain similarities and differences between the considered FDM schemes. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: Large-eddy simulations (LESS) are frequently used to model the planetary boundary layer, and the choice of the grid cell size, numerical schemes and sub grid model can significantly influence the simulation results. In the present paper the impact of grid spacing on LES of an idealized atmospheric convective boundary layer (CBL), for which the statistics and flow structures are well understood, is assessed for two mesoscale models: the Regional Atmospheric Modeling System (RAMS) and the Weather Research and Forecasting model (WRF). Nine simulations are performed on a fixed computational domain (6 x 6 x 2 km), combining three different horizontal (120, 60, 30 m) and vertical (20, 10, 5 m) spacings. The impact of the cell size on the CBL is investigated by comparing turbulence statistics and velocity spectra. The results demonstrate that both WRF and RAMS can perform LES of the CBL under consideration without requiring extremely high computational loads, but they also indicate the importance of adopting a computational grid that is adequate for the numerical schemes and subgrid models used. In both RAMS and WRF a horizontal cell size of 30 m is required to obtain a suitable turbulence reproduction throughout the CBL height. Considering the vertical grid spacing, WRF produced similar results for all the three tested values, while in RAMS it should be ensured that the aspect ratio of the cells does not exceed a value of 3. The two models were found to behave differently in function of the grid resolution, and they have different shortcomings in their prediction of CBL turbulence. WRF exhibits enhanced damping at the smallest scales, while RAMS is prone to the appearance of spurious fluctuations in the flow when the grid aspect ratio is too high. (C) 2015 Elsevier Ltd. All rights reserved.

Abstract: Reactive oxygen and nitrogen species (RONS) are involved in many biochemical processes, including nitrooxidative stress that causes cancer cell death, observed in cancer therapies such as photodynamic therapy and cold atmospheric plasma. However, their mechanisms of action and selectivity still remain elusive due to the complexity of biological cells. For example, it is not well known how RONS generated by cancer therapies permeate the cell membrane to cause nitro-oxidative damage. There are many studies dedicated to the perme­ation of RONS across native and oxidized membranes, but not across nitrated membranes, another lipid product also generated during nitro-oxidative stress. Herein, we performed molecular dynamics (MD) simulations to calculate the free energy barrier of RONS permeation across nitrated membranes. Our results show that hy­drophilic RONS, such as hydroperoxyl radical (HO2) and peroxynitrous acid (ONOOH), have relatively low barriers compared to hydrogen peroxide (H2O2) and hydroxyl radical (HO), and are more prone to permeate the membrane than for the native or peroxidized membranes, and similar to aldehyde-oxidized membranes. Hy­drophobic RONS like molecular oxygen (O2), nitrogen dioxide (NO2) and nitric oxide (NO) even have insignif­icant barriers for permeation. Compared to native and peroxidized membranes, nitrated membranes are more permeable, suggesting that we must not only consider oxidized membranes during nitro-oxidative stress, but also nitrated membranes, and their role in cancer therapies.

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: Structural, mechanical and electronic properties of two-dimensional single-layer hexagonal structures in the (111) crystal plane of IIIAs-ZnS systems (III = B, Ga and In) are studied by first-principles calculations based on density functional theory (DFT). Elastic and phonon dispersion relation display that 2D h-IIIAs systems (III = B, Ga and In) are both mechanical and dynamically stable. Electronic structures analysis show that the semiconducting nature of the 3D-IIIAs compounds is retained by their 2D single layer counterpart. Furthermore, density of states reveals the influence of sigma and pi bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Calculations of elastic constants show that the Young's modulus, bulk modulus and shear modulus decrease for 2D h-IIIAs binary compounds as we move down on the group of elements of the periodic table. In addition, as the bond length between the neighboring cation-anion atoms increases, the 2D h-IIIAs binary compounds display less stiffness and more plasticity. Our findings can be used to understand the contribution of the r and p bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Structural and electronic properties of h-IIIAs systems as a function of the number of layers have been also studied. It is shown that h-BAs keeps its planar geometry while both h-GAs and h-InAs retained their buckled ones obtained by their single layers. Bilayer h-IIIAs present the same bandgap nature of their counterpart in 3D. As the number of layers increase from 2 to 4, the bandgap width for layered h-IIIAs decreases until they become semimetal or metal. Interestingly, these results are different to those found for layered h-GaN. The results presented in this study for single and few-layer h-IIIAs structures could give some physical insights for further theoretical and experimental studies of 2D h-IIIV-like systems. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.

Abstract: Topological superconductivity is of great contemporary interest and has been proposed in doped Bi₂Se₃, in which electron-donating atoms such as Cu, Sr or Nb have been intercalated into the Bi₂Se₃structure. For Nb_{<italic>x</italic>}Bi₂Se₃, with<italic>T</italic>_c ~ 3 K, it is assumed in the literature that Nb is inserted in the van der Waals gap. However, in this work an alternative origin for the superconductivity in Nb-doped Bi₂Se₃is established. In contrast to previous reports, it is deduced that Nb intercalation in Bi₂Se₃does not take place. Instead, the superconducting behaviour in samples of nominal composition Nb_{<italic>x</italic>}Bi₂Se₃results from the (BiSe)_1.10NbSe₂misfit phase that is present in the sample as an impurity phase for small<italic>x</italic>(0.01 ≤ <italic>x</italic> ≤ 0.10) and as a main phase for large<italic>x</italic>(<italic>x</italic> = 0.50). The structure of this misfit phase is studied in detail using a combination of X-ray diffraction and transmission electron microscopy techniques.

Abstract: Colors in many paintings of great art historical value have changed over time, due to the combined effects of natural ageing, accumulated surface grime, and materials added during later conservation treatments. The physical restoration of the colors in such paintings is not possible. This article describes one part of work done to digitally restore the colors of Van Gogh's painting Field with Irises near Arles, dating from May 1888. We have used multispectral reflectance data to estimate absorption K and backscattering S parameters of Kubelka-Munk 2-constant theory. This was done for all 13 pigments known to have been used by Van Gogh in this painting, and based on this the concentration maps for each of these pigments were calculated. We validated the calculated concentration maps in several ways. For some pigments, we were able to predict spots on the painting where the pigment is expected to occur in unmixed form based on visual examination. For several other pigments, the concentration maps could be shown to agree with XRF data. Finally, for some other pigments the concentration maps were supported by additional evidence from microscopic examinations, remarks in Van Gogh's letters and from early color reproductions. For the 1.7 million pixels for which multispectral data is available, the average color difference between the calculated and measured spectral reflectance curves is CIEDE2000 = 1.05. This further confirms that the Kubelka-Munk calculations are well suited to describe the variety of spectral reflectance on the painting.

Abstract: Superparamagnetic iron oxide magnetic nanoparticles (MNPs) are successfully used as contrast agents in magnetic-resonance imaging. They can be easily functionalized for drug delivery functions, demonstrating great potential for both imaging and therapeutic applications. Here we developed new pH-responsive theranostic core-shell-corona nanoparticles consisting of superparamagentic Fe3O4 core that displays high T2 relaxivity, bovine serum albumin (BSA) shell that binds anticancer drug, doxorubicin (Dox) and poly(ethylene glycol) (PEG) corona that increases stability and biocompatibility. The nanoparticles were produced by adsorption of the BSA shell onto the Fe3O4 core followed by crosslinking of the protein layer and subsequent grafting of the PEG corona using monoamino-terminated PEG via carbodiimide chemistry. The hydrodynamic diameter, zeta-potential, composition and T2 relaxivity of the resulting nanoparticles were characterized using transmission electron microscopy, dynamic light scattering, thermogravimetric analysis and T2-relaxometry. Nanoparticles were shown to absorb Dox molecules, possibly through a combination of electrostatic and hydrophobic interactions. The loading capacity (LC) of the nanoparticles was 8 wt.%. The Dox loaded nanoparticles release the drug at a higher rate at pH 5.5 compared to pH 7.4 and display similar cytotoxicity against C6 and HEK293 cells as the free Dox. (C) 2015 Elsevier B.V. All rights reserved.

Abstract: The development of coatings with tunable performances is critical to meet a wide range of technological applications each one with different requirements. Using the plasma-enhanced chemical vapor deposition (PECVD) process, scientists can create hydrogenated amorphous carbon coatings doped with metal (a-C:H:Me) with a broad range of mechanical properties, varying from those resembling polymers to ones resembling diamond. These diverse properties, without clear relations between the different families, make the material selection and optimization difficult but also very rich. An innovative approach is proposed here based on projected performance indices related to fracture energy, strength, and stiffness in order to classify and optimize a-C:H:Me coatings. Four different a-C:H:Cr coatings deposited by PECVD with Ar/C2H2 discharge under different bias voltage and pressures are investigated. A path is found to produce coatings with a selective critical energy release rate between 5–125 J/m2 without compromising yield strength (1.6–2.7 GPa) and elastic limit (≈0.05). Finally, fine-tuned coatings are categorized to meet desired applications under different testing conditions.

Abstract: The Ricardian model is a widely used approach based on cross-sectional regression analysis to estimate climate change impacts on agricultural productivity. Up until now, researchers have focused on the impacts of gradual changes in temperature and precipitation, even though climate change is known to encompass also changes in the severity and frequency of extreme weather events. This research investigates the impact of heatwaves on European agriculture, additional to the impact of average climate change. Using a dataset of more than 60,000 European farms, the study examines whether adding a measure for heatwaves to the Ricardian model influences its results. We find that heatwaves have a minor impact on agricultural productivity and that this impact is moderated by average temperature. In colder regions, farm productivity increases with the number of heatwave days. For warmer regions, land values decrease with heatwave frequency. Despite the moderating effect, the marginal effect of heatwave frequency, i.e. the percentage change in agricultural land values caused by one more heatwave day per year, is small in comparison to the effect of average temperature increases. Non-marginal effects are found to be relevant, but only in the case of increased heatwave frequency. According to our results, farms are not expected to suffer more from extreme weather than from mean climate change, as was claimed by several previous studies.

Abstract: Economic and technological assessments have identified difficulties with the commercialization of bulk products from microalgae, like biofuels. To overcome these problems, a multi-product algal-based biorefinery has been proposed. This paper performs a techno-economic assessment of such a biorefinery. Four production pathways, ranging from a base case with commercial technologies to an improved case with innovative technologies, are analyzed. All region-specific parameters were adapted to Belgian conditions. Three scenarios result in techno-economically viable production plants. The most profitable scenario is the scenario which uses a specialized membrane for medium recycling and an open pond algae cultivation. Although the inclusion of a photobioreactor decreases the culture medium costs, the higher investment costs result in lower economic profits. The carotenoid content and price are identified as critical parameters. Furthermore, the economies of scale assumption for the photobioreactor is critical for the feasibility of this cultivation technology. The techno-economic assessment is an important methodology to guide and evaluate further improvements in research and shorten the time-to-market for innovative technologies in this field.

Abstract: The value of phytoextracting crops (plants cultivated for soil remediation) depends on the profitability of the sequential investment in a conversion technology aimed at the economic valorization of the plants. However, the net present value (NPV) of an investment in such an innovative technology is risky due to technical and economic uncertainties. Therefore, decision makers want to dispose of information about the probability of a positive NPV, the largest possible loss, and the crucial economic and technical parameters influencing the NPV. This paper maps the total uncertainty in the NPV of an investment in fast pyrolysis for the production of combined heat and power from willow cultivated for phytoextraction in the Belgian Campine. The probability of a positive NPV has been calculated by performing Monte Carlo simulations. Information about possible losses has been provided by means of experimental design. Both methods are then combined in order to identify the key economic and technical parameters influencing the project's profitability. It appears that the case study has a chance of 87% of generating a positive NPV with an expected value of 3 million euro (MEUR), while worst-case scenarios predict possible losses of 7 MEUR. The amount of arable land, the biomass yield, the purchase price of the crop, the policy support, and the product yield of fast pyrolysis are identified as the most influential parameters. It is concluded that both methods, i.e., Monte Carlo simulations and experimental design, provide decision makers with complementary information with regard to economic risk.

Abstract: MXenes, a new family of functional two-dimensional (2D) materials, have shown great potential for an extensive variety of applications within the last decade. Atomic defects and functional groups in MXenes are known to have a tremendous influence on the functional properties. In this review, we focus on recent progress in the characterization of atomic defects and functional group chemistry in MXenes, and how to control them to directly influence various properties (e.g., electron transport, Li' adsorption, hydrogen evolution reaction (HER) activity, and magnetism) of 2D MXenes materials. Dynamic structural transformations such as oxidation and growth induced by atomic defects in MXenes are also discussed. The review thus provides perspectives on property optimization through atomic defect engineering, and bottom-up synthesis methods based on defect-assisted homoepitaxial growth of MXenes. (C) 2020 Chinese Chemical Society and Institute of Materia Medica, Chinese Academy of Medical Sciences. Published by Elsevier B.V. All rights reserved.