Publications

2011
Katono M, Bessho T, Meng S, Humphry-Baker R, Rothenberger G, Zakeeruddin SM, Kaxiras E, Graetzel M. D-π-A Dye System Containing Cyano-Benzoic Acid as Anchoring Group for Dye-Sensitized Solar Cells. Langmuir. 2011;27 :14248-14252.Abstract

A D-π-A dye (KM-1) incorporating cyanobenzoic acid as a new acceptor/anchoring group has been synthesized for dye-sensitized solar cells (DSCs) with a high molar extinction coefficient of 66 700 M-1 cm-1 at 437 nm. Theoretical calculations show that the hydrogen bond between CN and surface hydroxyl leads to the most stable configuration on the surface of TiO2. In addition, the adsorption of the dye on TiO2 follows a BrunauerEmmettTeller (BET) isotherm. Multilayer adsorption of KM-1 on TiO2 seems to take place particularly at higher dye concentrations. DSC device using KM-1 reached a maximum incident photon-to-current conversion efficiency (IPCE) of 84%, with a solar to electric power conversion efficiency (PCE) of 3.3% at AM1.5 G illumination (100 mW cm-2 ). This new type of anchoring group paves a way to design new dyes that combine good visible light harvesting with strong binding to the metal oxide surface

pdf
Melchionna S, Kaxiras E, Bernaschi M, Succi S. Endothelial shear stress from large-scale blood flow simulations. Philosophical Transactions of the Royal Society A. 2011;369 :2354-2361.Abstract

We discuss the optimal evaluation of endothelial shear stress for real-life case studies based on anatomic data acquisition. The fluid dynamic simulations require smoothing of the geometric dataset to avoid major artefacts in the flow patterns, especially in the proximity of bifurcations. A systematic series of simulations at different corrugation levels shows that, below a smoothing length of about 0.5 mm, the numerical data are insensitive to further smoothing.

pdf
Mao Y, Wang WL, Wei D, Kaxiras E, Sodroski JG. Graphene Structures at an Extreme Degree of Buckling. ACS Nano. 2011;5 (2) :1395-1400.Abstract

T The distinctive properties of graphene sheets may be significantly influenced by the presence of corrugation structures. Our understanding of these graphene structures has been limited to the mesoscopic scale. Here we characterize angstrom-scale periodic buckling structures in freestanding graphene bilayers produced by liquid-phase processing in the absence of specific substrates. Monochromated, aberration-corrected transmission electron microscopy with sub-angstrom resolution revealed that the unit structures in the major buckling direction consist of only two and three unit cells of graphene's honeycomb lattice, resulting in buckling wavelengths of 3.6±0.5 and 6.4±0.8 Å, respectively. The buckling shows a strong preference of chiral direction and spontaneously chooses the orientation of the lowest deformation energy, governed by simple geometry rules agreeing with Euler buckling theory. Unexpectedly, the overall buckled structures demonstrate geometric complexity with cascaded features. First-principles calculations suggest that significant anisotropic changes in the electronic structure of graphene are induced by the buckling.

pdf
Knaup JM, Li H, Vlassak JJ, Kaxiras E. Influence of CH2 content and network defects on the elastic properties of organosilicate glasses. Physical Review B. 2011;83 :054204.Abstract

Organosilicate glasses (OSGs) are used as low-k intermetal dielectrics for advanced integrated circuits. In this application, the material must fulfill two conflicting requirements: It has to have low density to reduce the dielectric constant while being sufficiently mechanically stable to withstand thermomechanical and other stresses during subsequent steps of integrated circuit manufacture. Recent experimental advances in improving the mechanical and electrical properties of these materials have not yet been systematically studied theoretically at the ab initio level due to the large model sizes necessary to realistically describe amorphous materials. In this paper we employ the density-functional-based tight-binding method to achieve an accurate description of OSG properties at different compositions. We analyze the influence of composition and local network defects on the density and bulk modulus of nonporous OSG. We find that the dependence of density and that of mechanical stiffness on chemical composition are of different natures. This difference is traced to a transition between mechanisms of elastic deformation in silica glass and in silicon hydrocarbide, which is also the reason for the two materials’ different sensitivities to network defects.

pdf
Zhao K, Wang WL, Gregoire J, Pharr M, Suo Z, Vlassak JJ, Kaxiras E. Lithium-Assisted Plastic Deformation of Silicon Electrodes in Lithium-Ion Batteries: A First-Principles Theoretical Study. NanoLetters. 2011;11 :2962-2967.Abstract

Silicon can host a large amount of lithium, making it a promising electrode for high-capacity lithium-ion batteries. Recent experiments indicate that silicon experiences large plastic deformation upon Li absorption, which can significantly decrease the stresses induced by lithiation and thus mitigate fracture failure of electrodes. These issues become especially relevant in nanostructured electrodes with confined geometries. On the basis of first-principles calculations, we present a study of the microscopic deformation mechanism of lithiated silicon at relatively low Li concentration, which captures the onset of plasticity induced by lithiation. We find that lithium insertion leads to breaking of SiSi bonds and formation of weaker bonds between neighboring Si and Li atoms, which results in a decrease in Young’s modulus, a reduction in strength, and a brittle-to-ductile transition with increasing Li concentration. The microscopic mechanism of large plastic deformation is attributed to continuous lithium-assisted breaking and re-forming of SiSi bonds and the creation of nanopores.

pdf
Kuehne TD, Pascal TA, Kaxiras E, Jung Y. New Insights into the Structure of the Vapor/Water Interface from Large-Scale First-Principles Simulations. Journal of Physical Chemistry Letters. 2011;2 :105-113.Abstract

We present extensive ab initio simulations of the molecular arrangements at the vapor/water interface, which provide valuable insights into the interface structure. In particular, the simulations address the controversy of whether there is a significant amount of nondonor configurations at this prototypical interface, using a novel Car-Parrinello-like ab initio molecular dynamics approach. The interface is modeled by a system of 384 water molecules for 125 ps in a two-dimensional periodic slab, the most extensive ab initio molecular dynamics simulation to date. In contrast to previous theoretical simulations and X-ray absorption spectroscopy, but consistent with sum-frequency generation experiments, we observe no evidence for a significant occurrence of acceptoronly species at the vapor/water interface. Besides a distinct surface relaxation effect, we find that only the topmost layers of the interface obey structural order.

pdf
Ren J, Meng S, Wang Y-L, Ma X-C, Xue Q-K, Kaxiras E. Properties of copper (fluoro-)phthalocyanine layers deposited on epitaxial graphene. Journal of Chemical Physics. 2011;134 :194706.Abstract

We investigate the atomic structure and electronic properties of monolayers of copper phthalocyanines (CuPc) deposited on epitaxial graphene substrate. We focus in particular on hexadecafluorophthalocyanine (F16CuPc), using both theoretical and experimental (scanning tunneling microscopy – STM) studies. For the individual CuPc and F16CuPc molecules, we calculated the electronic and optical properties using density functional theory (DFT) and time-dependent DFT and found a red-shift in the absorption peaks of F16CuPc relative to those of CuPc. In F16CuPc, the electronic wavefunctions are more polarized toward the electronegative fluorine atoms and away from the Cu atom at the center of the molecule. When adsorbed on graphene, the molecules lie flat and form closely packed patterns: F16CuPc forms a hexagonal pattern with two well-ordered alternating α and β stripes while CuPc arranges into a square lattice. The competition between molecule-substrate and intermolecular van der Waals interactions plays a crucial role in establishing the molecular patterns leading to tunable electron transfer from graphene to the molecules. This transfer is controlled by the layer thickness of, or the applied voltage on, epitaxial graphene resulting in selective F16CuPc adsorption, as observed in STM experiments. In addition, phthalocyanine adsorption modifies the electronic structure of the underlying graphene substrate introducing intensity smoothing in the range of 2–3 eV below the Dirac point (ED) and a small peak in the density of states at ∼0.4 eV above ED.

pdf
Maze JR, Gali A, Togan E, Chu Y, Trifonov A, Kaxiras E, Lukin MD. Properties of nitrogen-vacancy centers in diamond: the group theoretic approach. New Journal of Physics. 2011;13 :025025.Abstract

We present a procedure that makes use of group theory to analyze and predict the main properties of the negatively charged nitrogen-vacancy (NV) center in diamond. We focus on the relatively low temperature limit where both the spin–spin and spin–orbit effects are important to consider. We demonstrate that group theory may be used to clarify several aspects of the NV structure, such as ordering of the singlets in the (e 2 ) electronic configuration and the spin–spin and spin–orbit interactions in the (ae) electronic configuration. We also discuss how the optical selection rules and the response of the center to electric field can be used for spin–photon entanglement schemes. Our general formalism is applicable to a broad class of local defects in solids. The present results have important implications for applications in quantum information science and nanomagnetometry.

newjphys_13_025025_2011.pdf
Baker TA, Xu B, Jensen SC, Friend CM, Kaxiras E. Role of defects in propene adsorption and reaction on a partially O-covered Au(111) surface. Catalysis Science & Technology. 2011;1 :1166-1174.Abstract

We investigate the role of defects—adatoms, vacancies, and steps—in the bonding and reaction of propene on Au(111) containing atomic oxygen, using density functional theory (DFT) calculations. The adsorption of propene is stronger on a surface containing defects compared to the flat, bulk-terminated surface, with the largest gain in binding (~0.7 eV) on a surface with a 1/9 monolayer (ML) of Au adatoms. Charge-density difference plots reveal that the difference between defective surfaces and the bulk-terminated surface is a more pronounced depletion of electron density from the carbon–carbon p bond and a charge accumulation between the double bond and the gold atom to which the propene is bound. We calculate the energy barriers for two competing reactions that are important in determining the selectivity for propene oxidation. Allylic H abstraction by adsorbed O leads to combustion, whereas O addition to form an oxymetallacycle is the first step in propene epoxide formation. A comparison of the energetics of these two pathways on flat and defect-containing Au surfaces indicates that the reactivity depends on the nature and prevalence of surface defects. Both electronic and geometric factors, such as the path and the distance the oxygen must travel to meet the allylic hydrogen for abstraction, are important in explaining the reaction trends.

pdf
Haubrich J, Kaxiras E, Friend CM. The Role of Surface and Subsurface Point Defects for Chemical Model Studies on TiO2: A First-Principles Theoretical Study of Formaldehyde Bonding on Rutile TiO2(110). Chemistry - a European Journal. 2011;17 :4496-4506.Abstract

We report a systematic investigation of the effects of different surface and subsurface point defects on the adsorption of formaldehyde on rutile TiO2(110) surfaces using density functional theory (DFT). All point defects investigated—including surface bridging oxygen vacancies, titanium interstitials, and subsurface oxygen vacancies—stabilize the adsorption significantly by up to 56 kJ mol-1 at a coverage of 0.1 monolayer (ML). The stabilization is due to a decrease of the coordination (covalent saturation) of the surface Ti adsorption sites adjacent to the defects, which leads to a stronger molecule–surface interaction. This change in the Ti is caused by the removal of a neighboring atom (oxygen vacancies) or substantial lattice relaxations induced by the subsurface defects. On the stoichiometric reference surface, the most stable adsorption geometry of formaldehyde is a tilted η2 -dioxymethylene (with an adsorption energy Eads=-125 kJ mol-1 ), in which a bond forms to a nearby bridging O atom and the carbonyl-O atom in the formaldehyde binds to a Ti atom in the adjacent fivefold coordinated lattice site. The η1 -top configuration on fivecoordinate Ti4+ is much less favorable (Eads=-69 kJ mol-1 ). The largest stabilization is exerted by subsurface Ti interstitials between the first and second layers. These defects stabilize the η2 -dioxymethylene structure by nearly 40 kJ mol-1 to an adsorption energy of 164 kJ mol-1 . Contrary to popular belief, adsorption in a bridging oxygen vacancy (Eads=-86 kJ mol-1 ) is much less favorable for formaldehyde compared to the η2 - dioxymethylene structures. From these results we conclude that formaldehyde will bind in the η2 - dioxymethylene structure on the stoichiometric surface as well as in the presence of Ti interstitials and bridging oxygen vacancies. In the light of these substantial effects, we conclude that it is essential to include all the types of point defects present in typical, reduced rutile samples used for model studies, at realistic concentrations to obtain correct adsorption sites, structures, energetic, and chemi-physical properties. 

pdf
Li H, Knaup JM, Kaxiras E, Vlassak JJ. Stiffening of organosilicate glasses by organic cross-linking. Acta Materialia. 2011;59 :44-52.Abstract

Atomistic simulations show that organosilicates, used as low permittivity dielectric materials in advanced integrated circuits, can be made substantially stiffer than amorphous silica, while maintaining a lower mass density. The enhanced stiffness is achieved by incorporating organic cross-links to replace bridging oxygen atoms in the silica network. To elucidate the mechanism responsible for the enhanced stiffness, the conformational changes in the network upon hydrostatic and shear loading are examined. The structural and mechanical impact of terminal methyl groups is also assessed quantitatively and compared with continuous random network theory.

actamater_59_44-52_2011.pdf
Xu B, Haubrich J, Baker TA, Kaxiras E, Friend CM. Theoretical Study of O-Assisted Selective Coupling of Methanol on Au(111). Journal of Physical Chemistry C. 2011;115 :3703-3708.Abstract

We report the first systematic theoretical study of the oxidative self-coupling of methanol to form the ester, methylformate, on atomic-oxygen-covered Au(111) using density functional theory calculations. The first step in the process —dissociation of the O-H bond in methanol—has a lower barrier for transfer of the proton to adsorbed oxygen than for transfer of H to gold, consistent with experimental observations that O is necessary to initiate the reaction. The computed barrier for formation of methoxy (CH3O) and OH is 0.41 eV, compared with 1.58 eV calculated for the transfer of H to the clean Au surface. Several different pathways for the ensuing β-H elimination in CH3O(ads) to form formaldehyde have been considered, namely, attack by adsorbed O, OH, or a second CH3O, and transfer to the Au metal. Methoxy attacked by surface oxygen has the lowest calculated barrier, 0.49 eV, and leads to adsorbed H2C=O and OH. Subsequent coupling of methoxy and formaldehyde has no apparent barrier in the calculation, consistent with the experimental conclusion that β-H elimination is the rate-limiting step for the overall reaction. With the exception of surface oxygen, all other surface species have low diffusion barriers, suggesting that rearrangement and movement of these species from the preferred adsorption sites to configurations necessary for reactions occur readily, thus contributing to the activity for coupling on gold.

pdf
Lin M-W, Ling C, Agapito LA, Kioussis N, Zhang Y, Cheng MM-C, Wang WL, Kaxiras E, Zhou Z. Approaching the intrinsic band gap in suspended high-mobility graphene nanoribbons. Physical Review B. 2011;84 :125411.Abstract

We report electrical transport measurements on a suspended ultra-low-disorder graphene nanoribbon (GNR) with nearly atomically smooth edges that reveal a high mobility exceeding 3000 cm2 V−1 s−1 and an intrinsic band gap. The experimentally derived band gap is in quantitative agreement with the results of our electronic structure calculations on chiral GNRs with comparable width taking into account the electron-electron interactions, indicating that the origin of the band gap in nonarmchair GNRs is partially due to the magnetic zigzag edges

pdf
Gali A, Kaxiras E, Zimanyi GT, Meng S. Effect of symmetry breaking on the optical absorption of semiconductor nanoparticles. Physical Review B. 2011;84 :035325.Abstract

The shape of semiconductor nanoparticles (NPs) can have an important influence on their optical absorption spectra with enhanced absorption resulting from breaking their symmetries. To illustrate this effect, we present broad energy optical excitation spectra of ∼2 nm Si and PbSe nanoparticles, obtained from extensive timedependent density functional theory calculations. We considered both highly symmetric spherical shapes and low-symmetry rodlike and disklike shapes. The low-symmetry shapes exhibit an increase in absorption at low and relatively high energies compared to the absorption of spherical NPs of similar volume, independent of their chemical compositions and surface structures. Our results elucidate the mechanism of enhanced multiexciton generation in semiconductor NPs, which is important in the quest to improve their photovoltaic applications.

pdf
2010
Chang CL, Sankaranarayanan SKRS, Ruzmetov D, Engelhard MH, Kaxiras E, Ramanathan S. Compositional tuning of ultrathin surface oxides on metal and alloysubstrates using photons: Dynamic simulations and experiments. Phys. Rev. B. 2010;81 :085406.
Meng S, Kaxiras E. Electron and Hole Dynamics in Dye-Sensitized Solar Cells:Influencing Factors and Systematic Trends. NanoLett. 2010;10 :1238.
Chen HP, Kalia RK, Kaxiras E, Lu G, Nakano A, Nomura K, van Duin ACT, Vashishta P, Yuan ZS. Embrittlement of Metal by Solute Segregation-Induced Amorphization. Phys. Rev. Lett. 2010;104 :155502.
Bernaschi M, Fatica M, Melchionna S, Succi S, Kaxiras E. A flexible high-performance Lattice Boltzmann GPU code for thesimulations of fluid flows in complex geometries. Concurrency Computat.: Pract. Exper. 2010;22 :1.
Melchionna S, Bernaschi M, Succi S, Kaxiras E, Rybicki FJ, Mitsouras D, Coskun AU, Feldman CL. Hydrokinetic approach to large-scale cardiovascular blood flow. Comp. Phys. Comm. 2010;181 :462.
Baker TA, Friend CM, Kaxiras E. Local Bonding Effects in the Oxidation of CO on Oxygen-CoveredAu(111) from Ab Initio Molecular Dynamics Simulations. J. Chem. Theory Comput. 2010;6 :279.

Pages