Search

In this work, a set of density-functional tight-binding (DFTB) parameters for the Zr–Zr, Zr–O, Y–Y, Y–O, and Zr–Y interactions was developed for bulk and surface simulations of ZrO2 (zirconia), Y2O3 (yttria), and yttria-stabilized zirconia (YSZ) materials. The parameterization lays the ground work for realistic simulations of zirconia-, yttria-, and YSZ-based electrolytes in solid oxide fuel cells and YSZ-based catalysts on long timescales and relevant size scales. The parameterization was validated for the zirconia and yttria polymorphs observed under standard conditions based on density functional theory calculations and experimental data. Additionally, we performed DFTB-based molecular dynamics (MD) simulations to compute structural and vibrational properties of these materials. The results show that the parameters can give a qualitatively correct phase ordering of zirconia, where the tetragonal phase is more stable than the cubic phase at a lower temperature. The lattice parameters are only slightly overestimated by 0.05–0.1 Å (2% error), still within the typical accuracy of first-principles methods. Additionally, the MD results confirm that zirconia and yttria phases are stable against transformations under standard conditions. The parameterization also predicts that vibrational spectra are within the range of 100–1000 cm–1 for zirconia and 100–800 cm–1 for yttria, which is in good agreement with predictions both from full quantum mechanics and a recently developed classical force field. To further demonstrate the advantage of the developed DFTB parameters in terms of computational resources, we conducted DFTB/MD simulations of the YSZ4 and YS12 models containing approximately 750 atoms.

Ceria (CeO2) is a promising metal-oxide support that is used in three-way catalysis (TWC). The activity of ceria-supported TWC depends on the location and concentration of oxygen vacancies. Oxygen ...

Carbonaceous or oxy-carbon species are intermediates formed during C x H y combustion on a Pt n /Al2O3 catalyst, which contain carbon, hydrogen, and oxygen atoms. The accumulation of the carbonaceous species, arguably, leads to catalytic deactivation; therefore, their removal is of importance. As the diffusion process is occasionally the rate-determining step in the growth of carbonaceous species, the present study aims to reveal the diffusion mechanisms. The free energy barriers of acetate, formate, and methoxy diffusion on the (100)-γ-Al2O3 surface were evaluated through extensive metadynamics simulations at the density-functional tight-binding level. The present work deduces that each adopted carbonaceous species exhibits different diffusion mechanisms and supports experimental evidence that the acetate species exhibits the slowest diffusivity among the adopted carbonaceous species.

Water confined within one-dimensional (1D) hydrophobic nanochannels has attracted significant interest due to its unusual structure and dynamic properties. As a representative system, water-filled carbon nanotubes (CNTs) are generally studied, but direct observation of the crystal structure and proton transport is difficult for CNTs due to their poor crystallinity and high electron conduction. Here, we report the direct observation of a unique water-cluster structure and high proton conduction realized in a metal-organic nanotube, [Pt(dach)(bpy)Br]4(SO4)4·32H2O (dach: (1R, 2R)-(–)-1, 2-diaminocyclohexane; bpy: 4, 4’-bipyridine). In the crystalline state, a hydrogen-bonded ice nanotube composed of water tetramers and octamers is found within the hydrophobic nanochannel. Single-crystal impedance measurements along the channel direction reveal a high proton conduction of 10−2 Scm−1. Moreover, fast proton diffusion and continuous liquid-to-solid transition are confirmed using solid-state 1H-NMR measurements. Our study provides valuable insight into the structural and dynamical properties of confined water within 1D hydrophobic nanochannels., ナノサイズの空間に閉じ込められた水が示す不思議な性質を実証 --高速で水素イオンを運ぶ水の状態を解明--. 京都大学プレスリリース. 2020-02-18.

Density-functional tight-binding (DFTB) parameters are presented for the simulation of the bulk phases of zirconium. Electronic parameters were obtained using a band structure fitting strategy, while two-center repulsive potentials were created by particle swarm optimization. As objective functions for the repulsive potential fitting, we employed the Birch-Murnaghan equations of state for hexagonal close-packed (HCP), body-centered cubic (BCC) and ω phases of Zr from density-functional theory (DFT). When fractional atomic coordinates are not allowed to change in the generation of the equation-of-state curves, long-range repulsive DFTB potentials are able to almost perfectly reproduce equilibrium structures, relative DFT energies of the bulk phases, and bulk moduli. However, the same potentials lead to artifacts in the DFTB potential energy surfaces when atom positions in the unit cell are allowed to fully relax during the change of unit cell parameters. Conventional short-range repulsive DFTB potentials, while inferior in their ability to reproduce DFT bulk energetics, are able to correctly reproduce the qualitative shape of the DFT potential energy surfaces, including the location of global minima, and can therefore be considered more transferable.

We report a closed-form formula for the Zhang&ndash, Zhang polynomial (also known as ZZ polynomial or Clar covering polynomial) of an important class of elementary peri-condensed benzenoids Rbn1,n2,m1,m2, usually referred to as ribbons. A straightforward derivation is based on the recently developed interface theory of benzenoids [Langner and Witek, MATCH Commun. Math. Comput. Chem.2020, 84, 143&ndash, 176 ]. The discovered formula provides compact expressions for various topological invariants of Rbn1,n2,m1,m2: the number of Kekul&eacute, structures, the number of Clar covers, its Clar number, and the number of Clar structures. The last two classes of elementary benzenoids, for which closed-form ZZ polynomial formulas remain to be found, are hexagonal flakes Ok,m,n and oblate rectangles Orm,n.

The density-functional tight-binding (DFTB) method is one of the useful quantum chemical methods, which provides a good balance between accuracy and computational efficiency. In this account, we reviewed the basis of the DFTB method, the linear-scaling divide-and-conquer (DC) technique, as well as the parameterization process. We also provide some refinement, modifications, and extension of the existing parameters that can be applicable for lithium-ion battery systems. The diffusion constants of common electrolyte molecules and LiTFSA salt in solution have been estimated using DC-DFTB molecular dynamics simulation with our new parameters. The resulting diffusion constants have good agreement to the experimental diffusion constants.

We investigated the atomistic and dynamical mechanism of polaron formation in methylammonium lead iodide perovskite (MAPbI3), which is a representative perovskite solar cell absorber, through the quantum mechanical molecular dynamics simulations. The simulations were conducted on the spatial scale of several nanometres, which can describe charge localization and the associated structural deformation, using the divide-and-conquer-type density-functional tight-binding method, which enables a quantum chemical treatment of systems comprising thousands of atoms. We found that both the structural parts of MAPbI3, namely, the inorganic framework (PbI3−) and the MA cations, involve the structural deformation associated with polaron formation. We elucidated that in the process of polaron formation, charge localization is invoked by thermal structural fluctuation, and a further structural deformation is caused by the relaxation of the charge carrier. Finally, importance of the two structural parts, PbI3− and MA, was examined from the energetical viewpoint.

Superconcentrated electrolyte solutions are receiving increasing attention as a novel class of liquid electrolyte for secondary batteries because of their unusual and favorable characteristics, which arise from a unique solution structure with a very small number of free solvent molecules. The present theoretical study investigates the concentration dependence of the structural and dynamical properties of these electrolyte solutions for Na-ion batteries using large-scale quantum molecular dynamics simulations. Microscopic analysis of the dynamical properties of Na+ ions reveals that ligand (solvent/anion) exchange reactions, an alternative diffusion pathway for Na+ ions, are responsible for carrier ion diffusion in the superconcentrated conditions.

The structural, dynamical, and energetic properties of the excess proton in ice were studied using density-functional tight-binding molecular dynamics simulations. The ice systems investigated herein consisted of low-density hexagonal and cubic crystalline variants (ice Ih and Ic) and high-density structures (ice III and melted ice VI). Analysis of the temperature dependence of radial distribution function and bond order parameters served to characterize the distribution and configuration of hundreds of water molecules in a unit cell. We confirmed that ice Ih and Ic possess higher hexagonal symmetries than ice III and melted ice VI. The estimated Grotthuss shuttling diffusion coefficients in ice were larger than that of liquid water, indicating a slower proton diffusion process in high-density structures than in low-density systems. The energy barriers calculated on the basis of the Arrhenius plot of diffusion coefficients were in reasonable agreement with experimental measurement for ice Ih. Furthermore,...

Aqueous Na- or K-ion batteries could virtually eliminate the safety and cost concerns raised from Li-ion batteries, but their widespread applications have generally suffered from narrow electrochemical potential window (ca. 1.23 V) of aqueous electrolytes that leads to low energy density. Herein, by exploring optimized eutectic systems of Na and K salts with asymmetric imide anions, we discovered, for the first time, room-temperature hydrate melts for Na and K systems, which are the second and third alkali metal hydrate melts reported since the first discovery of Li hydrate melt by our group in 2016. The newly discovered Na- and K- hydrate melts could significantly extend the potential window up to 2.7 and 2.5 V (at Pt electrode), respectively, owing to the merit that almost all water molecules participate in the Na+ or K+ hydration shells. As a proof-of-concept, a prototype Na3 V2 (PO4 )2 F3 |NaTi2 (PO4 )3 aqueous Na-ion full-cell with the Na-hydrate-melt electrolyte delivers an average discharge voltage of 1.75 V, that is among the highest value ever reported for all aqueous Na-ion batteries.

Alkaline metals are an ideal negative electrode for rechargeable batteries. Forming a fluorine‐rich interphase by a fluorinated electrolyte is recognized as key to utilizing lithium metal electrodes, and the same strategy is being applied to sodium metal electrodes. However, their reversible plating/stripping reactions have yet to be achieved. Herein, we report a contrary concept of fluorine‐free electrolytes for sodium metal batteries. A sodium tetraphenylborate/monoglyme electrolyte enables reversible sodium plating/stripping at an average Coulombic efficiency of 99.85 % over 300 cycles. Importantly, the interphase is composed mainly of carbon, oxygen, and sodium elements with a negligible presence of fluorine, but it has both high stability and extremely low resistance. This work suggests a new direction for stabilizing sodium metal electrodes via fluorine‐free interphases.

We show that the Zhang-Zhang (ZZ) polynomial of a benzenoid obtained by fusing a parallelogram M ( m , n ) with an arbitrary benzenoid structure A B C can be simply computed as a product of the ZZ polynomials of both fragments. It seems possible to extend this important result also to cases where both fused structures are arbitrary Kekulean benzenoids. Formal proofs of explicit forms of the ZZ polynomials for prolate rectangles P r ( m , n ) and generalized prolate rectangles P r ( m 1 , m 2 , ? , m n , n ) follow as a straightforward application of the general theory, giving Z Z ( P r ( m , n ) , x ) = ( 1 + ( 1 + x ) ? m ) n and Z Z ( P r ( m 1 , m 2 , ? , m n , n ) , x ) = ? k = 1 n ( 1 + ( 1 + x ) ? m k ) .

We present a computational study of convergence properties of vibrational IR and Raman spectra for a series of increasingly long units of polyethylene, cis- and trans-polyacetylenes, and polyynes. Convergent behavior to the spectra of infinitely long polymers was observed in all cases when chains reached lengths of approximately 60 carbon atoms, both with respect to the positions of the bands and to their intensities. The vibrational spectra of longer chains are practically indistinguishable. The convergence rate depends on the degree of the π conjugation in a studied system: Vibrational spectra for oligoethylenes converge noticeably faster than the spectra for the conjugated systems. The slowest convergence is observed for skeletal motions of the oligomer chains, which may require more than a hundred carbon atoms in the chain to show deviations smaller than 1 cm−1 to the corresponding solid-state calculations. The results suggest that the boundary between the properties of finite and infinite molecular systems fades away for a surprisingly small number of atoms.

The structural, dynamical, and energetic properties of the excess proton in ice were studied using density-functional tight-binding molecular dynamics simulations. The ice systems investigated herein consisted of low-density hexagonal and cubic crystalline variants (ice I

We report the preliminary results of the development of density-functional tight-binding (DFTB) repulsive potentials for the Zr – O element pair. The repulsive potentials were created using a computer code based on the particle swarm optimization. The potentials were tested on a set of systems for high temperature phases of bulk ZrO2, namely cubic and tetragonal. The potential sets were primarily developed for simulation of zirconia phase transitions at elevated temperatures.

We present an efficient numerical integration scheme (TWOCENT) to be used in the context of automatized parameterization of the density-functional tight-binding (DFTB) method. The accuracy of the integration process is assessed and its range of applicability is discussed. The functionality of the developed code is tested by reproducing the electronic portion of the existing mio parameter sets and by reproducing a series of reference DFT band structures of elemental solids.

In this work, we derive and test a new automatized strategy to construct repulsive potentials for the self-consistent charge density functional tight-binding (SCC-DFTB) method. This approach allows one to explore the parameter space in a systematic fashion in order to find optimal solutions. We find that due to the limited flexibility of the SCC-DFTB electronic part, not all properties can be optimized simultaneously. For example, the optimization of heats of formation is in conflict with the optimization of vibrational frequencies. Therefore, a special parametrization for vibrational frequencies is derived. It is shown that the performance of SCC-DFTB can be significantly improved using a more elaborate fitting strategy. A new fit for C and H is presented, which results in an average error of 2.6 kcal/mol for heats of formations for a large set of hydrocarbons, indicating that the performance of SCC-DFTB can be systematically improved also for other elements.

Variational Transition State Theory with Multidimensional Tunneling (VTST/MT) has been successfully used for calculating rate constants of reactions in gas and condensed phases. The current software implementation of VTST/MT is, however, based on the assumption of a fast, serial evaluation of the energetic information of a given molecular structure. We propose a simple and effective parallel method for performing VTST/MT calculations utilizing a cost effective Linux based PC cluster. Five different parallel computing schemes for choosing structures and computing their Hessians along a pre-defined Minimum Energy Path were developed. We found that the Energy Block and Asymmetric Cyclic Execution (EBACE) scheme, which is also most physically intuitive, results in converged rate constants with the least number of Hessians computed. We believe that carrying out the VTST/MT calculation in parallel makes it more attractive for calculating the rate constants of complex chemical systems.

We present a novel density-functional tight-binding (DFTB) parametrization toolkit developed to optimize the parameters of various DFTB models in a fully automatized fashion. The main features of the algorithm, based on the particle swarm optimization technique, are discussed, and a number of initial pilot applications of the developed methodology to molecular and solid systems are presented.

Explicit formulas for computing the Zhang–Zhang polynomial of cyclo-polyphenacenes, presented by Guo et al. (J. Math. Chem. 46:347, 2009) are found to be erroneous. In the present comment, a corrected version of Theorem 4.2 is given. The new formulation has been extensively tested by comparison with the Zhang–Zhang polynomials computed by brute force using a new, completely automatized computer code.

We present a detailed analysis of the factors influencing the formation of epoxide and ether groups in graphene nanoflakes using conventional density functional theory (DFT), the density-functional tight-binding (DFTB) method, π-Hückel theory, and graph theoretical invariants. The relative thermodynamic stability associated with the chemisorption of oxygen atoms at various positions on hexagonal graphene flakes (HGFs) of D(6h)-symmetry is determined by two factors - viz. the disruption of the π-conjugation of the HGF and the geometrical deformation of the HGF structure. The thermodynamically most stable structure is achieved when the former factor is minimized, and the latter factor is simultaneously maximized. Infrared (IR) spectra computed using DFT and DFTB reveal a close correlation between the relative thermodynamic stabilities of the oxidized HGF structures and their IR spectral activities. The most stable oxidized structures exhibit significant IR activity between 600 and 1800 cm(-1), whereas less stable oxidized structures exhibit little to no activity in this region. In contrast, Raman spectra are found to be less informative in this respect.