Abstract: The adsorption of CO2 on surfaces of graphene and carbon nanotubes (CNTs), decorated with Fe atoms, are investigated using the self-consistent-charge density-functional tight-binding (SCC-DFTB) method, neglecting the heat effects. Fe ad-atoms are more stable when they are dispersed on hollow sites. They introduce a large density of states at the Fermi level (N-F); where keeping such density low would help in gas sensing. Furthermore, the Fe ad-atom can weaken the C=O double bonds of the chemisorbed CO2 molecule, paving the way for oxygen atoms to drain more charges from Fe. Consequently, chemisorption of CO2 molecules reduces both N-F and the conductance while it enhances the sensitivity with the increasing gas dose. Conducting armchair CNTs (ac-CNTs) have higher sensitivity than graphene and semiconducting zigzag CNTs (zz-CNT5). Comparative study of sensitivity of ac-CNT-Fe composite towards various gases (e.g., O-2, N-2, H-2, H2O, CO and CO2) has shown high sensitivity and selectivity towards CO, CO2 and H2O gases. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: Bacterial and mammalian proteins, such as lysozyme, are gaining increasing interest as anticancer drugs. This study aims to modify the lysozyme structure using cold atmospheric plasma to boost its cancer cell killing effect. We investigated the structure at acidic and neutral pH using various experimental techniques (circular dichroism, fluorescence, and mass spectrometry) and molecular dynamics simulations. The controlled structural modification of lysozyme at neutral pH enhances its activity, while the activity was lost at acidic pH at the same treatment conditions. Indeed, a larger number of amino acids were oxidized at acidic pH after plasma treatment, which results in a greater distortion of the lysozyme structure, whereas only limited structural changes were observed in lysozyme after plasma treatment at neutral pH. We found that the plasma-treated lysozyme significantly induced apoptosis to the cancer cells. Our results reveal that plasma-treated lysozyme could have potential as a new cancer cell killing drug.

Abstract: The free energy surface (FES) for carbon segregation from nickel nanoparticles is obtained from advanced molecular dynamics simulations. A suitable reaction coordinate is developed that can distinguish dissolved carbon atoms from segregated dimers, chains and junctions on the nanoparticle surface. Because of the typically long segregation time scale (up to ms), metadynamics simulations along the developed reaction coordinate are used to construct FES over a wide range of temperatures and carbon concentrations. The FES revealed the relative stability of different stages in the segregation process, and free energy barriers and rates of the individual steps could then be calculated and decomposed into enthalpic and entropic contributions. As the carbon concentration in the nickel nanoparticle increases, segregated carbon becomes more stable in terms of both enthalpy and entropy. The activation free energy of the reaction also decreases with the increase of carbon concentration, which can be mainly attributed to entropic effects. These insights and the methodology developed to obtain them improve our understanding of carbon segregation process across materials science in general, and the nucleation and growth of carbon nanotube in particular.

Abstract: Despite modern therapeutic advances, the survival prospects of pancreatic cancer patients remain poor, due to chemoresistance and dysregulated oncogenic kinase signaling networks. We applied a novel kinome activitymapping approach using biological peptide targets as phospho-sensors to identify vulnerable kinase de­pendencies for therapy sensitization by physical plasma. Ser/Thr-kinome specific activity changes were mapped upon induction of ferroptotic cell death in pancreatic tumor cells exposed to reactive oxygen and nitrogen species of plasma-treated water (PTW). This revealed a broad kinome activity response involving the CAMK, the AGC and CMGC family of kinases. This systems-level kinome network response supports stress adaptive switches between chemoresistant anti-oxidant responses of Kelch-like ECH-associated protein 1 (KEAP1)/Heme Oxy­genase 1 (HMOX1) and ferroptotic cell death sensitization upon suppression of Nuclear factor (erythroid derived 2)-like 2 (NRF2) and Glutathione peroxidase 4 (GPX4). This is further supported by ex vivo experiments in the chicken chorioallantoic membrane assay, showing decreased GPX4 and Glutathione (GSH) expression as well as increased lipid peroxidation, along with suppressed BxPC-3 tumor growth in response to PTW. Taken all together, we demonstrate that plasma treated water-derived oxidants sensitize pancreatic cancer cells to fer­roptotic cell death by targeting a NRF2-HMOX1-GPX4 specific kinase signaling network.

Abstract: Amorphous platinum clusters supported on porous carbon have been envisaged for high-performance fuel cell electrodes. For this application, it is crucial to control the morphology of the Pt layer and the Ptsubstrate interaction to maximize activity and stability. We thus investigate the morphology evolution during Pt cluster growth on a porous carbon substrate employing atomic scale molecular dynamics simulations. The simulations are based on the Pt-C interaction potential using parameters derived from density functional theory and are found to yield a Pt cluster morphology similar to that observed in low loaded fuel cell electrodes prepared by plasma sputtering. Moreover, the simulations show amorphous Pt cluster growth in agreement with X-ray diffraction and transmission electron microscopy experiments on high performance low Pt content (10 μgPt cm−2) loaded fuel cell electrodes and provide a fundamental insight in the cluster growth mechanism.

Abstract: In this paper we employ first-principles calculations to investigate the effect of substitutional Si doping in the amorphous calcium-phosphate (a-HAP) structure on the work of adhesion, integral charge transfer, charge density difference and theoretical tensile strengths between an a-HAP coating and amorphous titanium dioxide (a-TiO2) substrate systemically. Our calculations demonstrate that substitution of a P atom by a Si atom in a-HAP (a-Si-HAP) with the creation of OH-vacancies as charge compensation results in a significant increase of the bonding strength of the coating to the substrate. The work of adhesion of the optimized Si-doped interfaces reaches a value of up to -2.52 J m(-2), which is significantly higher than for the stoichiometric a-HAP/a-TiO2. Charge density difference analysis indicates that the dominant interactions at the interface have significant covalent character, and in particular two Ti-O and three Ca-O bonds are formed for a-Si-HAP/a-TiO2 and one Ti-O and three Ca-O bonds for a-HAP/a-TiO2. From the stress-strain curve, the Young's modulus of a-Si-HAP/a-TiO2 is calculated to be about 25% higher than that of the a-HAP/a-TiO2, and the yielding stress is about 2 times greater than that of the undoped model. Our calculations therefore demonstrate that the presence of Si in the a-HAP structure strongly alters not only the bioactivity and resorption rates, but also the mechanical properties of the a-HAP/a-TiO2 interface. The results presented here provide an important theoretical insight into the nature of the chemical bonding at the a-HAP/a-TiO2 interface, and are particularly significant for the practical medical applications of HAP-based biomaterials.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: Methane reforming by plasma catalysis is a complex process that is far from understood. It requires a multidisciplinary approach which ideally takes into account all effects from the plasma on the catalyst, and vice versa. In this contribution, we focus on the interactions of CHx (x = {1,2,3}) radicals that are created in the plasma with several nickel catalyst surfaces. To this end, we perform reactive molecular dynamics simulations, based on the ReaxFF potential, in a wide temperature range of 4001600 K. First, we focus on the H2 formation as a function of temperature and surface structure. We observe that substantial H2 formation is obtained at 1400 K and above, while the role of the surface structure seems limited. Indeed, in the initial stage, the type of nickel surface influences the CH bond breaking efficiency of adsorbed radicals; however, the continuous carbon diffusion into the surface gradually diminishes the surface crystallinity and therefore reduces the effect of surface structure on the H2 formation probability. Furthermore, we have also investigated to what extent the species adsorbed on the catalyst surface can participate in surface reactions more in general, for the various surface structures and as a function of temperature. These results are part of the ongoing research on the methane reforming by plasma catalysis, a highly interesting yet complex alternative to conventional reforming processes.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.