Researcher profile

Takeo Hoshi

Takeo Hoshi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
7topics
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

11 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.

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.

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.

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.

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.

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.

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.

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.

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/.

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.

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.