Theoretical Study of Heteroatom Doping in Tuning the Catalytic Activity of Graphene for Triiodide Reduction

Graphene with heteroatom doping has found increasing applications in a broad range of catalytic reactions. However, the doping effects accounting for the enhanced catalytic activity still remain elusive. In this work, taking the triiodide electroreduction reaction as an example, we study systematically the intrinsic activity of graphene and explore the origin of doping-induced activity variation using first-principles calculations, in which two typical N and S dopants are tested. The most common graphene structures, basal plane, armchair edge, and zigzag edge, are considered, and it is found that the former two structures show a weak adsorption ability for the iodine atom (the key intermediate in the triiodide electroreduction reaction), corresponding to a low catalytic activity. Doping either N or S can strengthen the adsorption and thus increase the activity, and the codoping of N and S (NS-G) exhibits a synergistic effect. A detailed investigation into the whole process of the triiodide electroreduction reaction at the CH3CN/NS-G interface is also carried out to verify these activity trends. It is found that the zigzag edges which contain spin electrons show a relatively stronger adsorption strength compared with the basal plane and armchair edge, and initial doping would result in the spin disappearance that evidently weakens the adsorption; with the disappearance of spin, however, further doping can increase the adsorption again, suggesting that the spin electrons may play a preliminary role in affecting the intrinsic activity of graphene. We also analyzed extensively the origin of doping-induced adsorption enhancement of graphene in the absence of spin; it can be rationalized from the electronic and geometric factors. Specifically, N doping can result in a more delocalized “electron-donating area” to enhance I adsorption, while S doping provides a localized structural distortion, which activates the nearest sp2-C into coordinatively unsaturated sp3-C. These results explain well the improved activity of the doping and the synergistic effect of the codoping. The understandings are generalized to provide insight into the enhanced activity of the oxygen reduction reaction on heteroatom doped graphene. This work may be of importance toward the design of high-activity graphene based material.

Multiscale Analysis for Field-Effect Penetration through Two-Dimensional Materials

Gate-tunable two-dimensional (2D) materials-based quantum capacitors (QCs) and van der Waals heterostructures involve tuning transport or optoelectronic characteristics by the field effect. Recent studies have attributed the observed gate-tunable characteristics to the change of the Fermi level in the first 2D layer adjacent to the dielectrics, whereas the penetration of the field effect through the one-molecule-thick material is often ignored or oversimplified. Here, we present a multiscale theoretical approach that combines first-principles electronic structure calculations and the Poisson–Boltzmann equation methods to model penetration of the field effect through graphene in a metal–oxide–graphene–semiconductor (MOGS) QC, including quantifying the degree of “transparency” for graphene two-dimensional electron gas (2DEG) to an electric displacement field. We find that the space charge density in the semiconductor layer can be modulated by gating in a nonlinear manner, forming an accumulation or inversion layer at the semiconductor/graphene interface. The degree of transparency is determined by the combined effect of graphene quantum capacitance and the semiconductor capacitance, which allows us to predict the ranking for a variety of monolayer 2D materials according to their transparency to an electric displacement field as follows: graphene > silicene > germanene > WS2 > WTe2 > WSe2 > MoS2 > phosphorene > MoSe2 > MoTe2, when the majority carrier is electron. Our findings reveal a general picture of operation modes and design rules for the 2D-materials-based QCs.

Dielectric screening in atomically thin boron nitride nanosheets

Two-dimensional (2D) hexagonal boron nitride (BN) nanosheets are excellent dielectric substrate for graphene, molybdenum disulfide, and many other 2D nanomaterial-based electronic and photonic devices. To optimize the performance of these 2D devices, it is essential to understand the dielectric screening properties of BN nanosheets as a function of the thickness. Here, electric force microscopy along with theoretical calculations based on both state-of-the-art first-principles calculations with van der Waals interactions under consideration, and nonlinear Thomas-Fermi theory models are used to investigate the dielectric screening in high-quality BN nanosheets of different thicknesses. It is found that atomically thin BN nanosheets are less effective in electric field screening, but the screening capability of BN shows a relatively weak dependence on the layer thickness.

Direct observation of a long-lived single-atom catalyst chiseling atomic structures in graphene

Fabricating stable functional devices at the atomic scale is an ultimate goal of nanotechnology. In biological processes, such high-precision operations are accomplished by enzymes. A counterpart molecular catalyst that binds to a solid-state substrate would be highly desirable. Here, we report the direct observation of single Si adatoms catalyzing the dissociation of carbon atoms from graphene in an aberration-corrected high-resolution transmission electron microscope (HRTEM). The single Si atom provides a catalytic wedge for energetic electrons to chisel off the graphene lattice, atom by atom, while the Si atom itself is not consumed. The products of the chiseling process are atomic-scale features including graphene pores and clean edges. Our experimental observations and first-principles calculations demonstrated the dynamics, stability, and selectivity of such a single-atom chisel, which opens up the possibility of fabricating certain stable molecular devices by precise modification of materials at the atomic scale.

Interplay of Defects, Magnetism, Ripples and Strain in Graphene

We present a comprehensive study based on first-principles calculations about the interplay of four important ingredients on the electronic structure of graphene: defects + magnetism + ripples + strain. So far they have not been taken into account simultaneously in a set of ab initio calculations. Furthermore, we focus on the strain dependence of the properties of carbon monovacancies, with special attention to magnetic spin moments. We demonstrated that such defects show a very rich structural and spin phase-diagram with many spin solutions as function of strain. At zero strain the vacancy shows a spin moment of 1.5 Bohrs that increases up to 2 Bohrs with stretching. Changes are more dramatic under compression: the vacancy becomes non-magnetic under a compression larger than 2%. This transition is linked to the structural modifications associated with the formation of ripples in the graphene layer. Our results suggest that such interplay could have important implications for the design of future spintronics devices based on graphene derivatives, as for example a spin-strain switch based on vacancies.

Investigating routes toward atomic layer deposition of silicon carbide: Ab initio screening of potential silicon and carbon precursors

Silicon carbide (SiC) is a promising material for electronics due to its hardness, and ability to carry high currents and high operating temperature. SiC films are currently deposited using chemical vapor deposition (CVD) at high temperatures 1500–1600 °C. However, there is a need to deposit SiC-based films on the surface of high aspect ratio features at low temperatures. One of the most precise thin film deposition techniques on high-aspect-ratio surfaces that operates at low temperatures is atomic layer deposition (ALD). However, there are currently no known methods for ALD of SiC. Herein, the authors present a first-principles thermodynamic analysis so as to screen different precursor combinations for SiC thin films. The authors do this by calculating the Gibbs energy ΔGΔG of the reaction using density functional theory and including the effects of pressure and temperature. This theoretical model was validated for existing chemical reactions in CVD of SiC at 1000 °C. The precursors disilane (Si2H6), silane (SiH4), or monochlorosilane (SiH3Cl) with ethyne (C2H2), carbontetrachloride (CCl4), or trichloromethane (CHCl3) were predicted to be the most promising for ALD of SiC at 400 °C.

3D Ionospheric Effects on HF Propagation and Heating

The thesis uses a three-dimensional, first-principles model of the ionosphere in combination with High Frequency (HF) raytracing model to address key topics related to the physics of HF propagation and artificial ionospheric heating. In particular: 1. Explores the effect of the ubiquitous electron density gradients caused by Medium Scale Traveling Ionospheric Disturbances (MSTIDs) on high-angle of incidence HF radio wave propagation. Previous studies neglected the all-important presence of horizontal gradients in both the cross- and down-range directions, which refract the HF waves, significantly changing their path through the ionosphere. The physics-based ionosphere model SAMI3/ESF is used to generate a self-consistently evolving MSTID that allows for the examination of the spatio-temporal progression of the HF radio waves in the ionosphere. 2. Tests the potential and determines engineering requirements for ground- based high power HF heaters to trigger and control the evolution of Equatorial Spread F (ESF). Interference from ESF on radio wave propagation through the ionosphere remains a critical issue on HF systems reliability. Artificial HF heating has been shown to create plasma density cavities in the ionosphere similar to those that may trigger ESF bubbles. The work explores whether HF heating may trigger or control ESF bubbles. 3. Uses the combined ionosphere and HF raytracing models to create the first self-consistent HF Heating model. This model is utilized to simulate results from an Arecibo experiment and to provide understanding of the physical mechanism behind observed phenomena. The insights gained provide engineering guidance for new artificial heaters that are being built for use in low to middle latitude regions. In accomplishing the above topics: (i) I generated a model MSTID using the SAMI3/ESF code, and used a raytrace model to examine the effects of the MSTID gradients on radio wave propagation observables; (ii) I implemented a three- dimensional HF heating model in SAMI3/ESF and used the model to determine whether HF heating could artificially generate an ESF bubble; (iii) I created the first self-consistent model for artificial HF heating using the SAMI3/ESF ionosphere model and the MoJo raytrace model and ran a series of simulations that successfully modeled the results of early artificial heating experiments at Arecibo.

Vector Graphics: From PostScript and Flash to SVG

The XML-based specification for Scalable Vector Graphics (SVG), sponsored by the World Wide Web consortium, allows for compact and descriptive vector graphics for the Web. SVG s domain of discourse is that of graphic primitives whose optional attributes express line thickness, fill patterns, text size and so on. These primitives have very different properties from those of traditional document components (e.g. sections, paragraphs etc.) that XML is normally called upon to express. This paper describes a set of three tools for creating SVG, either from first principles or via the conversion of existing formats. The ab initio generation of SVG is effected from a server-side CGI script, using a PERL library of drawing functions; later sections highlight the problems of converting Adobe PostScript and Macromedia s Shockwave format (SWF) into SVG.

A new age of fuel performance code criteria studied through advanced atomistic simulation techniques

A fundamental step in understanding the effects of irradiation on metallic uranium and uranium dioxide ceramic fuels, or any material, must start with the nature of radiation damage on the atomic level. The atomic damage displacement results in a multitude of defects that influence the fuel performance. Nuclear reactions are coupled, in that changing one variable will alter others through feedback. In the field of fuel performance modeling, these difficulties are addressed through the use of empirical models rather than models based on first principles. Empirical models can be used as a predictive code through the careful manipulation of input variables for the limited circumstances that are closely tied to the data used to create the model. While empirical models are efficient and give acceptable results, these results are only applicable within the range of the existing data. This narrow window prevents modeling changes in operating conditions that would invalidate the model as the new operating conditions would not be within the calibration data set. This work is part of a larger effort to correct for this modeling deficiency. Uranium dioxide and metallic uranium fuels are analyzed through a kinetic Monte Carlo code (kMC) as part of an overall effort to generate a stochastic and predictive fuel code. The kMC investigations include sensitivity analysis of point defect concentrations, thermal gradients implemented through a temperature variation mesh-grid, and migration energy values. In this work, fission damage is primarily represented through defects on the oxygen anion sublattice. Results were also compared between the various models. Past studies of kMC point defect migration have not adequately addressed non-standard migration events such as clustering and dissociation of vacancies. As such, the General Utility Lattice Program (GULP) code was utilized to generate new migration energies so that additional non-migration events could be included into kMC code in the future for more comprehensive studies. Defect energies were calculated to generate barrier heights for single vacancy migration, clustering and dissociation of two vacancies, and vacancy migration while under the influence of both an additional oxygen and uranium vacancy.

Improving the use of G-CSF during chemotherapy using physiological mathematical modelling : a quantitative systems pharmacology approach

La diminution des doses administrées ou même la cessation complète d'un traitement chimiothérapeutique est souvent la conséquence de la réduction du nombre de neutrophiles, qui sont les globules blancs les plus fréquents dans le sang. Cette réduction dans le nombre absolu des neutrophiles, aussi connue sous le nom de myélosuppression, est précipitée par les effets létaux non spécifiques des médicaments anti-cancéreux, qui, parallèlement à leur effet thérapeutique, produisent aussi des effets toxiques sur les cellules saines. Dans le but d'atténuer cet impact myélosuppresseur, on administre aux patients un facteur de stimulation des colonies de granulocytes recombinant humain (rhG-CSF), une forme exogène du G-CSF, l'hormone responsable de la stimulation de la production des neutrophiles et de leurs libération dans la circulation sanguine. Bien que les bienfaits d'un traitement prophylactique avec le G-CSF pendant la chimiothérapie soient bien établis, les protocoles d'administration demeurent mal définis et sont fréquemment déterminés ad libitum par les cliniciens. Avec l'optique d'améliorer le dosage thérapeutique et rationaliser l'utilisation du rhG-CSF pendant le traitement chimiothérapeutique, nous avons développé un modèle physiologique du processus de granulopoïèse, qui incorpore les connaissances actuelles de pointe relatives à la production des neutrophiles des cellules souches hématopoïétiques dans la moelle osseuse. À ce modèle physiologique, nous avons intégré des modèles pharmacocinétiques/pharmacodynamiques (PK/PD) de deux médicaments: le PM00104 (Zalypsis®), un médicament anti-cancéreux, et le rhG-CSF (filgrastim). En se servant des principes fondamentaux sous-jacents à la physiologie, nous avons estimé les paramètres de manière exhaustive sans devoir recourir à l'ajustement des données, ce qui nous a permis de prédire des données cliniques provenant de 172 patients soumis au protocol CHOP14 (6 cycles de chimiothérapie avec une période de 14 jours où l'administration du rhG-CSF se fait du jour 4 au jour 13 post-chimiothérapie). En utilisant ce modèle physio-PK/PD, nous avons démontré que le nombre d'administrations du rhG-CSF pourrait être réduit de dix (pratique actuelle) à quatre ou même trois administrations, à condition de retarder le début du traitement prophylactique par le rhG-CSF. Dans un souci d'applicabilité clinique de notre approche de modélisation, nous avons investigué l'impact de la variabilité PK présente dans une population de patients, sur les prédictions du modèle, en intégrant des modèles PK de population (Pop-PK) des deux médicaments. En considérant des cohortes de 500 patients in silico pour chacun des cinq scénarios de variabilité plausibles et en utilisant trois marqueurs cliniques, soient le temps au nadir des neutrophiles, la valeur du nadir, ainsi que l'aire sous la courbe concentration-effet, nous avons établi qu'il n'y avait aucune différence significative dans les prédictions du modèle entre le patient-type et la population. Ceci démontre la robustesse de l'approche que nous avons développée et qui s'apparente à une approche de pharmacologie quantitative des systèmes (QSP). Motivés par l'utilisation du rhG-CSF dans le traitement d'autres maladies, comme des pathologies périodiques telles que la neutropénie cyclique, nous avons ensuite soumis l'étude du modèle au contexte des maladies dynamiques. En mettant en évidence la non validité du paradigme de la rétroaction des cytokines pour l'administration exogène des mimétiques du G-CSF, nous avons développé un modèle physiologique PK/PD novateur comprenant les concentrations libres et liées du G-CSF. Ce nouveau modèle PK a aussi nécessité des changements dans le modèle PD puisqu’il nous a permis de retracer les concentrations du G-CSF lié aux neutrophiles. Nous avons démontré que l'hypothèse sous-jacente de l'équilibre entre la concentration libre et liée, selon la loi d'action de masse, n'est plus valide pour le G-CSF aux concentrations endogènes et mènerait en fait à la surestimation de la clairance rénale du médicament. En procédant ainsi, nous avons réussi à reproduire des données cliniques obtenues dans diverses conditions (l'administration exogène du G-CSF, l'administration du PM00104, CHOP14). Nous avons aussi fourni une explication logique des mécanismes responsables de la réponse physiologique aux deux médicaments. Finalement, afin de mettre en exergue l’approche intégrative en pharmacologie adoptée dans cette thèse, nous avons démontré sa valeur inestimable pour la mise en lumière et la reconstruction des systèmes vivants complexes, en faisant le parallèle avec d’autres disciplines scientifiques telles que la paléontologie et la forensique, où une approche semblable a largement fait ses preuves. Nous avons aussi discuté du potentiel de la pharmacologie quantitative des systèmes appliquées au développement du médicament et à la médecine translationnelle, en se servant du modèle physio-PK/PD que nous avons mis au point.

On the structure of the Au18(SR)14 cluster

First principles calculations are used for a systematic search of the lowest-energy (most-stable) structure of the recently synthesized Au18(SR)14 cluster. A comparison of the calculated optical absorption and electronic circular dichroism spectra, which are highly sensitive to the cluster structure and chirality, with the experimental spectra of the glutathione-protected gold cluster, Au18(SG)14, is used to discriminate between low-energy isomers of the Au18(SR)14 (R = CH3) cluster. From the good agreement between calculated and measured spectra, it is predicted that the structure of the Au18(SR)14 cluster consists of a prolate Au8 core covered with two dimer (SR–Au–SR–Au–SR) and two trimer (SR–Au–SR–Au–SR–Au–SR) motifs. These results provide additional evidence on the existence of longer trimer motifs as protecting units of small thiolated gold clusters.

A THEORETICAL STUDY OF INTERACTION OF NANOPARTICLES WITH BIOMOLECULE

Many types of materials at nanoscale are currently being used in everyday life. The production and use of such products based on engineered nanomaterials have raised concerns of the possible risks and hazards associated with these nanomaterials. In order to evaluate and gain a better understanding of their effects on living organisms, we have performed first-principles quantum mechanical calculations and molecular dynamics simulations. Specifically, we will investigate the interaction of nanomaterials including semiconducting quantum dots and metallic nanoparticles with various biological molecules, such as dopamine, DNA nucleobases and lipid membranes. Firstly, interactions of semiconducting CdSe/CdS quantum dots (QDs) with the dopamine and the DNA nucleobase molecules are investigated using similar quantum mechanical approach to the one used for the metallic nanoparticles. A variety of interaction sites are explored. Our results show that small-sized Cd4Se4 and Cd4S4 QDs interact strongly with the DNA nucleobase if a DNA nucleobase has the amide or hydroxyl chemical group. These results indicate that these QDs are suitable for detecting subcellular structures, as also reported by experiments. The next two chapters describe a preparation required for the simulation of nanoparticles interacting with membranes leading to accurate structure models for the membranes. We develop a method for the molecular crystalline structure prediction of 1,2-Dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC), 1,2-Dimyristoyl-sn-glycero-3-phosphorylethanolamine (DMPE) and cyclic di-amino acid peptide using first-principles methods. Since an accurate determination of the structure of an organic crystal is usually an extremely difficult task due to availability of the large number of its conformers, we propose a new computational scheme by applying knowledge of symmetry, structural chemistry and chemical bonding to reduce the sampling size of the conformation space. The interaction of metal nanoparticles with cell membranes is finally carried out by molecular dynamics simulations, and the results are reported in the last chapter. A new force field is developed which accurately describes the interaction forces between the clusters representing small-sized metal nanoparticles and the lipid bilayer molecules. The permeation of nanoparticles into the cell membrane is analyzed together with the RMSD values of the membrane modeled by a lipid bilayer. The simulation results suggest that the AgNPs could cause the same amount of deformation as the AuNPs for the dysfunction of the membrane.

Microwave-assisted synthesis and local analyses of positive insertion electrodes for Li-ion batteries

Efficient energy storage holds the key to reducing waste energy and enabling the use of advanced handheld electronic devices, hydrid electric vehicles and residential energy storage. Recently, Li-ion batteries have been identified and employed as energy storage devices due to their high gravimetric and volumetric energy densities, in comparison to previous technologies. However, more research is required to enhance the efficiency of Li-ion batteries by discovering electrodes with larger electrochemical discharge capacities, while maintaining electrochemical stability. The aims of this study are to develop new microwave-assisted synthesis routes to nanostructured insertion cathodes, which harbor a greater affinity for lithium extraction and insertion than bulk materials. Subsequent to this, state-of-the-art synchrotron based techniques have been employed to understand structural and dynamic behaviour of nanostructured cathode materials during battery cell operation. In this study, microwave-assisted routes to a-LiFePO4, VO2(B), V3O7, H2V3O8 and V4O6(OH)4 have all been developed. Muon spin relaxation has shown that the presence of b-LiFePO4 has a detrimental effect on the lithium diffusion properties of a-LiFePO4, in agreement with first principles calculations. For the first time, a-LiFePO4 nanostructures have been obtained by employing a deep eutectic solvent reaction media showing near theoretical capacity (162 mAh g–1). Studies on VO2(B) have shown that the discharge capacity obtained is linked to the synthesis method. Electrochemical studies of H2V3O8 nanowires have shown outstanding discharge capacities (323 mAh g–1 at 100 mA g–1) and rate capability (180 mAh g–1 at 1 A g–1). The electrochemcial properties of V4O6(OH)4 have been investigated for the first time and show a promising discharge capacity of (180 mAh g–1). Lastly, in situ X-ray absorption spectroscopy has been utilised to track the evolution of the oxidation states in a-LiFePO4, VO2(B) and H2V3O8, and has shown these can all be observed dynamically.

Electronic structure and magnetism of Mn-doped GaSb for spintronic applications: A DFT study

We have carried out first-principles spin polarized calculations to obtain comprehensive information regarding the structural, magnetic, and electronic properties of the Mn-doped GaSb compound with dopant concentrations: x¼0.062, 0.083, 0.125, 0.25, and 0.50. The plane-wave pseudopotential method was used in order to calculate total energies and electronic structures. It was found that the MnGa substitution is the most stable configuration with a formation energy of 1.60 eV/Mn-atom. The calculated density of states shows that the half-metallic ferromagnetism is energetically stable for all dopant concentrations with a total magnetization of about 4.0 lB/Mn-atom. The results indicate that the magnetic ground state originates from the strong hybridization between Mn-d and Sb-p states, which agree with previous studies on Mn-doped wide gap semiconductors. This study gives new clues to the fabrication of diluted magnetic semiconductors

Computing the structural and vibrational properties of polymorphic organic molecular crystals through van der Waals corrected density functional theory and the electronic properties of organic thin films through microelectrostatic calculations.

The present Thesis reports on the various research projects to which I have contributed during my PhD period, working with several research groups, and whose results have been communicated in a number of scientific publications. The main focus of my research activity was to learn, test, exploit and extend the recently developed vdW-DFT (van der Waals corrected Density Functional Theory) methods for computing the structural, vibrational and electronic properties of ordered molecular crystals from first principles. A secondary, and more recent, research activity has been the analysis with microelectrostatic methods of Molecular Dynamics (MD) simulations of disordered molecular systems. While only very unreliable methods based on empirical models were practically usable until a few years ago, accurate calculations of the crystal energy are now possible, thanks to very fast modern computers and to the excellent performance of the best vdW-DFT methods. Accurate energies are particularly important for describing organic molecular solids, since they often exhibit several alternative crystal structures (polymorphs), with very different packing arrangements but very small energy differences. Standard DFT methods do not describe the long-range electron correlations which give rise to the vdW interactions. Although weak, these interactions are extremely sensitive to the packing arrangement, and neglecting them used to be a problem. The calculations of reliable crystal structures and vibrational frequencies has been made possible only recently, thanks to development of some good representations of the vdW contribution to the energy (known as “vdW corrections”).