Search

We present X-ray observations of the upper atmospheric density disturbance caused by the explosive eruption of the Hunga Tonga-Hunga Ha'apai (HTHH) volcano on 15 January 2022. From 14 January to 16 January, the Chinese X-ray astronomy satellite, Insight-HXMT, was observing the supernova remnant Cassiopeia A. The X-ray data obtained during Earth's atmospheric occultations allowed us to measure neutral densities in the altitude range of ~90-150 km. The density profiles above 110 km altitude obtained before the major eruption are in reasonable agreement with expectations by both GAIA and NRLMSIS 2.0 models. In contrast, after the HTHH eruption, a severe density depletion was found up to ~1,000 km away from the epicenter, and a relatively weak depletion extending up to ~7,000 km for over 8 hr after the eruption. In addition, density profiles showed wavy structures with a typical length scale of either ~20 km (vertical) or ~1,000 km (horizontal). This may be caused by Lamb waves or gravity waves triggered by the volcanic eruption., Comment: 14 pages, 4 figures, accepted for publication in Geophysical Research Letters

We investigate if and how the valence-bond-solid (VBS) state emerges in the Hubbard model on the honeycomb lattice when the Peierls-type electron-lattice coupling is introduced. We consider all possible lattice-distortion patterns allowed for this lattice model for graphene which preserve the reflection symmetry and determine the most stable configuration in the adiabatic limit by using an unbiased quantum Monte Carlo method. The VBS phase with Kekul\'e dimerization is found to appear as an intermediate phase between a semimetal and an antiferromagnetic Mott insulator for a moderately rigid lattice. This implies that the undistorted semimetallic graphene can be driven into the VBS phase by applying strain, accompanied by the single-particle excitation gap opening., Comment: 8 pages, 9 figures

We propose a quantum-classical hybrid scheme for implementing the nonunitary Gutzwiller factor using a discrete Hubbard-Stratonovich transformation, which allows us to express the Gutzwiller factor as a linear combination of unitary operators involving only single-qubit rotations, at the cost of the sum over the auxiliary fields. To perform the sum over the auxiliary fields, we introduce two approaches that have complementary features. The first approach employs a linear-combination-of-unitaries circuit, which enables one to probabilistically prepare the Gutzwiller wave function on a quantum computer, while the second approach uses importance sampling to estimate observables stochastically, similar to a quantum Monte Carlo method in classical computation. The proposed scheme is demonstrated with numerical simulations for the half-filled Fermi-Hubbard model. Furthermore, we perform quantum simulations using a real quantum device, demonstrating that the proposed scheme can reproduce the exact ground-state energy of the two-site Fermi-Hubbard model within error bars., Comment: 19 pages, 8 figures

We investigate a quantum criticality of an antiferromagnetic phase transition in the Hubbard model on a square lattice with a $d$-wave pairing field by large-scale auxiliary-field quantum Monte Carlo simulations. Since the $d$-wave pairing filed induces Dirac cones in the non-interacting single-particle spectrum, the quantum criticality should correspond to the chiral Heisenberg universality class in terms of the Gross-Neveu theory, which is the same as those expected in the Hubbard model on the honeycomb lattice, despite the unit cells being different (e.g., they contain one and two sites, respectively). We show that both the two phase transitions, expected to occur on the square and on the honeycomb lattices, indeed have the same quantum criticality. We also argue that details of the models, i.e., the way of counting the total number $N$ of fermion components and the anisotropy of the Dirac cones, do not change the critical exponents., Comment: 5 pages, 3 figures, proceeding for the XXXII IUPAP Conference on Computational Physics (CCP2021)

We provide the first unbiased evidence for a higher-order topological Mott insulator in three dimensions by numerically exact quantum Monte Carlo simulations. This insulating phase is adiabatically connected to a third-order topological insulator in the noninteracting limit, which features gapless modes around the corners of the pyrochlore lattice and is characterized by a $\mathbb{Z}_{4}$ spin-Berry phase. The difference between the correlated and non-correlated topological phases is that in the former phase the gapless corner modes emerge only in spin excitations being Mott-like. We also show that the topological phase transition from the third-order topological Mott insulator to the usual Mott insulator occurs when the bulk spin gap solely closes., Comment: 10 pages, 9 figures

We develop a quantum-inspired numerical procedure for searching low-energy states of a classical Hamiltonian composed of two-body fully-connected random Ising interactions and a random local longitudinal magnetic field. In this method, we introduce infinitesimal quantum interactions that do not commute with the original Ising Hamiltonian, and repeatedly generate and truncate direct product states, inspired by the Krylov subspace method, to obtain the low-energy states of the original classical Ising Hamiltonian. The computational cost is controlled by the form of infinitesimal quantum interactions (e.g., one-body or two-body interactions) and the numbers of infinitesimal interaction terms introduced, different initial states considered, and low-energy states kept during the iteration. For a demonstrate of the method, here we introduce as the infinitesimal quantum interactions pair products of Pauli $X$ operators acting on different sites and on-site Pauli $X$ operators into the random Ising Hamiltonian, in which the numerical cost is $O(N^3)$ per iteration with the system size $N$. We consider 120 instances of the random coupling realizations for the random Ising Hamiltonian with $N$ up to 600 and search the 120 lowest-energy states for each instance. We find that the time-to-solution by the quantum-inspired method proposed here, with parallelization in terms of the different initial states, for searching the ground state of the random Ising Hamiltonian scales approximately as $N^5$ for $N$ up to 600. We also examine the basic physical properties such as the ensemble-averaged ground-state and first-excited energies and the ensemble-averaged number of states in the low-energy region of the random Ising Hamiltonian., Comment: 13 pages, 11 figures

We numerically investigate the quantum criticality of the chiral Heisenberg universality class with the total number of fermion components $N$=8 in terms of the Gross-Neveu theory. Auxiliary-field quantum Monte Carlo simulations are performed for the square lattice Hubbard model in the presence of a $d$-wave pairing field, inducing Dirac cones in the single particle spectrum. This property makes the model particularly interesting because it turns out to belong to the same universality class of the Hubbard model on the honeycomb lattice, that is the canonical model for graphene, despite the unit cells being apparently different (e.g. they contain one and two sites, respectively). We indeed show that the two phase transitions, expected to occur on the square and on the honeycomb lattices, have the same quantum criticality. We also argue that details of the models, i.e., the way of counting $N$ and the anisotropy of the Dirac cones, do not change the critical exponents. The present estimates of the exponents for the $N$=8 chiral Heisenberg universality class are $\nu$=1.05(5), $\eta_{\phi}$=0.75(4), and $\eta_{\psi}$=0.23(4), which are compared with the previous numerical estimations., Comment: 12 pages, 12 figures

The solar radiation focusing effect is related to the specific phenomenon of propagation of the Sun-emitted HF and VHF waves through terrestrial ionosphere. This natural effect is observed with ground-based radio instruments running within 10-200 MHz range, as distinctive patterns - the Spectral Caustics (SCs) - on the solar dynamic spectra. It has been suggested that SCs are associated with medium-scale traveling ionospheric disturbances (MSTIDs). In this paper, we present the first direct observations of SCs induced by MSTIDs, using solar dynamic spectra with SCs obtained by different European radio telescopes on January 8, 2014 and simultaneous two-dimensional detrended total electron content (dTEC) maps over Europe. Spatial examination of dTEC maps as well as precise timing analysis of the maps and the dynamic spectra have been performed. First, we found several pairs of one-to-one (TID-SC) correspondences. The study provides strong observational evidence supporting the suggestion that MSTIDs are the cause of SCs., Comment: 10 pages, 5 figures

An unbiased zero-temperature auxiliary-field quantum Monte Carlo method is employed to analyze the nature of the semimetallic phase of the two-dimensional Hubbard model on the honeycomb lattice at half filling. It is shown that the quasiparticle weight $Z$ of the massless Dirac fermions at the Fermi level, which characterizes the coherence of zero-energy single-particle excitations, can be evaluated in terms of the long-distance equal-time single-particle Green's function. If this quantity remains finite in the thermodynamic limit, the low-energy single-particle excitations of the correlated semimetallic phase are described by a Fermi-liquid-type single-particle Green's function. Based on the unprecedentedly large-scale numerical simulations on finite-size clusters containing more than ten thousands sites, we show that the quasiparticle weight remains finite in the semimetallic phase below a critical interaction strength. This is also supported by the long-distance algebraic behavior ($\sim r^{-2}$, where $r$ is distance) of the equal-time single-particle Green's function that is expected for the Fermi liquid. Our result thus provides a numerical confirmation of Fermi-liquid theory in two-dimensional correlated metals., Comment: 10 pages, 6 figures

We investigate a semimetal-superconductor phase transition of two-dimensional Dirac electrons at zero temperature by large-scale and essentially unbiased quantum Monte Carlo simulations for the half-filled attractive Hubbard model on the triangular lattice, in the presence of alternating magnetic $\pi$-flux, that is introduced to construct two Dirac points in the one-particle bands at the Fermi level. This phase transition is expected to describe quantum criticality of the chiral XY class in the framework of the Gross-Neveu model, where, in the ordered phase, the $U(1)$ symmetry is spontaneously broken and a mass gap opens in the excitation spectrum. We compute the order parameter of the s-wave superconductivity and estimate the quasiparticle weight from the long-distance behavior of the single-particle Green's function. These calculations allow us to obtain the critical exponents of this transition in a reliable and accurate way. Our estimate for the critical exponents is in good agreement with those obtained for a transition to a Kekul\'{e} valence bond solid, where an emergent $U(1)$ symmetry is proposed [Z.-X. Li et al., Nat. Commun. 8, 314 (2017)]., Comment: 9 pages, 7 figures

The metal-insulator transition has been a subject of intense research since Nevil Mott has first proposed that the metallic behavior of interacting electrons could turn to the insulating one as electron correlations increase. Here, we consider electrons with massless Dirac-like dispersion in two spatial dimensions, described by the Hubbard models on two geometrically different lattices, and perform numerically exact calculations on unprecedentedly large systems that, combined with a careful finite size scaling analysis, allow us to explore the quantum critical behavior in the vicinity of the interaction-driven metal-insulator transition. We find thereby that the transition is continuous and determine the quantum criticality for the corresponding universality class, which is described in the continuous limit by the Gross-Neveu model, a model extensively studied in quantum field theory. We furthermore discuss a fluctuation-driven scenario for the metal-insulator transition in the interacting Dirac electrons: the metal-insulator transition is triggered only by the vanishing of the quasiparticle weight but not the Dirac Fermi velocity, which instead remains finite near the transition. This important feature cannot be captured by a simple mean-field or Gutzwiller-type approximate picture, but is rather consistent with the low energy behavior of the Gross-Neveu model., Comment: 18 pages, 17 figures; to be published in Phys. Rev. X