Source author record

Kazuhide Ichikawa

Kazuhide Ichikawa 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

18works
11topics
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

18 published item(s)

preprint2021arXiv

Deep variational quantum eigensolver for excited states and its application to quantum chemistry calculation of periodic materials

A programmable quantum device that has a large number of qubits without fault-tolerance has emerged recently. Variational Quantum Eigensolver (VQE) is one of the most promising ways to utilize the computational power of such devices to solve problems in condensed matter physics and quantum chemistry. As the size of the current quantum devices is still not large for rivaling classical computers at solving practical problems, Fujii et al. proposed a method called "Deep VQE" which can provide the ground state of a given quantum system with the smaller number of qubits by combining the VQE and the technique of coarse-graining [K. Fujii, et al, arXiv:2007.10917]. In this paper, we extend the original proposal of Deep VQE to obtain the excited states and apply it to quantum chemistry calculation of a periodic material, which is one of the most impactful applications of the VQE. We first propose a modified scheme to construct quantum states for coarse-graining in Deep VQE to obtain the excited states. We also present a method to avoid a problem of meaningless eigenvalues in the original Deep VQE without restricting variational quantum states. Finally, we classically simulate our modified Deep VQE for quantum chemistry calculation of a periodic hydrogen chain as a typical periodic material. Our method reproduces the ground-state energy and the first-excited-state energy with the errors up to O(1)% despite the decrease in the number of qubits required for the calculation by two or four compared with the naive VQE. Our result will serve as a beacon for tackling quantum chemistry problems with classically-intractable sizes by smaller quantum devices in the near future.

preprint2016arXiv

Computational Method for the Retarded Potential in the Real-Time Simulation of Quantum Electrodynamics

We discuss the method to compute the integrals which appear in the retarded potential term for a real-time simulation based on QED (Quantum Electrodynamics). We show that the oscillatory integrals over the infinite interval involved in them can be efficiently performed by the method developed by Ooura and Mori based on the double exponential (DE) formula.

preprint2016arXiv

Theoretical study of atoms by the electronic kinetic energy density and stress tensor density

We analyze the electronic structure of atoms in the first, second and third periods using the electronic kinetic energy density and stress tensor density, which are local quantities motivated by quantum field theoretic consideration, specifically the rigged quantum electrodynamics. We compute the zero surfaces of the electronic kinetic energy density, which we call the electronic interfaces, of the atoms. We find that their sizes exhibit clear periodicity and are comparable to the conventional atomic and ionic radii. We also compute the electronic stress tensor density and its divergence, tension density, of the atoms, and discuss how their electronic structures are characterized by them.

preprint2016arXiv

Theoretical Study of Lithium Ionic Conductors by Electronic Stress Tensor Density and Electronic Kinetic Energy Density

We analyze the electronic structure of lithium ionic conductors, ${\rm Li_3PO_4}$ and ${\rm Li_3PS_4}$, using the electronic stress tensor density and kinetic energy density with special focus on the ionic bonds among them. We find that, as long as we examine the pattern of the eigenvalues of the electronic stress tensor density, we cannot distinguish between the ionic bonds and bonds among metalloid atoms. We then show that they can be distinguished by looking at the morphology of the electronic interface, the zero surface of the electronic kinetic energy density.

preprint2015arXiv

Electronic stress tensor analysis of molecules in gas phase of CVD process for GeSbTe alloy

We analyze the electronic structure of molecules which may exist in gas phase of chemical vapor deposition process for GeSbTe alloy using the electronic stress tensor, with special focus on the chemical bonds between Ge, Sb and Te atoms. We find that, from the viewpoint of the electronic stress tensor, they have intermediate properties between alkali metals and hydrocarbon molecules. We also study the correlation between the bond order which is defined based on the electronic stress tensor, and energy-related quantities. We find that the correlation with the bond dissociation energy is not so strong while one with the force constant is very strong. We interpret these results in terms of the energy density on the "Lagrange surface", which is considered to define the boundary surface of atoms in a molecule in the framework of the electronic stress tensor analysis.

preprint2014arXiv

Study of Simulation Method of Time Evolution of Atomic and Molecular Systems by Quantum Electrodynamics

We discuss a method to follow step-by-step time evolution of atomic and molecular systems based on QED (Quantum Electrodynamics). Our strategy includes expanding the electron field operator by localized wavepackets to define creation and annihilation operators and following the time evolution using the equations of motion of the field operator in the Heisenberg picture. We first derive a time evolution equation for the excitation operator, the product of two creation or annihilation operators, which is necessary for constructing operators of physical quantities such as the electronic charge density operator. We then describe our approximation methods to obtain time differential equations of the electronic density matrix, which is defined as the expectation value of the excitation operator. By solving the equations numerically, we show "electron-positron oscillations", the fluctuations originated from virtual electron-positron pair creations and annihilations, appear in the charge density of a hydrogen atom and molecule. We also show that the period of the electron-positron oscillations becomes shorter by including the self-energy process, in which the electron emits a photon and then absorbs it again, and it can be interpreted as the increase in the electron mass due to the self-energy.

preprint2013arXiv

Theoretical study of lithium clusters by electronic stress tensor

We study the electronic structure of small lithium clusters Li_n (n=2-8) using the electronic stress tensor. We find that the three eigenvalues of the electronic stress tensor of the Li clusters are negative and degenerate, just like the stress tensor of liquid. This leads us to propose that we may characterize a metallic bond in terms of the electronic stress tensor. Our proposal is that in addition to the negativity of the three eigenvalues of the electronic stress tensor, their degeneracy characterizes some aspects of the metallic nature of chemical bonding. To quantify the degree of degeneracy, we use the differential eigenvalues of the electronic stress tensor. By comparing the Li clusters and hydrocarbon molecules, we show that the sign of the largest eigenvalue and the differential eigenvalues could be useful indices to evaluate the metallicity or covalency of a chemical bond.

preprint2012arXiv

Study of Simulation Method of Time Evolution in Rigged QED

We discuss how we formulate time evolution of physical quantities in the framework of the Rigged QED (Quantum Electrodynamics). The Rigged QED is a theory which has been proposed to treat dynamics of electrons, photons and atomic nuclei in atomic and molecular systems in a quantum field theoretic way. To solve the dynamics in the Rigged QED, we need different techniques from those developed for the conventional QED. As a first step toward this issue, we propose a procedure to expand the Dirac field operator, which represents electrons, by the electron annihilation/creation operators and solutions of the Dirac equation for electrons in nuclear potential. Similarly, the Schrodinger field operators, which represent atomic nuclei, are expanded by nucleus annihilation/creation operators. Then we derive time evolution equations for these annihilation and creation operators and discuss how time evolution of the operators for physical quantities can be calculated. In the end, we propose a method to approximate the evolution equations of the operators by the evolution equations for the density matrices of electrons and atomic nuclei. Under this approximation, we carry out numerical simulation of the time evolution of electron charge density of a hydrogen atom.

preprint2011arXiv

Electronic stress tensor analysis of hydrogenated palladium clusters

We study the chemical bonds of small palladium clusters Pd_n (n=2-9) saturated by hydrogen atoms using electronic stress tensor. Our calculation includes bond orders which are recently proposed based on the stress tensor. It is shown that our bond orders can classify the different types of chemical bonds in those clusters. In particular, we discuss Pd-H bonds associated with the H atoms with high coordination numbers and the difference of H-H bonds in the different Pd clusters from viewpoint of the electronic stress tensor. The notion of "pseudo-spindle structure" is proposed as the region between two atoms where the largest eigenvalue of the electronic stress tensor is negative and corresponding eigenvectors forming a pattern which connects them.

preprint2011arXiv

Inverted-sandwich-type and open-lantern-type dinuclear transition metal complexes: theoretical study of chemical bonds by electronic stress tensor

We study the electronic structure of two types of transition metal complexes, the inverted-sandwich-type and open-lantern-type, by the electronic stress tensor. In particular, the bond order b_e measured by the energy density which is defined from the electronic stress tensor is studied and compared with the conventional MO based bond order. We also examine the patterns found in the largest eigenvalue of the stress tensor and corresponding eigenvector field, the "spindle structure" and "pseudo-spindle structure". As for the inverted-sandwich-type complex, our bond order b_e calculation shows that relative strength of the metal-benzene bond among V, Cr and Mn complexes is V > Cr > Mn which is consistent with the MO based bond order. As for the open-lantern-type complex, we find that our energy density based bond order can properly describe the relative strength of Cr--Cr and Mo--Mo bonds by the surface integration of the energy density over the "Lagrange surface" which can take into account the spatial extent of the orbitals.

preprint2011arXiv

Microwave Emission from the Edgeworth-Kuiper Belt and the Asteroid Belt Constrained from WMAP

Objects in the Edgeworth-Kuiper belt and the main asteroid belt should emit microwaves that may give rise to extra anisotropy signals in the multipole of the cosmic microwave background (CMB) experiment. Constraints are derived from the absence of positive detection of such anisotropies for ell < 50, giving the total mass of Edgeworth-Kuiper belt objects to be smaller than 0.2 earth mass. This limit is consistent with the mass extrapolated from the observable population with the size of a > 15 km, assuming that the small-object population follows the power law in size dN/da ~ a^{-q} with the canonical index expected for collisional equilibrium, q ~ 3.5, with which 23% of the mass is ascribed to objects smaller than are observationally accessible down to grains. A similar argument applied to the main asteroid belt indicates that the grain population should not increase faster than q ~ 3.6 towards smaller radii, if it follows the power law continued to observed asteroids with larger radii. It is underlined that both cases are at or only slightly above the limit that can be physically significant, implying the importance of tightening further the CMB anisotropy limit, which may be attained with the observation at higher radio frequencies.

preprint2011arXiv

Theoretical Study of Hydrogenated Tetrahedral Aluminum Clusters

We report on the structures of aluminum hydrides derived from a tetrahedral aluminum Al4 cluster using ab initio quantum chemical calculation. Our calculation of binding energies of the aluminum hydrides reveals that stability of these hydrides increases as more hydrogen atoms are adsorbed, while stability of Al-H bonds decreases. We also analyze and discuss the chemical bonds of those clusters by using recently developed method based on the electronic stress tensor.

preprint2009arXiv

Constraining the time variation of the coupling constants from cosmic microwave background: effect of Λ_{QCD}

We investigate constraints on the time variation of the fine structure constant between the recombination epoch and the present epoch, Δα/α\equiv (α_{rec} - α_{now})/α_{now}, from cosmic microwave background (CMB) taking into account simultaneous variation of other physical constants, namely the electron mass m_{e} and the proton mass m_{p}. In other words, we consider the variation of Yukawa coupling and the QCD scale Λ_{QCD} in addition to the electromagnetic coupling. We clarify which parameters can be determined from CMB temperature anisotropy in terms of singular value decomposition. Assuming a relation among variations of coupling constants governed by a single scalar field (the dilaton), the 95% confidence level (C.L.) constraint on Δα/αis found to be -8.28 \times 10^{-3} < Δα/α< 1.81 \times 10^{-3}, which is tighter than the one obtained by considering only the change of αand m_{e}. We also obtain the constraint on the time variation of the proton-to-electron mass ratio μ\equiv m_{p}/m_{e} to be -0.52 < Δμ/μ< 0.17 (95% C.L.) under the same assumption. Finally, we also implement a forecast for constraints from the PLANCK survey.

preprint2009arXiv

Measuring the History of Cosmic Reionization using the 21-cm PDF from Simulations

The 21-cm PDF (i.e., distribution of pixel brightness temperatures) is expected to be highly non-Gaussian during reionization and to provide important information on the distribution of density and ionization. We measure the 21-cm PDF as a function of redshift in a large simulation of cosmic reionization and propose a simple empirical fit. Guided by the simulated PDF, we then carry out a maximum likelihood analysis of the ability of upcoming experiments to measure the shape of the 21-cm PDF and derive from it the cosmic reionization history. Under the strongest assumptions, we find that upcoming experiments can measure the reionization history in the mid to late stages of reionization to 1-10% accuracy. Under a more flexible approach that allows for four free parameters at each redshift, a similar accuracy requires the lower noise levels of second-generation 21-cm experiments.

preprint2009arXiv

Neutrino mass from cosmology: Impact of high-accuracy measurement of the Hubble constant

Non-zero neutrino mass would affect the evolution of the Universe in observable ways, and a strong constraint on the mass can be achieved using combinations of cosmological data sets. We focus on the power spectrum of cosmic microwave background (CMB) anisotropies, the Hubble constant H_0, and the length scale for baryon acoustic oscillations (BAO) to investigate the constraint on the neutrino mass, m_nu. We analyze data from multiple existing CMB studies (WMAP5, ACBAR, CBI, BOOMERANG, and QUAD), recent measurement of H_0 (SHOES), with about two times lower uncertainty (5%) than previous estimates, and recent treatments of BAO from the Sloan Digital Sky Survey (SDSS). We obtained an upper limit of m_nu < 0.2eV (95% C.L.), for a flat LambdaCDM model. This is a 40% reduction in the limit derived from previous H_0 estimates and one-third lower than can be achieved with extant CMB and BAO data. We also analyze the impact of smaller uncertainty on measurements of H_0 as may be anticipated in the near term, in combination with CMB data from the Planck mission, and BAO data from the SDSS/BOSS program. We demonstrate the possibility of a 5 sigma detection for a fiducial neutrino mass of 0.1eV or a 95% upper limit of 0.04eV for a fiducial of m_nu = 0eV. These constraints are about 50% better than those achieved without external constraint. We further investigate the impact on modeling where the dark-energy equation of state is constant but not necessarily -1, or where a non-flat universe is allowed. In these cases, the next-generation accuracies of Planck, BOSS, and 1% measurement of H_0 would all be required to obtain the limit m_nu < 0.05 - 0.06eV (95% C.L.) for the fiducial of m_nu = 0eV. The independence of systematics argues for pursuit of both BAO and H_0 measurements.

preprint2007arXiv

Dark Energy Parametrizations and the Curvature of the Universe

We investigate observational constraints on the curvature of the universe not restricting ourselves to a cosmological constant as dark energy, in particular allowing a dark energy equation of state to evolve with time in several ways. We use type Ia supernovae (SNeIa) data from the latest gold data set which includes 182 SNeIa, along with the CMB shift parameter and the baryon acoustic oscillation peak. We show quantitatively that the constraint on the curvature of the universe depends on dark energy model: some popular parametrizations give constraints closely around the flat universe at 5% level (2 sigma C.L.) whereas some parametrizations allow the universe to be as open as Omega_k ~ 0.2.

preprint2007arXiv

Increasing the effective number of neutrinos with decaying particles

We present models of decaying particles for increasing the effective number of neutrinos N_nu after big bang nucleosynthesis but before the structure formation begins. We point out that our scenario not only solves the discrepancy between the constraints on N_nu from these two epochs, but also provides a possible answer to deeper inconsistency in the estimation of the matter power spectrum amplitude at small scales, represented by sigma_8, between the WMAP and some small scale matter power measurements such as the Lyman-alpha forest and weak lensing. We consider (a) saxion decay into two axions; (b) gravitino decay into axino and axion; (c) Dirac right-handed sneutrino decay into gravitino and right-handed neutrino.

preprint2005arXiv

Constraining Neutrino Masses by CMB Experiments Alone

It is shown that a subelectronvolt upper limit can be derived on the neutrino mass from the CMB data alone in the Lambda CDM model with the power-law adiabatic perturbations, without the aid of any other cosmological data. Assuming the flatness of the universe, the constraint we can derive from the current WMAP observations is \sum m_nu < 2.0 eV at the 95% confidence level for the sum over three species of neutrinos (m_nu < 0.66 eV for the degenerate neutrinos). This constraint modifies little even if we abandon the flatness assumption for the spatial curvature. We argue that it would be difficult to improve the limit much beyond \sum m_nu \lesssim 1.5 eV using only the CMB data, even if their statistics are substantially improved. However, a significant improvement of the limit is possible if an external input is introduced that constrains the Hubble constant from below. The parameter correlation and the mechanism of CMB perturbations that give rise to the limit on the neutrino mass are also elucidated.