BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Takeo Hoshi

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

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.comp-phmath.NANumerical AnalysisDistributed, Parallel, and Cluster Computingcond-mat.str-elphysics.plasm-ph

Catalog footprint

What is connected

21works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

21 published item(s)

preprint2025arXiv

Monte Carlo simulation method for incoherent Thomson scattering spectra from arbitrary electron distribution functions

We developed a Monte Carlo simulation method to calculate incoherent Thomson scattering spectra in high temperature plasmas. The basic idea is to treat the entire scattering process as the superposition of individual photon-electron interactions. We introduce macro-particles, referred from particle-in-cell simulations, to reduce the computational cost, and obtain scattered spectra within a reasonable computational time. Since the velocity of the interacting electron is randomly sampled from an electron distribution function, the method can be applied to arbitrary electron distribution functions provided an appropriate sampling scheme is available. We present simulation results for relativistic Maxwellian and kappa distribution functions, and compare them with both analytical and numerical spectra for validation. The simulated spectra show good agreement with both analytical and numerical results, demonstrating that the Monte Carlo simulation method can reliably reproduce incoherent Thomson scattering spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Data-analysis software framework 2DMAT and its application to experimental measurements for two-dimensional material structures

An open-source data-analysis framework 2DMAT has been developed for experimental measurements of two-dimensional material structures. 2DMAT offers five analysis methods: (i) Nelder-Mead optimization, (ii) grid search, (iii) Bayesian optimization, (iv) replica exchange Monte Carlo method, and (v) population-annealing Monte Carlo method. Methods (ii) through (v) are implemented by parallel computation,which is efficient not only for personal computers but also for supercomputers.The current version of 2DMAT is applicable to total-reflection high-energy positron diffraction (TRHEPD), surface X-ray diffraction (SXRD), and low-energy electron diffraction (LEED) experiments by installing corresponding forward problem solvers that generate diffraction intensity data from a given dataset of the atomic positions. The analysis methods are general and can be applied also to other experiments and problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

sim-trhepd-rheed -- Open-source simulator of total-reflection high-energy positron diffraction (TRHEPD) and reflection high-energy electron diffraction (RHEED)

The present paper reports sim-trhepd-rheed (STR), an open-source simulator of total-reflection high-energy positron diffraction (TRHEPD) and reflection high-energy electron diffraction (RHEED) experiments which are used for atom-scale surface structure determination of a material. The STR simulator is used for the analysis of experimental diffraction data by simulating the rocking curve from a given trial surface structure by solving the partial differential equation of the dynamical quantum diffraction theory for positron or electron wavefunctions. Using the obtained surface structure, electronic structure, and other physical quantities can be evaluated through first-principles calculations. For this purpose, a utility software was also developed in order to realize a first principles calculation with the Quantum ESPRESSO suite.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Data-driven sensitivity analysis in a total-reflection high-energy positron diffraction (TRHEPD) experiment of the Si$_4$O$_5$N$_3$ / 6H-SiC (0001)-($\sqrt{3} \times \sqrt{3}$) R30$^\circ$

The present article proposes a data analysis method for experimentally-derived measurements, which consists of an auto-optimization procedure and a sensitivity analysis. The method was applied to the results of a total-reflection high-energy positron diffraction (TRHEPD) experiment, a novel technique of determining surface structures or the position of the atoms near the material surface. This method solves numerically the partial differential equation in the fully-dynamical quantum diffraction theory with many trial surface structures. In the sensitivity analysis, we focused on the experimental uncertainties and the variation over individual fitting parameters, which was analyzed by solving the eigenvalue problem of the variance-covariance matrix. A modern massively parallel supercomputer was used to complete the analysis within a moderate computational time. The sensitivity analysis provides a basis for the choice of variables in the data analysis for practical reliability. The effectiveness of the present analysis method was demonstrated in the structure determination of a Si$_4$O$_5$N$_3$ / 6H-SiC(0001)-($\sqrt{3} \times \sqrt{3}$) R30$^\circ$ surface. Furthermore, this analysis method is applicable to many experiments other than TRHEPD.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Performance prediction of massively parallel computation by Bayesian inference

A performance prediction method for massively parallel computation is proposed. The method is based on performance modeling and Bayesian inference to predict elapsed time T as a function of the number of used nodes P (T=T(P)). The focus is on extrapolation for larger values of P from the perspective of application researchers. The proposed method has several improvements over the method developed in a previous paper, and application to real-symmetric generalized eigenvalue problem shows promising prediction results. The method is generalizable and applicable to many other computations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Parallelizable Energy-Preserving Integrator MB4 and Its Application to Quantum-Mechanical Wavepacket Dynamics

In simulating physical systems, conservation of the total energy is often essential, especially when energy conversion between different forms of energy occurs frequently. Recently, a new fourth order energy-preserving integrator named MB4 was proposed based on the so-called continuous stage Runge--Kutta methods (Y.~Miyatake and J.~C.~Butcher, SIAM J.~Numer.~Anal., 54(3), 1993-2013). A salient feature of this method is that it is parallelizable, which makes its computational time for one time step comparable to that of second order methods. In this paper, we illustrate how to apply the MB4 method to a concrete ordinary differential equation using the nonlinear Schrödinger-type equation on a two-dimensional grid as an example. This system is a prototypical model of two-dimensional disordered organic material and is difficult to solve with standard methods like the classical Runge--Kutta methods due to the nonlinearity and the $δ$-function like potential coming from defects. Numerical tests show that the method can solve the equation stably and preserves the total energy to 16-digit accuracy throughout the simulation. It is also shown that parallelization of the method yields up to 2.8 times speedup using 3 computational nodes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An a posteriori verification method for generalized real-symmetric eigenvalue problems in large-scale electronic state calculations

An a posteriori verification method is proposed for the generalized real-symmetric eigenvalue problem and is applied to densely clustered eigenvalue problems in large-scale electronic state calculations. The proposed method is realized by a two-stage process in which the approximate solution is computed by existing numerical libraries and is then verified in a moderate computational time. The procedure returns intervals containing one exact eigenvalue in each interval. Test calculations were carried out for organic device materials, and the verification method confirms that all exact eigenvalues are well separated in the obtained intervals. This verification method will be integrated into EigenKernel (https://github.com/eigenkernel/), which is middleware for various parallel solvers for the generalized eigenvalue problem. Such an a posteriori verification method will be important in future computational science.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

K$ω$ -- Open-source library for the shifted Krylov subspace method of the form $(zI-H)x=b$

We develop K$ω$, an open-source linear algebra library for the shifted Krylov subspace methods. The methods solve a set of shifted linear equations $(z_k I-H)x^{(k)}=b\, (k=0,1,2,...)$ for a given matrix $H$ and a vector $b$, simultaneously. The leading order of the operational cost is the same as that for a single equation. The shift invariance of the Krylov subspace is the mathematical foundation of the shifted Krylov subspace methods. Applications in materials science are presented to demonstrate the advantages of the algorithm over the standard Krylov subspace methods such as the Lanczos method. We introduce benchmark calculations of (i) an excited (optical) spectrum and (ii) intermediate eigenvalues by the contour integral on the complex plane. In combination with the quantum lattice solver $\mathcal{H} Φ$, K$ω$ can realize parallel computation of excitation spectra and intermediate eigenvalues for various quantum lattice models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

White Paper from Workshop on Large-scale Parallel Numerical Computing Technology (LSPANC 2020): HPC and Computer Arithmetic toward Minimal-Precision Computing

In numerical computations, precision of floating-point computations is a key factor to determine the performance (speed and energy-efficiency) as well as the reliability (accuracy and reproducibility). However, precision generally plays a contrary role for both. Therefore, the ultimate concept for maximizing both at the same time is the minimal-precision computing through precision-tuning, which adjusts the optimal precision for each operation and data. Several studies have been already conducted for it so far (e.g. Precimoniuos and Verrou), but the scope of those studies is limited to the precision-tuning alone. Hence, we aim to propose a broader concept of the minimal-precision computing system with precision-tuning, involving both hardware and software stack. In 2019, we have started the Minimal-Precision Computing project to propose a more broad concept of the minimal-precision computing system with precision-tuning, involving both hardware and software stack. Specifically, our system combines (1) a precision-tuning method based on Discrete Stochastic Arithmetic (DSA), (2) arbitrary-precision arithmetic libraries, (3) fast and accurate numerical libraries, and (4) Field-Programmable Gate Array (FPGA) with High-Level Synthesis (HLS). In this white paper, we aim to provide an overview of various technologies related to minimal- and mixed-precision, to outline the future direction of the project, as well as to discuss current challenges together with our project members and guest speakers at the LSPANC 2020 workshop; https://www.r-ccs.riken.jp/labs/lpnctrt/lspanc2020jan/.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Development of data-analysis software for total-reflection high-energy positron diffraction (TRHEPD)

The present paper reports on the recent activity of the data analysis software development for total-reflection high-energy positron diffraction (TRHEPD), a novel experimental technique for surface structure determination. Experiments using TRHEPD are being conducted intensively at the Slow Positron Facility, Institute of Materials Structure Science, High Energy Accelerator Research Organization, revealing surface structure of interest. The data analysis software provides a solution to the inverse problem in which the atomic positions of a surface structure are determined from the experimental diffraction data (rocking curve). The forward problem is solved by the numerical solution of the partial differential equation in the quantum scattering problem. A technical demonstration with a test problem was carried out to confirm the software functioned as expected. Since the analysis method has a general mathematical foundation, it is also applicable to other experiments, such as X-ray or electron diffraction experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Numerical aspect of large-scale electronic state calculation for flexible device material

Numerical aspects of large-scale electronic state calculation are explored on flexible organic device materials. Physical theory, numerical method and real application studies are discussed in the context of application-algorithm-architecture co-design. An application study was carried out for disordered organic thin film. Participation ratio, a measure for the spatial extension of electronic wavefunction is focused on, since it is crucial for device property. A data scientific research is reported for a classification problem of disordered organic polymers, in which participation ratio is used as descriptor. These application studies indicate the potential need of purpose-specific solvers for internal eigenpairs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extremely scalable algorithm for 10$^8$-atom quantum material simulation on the full system of the K computer

An extremely scalable linear-algebraic algorithm was developed for quantum material simulation (electronic state calculation) with 10$^8$ atoms or 100-nm-scale materials. The mathematical foundation is generalized shifted linear equations ((zB - A) x = b), instead of conventional generalized eigenvalue equations. The method has a highly parallelizable mathematical structure. The fundamental theory is mathematical and is applicable also to other scientific fields. The benchmark shows an extreme strong scaling and a qualified time-to-solution on the full system of the K computer. The method was demonstrated in a real material research for ultra-flexible (organic) devices, key devices of next-generation IoT products. The present paper shows that an innovative scalable algorithm for a real research can appear by the co-design among application, algorithm and architecture.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

One-hundred-nm-scale electronic structure and transport calculations of organic polymers on the K computer

One-hundred-nm-scale electronic structure calculations were carried out on the K supercomputer by our original simulation code ELSES (http://www.elses.jp/) The present paper reports preliminary results of transport calculations for condensed organic polymers. Large-scale calculations are realized by novel massively parallel order-N algorithms. The transport calculations were carried out as a theoretical extension for the quantum wavepacket dynamics simulation. The method was applied to a single polymer chain and condensed polymers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient numerical solver for first-principles transport calculation based on real-space finite-difference method

We propose an efficient procedure to obtain Green's functions by combining the shifted conjugate orthogonal conjugate gradient (shifted COCG) method with the nonequilibrium Green's function (NEGF) method based on a real-space finite-difference (RSFD) approach. The bottleneck of the computation in the NEGF scheme is matrix inversion of the Hamiltonian including the self-energy terms of electrodes to obtain perturbed Green's function in the transition region. This procedure first computes unperturbed Green's functions and calculates perturbed Green's functions from the unperturbed ones using a mathematically strict relation. Since the matrices to be inverted to obtain the unperturbed Green's functions are sparse, complex-symmetric and shifted for a given set of sampling energy points, we can use the shifted COCG method, in which once the Green's function for a reference energy point has been calculated, the Green's functions for the other energy points can be obtained with a moderate computational cost. We calculate the transport properties of a C$_{60}$@(10,10) carbon nanotube (CNT) peapod suspended by (10,10)CNTs as an example of a large-scale transport calculation. The proposed scheme opens the possibility of performing large-scale RSFD-NEGF transport calculations using massively parallel computers without the loss of accuracy originating from the incompleteness of the localized basis set.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hybrid Numerical Solvers for Massively Parallel Eigenvalue Computation and Their Benchmark with Electronic Structure Calculations

Optimally hybrid numerical solvers were constructed for massively parallel generalized eigenvalue problem (GEP).The strong scaling benchmark was carried out on the K computer and other supercomputers for electronic structure calculation problems in the matrix sizes of M = 10^4-10^6 with upto 105 cores. The procedure of GEP is decomposed into the two subprocedures of the reducer to the standard eigenvalue problem (SEP) and the solver of SEP. A hybrid solver is constructed, when a routine is chosen for each subprocedure from the three parallel solver libraries of ScaLAPACK, ELPA and EigenExa. The hybrid solvers with the two newer libraries, ELPA and EigenExa, give better benchmark results than the conventional ScaLAPACK library. The detailed analysis on the results implies that the reducer can be a bottleneck in next-generation (exa-scale) supercomputers, which indicates the guidance for future research. The code was developed as a middleware and a mini-application and will appear online.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Novel linear algebraic theory and one-hundred-million-atom quantum material simulations on the K computer

The present paper gives a review of our recent progress and latest results for novel linear-algebraic algorithms and its application to large-scale quantum material simulations or electronic structure calculations. The algorithms are Krylov-subspace (iterative) solvers for generalized shifted linear equations, in the form of (zS-H)x=b,in stead of conventional generalized eigen-value equation. The method was implemented in our order-$N$ calculation code ELSES (http://www.elses.jp/) with modelled systems based on ab initio calculations. The code realized one-hundred-million-atom, or 100-nm-scale, quantum material simulations on the K computer in a high parallel efficiency with up to all the built-in processor cores. The present paper also explains several methodological aspects, such as use of XML files and 'novice' mode for general users. A sparse matrix data library in our real problems (http://www.elses.jp/matrix/) was prepared. Internal eigen-value problem is discussed as a general need from the quantum material simulation. The present study is a interdisciplinary one and is sometimes called 'Application-Algorithm-Architecture co-design'. The co-design will play a crucial role in exa-scale scientific computations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Novel linear algebraic theory and one-hundred-million-atom electronic structure calculation on the K computer

A novel linear-algebraic algorithm, multiple Arnoldi method, was developed in an interdisciplinary study between physics and applied mathematics and realized one-hundred-million-atom (100-nm-scale) electronic state calculations on the K computer. The algorithms are Krylov-subspace solvers for generalized shifted linear equations and were implemented in our order-N calculation code ELSES (http://www.elses.jp/). Moreover, a method for calculating eigen states is presented as a theoretical extension.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ten-million-atom electronic structure calculations on the K computer with a massively parallel order-N theory

A massively parallel order-N electronic structure theory was constructed by an interdisciplinary research between physics, applied mathematics and computer science. (1) A high parallel efficiency with ten-million-atom nanomaterials was realized on the K computer with upto 98,304 processor cores. The mathematical foundation is a novel linear algebraic algorithm for the generalized shifted linear equation. The calculation was carried out by our code ' ELSES ' (www.elses.jp) with modelled (tight-binding-form) systems based on ab initio calculations. (2) A post-calculation analysis method, called pi-orbital crystalline orbital Hamiltonian population (pi-COHP) method, is presented, since the method is ideal for huge electronic structure data distributed among massive nodes. The analysis method is demonstrated in an sp2-sp3 nano-composite carbon solid, with an original visualization software 'VisBAR'. The present research indicates general aspects of computational physics with current or next-generation supercomputers.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A hierarchical research by large-scale and ab initio electronic structure theories -- Si and Ge cleavage and stepped (111)-2x1 surfaces --

The ab initio calculation with the density functional theory and plane-wave bases is carried out for stepped Si(111)-2x1 surfaces that were predicted in a cleavage simulation by the large-scale (order-N) electronic structure theory (T. Hoshi, Y. Iguchi and T. Fujiwara, Phys. Rev. B72 (2005) 075323). The present ab initio calculation confirms the predicted stepped structure and its bias-dependent STM image. Moreover, two (meta)stable step-edge structures are found and compared. The investigation is carried out also for Ge(111)-2x1 surfaces, so as to construct a common understanding among elements. The present study demonstrates the general importance of the hierarchical research between large-scale and ab initio electronic structure theories.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Large scale simulation of quantum-mechanical molecular dynamics for nano-polycrystalline diamond

Quantum-mechanical molecular-dynamics simulations are carried out to explore possible precursor states of nano-polycrystalline diamond, a novel ultra-hard material produced directly from graphite. Large-scale simulation with 10^5 atoms is realized by using the ' order-N' simulation code 'ELSES' (http://www.elses.jp). The simulation starts with a diamond structure that contains initial structural defects and results in a mixture of graphite(sp^2)-like and diamond(sp^3)-like regions as nano-meter-scale domains. We speculate that the domains are metastable and are possible candidates of the precursor structures.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Accuracy control in ultra-large-scale electronic structure calculation

Numerical aspects are investigated in ultra-large-scale electronic structure calculation. Accuracy control methods in process (molecular-dynamics) calculation are focused. Flexible control methods are proposed so as to control variational freedoms, automatically at each time step, within the framework of generalized Wannier state theory. The method is demonstrated in silicon cleavage simulation with 10^2-10^5 atoms. The idea is of general importance among process calculations and is also used in Krylov subspace theory, another large-scale-calculation theory.

0 citations0 reviewsNo code linkedOpen access