Search

With the evolution of numerical methods, we are now aiming at not only qualitative understanding but also quantitative prediction and design of quantum many-body phenomena. As a novel numerical approach, machine learning techniques have been introduced in 2017 to analyze quantum many-body problems. Since then, proposed various novel approaches have opened a new era, in which challenging and fundamental problems in physics can be solved by machine learning methods. Especially, quantitative and accurate estimates of material-dependent physical properties of strongly correlated matter have now become realized by combining first-principles calculations with highly accurate quantum many-body solvers developed with the help of machine learning methods. Thus developed quantitative description of electron correlations will constitute a key element of materials science in the next generation., Comment: Review paper, 12 pages, 12 figures

Solar flares exhibit complex variations in elemental abundances compared to photospheric values. We examine the spatial and temporal evolution of coronal abundances in the X8.2 flare on 2017 September 10, aiming to interpret the often observed high first ionization potential (FIP) bias at loop tops and provide insights into differences between spatially resolved and Sun-as-a-star flare composition measurements. We analyze 12 Hinode/EIS raster scans spanning 3.5 hours, employing Ca XIV 193.87 A/Ar XIV 194.40 A and Fe XVI 262.98 A/S XIII 256.69 A composition diagnostics to derive FIP bias values. Both diagnostics consistently show that flare loop tops maintain high FIP bias values of >2-6, with peak phase values exceeding 4, over the extended duration, while footpoints exhibit photospheric FIP bias of ~1. We propose that this variation arises from a combination of two distinct processes: high FIP bias plasma downflows from the plasma sheet confined to loop tops, and chromospheric evaporation filling the loop footpoints with low FIP bias plasma. Mixing between these two sources produces the observed gradient. Our observations show that the localized high FIP bias signature at loop tops is likely diluted by the bright footpoint emission in spatially averaged measurements. The spatially resolved spectroscopic observations enabled by EIS prove critical for revealing this complex abundance variation in loops. Furthermore, our observations show clear evidence that the origin of hot flare plasma in flaring loops consists of a combination of both directly heated plasma in the corona and from ablated chromospheric material; and our results provide valuable insights into the formation and composition of loop top brightenings, also known as EUV knots, which are a common feature at the tops of flare loops., Comment: 13 pages, 7 figures, 1 table. Accepted in A&A. Comments and criticisms are welcomed!

A plasma scenario with negative triangularity (NT) shaping is achieved on MAST-U for the first time. While edge localized modes (ELMs) are eventually suppressed as the triangularity is decreased below $\delta$ < -0.06, an extended period of H-mode operation with Type-III ELMs is sustained at less negative $\delta$ even through access to the second stability region for ideal ballooning modes is closed. This documents a qualitative difference from the ELM-free access conditions documented in NT scenarios on conventional aspect ratio machines. The electron temperature at the pedestal top drops across the transition to ELM-free operation, but a steady rise in core temperature as $\delta$ is decreased allows for similar normalized beta in the ELM-free NT and H-mode positive triangularity shapes.

Stimulated by the successful descriptions of strongly correlated electron systems by fractionalized fermions, correspondence between interacting fermions and non-interacting multi-component fermions is formulated in examples of the Hubbard model. The formalism enables constructions of the neural network for a quantum many-body solver represented by coupled noninteracting fermions. After showing the exact correspondence of 1- and 2-site Hubbard model to two-component noninteracting fermions, numerical algorithm of the quantum machine learning for the Hubbard model is proposed. Benchmark for the 4-site systems is successfully presented and promising future directions as well as implications are discussed., Comment: 11 pages, 2 figures

We present a study of spectral line width measurements from the Extreme Ultraviolet Imaging Spectrometer (EIS) on {\it Hinode}. We used spectral line profiles of Fe {\sc xvi} 262.984 {\AA}, Fe {\sc xiv} 264.787 {\AA}, Fe {\sc xiv} 270.519 {\AA}, Fe {\sc xiv} 274.203 {\AA}, and Fe {\sc xv} 284.160 {\AA}, and studied 11 active regions. Previous studies of spectral line widths have shown that in hot loops in the cores of active regions, the observed non-thermal velocities are smaller than predicted from models of reconnection jets in the corona or shock heating associated with Alfv\'{e}n waves. The observed line widths are also inconsistent with models of chromospheric evaporation due to coronal nanoflares. We show that recent advances in higher resolution Alfv\'{e}n wave turbulence modeling enables us to obtain non-thermal velocities similar to those measured in active regions. The observed non-thermal velocities for the 11 active regions in our study are in the range of 17$-$30 $\rm km ~ s^{-1}$, consistent with the spectral line non-thermal widths predicted from our model of 16 interacting flux tubes, which are in the range of ~15$-$37 $\rm km ~ s^{-1}$.

We present Ir $4f$ core-level hard-x-ray photoemission spectroscopy (HAXPES) experiments conducted across antiferromagnetic (AFM) ordering transition in Ruddlesden-Popper iridates Sr$_2$IrO$_4$ and Sr$_3$Ir$_2$O$_7$. The Ir $4f$ spectra exhibit distinct changes between the AFM and paramagnetic (PM) phases, with the spectral difference $I_\text{PM}-I_\text{AFM}$ showing a contrasting behavior in the two compounds. By employing computational simulations using the local-density approximation combined with the dynamical mean-field theory method, we elucidate that $I_\text{PM}-I_\text{AFM}$ primary reflects the Slater or Mott-Hubbard character of the AFM insulating state rather than material specific details. This sensitivity to fine low-energy electronic structure arises from the dependence of charge-transfer responses to the sudden creation of a localized core hole on both metal-insulator transitions and long-range AFM ordering. Our result broadens the applications of core-level HAXPES as a tool for characterization of electronic structure., Comment: Main: 7 pages, 3 figures ; SM: 4 pages, 4 figures

Numerically computing physical quantities of time-evolved states in quantum many-body systems is a challenging task in general. Here, we propose a method that utilizes reliable short-time data of physical quantities to accurately forecast long-time behavior. The method is based on the dynamic mode decomposition (DMD), which is commonly used in fluid dynamics. The effectiveness and applicability of the DMD in quantum many-body systems such as the Ising model in the transverse field at the critical point are studied, even when the input data exhibits complicated features such as multiple oscillatory components and a power-law decay with long-ranged quantum entanglements unlike fluid dynamics. It is demonstrated that the present method enables accurate forecasts at time as long as nearly an order of magnitude longer than that of the short-time training data. Effects of noise on the accuracy of the forecast are also investigated, because they are important especially when dealing with the experimental data. We find that a few percent of noise does not affect the prediction accuracy destructively., Comment: 18 pages, 24 figures

We detected a giant X-ray flare from the RS-CVn type binary star UX Ari using MAXI on 2020 August 17 and started a series of NICER observations 89 minutes later. For a week, the entire duration of the flare was covered with 32 snapshot observations including the rising phase. The X-ray luminosity reached 2$\times$10$^{33}$ erg s$^{-1}$ and the entire energy release was $\sim 10^{38}$ erg in the 0.5--8.0~keV band. X-ray spectra characterized by continuum emission with lines of Fe XXV He$\alpha$ and Fe XXVI Ly$\alpha$ were obtained. We found that the temperature peaks before that of the flux, which suggests that the period of plasma formation in the magnetic flare loop was captured. Using the continuum information (temperature, flux, and their delay time), we estimated the flare loop size to be $\sim 3 \times 10^{11}$ cm and the peak electron density to be $\sim 4\times10^{10}$ cm$^{-3}$. Furthermore, using the line ratio of Fe XXV and Fe XXVI, we investigated any potential indications of deviation from collisional ionization equilibrium (CIE). The X-ray spectra were consistent with CIE plasma throughout the flare, but the possibility of an ionizing plasma away from CIE was not rejected in the flux rising phase., Comment: 9 pages, 6 figures. Accepted by ApJ

Non-periodic folding of periodic crease patterns paves the way to novel nonlinear phenomena that cannot be feasible through periodic folding. This paper focuses on the non-periodic folding of recursive crease patterns generalized from Spidron. Although it is known that the Spidron has a $1$-DOF isotropic rigid folding motion, its general kinematics and dependence on the crease pattern remain unclear. Using the kinematics of a single unit cell of the Spidron and the recursive construction of the folded state of multiple unit cells, we consider the folding of the Spidron that is not necessarily isotropic. We found that as the number of unit cells increases, the non-periodic folding is restricted and the isotropic folding becomes dominant. Then, we analyze the three kinds of isotropic folding modes by constructing $1$-dimensional dynamical systems governing each of them. We show that these systems can possess different recursive natures depending on folding modes even in an identical crease pattern. Furthermore, we show their novel nonlinear nature, including the period-doubling cascade leading to the emergence of chaos., Comment: 16 pages, 8 figures. To Appear in Origami8: Proceedings of the 8th International Meeting on Origami in Science, Mathematics and Education (8OSME), Melbourne, August 16-18, 2024

The superconducting (SC) cuprate HgBa$_2$Ca$_2$Cu$_3$O$_8$ (Hg1223) has the highest $T_{c}^{\rm opt}\simeq 138$ K (the experimental SC transition temperature at optimal hole doping) among cuprates at ambient pressure $P_{\rm amb}$. $T_{c}^{\rm opt}$ increases under pressure $P$ and reaches $164$ K at $P_{\rm opt}\simeq 30$ GPa, then decreases with increasing $P>P_{\rm opt}$. To understand the microscopic origin of this dome-like $P$ dependence of $T_{c}^{\rm opt}$, we consider the $ab$ $initio$ low-energy effective Hamiltonian (LEH) for the antibonding (AB) Cu$3d_{x^2-y^2}$/O$2p_{\sigma}$ band. In the AB LEH for cuprates with $N_\ell \leq 2$ laminated CuO$_2$ planes between block layers, it was proposed that $T_{c}^{\rm opt}\simeq 0.16|t_1|F_{\rm SC}$, where $t_1$ is the nearest neighbor hopping and the SC order parameter $F_{\rm SC}$ at optimal hole doping mainly depends on $u=U/|t_1|$ ($U$ is the onsite effective Coulomb repulsion): $F_{\rm SC}$ is maximal at $u_{\rm opt}\simeq 8.5$ and decreases sharply with decreasing $u$ for $uP_{\rm opt}$, the decrease in $T_{c}^{\rm opt}$ is accounted for by (II): The rapid decrease in $F_{\rm SC}$ with decreasing $u$ for $u

We present a study of the impact of an uncertainty in the beam far side-lobe knowledge on the measurement of the Cosmic Microwave Background $B$-mode signal at large scale. It is expected to be one of the main source of systematic effects in future CMB observations. Because it is crucial for all-sky survey missions to take into account the interplays between beam systematic effects and all the data analysis steps, the primary goal of this paper is to provide the methodology to carry out the end-to-end study of their effect for a space-borne CMB polarization experiment, up to the cosmological results in the form of a bias $\delta r$ on the tensor-to-scalar ratio $r$. LiteBIRD is dedicated to target the measurement of CMB primordial $B$ modes by reaching a sensitivity of $\sigma \left( r \right) \leq 10^{-3}$ assuming $r=0$. As a demonstration of our framework, we derive the relationship between the knowledge of the beam far side-lobes and the tentatively allocated error budget under given assumptions on design, simulation and component separation method. We assume no mitigation of the far side-lobes effect at any stage of the analysis pipeline. We show that $\delta r$ is mostly due to the integrated fractional power difference between the estimated beams and the true beams in the far side-lobes region, with little dependence on the actual shape of the beams, for low enough $\delta r$. Under our set of assumptions, in particular considering the specific foreground cleaning method we used, we find that the integrated fractional power in the far side-lobes should be known at a level as tight as $\sim 10^{-4}$, to achieve the required limit on the bias $\delta r < 1.9 \times 10^{-5}$. The framework and tools developed for this study can be easily adapted to provide requirements under different design, data analysis frameworks and for other future space-borne experiments beyond LiteBIRD.

The quantum antiferromagnet on the pyrochlore lattice offers an archetypal frustrated system, which potentially realizes a quantum spin liquid characterized by the absence of standard spontaneous symmetry breaking even at zero temperature, unusually as an isotropic 3D system. Despite tremendous progress in the literature, however, the nature of the ground state of the fully quantum-mechanical spin Hamiltonian on the pyrochlore lattice still remains elusive. Here, we show that an unconventional type of quantum spin liquid is born out from the pyrochlore system after the self-organized dimensional reduction leading to confined states in 2D layers. This conclusion is obtained from state-of-the-art variational Monte Carlo (VMC) simulations at zero temperature. Quantum spin liquids triggered by the emergent dimensional reduction is an unexplored route of the spin-liquid formation. The dimensional reduction from 3D to 2D is a consequence of a conventional spontaneous symmetry breaking, while the resultant decoupling of layers enables the emergence of a 2D quantum spin liquid that is adiabatically disconnected from trivial product states and exhibits strong quantum entanglement. The stabilized quantum spin liquid exhibits an algebraic decay of correlations and vanishing excitation gap in the thermodynamic limit. The wave-function structure supports the fractionalization of the spin into spinons. This spin-liquid ground state persists in the presence of spin-orbit interactions, which expands the possibilities of realizing quantum spin liquids in real pyrochlore-structured materials., Comment: 34 pages, 24 figures, 6 Tables

We investigate the thermal expansion and superconducting properties of a CuAl2-type (tetragonal) superconductor CoZr2 under high pressures. We perform high-pressure synchrotron X-ray diffraction in a pressure range of 2.9 GPa < P < 10.4 GPa and discover that CoZr2 exhibits volumetric negative thermal expansion under high pressures. Although the uniaxial positive thermal expansion (PTE) along the a-axis is observed under ambient pressure, that is suppressed by pressure, while the large uniaxial negative thermal expansion (NTE) along the c-axis is maintained under the pressure regime. As a result of a combination of the suppressed uniaxial PTE along the a-axis and uniaxial NTE along the c-axis, volumetric negative thermal expansion is achieved under high pressure in CoZr2. The mechanisms of volumetric NTE would be based on the flexible crystal structure caused by the soft Co-Co bond as seen in the iso-structural compound FeZr2, which exhibits uniaxial NTE along the c-axis. We also perform high-pressure electrical resistance measurements of CoZr2 to confirm the presence of superconductivity under the examined pressure regime in the range of 0.03 GPa < P < 41.9 GPa. We confirm the presence of superconductivity under all pressures and observe dome-like shape pressure dependence of superconducting transition temperature. Because of the coexistence of two phenomena, which are volumetric NTE and superconductivity, in CoZr2 under high pressure, the coexistence would be achievable under ambient pressure by tuning chemical compositions after our present observation., Comment: 22 pages, 7 figures, supporting information

In modern radio astronomy, one of the key technologies is to widen the frequency coverage of an instrument. The effects of higher-order modes on an instrument associated with wider bandwidths have been reported, which may degrade observation precision. It is important to quantify the higher-order propagation modes, though their power is too small to measure directly. Instead of the direct measurement of modes, we make an attempt to deduce them based on measurable radiation patterns. Assuming a linear system, whose radiated field is determined as a superposition of the mode coefficients in an instrument, we obtain a coefficient matrix connecting the modes and the radiated field and calculate the pseudo-inverse matrix. To investigate the accuracy of the proposed method, we demonstrate two cases with numerical simulations, axially-corrugated horn case and offset Cassegrain antenna case, and the effect of random errors on the precision. Both cases showed the deduced mode coefficients with a precision of 10e-6 with respect to the maximum mode amplitude and 10e-3 degrees in phase, respectively. The calculation errors were observed when the random errors were smaller than 0.01 percent of the maximum radiated field amplitude, which was a much lower level compared with measurement precision. The demonstrated method works independently of the details of a system. The method can quantify the propagation modes inside an instrument and will be applied to most of linear components and antennas, which leads to various applications such as diagnosis of feed alignment and higher-performance feed design., Comment: 9 pages, 13 figures, IEEE Transactions on Antennas and Propagation, accepted for publication

We evaluate the effect of half-wave plate (HWP) imperfections inducing intensity leakage to the measurement of Cosmic Microwave Background (CMB) $B$-mode polarization signal with future satellite missions focusing on the tensor-to-scalar ratio $r$. The HWP is modeled with the Mueller formalism, and coefficients are decomposed for any incident angle into harmonics of the HWP rotation frequency due to azimuthal angle dependence. Although we use a general formalism, band-averaged matrix coefficients are calculated as an example for a 9-layer sapphire HWP using EM propagation simulations. We perform simulations of multi-detector observations in a band centered at 140\,GHz using \LB instrumental configuration. We show both theoretically and with the simulations that most of the artefacts on Stokes parameter maps are produced by the dipole leakage on $B$-modes induced by the fourth harmonics $M^{(4f)}_{QI}$ and $M^{(4f)}_{UI}$. The resulting effect is strongly linked to the spin-2 focal plane scanning cross linking parameters. We develop a maximum likelihood-based method to correct the IP leakage by joint fitting of the Mueller matrix coefficients as well as the Stokes parameter maps. % by modifying the standard map-making procedure. We show that the residual leakage after correction leads to an additional noise limited uncertainty on $r$ of the order of $10^{-7}$, independently of the value of the Mueller matrix coefficients. We discuss the impact of the monopole signal and the potential coupling with other systematic effects such as gain variations and detector nonlinearities., Comment: 36 pages, 10 figures, submitted to JCAP

We report the ground-state symmetry of the Ce $4f$ states in the heavy fermion superconductor CeNi$_{2}$Ge$_{2}$, yielding anisotropic $cf$-hybridization between the Ce $4f$ states and conducting electrons. By analyzing linear dichroism in soft x-ray absorption and core-level hard x-ray photoemission spectra, the $4f$ symmetry is determined as $\Sigma$-type $\Gamma_{7}$, promoting predominant hybridization with the conducting electrons originating from the Ge site. The band structures probed by the soft x-ray angle-resolved photoemission indicates that the Ge $4p$ components contribute to the band renormalization through the anisotropic hybridization effects, suggesting that the control of the electronic structures of Ge orbital gives an impact to achieve the exotic phenomena in CeNi$_{2}$Ge$_{2}$.

The physical connection between thermal convection in the solar interior and the solar wind remains unclear due to their significant scale separation. Using an extended version of the three-dimensional radiative magnetohydrodynamic code RAMENS, we perform the first comprehensive simulation of the solar wind formation, starting from the wave excitation and the small-scale dynamo below the photosphere. The simulation satisfies various observational constraints as a slow solar wind emanating from the coronal hole boundary. The magnetic energy is persistently released in the simulated corona, showing a hot upward flow at the interface between open and closed fields. To evaluate the energetic contributions from Alfv\'en wave and interchange reconnection, we develop a new method to quantify the cross-field energy transport in the simulated atmosphere. The measured energy transport from closed coronal loops to open field accounts for approximately half of the total. These findings suggest a significant role of the supergranular-scale interchange reconnection in solar wind formation., Comment: 10 pages, 5 figures, accepted for publication in ApJL

In this study, we investigate the bulk electronic structure of Pb-substituted LaO$_{0.5}$F$_{0.5}$BiS$_2$ single crystals, using two types of hard-x-ray spectroscopy. High-energy-resolution fluorescence-detected x-ray absorption spectroscopy revealed a spectral change at low temperatures. Using density functional theory (DFT) simulations, we find that the temperature-induced change originates from a structural phase transition, similar to the pressure-induced transition in LaO$_{0.5}$F$_{0.5}$BiS$_2$. This finding suggests that the mechanism of bulk superconductivity induced by Pb substitution is the same as that under high pressure. Furthermore, a novel low-valence state with a mixture of divalent and trivalent Bi ions is discovered using hard x-ray photoemission spectroscopy with the aid of DFT calculations., Comment: 8 pages, 6 figures, 1 table

We numerically analyze superconductivity (SC) in the cuprate superconductors by using ab initio effective Hamiltonians consisting of the antibonding combination of Cu $3d_{x^2-y^2}$ and O $2p_{\sigma}$ orbitals. We perform variational Monte Carlo calculations for the four carrier doped cuprates with diverse experimental optimal SC critical temperature $T_{c}^{\rm opt}$: CaCuO$_2$ ($T_{c}^{\rm opt} \sim 110$ K), Bi$_2$Sr$_2$CuO$_6$ ($T_{c}^{\rm opt} \sim 10$-$40$ K), Bi$_2$Sr$_2$CaCu$_2$O$_8$ ($T_{c}^{\rm opt} \sim 85$-$100$ K), and HgBa$_2$CuO$_4$ ($T_{c}^{\rm opt} \sim 90$ K). Materials and hole doping concentration ($\delta$) dependencies of the SC order parameter $F_{\rm SC}$ and the competition with spin/charge order show essential and quantitative agreements with the available experiments in the following points: (1) The ground state is commonly the SC state, which is severely competing with the charge/spin stripe and antiferromagnetic states. (2) $F_{\rm SC}$ shows amplitude consistent with the superfluid density measured in the muon spin resonance and its dome structure found in $\delta$ dependence shows consistency with that of the SC gap in the tunneling and photoemission measurements. We further find insights into the universal SC mechanism: (I) $F_{\rm SC}$ increases with the ratio $U/|t_1|$, indicating that $U/|t_1|$ is the principal component controlling the SC. Here, $U$ and $t_1$ are the onsite Coulomb repulsion and the nearest neighbor hopping, respectively, in the Hamiltonians. (II) A universal scaling $T_{c}^{\rm opt}\sim 0.16 \lvert t_1 \rvert F_{\rm SC}$ holds. (III) SC is enhanced and optimized if $U$ is increased beyond the real available materials. It is further enhanced by decreasing the offsite interaction. The present findings provide useful clues for the design of new SC materials with even higher $T_{c}^{\rm opt}$., Comment: 27 pages, 30 figures

Metallic states near the Mott insulator show a variety of quantum phases including various magnetic, charge ordered states and high-temperature superconductivity in various transition metal oxides and organic solids. The emergence of a variety of phases and their competitions are likely intimately associated with quantum transitions between the electron-correlation driven Mott insulator and metals characterized by its criticality, and is related to many central questions of condensed matter. The quantum criticality is, however, not well understood when the transition is controlled by the bandwidth through physical parameters such as pressure. Here, we quantitatively estimate the universality class of the transition characterized by a comprehensive set of critical exponents by using a variational Monte Carlo method implemented as an open-source innovated quantum many-body solver, with the help of established scaling laws at a typical bandwidth-controlled Mott transition. The criticality indicates a weaker charge and density instability in contrast to the filling-controlled transition realized by carrier doping, implying a weaker instability to superconductivity as well. The present comprehensive clarification opens up a number of routes for quantitative experimental studies for complete understanding of elusive quantum Mott transition and nearby strange metal that cultivate future design of functionality., Comment: 18 pages, 16 figures

LiteBIRD is a planned JAXA-led CMB B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, $\delta r$, down to $\delta r<0.001$. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust SED, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compare the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the HFT (High-Frequency Telescope) frequency range is shifted logarithmically towards higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measure the tensor-to-scalar ratio $r$ uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on $r$ after foreground cleaning may be reduced by as much as 30--50 % by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to better ability to discriminate between models through higher $\chi^2$ sensitivity. (abridged), Comment: 18 pages, 13 figures. Published in A&A

The continued development of computational approaches to many-body ground-state problems in physics and chemistry calls for a consistent way to assess its overall progress. In this work, we introduce a metric of variational accuracy, the V-score, obtained from the variational energy and its variance. We provide an extensive curated dataset of variational calculations of many-body quantum systems, identifying cases where state-of-the-art numerical approaches show limited accuracy, and future algorithms or computational platforms, such as quantum computing, could provide improved accuracy. The V-score can be used as a metric to assess the progress of quantum variational methods toward a quantum advantage for ground-state problems, especially in regimes where classical verifiability is impossible., Comment: 27 pages, 6 figures