Showing 1–50 of 371 results for author: Ishikawa, K

Strong Field Optical Hall Effect in 2D Weyl Semimetal

Authors: M. Umar Farooq, Arqum Hashmi, Mizuki Tani, Kazuhiro Yabana, Kenichi L. Ishikawa, Li Huang, Tomohito Otobe

Abstract: The study of interplay between the geometric nature of Bloch electrons and transverse responses under strong field offers new opportunities for optoelectronic applications. Here, we present a comprehensive study of the strong-field response of Weyl Dirac nodes in bilayer T'-WTe2 using time-dependent first-principles formalism. The electron dynamics is explored focusing on the mid-infrared frequenc… ▽ More The study of interplay between the geometric nature of Bloch electrons and transverse responses under strong field offers new opportunities for optoelectronic applications. Here, we present a comprehensive study of the strong-field response of Weyl Dirac nodes in bilayer T'-WTe2 using time-dependent first-principles formalism. The electron dynamics is explored focusing on the mid-infrared frequency, ranging from the perturbative to nonperturbative regime. In the nonperturbative regime, the high-harmonic generation (HHG) spectra under a strong field clearly exhibit a plateau and energy cutoffs for both longitudinal and anomalous Hall (transverse) currents, with the latter being due to the large interband Berry curvature of the Weyl-Dirac semimetal. For the longitudinal harmonics, the intraband contributions increase with intensity, resulting in a complex interplay between interband polarization and intraband motions. Remarkably, if we take a comprehensive all-band perspective enabled by time-dependent density functional calculations, the anomalous Hall responses are purely attributed to the interband processes, even in the nonperturbative regime, thus Hall HHG can be crucial to understand the carrier dynamics. Our findings suggest that HHG associated with the ultrafast strong-field driven electron dynamics holds immense potential for exploring the nonlinear high Hall responses in Weyl semimetal. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 19 pages, 5 figures

Non-perturbative determination of the ${\cal N} = 1$ SUSY Yang-Mills gluino condensate at large $N$

Authors: Claudio Bonanno, Pietro Butti, Margarita García Pérez, Antonio González-Arroyo, Ken-Ichi Ishikawa, Masanori Okawa

Abstract: We present the first non-perturbative large-$N$ calculation of the $\mathcal{N}=1$ supersymmetric (SUSY) $\mathrm{SU}(N)$ Yang-Mills gluino condensate, obtained by means of numerical simulations of the lattice-discretized theory, and exploiting large-$N$ twisted volume reduction. We present two different determinations based, respectively, on the Banks-Casher formula and on the Gell-Mann-Oakes-Ren… ▽ More We present the first non-perturbative large-$N$ calculation of the $\mathcal{N}=1$ supersymmetric (SUSY) $\mathrm{SU}(N)$ Yang-Mills gluino condensate, obtained by means of numerical simulations of the lattice-discretized theory, and exploiting large-$N$ twisted volume reduction. We present two different determinations based, respectively, on the Banks-Casher formula and on the Gell-Mann-Oakes-Renner relation, both giving perfectly consistent results. By matching the lattice and the Novikov-Shifman-Vainshtein-Zakharov (NSVZ) regularization schemes, we are able for the first time to compare numerical and analytic computations. Our most accurate determination gives $Σ_{\rm RGI} /Λ_{\rm NSVZ}^3 = 1.55(57)$, in very good agreement with the $N$-dependence and the value predicted by the weak coupling instanton-based approach $Σ_{\rm RGI} /Λ_{\rm NSVZ}^3 = 1$. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: Main text: 5 pages, 3 figures. Supplemental Material: 4 pages, 3 figures

Report number: IFT-UAM/CSIC-24-90, HUPD-2404

Overlap integral of stationary scattering states

Authors: Kenzo Ishikawa, Yuya Nishio

Abstract: The overlap integrals of scattering states in potentials of finite widths are expressed with their asymptotic behaviors and those of energies $E_1$ and $E_2$ consist of diagonal terms that are proportional to $δ(E_1-E_2)$ and nondiagonal terms. Owing to the composition of nondiagonal terms, superpositions of stationary states have time-dependent norms and finite probability currents. These do not… ▽ More The overlap integrals of scattering states in potentials of finite widths are expressed with their asymptotic behaviors and those of energies $E_1$ and $E_2$ consist of diagonal terms that are proportional to $δ(E_1-E_2)$ and nondiagonal terms. Owing to the composition of nondiagonal terms, superpositions of stationary states have time-dependent norms and finite probability currents. These do not represent isolate states. In various exceptional potentials and in free theory, nondiagonal terms do not exist, and the superpositions of states with different energies represent isolate particles that exactly describe scattering processes. △ Less

Submitted 30 June, 2024; v1 submitted 21 April, 2024; originally announced June 2024.

Comments: 47 pages. 4 figures

Magnetization without spin: effective Lagrangian of itinerant electrons

Authors: Kenzo Ishikawa

Abstract: Effective Lagrangian of itinerant electron system of finite density at finite magnetic field is found to include Chern-Simons term of electromagnetic potentials of lower scale dimension than those studied before. This term has an origin in many-body wave function and unique topological property that is independent of a spin degree of freedom. The coupling strength is proportional to $\fracρ{eB}$,… ▽ More Effective Lagrangian of itinerant electron system of finite density at finite magnetic field is found to include Chern-Simons term of electromagnetic potentials of lower scale dimension than those studied before. This term has an origin in many-body wave function and unique topological property that is independent of a spin degree of freedom. The coupling strength is proportional to $\fracρ{eB}$, which is singular at $B=0$ for a constant charge density. The effective Lagrangian at a finite $B$ represents physical effects at $ B \neq 0$ properly. A universal shift of the magnetic field known as Slater-Pauling curve is derived from the effective Lagrangian. △ Less

Submitted 10 June, 2024; v1 submitted 5 June, 2024; originally announced June 2024.

Comments: 25 pages, 2 figures

Computing the associated groups of quandles with tools from group homology theory

Authors: Katsumi Ishikawa

Abstract: We give a complete description of the associated group of any quandle as a central extension of the inner-automorphism group. As an application, we compute the second quandle homology groups of quandles of some families, including those of Alexander quandles. We give a complete description of the associated group of any quandle as a central extension of the inner-automorphism group. As an application, we compute the second quandle homology groups of quandles of some families, including those of Alexander quandles. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 16 pages

MSC Class: 20N02 (Primary) 20E22; 57K10 (Secondary)

Discretization effects on nucleon root-mean-square radii from lattice QCD at the physical point

Authors: Ryutaro Tsuji, Yasumichi Aoki, Ken-Ichi Ishikawa, Yoshinobu Kuramashi, Shoichi Sasaki, Eigo Shintani, Takeshi Yamazaki

Abstract: We present results for the axial-vector coupling and root-mean-square (RMS) radii of the nucleon obtained from 2+1 flavor lattice QCD at the physical point with a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared $O(a)$ improved Wilson-clover quark action and Iwasaki gauge action a… ▽ More We present results for the axial-vector coupling and root-mean-square (RMS) radii of the nucleon obtained from 2+1 flavor lattice QCD at the physical point with a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared $O(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $β$ = 1.82 and 2.00 corresponding to lattice spacings of 0.085 fm and 0.063 fm, respectively. We first evaluate the value of the axial-vector coupling of the nucleon ($g_A$). In addition, the isovector electric, magnetic and axial radii and magnetic moment from the corresponding form factors are also determined. Combining the results at $β=1.82$ and $2.00$, we finally discuss the finite lattice spacing effect. It was found that the effect on $g_A$ is kept smaller than the statistical error of 2% while the effect on the isovector radii was observed as a possible discretization error of about 10%, regardless of the channel. △ Less

Submitted 29 November, 2023; originally announced January 2024.

Comments: 7 pages, 1 figure, Proceeding for the 40th International Symposium on Lattice Field Theory, 31 July 2023 - 4 August 2023, Chicago, USA

$|V_{us}|$ from kaon semileptonic form factor in $N_f = 2+1$ QCD at the physical point on (10 fm)$^4$

Authors: Takeshi Yamazaki, Ken-ichi Ishikawa, Naruhito Ishizuka, Yoshinobu Kuramashi, Yusuke Namekawa, Yusuke Taniguchi, Naoya Ukita, Tomoteru Yoshié for PACS Collaboration

Abstract: We present a preliminary result of the kaon semileptonic form factor calculated at the smallest lattice spacing in the PACS10 configurations, whose physical volumes are more than (10 fm)$^4$ at the physical point. The configurations were generated using the Iwasaki gauge action and $N_f=2+1$ stout-smeared nonperturbatively $O(a)$ improved Wilson quark action at the three lattice spacings, 0.085, 0… ▽ More We present a preliminary result of the kaon semileptonic form factor calculated at the smallest lattice spacing in the PACS10 configurations, whose physical volumes are more than (10 fm)$^4$ at the physical point. The configurations were generated using the Iwasaki gauge action and $N_f=2+1$ stout-smeared nonperturbatively $O(a)$ improved Wilson quark action at the three lattice spacings, 0.085, 0.063, and 0.041 fm. The value of $|V_{us}|$ in the continuum limit is estimated from our results including the preliminary one. We compare our result of $|V_{us}|$ with the previous results and those through the kaon leptonic decay. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 7 pages, 5 figures, Proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023)

Multi-Objective Bayesian Optimization with Active Preference Learning

Authors: Ryota Ozaki, Kazuki Ishikawa, Youhei Kanzaki, Shinya Suzuki, Shion Takeno, Ichiro Takeuchi, Masayuki Karasuyama

Abstract: There are a lot of real-world black-box optimization problems that need to optimize multiple criteria simultaneously. However, in a multi-objective optimization (MOO) problem, identifying the whole Pareto front requires the prohibitive search cost, while in many practical scenarios, the decision maker (DM) only needs a specific solution among the set of the Pareto optimal solutions. We propose a B… ▽ More There are a lot of real-world black-box optimization problems that need to optimize multiple criteria simultaneously. However, in a multi-objective optimization (MOO) problem, identifying the whole Pareto front requires the prohibitive search cost, while in many practical scenarios, the decision maker (DM) only needs a specific solution among the set of the Pareto optimal solutions. We propose a Bayesian optimization (BO) approach to identifying the most preferred solution in the MOO with expensive objective functions, in which a Bayesian preference model of the DM is adaptively estimated by an interactive manner based on the two types of supervisions called the pairwise preference and improvement request. To explore the most preferred solution, we define an acquisition function in which the uncertainty both in the objective functions and the DM preference is incorporated. Further, to minimize the interaction cost with the DM, we also propose an active learning strategy for the preference estimation. We empirically demonstrate the effectiveness of our proposed method through the benchmark function optimization and the hyper-parameter optimization problems for machine learning models. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Nucleon form factors in $N_f=2+1$ lattice QCD at the physical point : finite lattice spacing effect on the root-mean-square radii

Authors: Ryutaro Tsuji, Yasumichi Aoki, Ken-Ichi Ishikawa, Yoshinobu Kuramashi, Shoichi Sasaki, Kohei Sato, Eigo Shintani, Hiromasa Watanabe, Takeshi Yamazaki

Abstract: We present results for the nucleon form factors: electric ($G_E$), magnetic ($G_M$), axial ($F_A$), induced pseudoscalar ($F_P$) and pseudoscalar ($G_P$) form factors, using the second PACS10 ensemble that is one of three sets of $2+1$ flavor lattice QCD configurations at physical quark masses in large spatial volumes (exceeding $(10\ \mathrm{fm})^3$). The second PACS10 gauge configurations are ge… ▽ More We present results for the nucleon form factors: electric ($G_E$), magnetic ($G_M$), axial ($F_A$), induced pseudoscalar ($F_P$) and pseudoscalar ($G_P$) form factors, using the second PACS10 ensemble that is one of three sets of $2+1$ flavor lattice QCD configurations at physical quark masses in large spatial volumes (exceeding $(10\ \mathrm{fm})^3$). The second PACS10 gauge configurations are generated by the PACS Collaboration with the six stout-smeared $O(a)$ improved Wilson quark action and Iwasaki gauge action at the second gauge coupling $β=2.00$ corresponding to the lattice spacing of $a=0.063$ fm. We determine the isovector electric, magnetic and axial radii and magnetic moment from the corresponding form factors, as well as the axial-vector coupling $g_A$. Combining our previous results for the coarser lattice spacing [E. Shintani et al., Phys. Rev. D99 (2019) 014510; Phys. Rev. D102 (2020) 019902 (erattum)], the finite lattice spacing effects on the isovector radii, magnetic moment and axial-vector coupling are investigated using the difference between the two results. It was found that the effect on $g_A$ is kept smaller than the statistical error of 2% while the effect on the isovector radii was observed as a possible discretization error of about 10%, regardless of the channel. We also report the partially conserved axial vector current (PCAC) relation using a set of nucleon three-point correlation functions in order to verify the effect by $O(a)$-improvement of the axial-vector current. △ Less

Submitted 5 June, 2024; v1 submitted 17 November, 2023; originally announced November 2023.

Comments: 100 pages, 48 figures

Report number: UTHEP-783, UTCCS-P-149, HUPD-2307, YITP-23-143

The chiral condensate at large $N$

Authors: Claudio Bonanno, Pietro Butti, Margarita García Peréz, Antonio González-Arroyo, Ken-Ichi Ishikawa, Masanori Okawa

Abstract: We present results for the large-$N$ limit of the chiral condensate computed from twisted reduced models. We followed a two-fold strategy, one constiting in extracting the condensate from the quark-mass dependence of the pion mass, the other consisting in extracting the condensate from the mode number of the Dirac operator. We present results for the large-$N$ limit of the chiral condensate computed from twisted reduced models. We followed a two-fold strategy, one constiting in extracting the condensate from the quark-mass dependence of the pion mass, the other consisting in extracting the condensate from the mode number of the Dirac operator. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: 11 pages, 5 figures, contribution for the 40th International Symposium on Lattice Field Theory (Lattice 2023), July 31st - August 4th, 2023, Fermi National Accelerator Laboratory

Acousto-optic reconstruction of exterior sound field based on concentric circle sampling with circular harmonic expansion

Authors: Phuc Duc Nguyen, Kenji Ishikawa, Noboru Harada, Takehiro Moriya

Abstract: Acousto-optic sensing provides an alternative approach to traditional microphone arrays by shedding light on the interaction of light with an acoustic field. Sound field reconstruction is a fascinating and advanced technique used in acousto-optics sensing. Current challenges in sound-field reconstruction methods pertain to scenarios in which the sound source is located within the reconstruction ar… ▽ More Acousto-optic sensing provides an alternative approach to traditional microphone arrays by shedding light on the interaction of light with an acoustic field. Sound field reconstruction is a fascinating and advanced technique used in acousto-optics sensing. Current challenges in sound-field reconstruction methods pertain to scenarios in which the sound source is located within the reconstruction area, known as the exterior problem. Existing reconstruction algorithms, primarily designed for interior scenarios, often exhibit suboptimal performance when applied to exterior cases. This paper introduces a novel technique for exterior sound-field reconstruction. The proposed method leverages concentric circle sampling and a two-dimensional exterior sound-field reconstruction approach based on circular harmonic extensions. To evaluate the efficacy of this approach, both numerical simulations and practical experiments are conducted. The results highlight the superior accuracy of the proposed method when compared to conventional reconstruction methods, all while utilizing a minimal amount of measured projection data. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

Twisted Alexander vanishing order of knots

Authors: Katsumi Ishikawa, Takayuki Morifuji, Masaaki Suzuki

Abstract: Based on a vanishing theorem for non-fibered knots due to Friedl and Vidussi, we define the twisted Alexander vanishing order of a knot to be the order of the smallest finite group such that the corresponding twisted Alexander polynomial is zero. In this paper, we show its basic properties, and provide several explicit values for knots with $10$ or fewer crossings. Moreover, we characterize a fini… ▽ More Based on a vanishing theorem for non-fibered knots due to Friedl and Vidussi, we define the twisted Alexander vanishing order of a knot to be the order of the smallest finite group such that the corresponding twisted Alexander polynomial is zero. In this paper, we show its basic properties, and provide several explicit values for knots with $10$ or fewer crossings. Moreover, we characterize a finite group admitting the zero-twisted Alexander polynomial. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 16 pages, 1 figure

MSC Class: Primary 57K14; Secondary 57K10

Control of ion-photoelectron entanglement and coherence via Rabi oscillations

Authors: Kenichi L. Ishikawa, Kevin C. Prince, Kiyoshi Ueda

Abstract: We report a theoretical investigation of photoionization by a pair of coherent, ultrashort fundamental and second-harmonic extreme-ultraviolet pulses, where the photon energies are selected to yield the same photoelectron energy for ionization of two different sub-shells. This choice implies that the fundamental energy is equal to the difference in energy of the ionic states, and they are therefor… ▽ More We report a theoretical investigation of photoionization by a pair of coherent, ultrashort fundamental and second-harmonic extreme-ultraviolet pulses, where the photon energies are selected to yield the same photoelectron energy for ionization of two different sub-shells. This choice implies that the fundamental energy is equal to the difference in energy of the ionic states, and they are therefore coupled by the fundamental photon. By deriving analytical expressions using the essential-states approach, we show that this Rabi coupling creates coherence between the two photoelectron wave packets, which would otherwise be incoherent. We analyze how the coupling is affected by the parameters such as relative phase, pulse width, delay between the two pulses, the Rabi coupling strength, and photoelectron energy. Our discussion mostly considers Ne 2p and 2s photoionization, but it is generally valid for many other quantum systems where photoionization from two different shells is observed. △ Less

Submitted 23 November, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: 28 pages, 16 figures

On the Parallel Complexity of Multilevel Monte Carlo in Stochastic Gradient Descent

Authors: Kei Ishikawa

Abstract: In the stochastic gradient descent (SGD) for sequential simulations such as the neural stochastic differential equations, the Multilevel Monte Carlo (MLMC) method is known to offer better theoretical computational complexity compared to the naive Monte Carlo approach. However, in practice, MLMC scales poorly on massively parallel computing platforms such as modern GPUs, because of its large parall… ▽ More In the stochastic gradient descent (SGD) for sequential simulations such as the neural stochastic differential equations, the Multilevel Monte Carlo (MLMC) method is known to offer better theoretical computational complexity compared to the naive Monte Carlo approach. However, in practice, MLMC scales poorly on massively parallel computing platforms such as modern GPUs, because of its large parallel complexity which is equivalent to that of the naive Monte Carlo method. To cope with this issue, we propose the delayed MLMC gradient estimator that drastically reduces the parallel complexity of MLMC by recycling previously computed gradient components from earlier steps of SGD. The proposed estimator provably reduces the average parallel complexity per iteration at the cost of a slightly worse per-iteration convergence rate. In our numerical experiments, we use an example of deep hedging to demonstrate the superior parallel complexity of our method compared to the standard MLMC in SGD. △ Less

Submitted 10 October, 2023; v1 submitted 3 October, 2023; originally announced October 2023.

Comments: Fixed a typo in the title and added acknowledgement

The large-$N$ limit of the chiral condensate from twisted reduced models

Authors: Claudio Bonanno, Pietro Butti, Margarita García Peréz, Antonio González-Arroyo, Ken-Ichi Ishikawa, Masanori Okawa

Abstract: We compute the large-$N$ limit of the QCD chiral condensate on the lattice using twisted reduced models, and performing controlled continuum and chiral extrapolations. We perform two different calculations: one consists in extracting the chiral condensate from the quark mass dependence of the pion mass, and the other consists in extracting the chiral condensate from the behaviour of the mode numbe… ▽ More We compute the large-$N$ limit of the QCD chiral condensate on the lattice using twisted reduced models, and performing controlled continuum and chiral extrapolations. We perform two different calculations: one consists in extracting the chiral condensate from the quark mass dependence of the pion mass, and the other consists in extracting the chiral condensate from the behaviour of the mode number of the Wilson-Dirac operator for small eigenvalues. We find consistency between the results of the two methods, giving a joint estimate of $\lim_{N\to\infty}Σ(N)/N=[184(13)$ MeV$]^3$ ($\overline{\mathrm{MS}}$, $μ=2$ GeV, taking the square root of the string tension $\sqrtσ=440$ MeV to set the scale), in remarkable agreement with the $\mathrm{SU}(3)$ 2-flavor FLAG result. △ Less

Submitted 4 December, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 25 pages, 11 figures. v2: minor changes, matches accepted version on JHEP

Report number: IFT-UAM/CSIC-23-118, HUPD-2303

Transformation of polar nematic phases in the presence of electric field

Authors: A. V. Emelyanenko, V. Yu. Rudyak, F. Araoka, H. Nishikawa, K. Ishikawa

Abstract: Only a few years have passed since discovery of polar nematics, and now they are becoming the most actively studied liquid crystal materials. Despite numerous breakthrough findings made recently, a theoretical systematization is still lacking. In the present paper we are making a step on the way of systematization. A powerful technique that molecular-statistical physics is has been applied to an a… ▽ More Only a few years have passed since discovery of polar nematics, and now they are becoming the most actively studied liquid crystal materials. Despite numerous breakthrough findings made recently, a theoretical systematization is still lacking. In the present paper we are making a step on the way of systematization. A powerful technique that molecular-statistical physics is has been applied to an assembly of polar molecules influenced by electric field. Totally, the three polar nematic phases were found to be stable at various conditions: the double-splay ferroelectric nematic $N_F^{2D}$ (observed in the lower-temperature range in the absence or at low electric field), the double-splay antiferroelectric nematic $N_{AF}$ (observed at intermediate temperature in the absence or at low electric field) and the single-splay ferroelectric nematic $N_F^{1D}$ (observed at moderate electric field at any temperature below transition into paraelectric nematic $N$ and in the higher-temperature range (also below $N$) at low electric field or without it. A paradoxal transition from $N_F^{1D}$ to $N$ induced by application of higher electric field has been found and explained. A transformation of the structure of polar nematic phases at application of electric field has also been investigated by Monte Carlo simulations and experimentally by observation of POM images. In particular, it has been realized that, at planar anchoring, $N_{AF}$ in the presence of moderate out-of-plane electric field exhibits the twofold splay modulation: antiferroelectric in the plane of the substrate and ferroelectric in the plane normal to the substrate. Several additional sub-transitions related to fitting confined geometry of the cell by the structure of polar phases were detected. △ Less

Submitted 26 September, 2023; originally announced September 2023.

A Convex Framework for Confounding Robust Inference

Authors: Kei Ishikawa, Niao He, Takafumi Kanamori

Abstract: We study policy evaluation of offline contextual bandits subject to unobserved confounders. Sensitivity analysis methods are commonly used to estimate the policy value under the worst-case confounding over a given uncertainty set. However, existing work often resorts to some coarse relaxation of the uncertainty set for the sake of tractability, leading to overly conservative estimation of the poli… ▽ More We study policy evaluation of offline contextual bandits subject to unobserved confounders. Sensitivity analysis methods are commonly used to estimate the policy value under the worst-case confounding over a given uncertainty set. However, existing work often resorts to some coarse relaxation of the uncertainty set for the sake of tractability, leading to overly conservative estimation of the policy value. In this paper, we propose a general estimator that provides a sharp lower bound of the policy value using convex programming. The generality of our estimator enables various extensions such as sensitivity analysis with f-divergence, model selection with cross validation and information criterion, and robust policy learning with the sharp lower bound. Furthermore, our estimation method can be reformulated as an empirical risk minimization problem thanks to the strong duality, which enables us to provide strong theoretical guarantees of the proposed estimator using techniques of the M-estimation. △ Less

Submitted 1 November, 2023; v1 submitted 21 September, 2023; originally announced September 2023.

Comments: This is an extended version of the following work https://proceedings.mlr.press/v206/ishikawa23a.html. arXiv admin note: text overlap with arXiv:2302.13348

Wave-Packet Effects: A Solution for Isospin Anomalies in Vector-Meson Decay

Authors: Kenzo Ishikawa, Osamu Jinnouchi, Kenji Nishiwaki, Kin-ya Oda

Abstract: There is a long-standing anomaly in the ratio of the decay width for $ψ(3770)\to D^0\overline{D^0}$ to that for $ψ(3770)\to D^+D^-$ at the level of $9.5\,σ$. A similar anomaly exists for the ratio of $φ(1020)\to K_\text{L}^0K_\text{S}^0$ to $φ(1020)\to K^+K^-$ at $2.1\,σ$. In this study, we reassess the anomaly through the lens of Gaussian wave-packet formalism. Our comprehensive calculations incl… ▽ More There is a long-standing anomaly in the ratio of the decay width for $ψ(3770)\to D^0\overline{D^0}$ to that for $ψ(3770)\to D^+D^-$ at the level of $9.5\,σ$. A similar anomaly exists for the ratio of $φ(1020)\to K_\text{L}^0K_\text{S}^0$ to $φ(1020)\to K^+K^-$ at $2.1\,σ$. In this study, we reassess the anomaly through the lens of Gaussian wave-packet formalism. Our comprehensive calculations include the localization of the overlap of the wave packets near the mass thresholds as well as the composite nature of the initial-state vector mesons. The results align within $\sim 1 σ$ confidence level with the Particle Data Group's central values for a physically reasonable value of the form-factor parameter, indicating a resolution to these anomalies. We also check the deviation of a wave-packet resonance from the Briet-Wigner shape and find that wide ranges of the wave-packet size are consistent with the experimental data. △ Less

Submitted 24 October, 2023; v1 submitted 19 August, 2023; originally announced August 2023.

Comments: 50 pages, 16 figures; footnote 23 added, minor typos corrected, accepted for publication in EPJC

Journal ref: The European Physical Journal C, 83, Article number: 978 (2023)

Development of pericardial fat count images using a combination of three different deep-learning models

Authors: Takaaki Matsunaga, Atsushi Kono, Hidetoshi Matsuo, Kaoru Kitagawa, Mizuho Nishio, Hiromi Hashimura, Yu Izawa, Takayoshi Toba, Kazuki Ishikawa, Akie Katsuki, Kazuyuki Ohmura, Takamichi Murakami

Abstract: Rationale and Objectives: Pericardial fat (PF), the thoracic visceral fat surrounding the heart, promotes the development of coronary artery disease by inducing inflammation of the coronary arteries. For evaluating PF, this study aimed to generate pericardial fat count images (PFCIs) from chest radiographs (CXRs) using a dedicated deep-learning model. Materials and Methods: The data of 269 conse… ▽ More Rationale and Objectives: Pericardial fat (PF), the thoracic visceral fat surrounding the heart, promotes the development of coronary artery disease by inducing inflammation of the coronary arteries. For evaluating PF, this study aimed to generate pericardial fat count images (PFCIs) from chest radiographs (CXRs) using a dedicated deep-learning model. Materials and Methods: The data of 269 consecutive patients who underwent coronary computed tomography (CT) were reviewed. Patients with metal implants, pleural effusion, history of thoracic surgery, or that of malignancy were excluded. Thus, the data of 191 patients were used. PFCIs were generated from the projection of three-dimensional CT images, where fat accumulation was represented by a high pixel value. Three different deep-learning models, including CycleGAN, were combined in the proposed method to generate PFCIs from CXRs. A single CycleGAN-based model was used to generate PFCIs from CXRs for comparison with the proposed method. To evaluate the image quality of the generated PFCIs, structural similarity index measure (SSIM), mean squared error (MSE), and mean absolute error (MAE) of (i) the PFCI generated using the proposed method and (ii) the PFCI generated using the single model were compared. Results: The mean SSIM, MSE, and MAE were as follows: 0.856, 0.0128, and 0.0357, respectively, for the proposed model; and 0.762, 0.0198, and 0.0504, respectively, for the single CycleGAN-based model. Conclusion: PFCIs generated from CXRs with the proposed model showed better performance than those with the single model. PFCI evaluation without CT may be possible with the proposed method. △ Less

Submitted 25 July, 2023; v1 submitted 23 July, 2023; originally announced July 2023.

Decay Pattern of Pygmy States Observed in Neutron-Rich 26 Ne

Authors: J. Gibelin, D. Beaumel, T. Motobayashi, Y. Blumenfeld, N. Aoi, H. Baba, Z. Elekes, S. Fortier, N. Frascaria, N. Fukuda, T. Gomi, K. Ishikawa, Y. Kondo, T. Kubo, V. Lima, T. Nakamura, A. Saito, Y. Satou, J. -A. Scarpaci, E. Takeshita, S. Takeuchi, T. Teranishi, Y. Togano, A. M. Vinodkumar, Y. Yanagisawa , et al. (1 additional authors not shown)

Abstract: Coulomb excitation of the exotic neutron-rich nucleus Ne26 on a Pb208 target was measured at 58 MeV/u in order to search for low-lying E1 strength above the neutron emission threshold. This radioactive beam experiment was carried out at the RIKEN Accelerator Research Facility. Using the invariant mass method in the Ne25+n channel, we observe a sizable amount of E1 strength between 6 and 10 MeV exc… ▽ More Coulomb excitation of the exotic neutron-rich nucleus Ne26 on a Pb208 target was measured at 58 MeV/u in order to search for low-lying E1 strength above the neutron emission threshold. This radioactive beam experiment was carried out at the RIKEN Accelerator Research Facility. Using the invariant mass method in the Ne25+n channel, we observe a sizable amount of E1 strength between 6 and 10 MeV excitation energy. By performing a multipole decomposition of the differential cross section, a reduced dipole transition probability of B(E1)=0.49+-0.16e2fm2 is deduced, corresponding to 4.9+-1.6% of the Thomas-Reiche-Kuhn sum rule. For the first time, the decay pattern of low-lying strength in a neutron-rich nucleus is measured. The extracted decay pattern is not consistent with several mean-field theory descriptions of the pygmy states. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: arXiv admin note: text overlap with arXiv:0705.1753

Robustness of Baryon Acoustic Oscillations Measurements with Photometric Redshift Uncertainties

Authors: Keitaro Ishikawa, Tomomi Sunayama, Atsushi J. Nishizawa, Hironao Miyatake, Takahiro Nishimichi

Abstract: We investigate the robustness of baryon acoustic oscillations (BAO) measurements with a photometric galaxy sample using mock galaxy catalogues with various sizes of photometric redshift (photo-$z$) uncertainties. We first conduct the robustness of BAO measurements, assuming we have a perfect knowledge of photo-$z$ uncertainties. We find that the BAO shift parameter $α$ can be constrained in an unb… ▽ More We investigate the robustness of baryon acoustic oscillations (BAO) measurements with a photometric galaxy sample using mock galaxy catalogues with various sizes of photometric redshift (photo-$z$) uncertainties. We first conduct the robustness of BAO measurements, assuming we have a perfect knowledge of photo-$z$ uncertainties. We find that the BAO shift parameter $α$ can be constrained in an unbiased manner even for 3% photometric redshift uncertainties up to $z\sim 1$. For instance, $α=1.006 \pm 0.078$ with 95% confidence level is obtained from 3% photo-$z$ uncertainty data at $z=1.03$ using the sample of $M_* \ge 10^{10.25} M_{\odot}/h^2$. We also find that a sparse galaxy sample, e.g. $<2\times10^{-4}$ [$h$ Mpc$^{-1}]^3$ causes additional noise in the covariance matrix calculation and can bias the constraint on $α$. Following this, we look into the scenario where incorrect photometric redshift uncertainties are assumed in the fitting model. We find that underestimating the photo-$z$ uncertainty leads to a degradation in the constraining power on $α$. However, the constrained value of $α$ is not biased. We also quantify the constraining power on $Ω_{\rm m0}$ assuming the LSST-like covariance and find that the 95% confidence level is $σ(Ω_{\rm m0})\sim0.03$-$0.05$ corresponding to the photo-$z$ uncertainties of 1% to 3% respectively. Finally, we examine whether the skewness in the photometric redshift can bias the constraint on $α$ and confirm that the constraint on $α$ is unbiased, even assuming a Gaussian photo-$z$ uncertainty in our model. △ Less

Submitted 13 July, 2023; v1 submitted 2 June, 2023; originally announced June 2023.

Comments: 11 pages, 11 figures, 3 tables, submitted to MNRAS

Potential scattering in $L^2$ space: (2) Rigorous scattering probability of wave packets

Authors: Kenzo Ishikawa

Abstract: In this study, potential scatterings are formulated in experimental setups with Gaussian wave packets in accordance with a probability principle and associativity of products. A breaking of an associativity is observed in scalar products with stationary scattering states in a majority of short-range potentials. Due to the breaking, states of different energies are not orthogonal and their superpos… ▽ More In this study, potential scatterings are formulated in experimental setups with Gaussian wave packets in accordance with a probability principle and associativity of products. A breaking of an associativity is observed in scalar products with stationary scattering states in a majority of short-range potentials. Due to the breaking, states of different energies are not orthogonal and their superposition is not suitable for representing a normalized isolate state. Free wave packets in perturbative expansions in coupling strengths keep the associativity, and give a rigorous amplitude that preserves manifest unitarity and other principles of the quantum mechanics. An absolute probability is finite and comprises cross sections and new terms of unique properties. The results also demonstrate an interference term displaying unique behavior at an extreme forward direction. △ Less

Submitted 11 December, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 54 pages

Journal ref: annals of physics, 460,january 2024,169571

Potential scatterings in $L^2$ space: (1) non-orthogonality of stationary states

Authors: Kenzo Ishikawa

Abstract: Orthogonality of eigenstates of different energies and its implications in potential scattering are unlabeled. Scalar products of scattering states of different energies are found to have finite non-orthogonal terms in potentials of finite widths. Their superpositions have time-dependent norms, and are not suitable for isolate states. In these systems, a perturbative method and a variational metho… ▽ More Orthogonality of eigenstates of different energies and its implications in potential scattering are unlabeled. Scalar products of scattering states of different energies are found to have finite non-orthogonal terms in potentials of finite widths. Their superpositions have time-dependent norms, and are not suitable for isolate states. In these systems, a perturbative method and a variational method are viable methods for finding a rigorous transition probability that describes phenomena completely. In various exceptional potentials, an orthogonality is satisfied. △ Less

Submitted 30 June, 2024; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: Replaced with arXiv 2406.03595

Deep sound-field denoiser: optically-measured sound-field denoising using deep neural network

Authors: Kenji Ishikawa, Daiki Takeuchi, Noboru Harada, Takehiro Moriya

Abstract: This paper proposes a deep sound-field denoiser, a deep neural network (DNN) based denoising of optically measured sound-field images. Sound-field imaging using optical methods has gained considerable attention due to its ability to achieve high-spatial-resolution imaging of acoustic phenomena that conventional acoustic sensors cannot accomplish. However, the optically measured sound-field images… ▽ More This paper proposes a deep sound-field denoiser, a deep neural network (DNN) based denoising of optically measured sound-field images. Sound-field imaging using optical methods has gained considerable attention due to its ability to achieve high-spatial-resolution imaging of acoustic phenomena that conventional acoustic sensors cannot accomplish. However, the optically measured sound-field images are often heavily contaminated by noise because of the low sensitivity of optical interferometric measurements to airborne sound. Here, we propose a DNN-based sound-field denoising method. Time-varying sound-field image sequences are decomposed into harmonic complex-amplitude images by using a time-directional Fourier transform. The complex images are converted into two-channel images consisting of real and imaginary parts and denoised by a nonlinear-activation-free network. The network is trained on a sound-field dataset obtained from numerical acoustic simulations with randomized parameters. We compared the method with conventional ones, such as image filters, a spatiotemporal filter, and other DNN architectures, on numerical and experimental data. The experimental data were measured by parallel phase-shifting interferometry and holographic speckle interferometry. The proposed deep sound-field denoiser significantly outperformed the conventional methods on both the numerical and experimental data. Code is available on GitHub: https://github.com/nttcslab/deep-sound-field-denoiser. △ Less

Submitted 21 September, 2023; v1 submitted 27 April, 2023; originally announced April 2023.

Comments: 16 pages, 10 figures, 2 tables

Hitomi/XRISM micro-calorimeter

Authors: Kosuke Sato, Yuusuke Uchida, Kumi Ishikawa

Abstract: We present an overview of the ASTRO-H (Hitomi) Soft X-Ray Spectrometer (SXS) and the X-Ray Imaging and Spectrometer Mission (XRISM) {\it Resolve} spectrometer. In each, a 36-pixel X-ray micro-calorimeter array operated at 50 mK covers a $3 \times 3$ arc-minute field of view. The instruments are designed to achieve an energy resolution of better than 7 eV over the 0.3 -- 12 keV energy range and ope… ▽ More We present an overview of the ASTRO-H (Hitomi) Soft X-Ray Spectrometer (SXS) and the X-Ray Imaging and Spectrometer Mission (XRISM) {\it Resolve} spectrometer. In each, a 36-pixel X-ray micro-calorimeter array operated at 50 mK covers a $3 \times 3$ arc-minute field of view. The instruments are designed to achieve an energy resolution of better than 7 eV over the 0.3 -- 12 keV energy range and operate for more than 3 years in orbit. Actually, the SXS achieved the energy resolution of $\sim$5 eV in orbit, but it was lost after only a month of operation due to the loss of spacecraft attitude control. For the recovery mission, XRISM will be equipped with the {\it Resolve} spectrometer which has mostly the same design as SXS and is expected to have the same in-flight performance. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: Invited chapter for the edited book {\it High-Resolution X-Ray Spectroscopy: Instrumentation, Data Analysis, and Science} (Eds. C. Bambi and J. Jiang, Springer Singapore, expected in 2023)

Two-tone colorings and surjective dihedral representations for links

Authors: Kazuhiro Ichihara, Katsumi Ishikawa, Eri Matsudo, Masaaki Suzuki

Abstract: It is well-known that a knot is Fox $n$-colorable for a prime $n$ if and only if the knot group admits a surjective homomorphism to the dihedral group of degree $n$. However, this is not the case for links with two or more components. In this paper, we introduce a two-tone coloring on a link diagram, and give a condition for links so that the link groups admit surjective representations to the dih… ▽ More It is well-known that a knot is Fox $n$-colorable for a prime $n$ if and only if the knot group admits a surjective homomorphism to the dihedral group of degree $n$. However, this is not the case for links with two or more components. In this paper, we introduce a two-tone coloring on a link diagram, and give a condition for links so that the link groups admit surjective representations to the dihedral groups. In particular, it is shown that the link group of any link with at least 3 components admits a surjective homomorphism to the dihedral group of arbitrary degree. △ Less

Submitted 28 April, 2024; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: 13 pages, 6 figures. Several improvements based on the referees' comments

MSC Class: 57K10

Kernel Conditional Moment Constraints for Confounding Robust Inference

Authors: Kei Ishikawa, Niao He

Abstract: We study policy evaluation of offline contextual bandits subject to unobserved confounders. Sensitivity analysis methods are commonly used to estimate the policy value under the worst-case confounding over a given uncertainty set. However, existing work often resorts to some coarse relaxation of the uncertainty set for the sake of tractability, leading to overly conservative estimation of the poli… ▽ More We study policy evaluation of offline contextual bandits subject to unobserved confounders. Sensitivity analysis methods are commonly used to estimate the policy value under the worst-case confounding over a given uncertainty set. However, existing work often resorts to some coarse relaxation of the uncertainty set for the sake of tractability, leading to overly conservative estimation of the policy value. In this paper, we propose a general estimator that provides a sharp lower bound of the policy value. It can be shown that our estimator contains the recently proposed sharp estimator by Dorn and Guo (2022) as a special case, and our method enables a novel extension of the classical marginal sensitivity model using f-divergence. To construct our estimator, we leverage the kernel method to obtain a tractable approximation to the conditional moment constraints, which traditional non-sharp estimators failed to take into account. In the theoretical analysis, we provide a condition for the choice of the kernel which guarantees no specification error that biases the lower bound estimation. Furthermore, we provide consistency guarantees of policy evaluation and learning. In the experiments with synthetic and real-world data, we demonstrate the effectiveness of the proposed method. △ Less

Submitted 14 September, 2023; v1 submitted 26 February, 2023; originally announced February 2023.

Journal ref: AISTATS 2023

Efficient simulation of multielectron dynamics in molecules under intense laser pulses: Implementation of the multiconfiguration time-dependent Hartree-Fock method based on the adaptive finite element method

Authors: Yuki Orimo, Takeshi Sato, Kenichi L. Ishikawa

Abstract: We present an implementation of the multiconfiguration time-dependent Hartree-Fock method based on the adaptive finite element method for molecules under intense laser pulses. For efficient simulations, orbital functions are propagated by a stable propagator using the short iterative Arnoldi scheme and our implementation is parallelized for distributed memory computing. This is demonstrated by sim… ▽ More We present an implementation of the multiconfiguration time-dependent Hartree-Fock method based on the adaptive finite element method for molecules under intense laser pulses. For efficient simulations, orbital functions are propagated by a stable propagator using the short iterative Arnoldi scheme and our implementation is parallelized for distributed memory computing. This is demonstrated by simulating high-harmonic generation from a water molecule and achieves a simulation of multielectron dynamics with overwhelmingly less computational time, compared to our previous work. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Adjoint fermions at large-N on the lattice

Authors: P. Butti, M. García Pérez, A. González-Arroyo, K. I. Ishikawa

Abstract: Lattice simulations of Yang-Mills theories coupled with $N_f$ flavours of fermions in the adjoint representation provide a way to probe the non-perturbative regime of a plethora of different physical scenarios, such as Supersymmetric Yang-Mills theory to BSM models. We are interested in the large-$N_c$ limit of these theories, for which standard lattice techniques are limited by high-computational… ▽ More Lattice simulations of Yang-Mills theories coupled with $N_f$ flavours of fermions in the adjoint representation provide a way to probe the non-perturbative regime of a plethora of different physical scenarios, such as Supersymmetric Yang-Mills theory to BSM models. We are interested in the large-$N_c$ limit of these theories, for which standard lattice techniques are limited by high-computational costs. Our approach makes use of the well-established twisted volume reduction, which allows one to simulate these theories on a $1^4$ lattice. In this talk, we are going to present a detailed study of the scale setting of these theories, performed with Wilson flow techniques, but endowed with a procedure that allowed us to reduce finite-volume (finite $N_c$) systematic effects. We will apply this procedure to our configurations for $N_f=0,\frac{1}{2},1,2$ and analyze the dependence of the scale on the coupling, the adjoint fermion mass and the number of flavours. For the cases in which enough couplings are available we compared the resulting $β$-function with the predictions of perturbation theory, finding very good agreement. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Momentum transfer dependence of kaon semileptonic form factor on (10 fm)$^4$ at the physical point

Authors: Takeshi Yamazaki, Ken-ichi Ishikawa, Naruhito Ishizuka, Yoshinobu Kuramashi, Yusuke Namekawa, Yusuke Taniguchi, Naoya Ukita, Tomoteru Yoshié for PACS Collaboration

Abstract: We calculate the kaon semileptonic form factors using the two sets of the PACS10 configuration, whose physical volumes are more than (10 fm)$^4$ at the physical point. The lattice spacings are 0.063 and 0.085 fm. The configurations were generated using the Iwasaki gauge action and $N_f=2+1$ stout-smeared nonperturbatively $O(a)$-improved Wilson quark action. From the momentum transfer dependence o… ▽ More We calculate the kaon semileptonic form factors using the two sets of the PACS10 configuration, whose physical volumes are more than (10 fm)$^4$ at the physical point. The lattice spacings are 0.063 and 0.085 fm. The configurations were generated using the Iwasaki gauge action and $N_f=2+1$ stout-smeared nonperturbatively $O(a)$-improved Wilson quark action. From the momentum transfer dependence of the form factors, we evaluate the slope and curvature for the form factors at the zero momentum transfer. Furthermore, we calculate the phase space factor, which is used to obtain $|V_{us}|$ through the kaon semileptonic decay. These results are compared with previous lattice results and experimental values. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: 10 pages, 8 figures, Proceedings of the 39th International Symposium on Lattice Field Theory (LATTICE2022), 8th-13th August 2022, Bonn, Germany

Informative Sample-Aware Proxy for Deep Metric Learning

Authors: Aoyu Li, Ikuro Sato, Kohta Ishikawa, Rei Kawakami, Rio Yokota

Abstract: Among various supervised deep metric learning methods proxy-based approaches have achieved high retrieval accuracies. Proxies, which are class-representative points in an embedding space, receive updates based on proxy-sample similarities in a similar manner to sample representations. In existing methods, a relatively small number of samples can produce large gradient magnitudes (ie, hard samples)… ▽ More Among various supervised deep metric learning methods proxy-based approaches have achieved high retrieval accuracies. Proxies, which are class-representative points in an embedding space, receive updates based on proxy-sample similarities in a similar manner to sample representations. In existing methods, a relatively small number of samples can produce large gradient magnitudes (ie, hard samples), and a relatively large number of samples can produce small gradient magnitudes (ie, easy samples); these can play a major part in updates. Assuming that acquiring too much sensitivity to such extreme sets of samples would deteriorate the generalizability of a method, we propose a novel proxy-based method called Informative Sample-Aware Proxy (Proxy-ISA), which directly modifies a gradient weighting factor for each sample using a scheduled threshold function, so that the model is more sensitive to the informative samples. Extensive experiments on the CUB-200-2011, Cars-196, Stanford Online Products and In-shop Clothes Retrieval datasets demonstrate the superiority of Proxy-ISA compared with the state-of-the-art methods. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: Accepted at ACM Multimedia Asia (MMAsia) 2022

Time-dependent multiconfiguration self-consistent-field and time-dependent optimized coupled-cluster methods for intense laser-driven multielectron dynamics

Authors: Takeshi Sato, Himadri Pathak, Yuki Orimo, Kenichi L. Ishikawa

Abstract: We review time-dependent multiconfiguration self-consistent-field (TD-MCSCF) method and time-dependent optimized coupled-cluster (TD-OCC) method for first-principles simulations of high-field phenomena such as tunneling ionization and high-order harmonic generation in atoms and molecules irradiated by a strong laser field. These methods provide a flexible and systematically improvable description… ▽ More We review time-dependent multiconfiguration self-consistent-field (TD-MCSCF) method and time-dependent optimized coupled-cluster (TD-OCC) method for first-principles simulations of high-field phenomena such as tunneling ionization and high-order harmonic generation in atoms and molecules irradiated by a strong laser field. These methods provide a flexible and systematically improvable description of the multielectron dynamics by expressing the all-electron wavefunction by configuration interaction expansion or coupled-cluster expansion, using time-dependent one-electron orbital functions. The time-dependent variational principle plays a key role to derive these methods satisfying gauge invariance and Ehrenfest theorem. The real-time/real-space implementation with an absorbing boundary condition enables the simulation of high-field processes involving multiple excitation and ionization. We present a detailed, comprehensive discussion of such features of TD-MCSCF and TD-OCC methods. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Empirical Study on Optimizer Selection for Out-of-Distribution Generalization

Authors: Hiroki Naganuma, Kartik Ahuja, Shiro Takagi, Tetsuya Motokawa, Rio Yokota, Kohta Ishikawa, Ikuro Sato, Ioannis Mitliagkas

Abstract: Modern deep learning systems do not generalize well when the test data distribution is slightly different to the training data distribution. While much promising work has been accomplished to address this fragility, a systematic study of the role of optimizers and their out-of-distribution generalization performance has not been undertaken. In this study, we examine the performance of popular firs… ▽ More Modern deep learning systems do not generalize well when the test data distribution is slightly different to the training data distribution. While much promising work has been accomplished to address this fragility, a systematic study of the role of optimizers and their out-of-distribution generalization performance has not been undertaken. In this study, we examine the performance of popular first-order optimizers for different classes of distributional shift under empirical risk minimization and invariant risk minimization. We address this question for image and text classification using DomainBed, WILDS, and Backgrounds Challenge as testbeds for studying different types of shifts -- namely correlation and diversity shift. We search over a wide range of hyperparameters and examine classification accuracy (in-distribution and out-of-distribution) for over 20,000 models. We arrive at the following findings, which we expect to be helpful for practitioners: i) adaptive optimizers (e.g., Adam) perform worse than non-adaptive optimizers (e.g., SGD, momentum SGD) on out-of-distribution performance. In particular, even though there is no significant difference in in-distribution performance, we show a measurable difference in out-of-distribution performance. ii) in-distribution performance and out-of-distribution performance exhibit three types of behavior depending on the dataset -- linear returns, increasing returns, and diminishing returns. For example, in the training of natural language data using Adam, fine-tuning the performance of in-distribution performance does not significantly contribute to the out-of-distribution generalization performance. △ Less

Submitted 5 June, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: Accepted to TMLR

Towards the continuum limit of nucleon form factors at the physical point using lattice QCD

Authors: Ryutaro Tsuji, Yasumichi Aoki, Ken-Ichi Ishikawa, Yoshinobu Kuramashi, Shoichi Sasaki, Eigo Shintani, Takeshi Yamazaki

Abstract: We present results for the axial charge and root-mean-square (RMS) radii of the nucleon obtained from 2+1 flavor lattice QCD at the physical point with a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared $O(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $β$ = 1… ▽ More We present results for the axial charge and root-mean-square (RMS) radii of the nucleon obtained from 2+1 flavor lattice QCD at the physical point with a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared $O(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $β$ = 1.82 and 2.00 corresponding to lattice spacings of 0.085 fm and 0.063 fm respectively. We first evaluate the value of $g_A/g_V$ , which is not renormalized in the continuum limit and thus ends up with the renormalized axial charge. Moreover, we also calculate the nucleon elastic form factors and determine three kinds of isovector RMS radii such as electric, magnetic and axial ones at the two lattice spacings. We finally discuss the discretization uncertainties on renormalized axial charge and isovector RMS radii towards the continuum limit. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 9 pages, 2 figures, Proceeding for the 39th International Symposium on Lattice Field Theory, 8th-13th August 2022, Bonn, Germany

Modeling of geocoronal solar wind charge exchange events detected with Suzaku

Authors: Daiki Ishi, Kumi Ishikawa, Yoshizumi Miyoshi, Naoki Terada, Yuichiro Ezoe

Abstract: A model of geocoronal solar wind charge exchange (SWCX) emission was built and compared to five Suzaku detections of bright geocoronal SWCX events. An exospheric neutral hydrogen distribution model, charge exchange cross sections, solar wind ion data taken with the ACE and WIND satellites, and magnetic field models of the Earth's magnetosphere are all combined in order to predict time-variable geo… ▽ More A model of geocoronal solar wind charge exchange (SWCX) emission was built and compared to five Suzaku detections of bright geocoronal SWCX events. An exospheric neutral hydrogen distribution model, charge exchange cross sections, solar wind ion data taken with the ACE and WIND satellites, and magnetic field models of the Earth's magnetosphere are all combined in order to predict time-variable geocoronal SWCX emission depending on line-of-sight directions of the Suzaku satellite. The modeled average intensities of O VII emission lines were consistent with the observed ones within a factor of three in four out of the five cases except for an event in which a line-of-sight direction was toward the night side of the high-latitude magnetosheath and a major geomagnetic storm was observed. Those of O VIII emission lines were underestimated by a factor of three or more in all the five cases. On the other hand, the modeled O VII and O VIII light curves reproduced the observed ones after being scaled by ratios between the observed and modeled average intensities. In particular, short-term variations due to line-of-sight directions traversing cusp regions during an orbital motion of the Suzaku satellite were reproduced. These results are discussed in the context of model uncertainties. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: 47 pages, 23 figures, accepted for publication in PASJ

Time-dependent optimized coupled-cluster method with doubles and perturbative triples [TD-OCCD(T)] for first principles simulation of multielectron dynamics

Authors: Himadri Pathak, Takeshi Sato, Kenichi L. Ishikawa

Abstract: We report the formulation of a new, cost-effective approximation method in the time-dependent optimized coupled-cluster (TD-OCC) framework [T. Sato et al., J. Chem. Phys. 148, 051101 (2018)] for first-principles simulations of multielectron dynamics in an intense laser field. The method, designated as TD-OCCD(T), is a time-dependent, orbital-optimized extension of the "gold-standard" CCSD(T) metho… ▽ More We report the formulation of a new, cost-effective approximation method in the time-dependent optimized coupled-cluster (TD-OCC) framework [T. Sato et al., J. Chem. Phys. 148, 051101 (2018)] for first-principles simulations of multielectron dynamics in an intense laser field. The method, designated as TD-OCCD(T), is a time-dependent, orbital-optimized extension of the "gold-standard" CCSD(T) method in the ground-state electronic structure theory. The equations of motion for the orbital functions and the coupled-cluster amplitudes are derived based on the real-valued time-dependent variational principle using the fourth-order Lagrangian. The TD-OCCD(T) is size extensive and gauge invariant, and scales as O(N^7) with respect to the number of active orbitals N. The pilot application of the TD-OCCD(T) method to the strong-field ionization and high-order harmonic generation from a Kr atom is reported in comparison with the results of the previously developed methods, such as the time-dependent complete-active-space self-consistent field (TD-CASSCF), TD-OCC with double and triple excitations (TD-OCCDT), TD-OCC with double excitations (TD-OCCD), and the time-dependent Hartree-Fock (TDHF) methods. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Energy polarization and energy pumping in Rice-Mele chains

Authors: Kiminori Hattori, Kodai Ishikawa, Yuma Kaneko

Abstract: The one-dimensional Rice-Mele lattice consisting of two-site unit cells enables topological charge pumping at low enough temperatures. We show here that energy pumping is feasible in this system. We analyze energy polarization using the second-quantization formalism and connect it to intercell energy current by the relevant continuity equation. Following this formulation, we numerically evaluate t… ▽ More The one-dimensional Rice-Mele lattice consisting of two-site unit cells enables topological charge pumping at low enough temperatures. We show here that energy pumping is feasible in this system. We analyze energy polarization using the second-quantization formalism and connect it to intercell energy current by the relevant continuity equation. Following this formulation, we numerically evaluate the energy current pumped in an adiabatic cycle when the system parameters vary periodically and slowly with time. Unlike the pumped charge, the pumped energy is entirely temperature-independent at half filling. The global Berry phase for all bands and the related topological phase transitions during the pumping cycle account for the energy pumping. The present study paves the way for manipulating energy current without electric and thermal biases. △ Less

Submitted 10 February, 2023; v1 submitted 10 August, 2022; originally announced August 2022.

Comments: 7 pages, 5 figures

Journal ref: Phys. Rev. B 107, 115401 (2023)

Perturbative study of large $N$ principal chiral model with twisted reduction

Authors: Antonio González-Arroyo, Ken-Ichi Ishikawa, Yingbo Ji, Masanori Okawa

Abstract: We compute the first four perturbative coefficients of the internal energy for the twisted reduced principal chiral model (TRPCM) using numerical stochastic perturbation theory (NSPT). This matrix model has the same large $N$ limit as the ordinary principal chiral model (PCM) at infinite volume. Indeed, we verify that the first three coefficients match the analytic result for the PCM coefficients… ▽ More We compute the first four perturbative coefficients of the internal energy for the twisted reduced principal chiral model (TRPCM) using numerical stochastic perturbation theory (NSPT). This matrix model has the same large $N$ limit as the ordinary principal chiral model (PCM) at infinite volume. Indeed, we verify that the first three coefficients match the analytic result for the PCM coefficients at large $N$ with a precision of three to four significant digits. The fourth coefficient also matches our own NSPT calculation of the corresponding PCM coefficient at large $N$. The finite-$N$ corrections to all coefficients beyond the leading order are smaller for TRPCM than for PCM. We analyze the variance to determine the feasibility of extending the calculations to higher orders. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 27 pages, 12 figures, using ws-ijmpa.cls style, using siunitx version 3

Report number: HUPD-2208, IFT-UAM/CSIC-22-78

Nucleon isovector couplings in Nf = 2 + 1 lattice QCD at the physical point

Authors: Ryutaro Tsuji, Natsuki Tsukamoto, Yasumichi Aoki, Ken-Ichi Ishikawa, Yoshinobu Kuramashi, Shoichi Sasaki, Eigo Shintani, Takeshi Yamazaki

Abstract: We present results for the scalar and tensor isovector-couplings ($g_S$ and $g_T$) of the nucleon measured at the physical point ($M_π=135$ MeV) with a single lattice spacing of $0.085\ \mathrm{fm}$ in 2+1 flavor QCD. Our calculations are carried out with two ensembles of gauge configurations generated by the PACS Collaboration with nonperturbatively ${\cal O}(a)$ improved Wilson quark action and… ▽ More We present results for the scalar and tensor isovector-couplings ($g_S$ and $g_T$) of the nucleon measured at the physical point ($M_π=135$ MeV) with a single lattice spacing of $0.085\ \mathrm{fm}$ in 2+1 flavor QCD. Our calculations are carried out with two ensembles of gauge configurations generated by the PACS Collaboration with nonperturbatively ${\cal O}(a)$ improved Wilson quark action and Iwasaki gauge action on $(10.9\ {\rm fm})^4$ and $(5.5\ {\rm fm})^4$ lattices, where the finite-size effect on the nucleon mass was not shown at the level of the statistical precision less than 0.5%. We compute the nucleon three-point correlation functions in the vector, axial, scalar, and tensor channels. We confirm that our previous result of the nucleon axial coupling on the large spatial volume of $(10.9\ {\rm fm})^4$ has no finite-size effect at the level of the statistical precision of 1.9%. For the renormalization, we first renormalize $g_S$ and $g_T$ nonperturbatively using the RI/SMOM$_{(γ_μ)}$ scheme, a variant of Rome-Southampton RI/MOM scheme with reduced systematic errors, as the intermediate scheme. We evaluate our final results at the renormalization scale of 2 GeV in the $\overline{\rm MS}$ scheme through matching procedure between the RI/SMOM$_{(γ_μ)}$ and $\overline{\rm MS}$ schemes with the help of perturbation theory, and then obtain $g_S=0.927(71)_{\rm stat}(22)_{\rm syst}$ and $g_T=1.036(6)_{\rm stat}(20)_{\rm syst}$. △ Less

Submitted 20 October, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 52 pages, 11 figures

Report number: UTHEP-771, UTCCS-P-145, HUPD-2209

$K_{\ell 3}$ form factors at the physical point: Toward the continuum limit

Authors: Ken-ichi Ishikawa, Naruhito Ishizuka, Yoshinobu Kuramashi, Yusuke Namekawa, Yusuke Taniguchi, Naoya Ukita, Takeshi Yamazaki, Tomoteru Yoshié for PACS Collaboration

Abstract: We present updated results for the form factors of the kaon semileptonic $(K_{\ell 3})$ decay process calculated with $N_f = 2 + 1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action at the physical point on large volumes of more than (10 fm)$^4$. In addition to our previous calculation at the lattice spacing $a = 0.085$ fm, we perform a calculation at the second lattic… ▽ More We present updated results for the form factors of the kaon semileptonic $(K_{\ell 3})$ decay process calculated with $N_f = 2 + 1$ nonperturbatively $O(a)$-improved Wilson quark action and Iwasaki gauge action at the physical point on large volumes of more than (10 fm)$^4$. In addition to our previous calculation at the lattice spacing $a = 0.085$ fm, we perform a calculation at the second lattice spacing of $0.063$ fm. Using the results for the form factors extracted from 3-point functions with the local and also conserved vector currents at the two lattice spacings, continuum extrapolation and interpolation of the momentum transfer are carried out simultaneously to obtain the value of the form factor $f_+(0)$ at the zero momentum transfer in the continuum limit. After investigation of stability of $f_+(0)$ against several fit forms and different data, we obtain $f_+(0) = 0.9615(10)(^{+47}_{\ -3})(5)$, where the first, second, and third errors are statistical, systematic errors from choice of the fit forms and isospin breaking effect, respectively. Combining our value of $f_+(0)$ and experimental input of the $K_{\ell 3}$ decay, one of the Cabibbo-Kobayashi-Maskawa matrix elements $|V_{us}|$ is determined as $|V_{us}| = 0.2252(^{\ +5}_{-12})$, whose error contains the experimental one as well as that in the lattice calculation. This value is reasonably consistent with the ones determined from recent lattice QCD results of $f_+(0)$ and also the one determined through the kaon leptonic decay process. We observe some tension between our value and $|V_{us}|$ evaluated from the unitarity of the CKM matrix with $|V_{ud}|$, while it depends on the size of the error of $|V_{ud}|$. It is also found that $|V_{us}|$ determined with our phase space integrals through six $K_{\ell 3}$ decay processes is consistent with the above one using $f_+(0)$. △ Less

Submitted 21 November, 2022; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: 54 pages, 23 figures, published version

Report number: UTHEP-770, UTCCS-P-144, HUPD-2207, KUNS-2930

Journal ref: Phys.Rev.D 106 (2022) 9, 094501

Feature Space Particle Inference for Neural Network Ensembles

Authors: Shingo Yashima, Teppei Suzuki, Kohta Ishikawa, Ikuro Sato, Rei Kawakami

Abstract: Ensembles of deep neural networks demonstrate improved performance over single models. For enhancing the diversity of ensemble members while keeping their performance, particle-based inference methods offer a promising approach from a Bayesian perspective. However, the best way to apply these methods to neural networks is still unclear: seeking samples from the weight-space posterior suffers from… ▽ More Ensembles of deep neural networks demonstrate improved performance over single models. For enhancing the diversity of ensemble members while keeping their performance, particle-based inference methods offer a promising approach from a Bayesian perspective. However, the best way to apply these methods to neural networks is still unclear: seeking samples from the weight-space posterior suffers from inefficiency due to the over-parameterization issues, while seeking samples directly from the function-space posterior often results in serious underfitting. In this study, we propose optimizing particles in the feature space where the activation of a specific intermediate layer lies to address the above-mentioned difficulties. Our method encourages each member to capture distinct features, which is expected to improve ensemble prediction robustness. Extensive evaluation on real-world datasets shows that our model significantly outperforms the gold-standard Deep Ensembles on various metrics, including accuracy, calibration, and robustness. Code is available at https://github.com/DensoITLab/featurePI . △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: ICML2022

Enhanced energy deposition and carrier generation in silicon induced by two-color intense femtosecond laser pulses

Authors: Mizuki Tani, Kakeru Sasaki, Yasushi Shinohara, Kenichi L. Ishikawa

Abstract: We theoretically investigate the optical energy absorption of crystalline silicon subject to dual-color femtosecond laser pulses, using the time-dependent density functional theory (TDDFT). We employ the modified Becke-Johnson (mBJ) exchange-correlation potential which reproduces the experimental direct bandgap energy $E_g$. We consider situations where the one color is in the ultraviolet (UV) ran… ▽ More We theoretically investigate the optical energy absorption of crystalline silicon subject to dual-color femtosecond laser pulses, using the time-dependent density functional theory (TDDFT). We employ the modified Becke-Johnson (mBJ) exchange-correlation potential which reproduces the experimental direct bandgap energy $E_g$. We consider situations where the one color is in the ultraviolet (UV) range above $E_g$ and the other in the infrared (IR) range below it. The energy deposition is examined as a function of mixing ratio $η$ of the two colors with the total pulse energy conserved. Energy transfer from the laser pulse to the electronic system in silicon is dramatically enhanced by simultaneous dual-color irradiation and maximized at $η\sim 0.5$. Increased is the number of generated carriers, not the absorbed energy per carrier. The effect is more efficient for lower IR photon energy, or equivalently, larger vector-potential amplitude. As the underlying mechanism is identified the interplay between intraband electron motion in the {\it valence} band (before excitation) driven by the IR component and resonant valence-to-conduction interband excitation (carrier injection) induced by the UV component. The former increases excitable electrons which pass through the $k$ points of resonant transitions. The effect of different multiphoton absorption paths or intraband motion of carriers generated in the conduction band play a minor role. △ Less

Submitted 1 November, 2022; v1 submitted 31 May, 2022; originally announced May 2022.

Comments: 14 pages, 20 figures

Scale setting for large-$N$ SUSY Yang-Mills on the lattice

Authors: Pietro Butti, Margarita Garcia Perez, Antonio Gonzalez-Arroyo, Ken-Ichi Ishikawa, Masanori Okawa

Abstract: In this paper we study the large $N$ limit of four-dimensional Supersymmetric Yang-Mills on the lattice using twisted reduced models. We have generated configurations with dynamical massive gluinos and various lattice 't Hooft couplings, and verified that the Pfaffian remains positive. We have determined the lattice spacing in terms of various observables obtaining compatible results. Extrapolatin… ▽ More In this paper we study the large $N$ limit of four-dimensional Supersymmetric Yang-Mills on the lattice using twisted reduced models. We have generated configurations with dynamical massive gluinos and various lattice 't Hooft couplings, and verified that the Pfaffian remains positive. We have determined the lattice spacing in terms of various observables obtaining compatible results. Extrapolating to the massless gluino limit we obtain the lattice spacing dependence on the bare couplings for the supersymmetric theory. The observed dependence goes along the expected behaviour predicted by perturbation theory. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: 55 pages

Report number: IFT-UAM/CSIC-22-39, HUPD-2206

Cr doping-induced ferromagnetism in the spin-glass Cd1-xMnxTe studied by x-ray magnetic circular dichroism

Authors: V. K. Verma, S. Sakamoto, K. Ishikawa, V. R. Singh, K. Ishigami, G. Shibata, T. Kadono, T. Koide, S. Kuroda, A. Fujimori

Abstract: The prototypical diluted magnetic semiconductor Cd1-xMnxTe is a spin glass (x<0.6) or an antiferromagnet (x>0.6), but becomes ferromagnetic upon doping with a small amount of Cr atoms substituting for Mn. In order to investigate the origin of the ferromagnetism in Cd1-x-yMnxCryTe, we have studied its element specific magnetic properties by x-ray absorption spectroscopy (XAS) and x-ray magnetic cir… ▽ More The prototypical diluted magnetic semiconductor Cd1-xMnxTe is a spin glass (x<0.6) or an antiferromagnet (x>0.6), but becomes ferromagnetic upon doping with a small amount of Cr atoms substituting for Mn. In order to investigate the origin of the ferromagnetism in Cd1-x-yMnxCryTe, we have studied its element specific magnetic properties by x-ray absorption spectroscopy (XAS) and x-ray magnetic circular dichroism (XMCD) at the Cr and Mn L2,3 edges. Thin films were grown by molecular beam epitaxy with a fixed Mn content of x = 0.2 and varying Cr content in the range of y = 0 - 0.04. Measured XAS and XMCD spectra indicate that both Cr and Mn atoms are divalent and that the ferromagnetic or superparamagnetic components of Cr and Mn are aligned in the same directions. The magnetization of Mn increases with increasing Cr content. These results can be explained if ferromagnetic interaction exists between neighboring Mn and Cr ions although interaction between Mn atoms is largely antiferromagnetic. We conclude that each ferromagnetic or superparamagnetic cluster consists of ferromagnetically coupled several Cr and a much larger number of Mn ions. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 13 pages, 5 figures

Nucleon isovector tensor charge from lattice QCD with physical light quarks

Authors: Ryutaro Tsuji, Yasumichi Aoki, Ken-Ichi Ishikawa, Yoshinobu Kuramashi, Shoichi Sasaki, Eigo Shintani, Takeshi Yamazaki

Abstract: We present preliminary results for the axial, scalar and tensor charges of the nucleon measured in 2+1 flavor QCD with the physical light quarks ($m_π=135$ MeV). Our simulations are carried out with gauge configurations generated by the PACS Collaboration with the stout-smeared $O(a)$ improved Wilson fermions and Iwasaki gauge action at a single lattice spacing of $0.085\ (\mathrm{fm})$. There are… ▽ More We present preliminary results for the axial, scalar and tensor charges of the nucleon measured in 2+1 flavor QCD with the physical light quarks ($m_π=135$ MeV). Our simulations are carried out with gauge configurations generated by the PACS Collaboration with the stout-smeared $O(a)$ improved Wilson fermions and Iwasaki gauge action at a single lattice spacing of $0.085\ (\mathrm{fm})$. There are two lattice ensembles of the PACS gauge configurations, which have physical lattice sizes over $(10\ \mathrm{fm})^4$ and $(5\ \mathrm{fm})^4$, respectively. We compute the nucleon three-point correlation functions in the axial, scalar, and tensor channels. For the renormalization, we use the Rome-Southampton method as the intermediate scheme in order to evaluate the renormalization constants for the scalar and tensor currents in fully nonperturbative manner. We then evaluate the renormalized values of the scalar and tensor charges ($g_S$ and $g_T$) in the $\overline{\rm MS}$ scheme at the renormalization scale of 2 GeV with a help of the continuum perturbation theory for the matching between two schemes. We compare our preliminary results of $g_S$ and $g_T$ with those of other collaboration results. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 5 pages, 5 figures, Proceedings of the 24th International Spin Symposium, October 18-22, 2021, Matsue, Japan

Report number: HUPD-2202

The lower moments of nucleon structure functions in lattice QCD with physical quark masses

Authors: Ryutaro Tsuji, Yasumichi Aoki, Ken-Ichi Ishikawa, Yoshinobu Kuramashi, Shoichi Sasaki, Eigo Shintani, Takeshi Yamazaki

Abstract: We present results for the nucleon structure functions and form factors obtained from 2+1 flavor lattice QCD with physical light quark masses ($m_π=135$ MeV) in a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared ${\mathscr{O}}(a)$ improved Wilson-clover quark action and Iwasaki ga… ▽ More We present results for the nucleon structure functions and form factors obtained from 2+1 flavor lattice QCD with physical light quark masses ($m_π=135$ MeV) in a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared ${\mathscr{O}}(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $β=1.82$ and $2.00$ corresponding to lattice spacings of $0.085$ fm and $0.064$ fm respectively. The lower moments of structure functions, $\langle x \rangle_{u-d}$ and $\langle x \rangle_{Δu - Δd}$ given by the twist-2 operators being properly renormalized, are evaluated in the $\overline{\rm MS}$ scheme at the renormalization scale of 2 GeV only at $β=1.82$, since the renormalization factors at $β=2.00$ have not yet determined nonperturbatively in the RI/MOM scheme. Instead, at two lattice spacings, we evaluate appropriate ratios of $g_{A}/g_{V}$ and $\langle x \rangle_{u-d}/\langle x \rangle_{Δu -Δd}$, which are not renormalized in the continuum limit. These quantities thus can be directly compared with the experimental data without the renormalization. △ Less

Submitted 30 December, 2021; originally announced December 2021.

Comments: 9 pages, 4 figures, , talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Report number: HUPD-2112

Multigrid Solver on Fugaku

Authors: Ken-Ichi Ishikawa, Issaku Kanamori, Hideo Matsufuru

Abstract: We report an implementation of a multigrid solver for the Clover fermion on supercomputer Fugaku, which uses A64FX CPU with Arm architecture. On Fugaku, a highly optimized implementation of BiCGStab solver with domain decomposed preconditioner for the Clover fermion, called QCD Wide SIMD library (QWS), is available. We use the preconditioner in QWS as a smoother of the multigrid solver. Our implem… ▽ More We report an implementation of a multigrid solver for the Clover fermion on supercomputer Fugaku, which uses A64FX CPU with Arm architecture. On Fugaku, a highly optimized implementation of BiCGStab solver with domain decomposed preconditioner for the Clover fermion, called QCD Wide SIMD library (QWS), is available. We use the preconditioner in QWS as a smoother of the multigrid solver. Our implementation shows reasonably good scaling and performance. The code is developed by using Bridge++ code framework and its extension. △ Less

Submitted 1 December, 2021; originally announced December 2021.

Comments: 8 pages, 3 figures, talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

Report number: HUPD-2114

Atomic real-space perspective of light-field-driven currents in graphene

Authors: Yuya Morimoto, Yasushi Shinohara, Kenichi L. Ishikawa, Peter Hommelhoff

Abstract: When graphene is exposed to a strong few-cycle optical field, a directional electric current can be induced depending on the carrier-envelope phase of the field. This phenomenon has successfully been explained by the charge dynamics in reciprocal space, namely an asymmetry in the conduction band population left after the laser excitation. However, the corresponding real-space perspective has not b… ▽ More When graphene is exposed to a strong few-cycle optical field, a directional electric current can be induced depending on the carrier-envelope phase of the field. This phenomenon has successfully been explained by the charge dynamics in reciprocal space, namely an asymmetry in the conduction band population left after the laser excitation. However, the corresponding real-space perspective has not been explored so far although it could yield knowledge about the atomic origin of the macroscopic currents. In this work, by adapting the nearest-neighbor tight-binding model including overlap integrals and the semiconductor Bloch equation, we reveal the spatial distributions of the light-field-driven currents on the atomic scale and show how they are related to the light-induced changes of charge densities. The atomic-scale currents flow dominantly through the network of the π bonds and are the strongest at the bonds parallel to the field polarization, where an increase of the charge density is observed. The real-space maps of the currents and changes in charge densities are elucidated using simple symmetries connecting real and reciprocal space. We also discuss the strong-field-driven Rabi oscillations appearing in the atomic-scale charge densities. This work highlights the importance of real-space measurements and stimulates future time-resolved atomic-scale experimental studies with high-energy electrons or X-rays, for examples. △ Less

Submitted 1 March, 2022; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 23 pages, 5 figures

Object-oriented implementation of algebraic multi-grid solver for lattice QCD on SIMD architectures and GPU clusters

Authors: Issaku Kanamori, Ken-Ichi Ishikawa, Hideo Matsufuru

Abstract: A portable implementation of elaborated algorithm is important to use variety of architectures in HPC applications. In this work we implement and benchmark an algebraic multi-grid solver for Lattice QCD on three different architectures, Intel Xeon Phi, Fujitsu A64FX, and NVIDIA Tesla V100, in keeping high performance and portability of the code based on the object-oriented paradigm. Some parts of… ▽ More A portable implementation of elaborated algorithm is important to use variety of architectures in HPC applications. In this work we implement and benchmark an algebraic multi-grid solver for Lattice QCD on three different architectures, Intel Xeon Phi, Fujitsu A64FX, and NVIDIA Tesla V100, in keeping high performance and portability of the code based on the object-oriented paradigm. Some parts of code are specific to an architecture employing appropriate data layout and tuned matrix-vector multiplication kernels, while the implementation of abstract solver algorithm is common to all architectures. Although the performance of the solver depends on tuning of the architecture-dependent part, we observe reasonable scaling behavior and better performance than the mixed precision BiCGSstab solvers. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: 17 pages, 7 figures. The 21st International Conference on Computational Science and Its Applications, ICCSA 2021

Report number: HUPD-2111

Journal ref: LNCS, volume 12953, pp. 1-16, 2021

Scale setting for $\mathcal{N} = 1$ SUSY Yang-Mills at large-$N$ through volume-reduced twisted matrix model

Authors: P. Butti, M. Garcia Perez, A. Gonzalez-Arroyo, K. I. Ishikawa, M. Okawa

Abstract: $\mathcal{N}=1$ SUSY Yang-Mills theory is an appealing theoretical framework that has been studied in the literature using different methods, including standard lattice simulations. Among these, the volume-reduced twisted Eguchi-Kawai model, endowed with one adjoint Majorana fermion, could play an important role in studying its large-$N$ limit via the Curci-Veneziano prescription. In this talk, we… ▽ More $\mathcal{N}=1$ SUSY Yang-Mills theory is an appealing theoretical framework that has been studied in the literature using different methods, including standard lattice simulations. Among these, the volume-reduced twisted Eguchi-Kawai model, endowed with one adjoint Majorana fermion, could play an important role in studying its large-$N$ limit via the Curci-Veneziano prescription. In this talk, we present our results on the analysis of the scale of the theory, performed via different methods based on purely gluonic observables as well as (quenched) fundamental mesons in the chiral limit. These lattice results will be used as a scale setting for the analysis of the spectrum of the theory. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 9 pages, 4 figures, 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology

Report number: FTUAM-21-4, HUPD-2109, IFT-UAM/CSIC-21-120