BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Kousuke Nakano

Kousuke Nakano contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phphysics.comp-phcond-mat.mtrl-scicond-mat.supr-concond-mat.str-elcond-mat.other

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2021arXiv

Diffusion Monte Carlo evaluation of disiloxane linearization barrier

The disiloxane molecule is a prime example of silicate compounds containing the Si-O-Si bridge. The molecule is of significant interest within the field of quantum chemistry, owing to the difficulty in theoretically predicting its properties. Herein, the linearisation barrier of disiloxane is investigated using a fixed-node diffusion Monte Carlo (FNDMC) approach, which is currently the most reliable {\it ab initio} method in accounting for an electronic correlation. Calculations utilizing the density functional theory (DFT) and the coupled cluster method with single and double substitutions, including noniterative triples (CCSD(T))are carried out alongside FNDMC for comparison. Two families of basis sets are used to investigate the disiloxane linearisation barrier - Dunning's correlation-consistent basis sets cc-pV$x$Z ($x = $ D, T, and Q) and their core-valence correlated counterparts, cc-pCV$x$Z. It is concluded that FNDMC successfully predicts the disiloxane linearisation barrier and does not depend on the completeness of the basis sets as much as DFT or CCSD(T), thus establishing its suitability.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Diffusion Monte Carlo Study on Relative Stabilities of Boron Nitride Polymorphs

Although Boron nitride (BN) is a well-known compound widely used for engineering and scientific purposes, the phase stability of its polymorphs, one of its most fundamental properties, is still under debate. The ab initio determination of the ground state of the BN polymorphs, such as hexagonal and zinc-blende, is difficult because of the elusive Van der Waals interaction, which plays a decisive role in some of the polymorphs, making quantitative prediction highly challenging. Hence, despite multiple theoretical studies, there has been no consensus on the ground state yet, primarily due to contradicting reports. In this study, we apply a state-of-the-art ab initio framework - fixed-node diffusion Monte Carlo (FNDMC), to four well known BN polymorphs, namely hexagonal, rhombohedral, wurtzite, and zinc-blende BNs. Our FNDMC calculations show that hBN is thermodynamically the most stable among the four polymorphs at 0 K as well as at 300K. This result agrees with the experimental data of Corrigan~{\it et al.} and Fukunaga. The conclusions are consistent with those obtained using other high-level methods, such as coupled cluster. We demonstrate that the FNDMC is a powerful method to address polymorphs that exhibit bonds of various forms. It also provides valuable information, like reliable reference energies, when reliable experimental data are missing or difficult to access. Our findings should promote the application of FNDMC for other van der Waals materials.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High-pressure Mg-Sc-H phase diagram and its superconductivity from first-principles calculations

In this work, global search for crystal structures of ternary Mg-Sc-H hydrides (Mg$_x$Sc$_y$H$_z$) under high pressure ($100 \le P \le 200$ GPa) were performed using the evolutionary algorithm and first-principles calculations. Based on them, we computed the thermodynamic convex hull and pressure-dependent phase diagram of Mg$_x$Sc$_y$H$_z$ for $z/(x+y) < 4$. We have identified the stable crystal structures of four thermodynamically stable compounds with the higher hydrogen content, i.e., $R\bar{3}m$-MgScH$_{6}$, $C2/m$-Mg$_{2}$ScH$_{10}$, $Immm$-MgSc$_{2}$H$_{9}$ and $Pm\bar{3}m$-Mg(ScH$_{4}$)$_{3}$. Their superconducting transition temperatures were computationally predicted by the McMillan-Allen-Dynes formula combined with first-principles phonon calculations. They were found to exhibit superconductivity; among them, $R\bar{3}m$-MgScH$_{6}$ was predicted to have the highest $T_{c}$ (i.e. 23.34 K) at 200 GPa.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Space-warp coordinate transformation for efficient ionic force calculations in quantum Monte Carlo

Ab-initio quantum Monte Carlo (QMC) methods are a state-of-the-art computational approach to obtaining highly accurate many-body wave functions. Although QMC methods are widely used in physics and chemistry to compute ground-state energies, calculation of atomic forces is still under technical/algorithmic development. Very recently, force evaluation has started to become of paramount importance for the generation of machine-learning force-field potentials. Nevertheless, there is no consensus regarding whether an efficient algorithm is available for the QMC force evaluation, namely one that scales well with the number of electrons and the atomic numbers. In this study, we benchmark the accuracy of all-electron variational Monte Carlo (VMC) and lattice-regularized diffusion Monte Carlo (LRDMC) forces for various mono- and heteronuclear dimers. The VMC and LRDMC forces were calculated with and without the so-called space-warp coordinate transformation (SWCT) and appropriate regularization techniques to remove the infinite variance problem. The LRDMC forces were computed with the Reynolds (RE) and variational-drift (VD) approximations. The potential energy surfaces obtained from the LRDMC energies give equilibrium bond lengths ($r_{\rm eq}$) and harmonic frequencies ($ω$) very close to the experimental values for all dimers, improving the corresponding VMC results. The LRDMC forces improve the VMC forces, implying that it is worth computing the DMC forces beyond VMC in spite of the higher computational cost. We find that the ratio of computational costs between QMC energy and forces scales as $Z^{\sim 2.5}$ without the SWCT. In contrast, the application of the SWCT makes the ratio {\it independent} of $Z$. As such, the accessible QMC system size is not affected by the evaluation of ionic forces but governed by the same scaling as the total energy one.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The systematic study on the stability and superconductivity of Y-Mg-H compounds under high pressure

Motivated by recent discovery of yttrium-based high-temperature ternary superconducting hydrides (e.g., CaYH$_{12}$, LaYH$_{12}$, and ScYH$_{6}$), we have employed evolutionary algorithm and first-principles calculations to comprehensively examine the structural stability and superconductivity of the YMgH$_{x}$ system at high pressure. The hydrogen content $x$ and the pressure are both important factors in the stability of these candidate structures. We find that the stability of hydrogen-rich materials frequently necessitates higher pressure. For instance, the pressures to stabilize $P4/mmm$-YMgH$_{8}$ and $Cmmm$-YMgH$_{12}$ are both more than 250 GPa. Hydrogen-less materials, such as $I4_{1}/amd$-YMgH$_{2}$ and $P6_{3}/mmc$-YMgH$_{3}$, can be stable at pressures as low as 100 GPa. In addition, we find a metastable structure for YMgH$_{6}$ with the same space group as the $P4/mmm$-YMgH$_{8}$. A metastable sodalite-like face-centered cubic (FCC) structure is also found in YMgH$_{12}$. These four clathrate structures of $P4/mmm$-YMgH$_{6}$, $P4/mmm$-YMgH$_{8}$, $Cmmm$-YMgH$_{12}$, and $Fd\bar{3}m$-YMgH$_{12}$ is made up of H14, H18, H24, and H24 cages, respectively, in which the H-H pair exhibits weak covalent bonding. According to phonon calculations, $P4/mmm$-YMgH$_{6}$ and $P4/mmm$-YMgH$_{8}$ require a pressure of 300 GPa to maintain dynamic stability, however $Cmmm$-YMgH$_{12}$ and $Fd\bar{3}m$-YMgH$_{12}$ can maintain dynamic stability at pressures of 200 GPa and 250 GPa, respectively. Electron-phonon coupling calculations indicate that they might be potential high-temperature superconductors, with superconductivity intimately linked to the H cage structure. The sodalite structure $Fd\bar{3}m$-YMgH$_{12}$ has a $T_\mathrm{c}$ value of 190 K and a strong electron-phonon coupling constant of 2.18.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

General correlated geminal ansatz for electronic structure calculations: exploiting Pfaffians in place of determinants

We propose here a single Pfaffian correlated variational ansatz, that dramatically improves the accuracy with respect to the single determinant one, while remaining at a similar computational cost. A much larger correlation energy is indeed determined by the most general two electron pairing function, including both singlet and triplet channels, combined with a many-body Jastrow factor, including all possible spin-spin spin-density and density-density terms. The main technical ingredient to exploit this accuracy is the use of the Pfaffian for antisymmetrizing an highly correlated pairing function, thus recovering the Fermi statistics for electrons with an affordable computational cost. Moreover the application of the Diffusion Monte Carlo, within the fixed node approximation, allows us to obtain very accurate binding energies for the first preliminary calculations reported in this study: C$_2$, N$_2$ and O$_2$ and the benzene molecule. This is promising and remarkable, considering that they represent extremely difficult molecules even for computationally demanding multi-determinant approaches, and opens therefore the way for realistic and accurate electronic simulations with an algorithm scaling at most as the fourth power of the number of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TurboRVB: a many-body toolkit for {\it ab initio} electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for {\it ab initio} Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC), and Diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions. The electronic wave function (WF) is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal product with spin-singlet pairing, or as a Pfaffian, including both singlet and triplet correlations. This wave function can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) ansatz, first proposed by L. Pauling and P. W. Anderson in quantum chemistry and condensed matter physics, respectively. The RVB ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization is made possible thanks to state-of-the-art stochastic algorithms for energy minimization. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, that is also an ideal environment for exploiting the computational power of modern GPU accelerators.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

All-electron quantum Monte Carlo with Jastrow single determinant Ansatz: application to the sodium dimer

In this work, we report potential energy surfaces (PESs) of the sodium dimer calculated by variational (VMC) and lattice regularized diffusion Monte Carlo (LRDMC). The VMC calculation is accurate for determining the equilibrium distance and the qualitative shape of the experimental PES. Remarkably, after the application of the LRDMC projection to this single determinant ansatz, namely the Jastrow Antisymmetrized Geminal Product (JAGP), chemical accuracy (~ 1kcal/mol) is reached, and the obtained dissociation energy, equilibrium internuclear distance, and harmonic vibrational frequency are in very good agreement with the experimental ones. This outcome crucially depends on the quality of the optimization used to determine the best possible trial function within the chosen ansatz. The strategy adopted in this work is to minimize the variational energy by initializing the trial function with the DFT single determinant ansatz expanded exactly in the same atomic basis used for the corresponding VMC and LRDMC calculations. This atomic basis is ad-hoc reshaped for QMC calculations. Indeed, we multiply the standard Gaussian type atomic orbitals by a one-body Jastrow factor, satisfying in this way the electron-ion cusp conditions. This allows us to use a very small basis almost converged in the complete basis set limit, by reducing the computational effort as well as the statistical fluctuations on the total energy. In order to achieve these important advantages, we have defined a very efficient DFT algorithm in the mentioned basis, by estimating the corresponding matrix elements on a mesh, and by using a much finer mesh grid in the vicinity of nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Speeding up the ab initio diffusion Monte Carlo by a smart lattice regularization

One of the most significant drawbacks of the all-electron ab initio diffusion Monte Carlo (DMC) is that its computational cost drastically increases with the atomic number ($Z$), which typically scales with $Z^{\sim 6}$. In this study, we introduce an algorithm based on a very efficient implementation of the Lattice Regularized Diffusion Monte Carlo (LRDMC), where the conventional time discretization is replaced by its lattice space counterpart. This scheme enables us to conveniently adopt a small lattice space in the vicinity of nuclei, and a large one in the valence region, by which a considerable speedup is achieved, especially for large atomic number $Z$. Indeed, the computational performances of our algorithm can be theoretically established by using the Thomas-Fermi model for heavy atoms, yielding an almost affordable scaling with the atomic number, i.e., $Z^{\sim 5}$. This opens the way for efficient and accurate all-electron ab initio DMC in electronic structure calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Exotic phase transition and superconductivity in layered titanium-oxypnictides implied by a computational phonon analysis

We applied $ab$ $initio$ phonon analysis to layered titanium-oxypnictides, Na$_2$Ti$_2Pn_2$O ($Pn$ = As, Sb), and found a clear contrast between the cases with lighter/heavier pnictogen in comparisons with experiments. The result completely explains the experimental structure, $C$2/$m$ for $Pn$ = As, within the conventional electron-phonon framework, while there arise discrepancies when the pnictogen gets heavier, being in the same trend for the BaTi$_2Pn_2$O ($Pn$ = As, Sb, Bi) case. The fact implies a systematic dependence on pnictogen to tune the mechanism of the phase transition and superconductivity from conventional to exotic.

0 citations0 reviewsNo code linkedOpen access