Seiji Tsuzuki

Bio: Seiji Tsuzuki is an academic researcher from National Institute of Advanced Industrial Science and Technology. The author has contributed to research in topics: Ionic liquid & Ab initio. The author has an hindex of 63, co-authored 300 publications receiving 16252 citations. Previous affiliations of Seiji Tsuzuki include University of Tokyo & Yokohama National University.

Physicochemical Properties and Structures of Room Temperature Ionic Liquids. 1. Variation of Anionic Species

Abstract: Room-temperature ionic liquids (RTILs) based on 1-butyl-3-methylimidazolium ([bmim]) with a variety of fluorinated anions were prepared, and the thermal behavior, density, viscosity, self-diffusion coefficients of the cations and anions, and ionic conductivity were measured over a wide temperature range. The temperature dependencies of the self-diffusion coefficient, viscosity, ionic conductivity, and molar conductivity have been fitted to the Vogel−Fulcher−Tamman equation, and the best-fit parameters for the self-diffusion coefficient, viscosity, ionic conductivity, and molar conductivity have been estimated, together with the linear fitting parameters for the density. The self-diffusion coefficients determined for the individual ions by pulsed-field-gradient spin−echo NMR method exhibit higher values for the cation compared with the anion over a wide temperature range, even if its radius is larger than that of the anionic radii. The summation of the cationic and anionic diffusion coefficients for the RT...

How Ionic Are Room-Temperature Ionic Liquids? An Indicator of the Physicochemical Properties

Abstract: Room-temperature ionic liquids (RTILs) are liquids consisting entirely of ions, and their important properties, e.g., negligible vapor pressure, are considered to result from the ionic nature. However, we do not know how ionic the RTILs are. The ionic nature of the RTILs is defined in this study as the molar conductivity ratio (Λimp/ΛNMR), calculated from the molar conductivity measured by the electrochemical impedance method (Λimp) and that estimated by use of pulse-field-gradient spin−echo NMR ionic self-diffusion coefficients and the Nernst−Einstein relation (ΛNMR). This ratio is compared with solvatochromic polarity scales: anionic donor ability (Lewis basicity), ET(30), hydrogen bond donor acidity (α), and dipolarity/polarizability (π*), as well as NMR chemical shifts. The Λimp/ΛNMR well illustrates the degree of cation−anion aggregation in the RTILs at equilibrium, which can be explained by the effects of anionic donor and cationic acceptor abilities for the RTILs having different anionic and catio...

Origin of Attraction and Directionality of the π/π Interaction: Model Chemistry Calculations of Benzene Dimer Interaction

Abstract: A model chemistry for the evaluation of intermolecular interaction between aromatic molecules (AIMI Model) has been developed. The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium size basis set. The calculated interaction energies of the parallel, T-shaped,and slipped-parallel benzene dimers are -1.48, -2.46, and -2.48 kcal/mol, respectively. The substantial attractive interaction in benzene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge-transfer but long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases attraction significantly. The dispersion interaction is found to be the major source of attraction in the benzene dimer. The orientation dependence of the dimer interaction is mainly controlled by long-range interactions. Although electrostatic interaction is considerably weaker than dispersion interaction, it is highly orientation dependent. Dispersion and electrostatic interactions are both important for the directionality of the benzene dimer interaction.

Interaction energies of van der Waals and hydrogen bonded systems calculated using density functional theory: Assessing the PW91 model

Abstract: The performance of density functional theory using the Perdew and Wang’s exchange and correlation functionals (PW91) functional for the prediction of intermolecular interactionenergies is evaluated based on calculations on the neon, argon, methane, ethylene, and benzene dimers, as well as on 12 hydrogen bonded complexes (water, methanol, formic acid, hydrogen fluoride, ammonia, formamide dimers and water–methanol, water–dimethyl ether, water–formaldehyde, hydrogen cyanide–hydrogen fluoride, water–ammonia, water–formamide complexes). The results were compared with those obtained from Becke’s exchange and Lee, Yang, and Parr’s correlation functionals (BLYP), Becke’s 3 parameter functional combined with Lee, Yang, and Parr’s correlation functional (B3LYP), second order Mo/ller–Plesset perturbation (MP2), and coupled cluster calculations with single and double substitutions and with non-iterative triple corrections [CCSD(T)] calculations. The calculated interactionenergies show that the PW91 functional performs much better than the BLYP or B3LYP functionals. The error in the computed binding energies of the hydrogen bonded complexes is 20% in the worst case. The most demanding cases are the systems with large dispersion contributions to the binding energy, such as the benzene dimer. In contrast to the BLYP and B3LYP functionals which fail to account for dispersion, the PW91 functional at least partly recovers the attraction. The basis set dependence of the PW91 functionals is relatively small in contrast to the MP2 and CCSD(T) methods. Despite its occasional difficulties with dispersion interaction, the PW91 functional may be a viable alternative to the ab initio methods, certainly in situations where large complexes are being studied.

Oxidative-stability enhancement and charge transport mechanism in glyme-lithium salt equimolar complexes.

Abstract: The oxidative stability of glyme molecules is enhanced by the complex formation with alkali metal cations. Clear liquid can be obtained by simply mixing glyme (triglyme or tetraglyme) with lithium bis(trifluoromethylsulfonyl)amide (Li[TFSA]) in a molar ratio of 1:1. The equimolar complex [Li(triglyme or tetraglyme)(1)][TFSA] maintains a stable liquid state over a wide temperature range and can be regarded as a room-temperature ionic liquid consisting of a [Li(glyme)(1)](+) complex cation and a [TFSA](-) anion, exhibiting high self-dissociativity (ionicity) at room temperature. The electrochemical oxidation of [Li(glyme)(1)][TFSA] takes place at the electrode potential of ~5 V vs Li/Li(+), while the oxidation of solutions containing excess glyme molecules ([Li(glyme)(x)][TFSA], x > 1) occurs at around 4 V vs Li/Li(+). This enhancement of oxidative stability is due to the donation of lone pairs of ether oxygen atoms to the Li(+) cation, resulting in the highest occupied molecular orbital (HOMO) energy level lowering of a glyme molecule, which is confirmed by ab initio molecular orbital calculations. The solvation state of a Li(+) cation and ion conduction mechanism in the [Li(glyme)(x)][TFSA] solutions is elucidated by means of nuclear magnetic resonance (NMR) and electrochemical methods. The experimental results strongly suggest that Li(+) cation conduction in the equimolar complex takes place by the migration of [Li(glyme)(1)](+) cations, whereas the ligand exchange mechanism is overlapped when interfacial electrochemical reactions of [Li(glyme)(1)](+) cations occur. The ligand exchange conduction mode is typically seen in a lithium battery with a configuration of [Li anode|[Li(glyme)(1)][TFSA]|LiCoO(2) cathode] when the discharge reaction of a LiCoO(2) cathode, that is, desolvation of [Li(glyme)(1)](+) and insertion of the resultant Li(+) into the cathode, occurs at the electrode-electrolyte interface. The battery can be operated for more than 200 charge-discharge cycles in the cell voltage range of 3.0-4.2 V, regardless of the use of ether-based electrolyte, because the ligand exchange rate is much faster than the electrode reaction rate.

A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu

Abstract: The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.

Fast parallel algorithms for short-range molecular dynamics

Abstract: Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Semiempirical GGA-type density functional constructed with a long-range dispersion correction.

Abstract: A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

Density functionals with broad applicability in chemistry.

Abstract: Although density functional theory is widely used in the computational chemistry community, the most popular density functional, B3LYP, has some serious shortcomings: (i) it is better for main-group chemistry than for transition metals; (ii) it systematically underestimates reaction barrier heights; (iii) it is inaccurate for interactions dominated by medium-range correlation energy, such as van der Waals attraction, aromatic−aromatic stacking, and alkane isomerization energies. We have developed a variety of databases for testing and designing new density functionals. We used these data to design new density functionals, called M06-class (and, earlier, M05-class) functionals, for which we enforced some fundamental exact constraints such as the uniform-electron-gas limit and the absence of self-correlation energy. Our M06-class functionals depend on spin-up and spin-down electron densities (i.e., spin densities), spin density gradients, spin kinetic energy densities, and, for nonlocal (also called hybrid)...

The hydrogen bond in the solid state.

Abstract: The hydrogen bond is the most important of all directional intermolecular interactions. It is operative in determining molecular conformation, molecular aggregation, and the function of a vast number of chemical systems ranging from inorganic to biological. Research into hydrogen bonds experienced a stagnant period in the 1980s, but re-opened around 1990, and has been in rapid development since then. In terms of modern concepts, the hydrogen bond is understood as a very broad phenomenon, and it is accepted that there are open borders to other effects. There are dozens of different types of X-H.A hydrogen bonds that occur commonly in the condensed phases, and in addition there are innumerable less common ones. Dissociation energies span more than two orders of magnitude (about 0.2-40 kcal mol(-1)). Within this range, the nature of the interaction is not constant, but its electrostatic, covalent, and dispersion contributions vary in their relative weights. The hydrogen bond has broad transition regions that merge continuously with the covalent bond, the van der Waals interaction, the ionic interaction, and also the cation-pi interaction. All hydrogen bonds can be considered as incipient proton transfer reactions, and for strong hydrogen bonds, this reaction can be in a very advanced state. In this review, a coherent survey is given on all these matters.