Source author record

Tsuyoshi Miyazaki

Tsuyoshi Miyazaki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Large-scale DFT methods for calculations of materials with complex structures

Large-scale density functional theory (DFT) calculations provide a powerful tool to investigate the atomic and electronic structure of materials with complex structures. This article reviews a large-scale DFT calculation method, the multi-site support function (MSSF) method, in the CONQUEST code. MSSFs are linear combinations of the basis functions which belong to a group of atoms in a local region. The method can reduce the computational time while preserving accuracy. The accuracy of MSSFs has been assessed for bulk Si, Al, Fe and NiO and hydrated DNA, which demonstrate the applicability of the MSSFs for varied materials. The applications of MSSFs on large systems with several thousand atoms, which have complex interfaces and non-periodic structures, indicate that the MSSF method is promising for precise investigations of materials with complex structures.

preprint2021arXiv

Unsupervised learning-based structural analysis: Search for a characteristic low-dimensional space by local structures in atomistic simulations

Owing to the advances in computational techniques and the increase in computational power, atomistic simulations of materials can simulate large systems with higher accuracy. Complex phenomena can be observed in such state-of-the-art atomistic simulations. However, it has become increasingly difficult to understand what is actually happening and mechanisms, for example, in molecular dynamics (MD) simulations. We propose an unsupervised machine learning method to analyze the local structure around a target atom. The proposed method, which uses the two-step locality preserving projections (TS-LPP), can find a low-dimensional space wherein the distributions of datapoints for each atom or groups of atoms can be properly captured. We demonstrate that the method is effective for analyzing the MD simulations of crystalline, liquid, and amorphous states and the melt-quench process from the perspective of local structures. The proposed method is demonstrated on a silicon single-component system, a silicon-germanium binary system, and a copper single-component system.

preprint2020arXiv

Blue moon ensemble simulation of aquation free energy profiles applied to mono and bifunctional platinum anticancer drugs

Aquation free energy profiles of neutral cisplatin and cationic monofunctional derivatives, including triaminochloroplatinum(II) and cis-diammine(pyridine)chloroplatinum(II), were computed using state of the art thermodynamic integration, for which temperature and solvent were accounted for explicitly using density functional theory based canonical molecular dynamics (DFT-MD). For all the systems the "inverse-hydration" where the metal center acts as an acceptor of hydrogen bond has been observed. This has motivated to consider the inversely bonded solvent molecule in the definition of the reaction coordinate required to initiate the constrained DFT-MD trajectories. We found that there exists little difference in free enthalpies of activations, such that these platinum-based anticancer drugs are likely to behave the same way in aqueous media. Detailed analysis of the microsolvation structure of the square-planar complexes, along with the key steps of the aquation mechanism are discussed.

preprint2020arXiv

First-principles study of the charge ordered phase in $κ$-D$_3$(Cat-EDT-TTF/ST)$_2$: Stability of $π$-electron deuterium coupled ordering in hydrogen-bonded molecular conductors

We study the electronic and structural properties of the low-temperature ordered phase of hydrogen-bonded molecular conductors, $κ$-D$_3$(Cat-EDT-TTF)$_2$ and its selenium-substituted analog $κ$-D$_3$(Cat-EDT-ST)$_2$, by means of first-principles density functional theory~(DFT) calculations. In these compounds, the charge ordering in the $π$-electron system is coupled with the ordering of the displacements in the deuteriums forming the hydrogen-bond, equally shared by two oxygens in the high-temperature phase. While the structural optimization within the standard DFT method based on the generalized gradient approximation fails to reproduce the structural stability of the charge-ordered (CO) phase, we show that a hybrid functional of Heyd, Scuseria, and Ernzerhof can reproduce structural characters of the CO phase, owing to the more localized nature of the wave functions. Furthermore, using the ability of the hybrid functional to predict the electronic and structural properties, we find a stable noncentrosymmetric CO phase with another pattern of deuterium ordering.

preprint2020arXiv

Large scale and linear scaling DFT with the CONQUEST code

We survey the underlying theory behind the large-scale and linear scaling DFT code, Conquest, which shows excellent parallel scaling and can be applied to thousands of atoms with exact solutions, and millions of atoms with linear scaling. We give details of the representation of the density matrix and the approach to finding the electronic ground state, and discuss the implementation of molecular dynamics with linear scaling. We give an overview of the performance of the code, focussing in particular on the parallel scaling, and provide examples of recent developments and applications.

preprint2015arXiv

A linear scaling DFT study of the growth of a new \{105\} facet layer on a Ge hut cluster

Ge hut clusters, which appear during heteroepitaxy of Ge on Si(001), are prototypical examples of islands formed through strain. Experimentally, complete facets are observed to form rapidly, though the mechanism is unknown. We model the growth of new faces on Ge hut clusters, using linear scaling DFT. We build realistic small huts on substrates, and show that the growth of \{105\} facets proceeds from top-to-bottom, even for these small huts. The growth of the facet is driven by the reconstruction on the \{105\} facet, and nucleates at the boundaries of the facet. We thus resolve any ambiguities in kinetic models of hut cluster growth.

preprint2015arXiv

First-principles study of hydrogen-bonded molecular conductor $κ$-H$_3$(Cat-EDT-TTF/ST)$_2$

We theoretically study hydrogen-bonded molecular conductors synthesized recently, $κ$-H$_3$(Cat-EDT-TTF)$_2$ and its diselena analog, $κ$-H$_3$(Cat-EDT-ST)$_2$, by first-principles density-functional theory calculations. In these crystals, two H(Cat-EDT-TTF/ST) units share a hydrogen atom with a short O--H--O hydrogen bond. The calculated band structure near the Fermi level shows a quasi-two-dimensional character, with a rather large interlayer dispersion due to the absence of insulating layers in contrast with conventional molecular conductors. We discuss effective low-energy models based on H(Cat-EDT-TTF/ST) units and its dimers, respectively, where the microscopic character of the orbitals composing them are analyzed. Furthermore, we find a stable structure which is different from the experimentally determined structure, where the shared hydrogen atom becomes localized to one of the oxygen atoms, in which charge disproportionation between the two types of H(Cat-EDT-TTF) units is associated. The calculated potential energy surface for the H atom is very shallow near the minimum points, therefore the probability of the H atom can be delocalized between the two O atoms.

preprint2015arXiv

Fragment Model Study of Molecular Multi-Orbital System $X$[Pd(dmit)$_2$]$_2$

Electronic properties of quasi-two-dimensional molecular conductors $X$[Pd(dmit)$_2$]$_2$ are studied theoretically. We construct an effective model based on the fragment molecular orbital scheme developed recently, which can describe the multi-orbital degree of freedom in this system. The tight-binding parameters for a series of $β'$-type compounds with different cations $X$ are evaluated by fitting to first-principles band calculations. We find that the transfer integrals within the dimers of Pd(dmit)$_2$ molecules, along the intramolecular and intermolecular bonds including the diagonal ones, are the same order, leading to hybridization between different molecular orbitals. This results in charge disproportionation within each molecule, as seen in our previous ab initio study [T. Tsumuraya et al, J. Phys. Soc. Jpn. 82, 033709 (2013)], and also to a revised picture of an effective dimer model. Furthermore, we discuss broken-symmetry insulating states triggered by interaction effects, which show characteristic features owing to the multi-orbital nature. The on-site Coulomb interaction induces antiferromagnetic states with intramolecular antiparallel spin pattern, while electron-lattice couplings stabilize non-magnetic charge-lattice ordered states where two kinds of dimers with different charge occupation arrange periodically. These states showing different spatial patterns compete with each other as well as with the paramagnetic metallic state.

preprint2015arXiv

Optimized multi-site local orbitals in the large-scale DFT program CONQUEST

We introduce numerical optimization of multi-site support functions in the linear-scaling DFT code CONQUEST. Multi-site support functions, which are linear combinations of pseudo-atomic orbitals on a target atom and those neighbours within a cutoff, have been recently proposed to reduce the number of support functions to the minimal basis while keeping the accuracy of a large basis [J. Chem. Theory Comput., 2014, 10, 4813]. The coefficients were determined by using the local filter diagonalization (LFD) method [Phys. Rev. B, 2009, 80, 205104]. We analyse the effect of numerical optimization of the coefficients produced by the LFD method. Tests on crystalline silicon, a benzene molecule and hydrated DNA systems show that the optimization improves the accuracy of the multi-site support functions with small cutoffs. It is also confirmed that the optimization guarantees the variational energy minimizations with multi-site support functions.

preprint2014arXiv

Efficient Calculations with Multisite Local Orbitals in the large-scale DFT Code CONQUEST

Multisite local orbitals, which are formed from linear combinations of pseudo-atomic orbitals from a target atom and its neighbor atoms, have been introduced in the large-scale density functional theory calculation code CONQUEST. Multisite local orbitals correspond to local molecular orbitals so that the number of required local orbitals can be minimal. The multisite support functions are determined by using the localized filter diagonalization (LFD) method [Phys. Rev. B 2009, 80, 205104]. Two new methods, the double cutoff method and the smoothing method, are introduced to the LFD method to improve efficiency and stability. The Hamiltonian and overlap matrices with multisite local orbitals are constructed by efficient sparse-matrix multiplications in CONQUEST. The investigation of the calculated energetic and geometrical properties and band structures of bulk Si, Al and DNA systems demonstrate the accuracy and the computational efficiency of the present method. The representability of both occupied and unoccupied band structures with the present method has been also confirmed.

preprint2014arXiv

Stable and Efficient Linear Scaling First-Principles Molecular Dynamics for 10,000+ atoms

The recent progress of linear-scaling or O(N) methods in the density functional theory (DFT) is remarkable. We expect that first-principles molecular dynamics (FPMD) simulations based on DFT can now treat more realistic and complex systems using the O(N) technique. However, very few examples of O(N) FPMD simulations exist so far and the information for the accuracy or reliability of the simulations is very limited. In this paper, we show that efficient and robust O(N) FPMD simulations are now possible by the combination of the extended Lagrangian Born-Oppenheimer molecular dynamics method, which was recently proposed by Niklasson et al (Phys. Rev. Lett. 100, 123004 (2008)), and the density matrix method as an O(N) technique. Using our linear-scaling DFT code Conquest, we investigate the reliable calculation conditions for the accurate O(N) FPMD and demonstrate that we are now able to do actual and reliable self-consistent FPMD simulation of a very large system containing 32,768 atoms.

preprint2013arXiv

Cation Dependence of the Electronic States in Molecular Triangular Lattice System β'-X[Pd(dmit)_2]_2: A First-principles study

The electronic structure of an isostructural series of molecular conductors, β'-X[Pd(dmit)_2]_2 is systematically studied by a first-principles method based on the density-functional theory. The calculated band structures are fitted to the tight-binding model based on Pd(dmit)_2 dimers on the triangular lattice. We find systematic variation in the anisotropy of the transfer integrals along the three directions of the triangular lattice taking different values. The transfer integral along the face-to-face stacking direction of Pd(dmit)_2 dimers is always the largest. Around the quantum spin liquid, X = EtMe_3Sb, the other two transfer integrals become comparable. We also report sensible differences in the distribution of wavefunctions near the Fermi level between the two dmit ligands of the Pd(dmit)_2 molecule.

preprint2013arXiv

Density-functional theory study of gramicidin A ion channel geometry and electronic properties

Understanding the mechanisms underlying ion channel function from the atomic-scale requires accurate ab initio modelling as well as careful experiments. Here, we present a density functional theory (DFT) study of the ion channel gramicidin A, whose inner pore conducts only monovalent cations and whose conductance has been shown to depend on the side chains of the amino acids in the channel. We investigate the ground-state geometry and electronic properties of the channel in vacuum, focusing on their dependence on the side chains of the amino acids. We find that the side chains affect the ground state geometry, while the electrostatic potential of the pore is independent of the side chains. This study is also in preparation for a full, linear scaling DFT study of gramicidin A in a lipid bilayer with surrounding water. We demonstrate that linear scaling DFT methods can accurately model the system with reasonable computational cost. Linear scaling DFT allows ab initio calculations with 10,000 to 100,000 atoms and beyond, and will be an important new tool for biomolecular simulations.

preprint2013arXiv

Effect of magnetic field on the superconducting phase in the electron-doped metallic double-chain compound Pr$_{2}$Ba$_{4}$Cu$_{7}$O$_{15-δ}$

We report the magnetotransport and $dc$ magnetic susceptibility of the polycrystalline samples of Pr$_{2}$Ba$_{4}$Cu$_{7}$O$_{15-δ}$, to examine the effect of magnetic field on the superconducting phase of the metallic CuO double chain. The resistive critical magnetic field is estimated to be about 21 T at low temperatures from the resistive transition data. On the other hand, the corresponding critical field determined from the magnetization measurements gives rise to a very low value of $\sim 0.3$T at 2 K. These discrepancies in the magnetic response between the resistivity and magnetization data are caused by disappearance of the magnetically shielding effect even in relatively lower fields. In spite of the observation of the resistive drop associated with the superconducting transport currents, the suppression of the diamagnetic signal is probably related to the superconductivity of quasi one-dimensional CuO double-chain. The behavior of Seebeck coefficient in the superconducting Pr$_{2}$Ba$_{4}$Cu$_{7}$O$_{15-δ}$ is discussed on the basis of the double chain model from the density functional band calculation.

preprint2012arXiv

Linear-scaling implementation of exact exchange using localized numerical orbitals and contraction reduction integrals

We present enhancements to the computational efficiency of exact exchange calculations using the density matrix and local support functions. We introduce a numerical method which avoids the explicit calculation the four-center two-electron repulsion integrals and reduces the prefactor scaling by a factor N, where N is the number of atoms within the range of the exact exchange Hamiltonian. This approach is based on a contraction-reduction scheme, and takes advantage of the discretization space which enables the direct summation over the support functions in a localized space. Using the sparsity property of the density matrix, the scaling of the prefactor can be further reduced to reach asymptotically O(N).