Source author record

Satoshi Morita

Satoshi Morita appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

TeNeS: Tensor Network Solver for Quantum Lattice Systems

TeNeS (Tensor Network Solver) is a free/libre open-source software program package for calculating two-dimensional many-body quantum states based on the tensor network method and the corner transfer matrix renormalization group (CTMRG) method. This package calculates ground-state wavefunctions for user-defined Hamiltonians and evaluates user-defined physical quantities such as magnetization and correlation functions. For certain predefined models and lattices, there is a tool that makes it easy to generate input files. TeNeS uses an OpenMP/MPI hybrid parallelized tensor operation library and thus can perform large-scale calculations using massively parallel machines.

preprint2021arXiv

A Three-Dimensional Velocity of an Erupting Prominence Prior to a Coronal Mass Ejection

We present a detailed three-dimensional (3D) view of a prominence eruption, coronal loop expansion, and coronal mass ejections (CMEs) associated with an M4.4 flare that occurred on 2011 March 8 in the active region NOAA 11165. Full-disk H$α$ images of the flare and filament ejection were successfully obtained by the Flare Monitoring Telescope (FMT) following its relocation to Ica University, Peru. Multiwavelength observation around the H$α$ line enabled us to derive the 3D velocity field of the H$α$ prominence eruption. Features in extreme ultraviolet were also obtained by the Atmospheric Imager Assembly onboard the {\it Solar Dynamic Observatory} and the Extreme Ultraviolet Imager on board the {\it Solar Terrestrial Relations Observatory - Ahead} satellite. We found that, following collision of the erupted filament with the coronal magnetic field, some coronal loops began to expand, leading to the growth of a clear CME. We also discuss the succeeding activities of CME driven by multiple interactions between the expanding loops and the surrounding coronal magnetic field.

preprint2021arXiv

Multiple magnetization plateaus induced by farther neighbor interaction in an $S = 1$ two-leg Heisenberg spin ladder

We study the magnetization process of the $S=1$ Heisenberg model on a two-leg ladder with farther neighbor spin-exchange interaction. We consider the interaction that couples up to the next-nearest neighbor rungs and find an exactly solvable regime where the ground states become product states. The next-nearest neighbor interaction tends to stabilize magnetization plateaus at multiples of 1/6. In most of the exactly solvable regime, a single magnetization curve shows two series of plateaus with different periodicities.

preprint2020arXiv

Boundary conformal spectrum and surface critical behaviors of the classical spin systems: a tensor network renormalization study

We numerically obtain the conformal spectrum of several classical spin models on a two-dimensional lattice with open boundaries, for every boundary fixed point obtained by the Cardy's derivation [J. L. Cardy, Nucl. Phys. B 324, 581 (1989)]. In order to extract accurate conformal data, we implement the tensor network renormalization algorithm [G. Evenbly and G. Vidal, Phys. Rev. Lett. 115, 180405 (2015)] extended so as to be applicable to a square lattice with open boundaries. We successfully compute the boundary conformal spectrum consistent with the underlying boundary conformal field theories (BCFTs) for the Ising, tri-critical Ising, and 3-state Potts models on the lattice, which allows us to confirm the validity of the BCFT analyses for the surface critical behaviors of those lattice models.

preprint2020arXiv

Continuous phase transition between Neel and valence bond solid phases in a J-Q-like spin ladder system

We investigate a quantum phase transition between a Neel phase and a valence bond solid (VBS) phase, in each of which a different Z2 symmetry is broken, in a spin-1/2 two-leg XXZ ladder with a four-spin interaction. The model can be viewed as a one-dimensional variant of the celebrated J-Q model on a square lattice. By means of variational uniform matrix product state calculations and an effective field theory, we determine the phase diagram of the model, and present evidences that the Neel-VBS transition is continuous and belongs to the Gaussian universality class with the central charge c=1. In particular, the critical exponents $β, η,$ and, $ν$ are found to satisfy the constraints expected for a Gaussian transition within numerical accuracy. These exponents do not detectably change along the phase boundary while they are in general allowed to do so for the Gaussian class.

preprint2020arXiv

Global optimization of tensor renormalization group using the corner transfer matrix

A tensor network renormalization algorithm with global optimization based on the corner transfer matrix is proposed. Since the environment is updated by the corner transfer matrix renormalization group method, the forward-backward iteration is unnecessary, which is a time-consuming part of other methods with global optimization. In addition, a further approximation reducing the order of the computational cost of contraction for the calculation of the coarse-grained tensor is proposed. The computational time of our algorithm in two dimensions scales as the sixth power of the bond dimension while the higher-order tensor renormalization group and the higher-order second renormalization group methods have the seventh power. We perform benchmark calculations in the Ising model on the square lattice and show that the time-to-solution of the proposed algorithm is faster than that of other methods.

preprint2018arXiv

Calculation of higher-order moments by higher-order tensor renormalization group

A calculation method for higher-order moments of physical quantities, including magnetization and energy, based on the higher-order tensor renormalization group is proposed. The physical observables are represented by impurity tensors. A systematic summation scheme provides coarse-grained tensors including multiple impurities. Our method is compared with the Monte Carlo method on the two-dimensional Potts model. While the nature of the transition of the $q$-state Potts model has been known for a long time owing to the analytical arguments, a clear numerical confirmation has been difficult due to extremely long correlation length in the weakly first-order transitions, e.g., for $q=5$. A jump of the Binder ratio precisely determines the transition temperature. The finite-size scaling analysis provides critical exponents and distinguishes the weakly first-order and the continuous transitions.

preprint2016arXiv

Pre-eruption Oscillations in Thin and Long Features in a Quiescent Filament

We investigate the eruption of a quiescent filament located close to an active region. Large-scale activation was observed in only half of the filament in the form of pre-eruption oscillations. Consequently only this half erupted nearly 30 hr after the oscillations commenced. Time-slice diagrams of 171 Å images from Atmospheric Imaging Assembly were used to study the oscillations. The oscillations were observed in several thin and long features connecting the filament spine to the chromosphere below. This study traces the origin of such features and proposes their possible interpretation. Small-scale magnetic flux cancellation accompanied by a brightening was observed at the footpoint of the features shortly before their appearance, in images recorded by the Helioseismic and Magnetic Imager. Slow rise of the filament was detected in addition to the oscillations, indicating a gradual loss of equilibrium. Our analysis indicates that a change in magnetic field connectivity between two neighbouring active regions and the quiescent filament resulted in weakening of the overlying arcade of the filament leading to its eruption. It is also suggested that the oscillating features are filament barbs, and the oscillations are a manifestation during the pre-eruption phase of filaments.

preprint2016arXiv

Temporal evolution and spatial distribution of white-light flare kernels in a solar flare

On 2011 September 6, we observed an X2.1-class flare in continuum and H$α$ with a frame rate of about 30~Hz. After processing images of the event by using a speckle-masking image reconstruction, we identified white-light (WL) flare ribbons on opposite sides of the magnetic neutral line. We derive the lightcurve decay times of the WL flare kernels at each resolution element by assuming that the kernels consist of one or two components that decay exponentially, starting from the peak time. As a result, 42% of the pixels have two decay-time components with average decay times of 15.6 and 587 s, whereas the average decay time is 254 s for WL kernels with only one decay-time component. The peak intensities of the shorter decay-time component exhibit good spatial correlation with the WL intensity, whereas the peak intensities of the long decay-time components tend to be larger in the early phase of the flare at the inner part of the flare ribbons, close to the magnetic neutral line. The average intensity of the longer decay-time components is 1.78 times higher than that of the shorter decay-time components. If the shorter decay time is determined by either the chromospheric cooling time or the nonthermal ionization timescale and the longer decay time is attributed to the coronal cooling time, this result suggests that WL sources from both regions appear in 42% of the WL kernels and that WL emission of the coronal origin is sometimes stronger than that of chromospheric origin.

preprint2015arXiv

Phase transition of two-dimensional Ising models on the honeycomb and related lattices with striped random impurities

Two-dimensional Ising models on the honeycomb lattice and the square lattice with striped random impurities are studied to obtain their phase diagrams. Assuming bimodal distributions of the random impurities where all the non-zero couplings have the same magnitude, exact critical values for the fraction p of ferromagnetic bonds at the zero-temperature (T=0) are obtained. The critical lines in the p-T plane are drawn by numerically evaluating the Lyapunov exponent of random matrix products.

preprint2015arXiv

Quantum Spin Liquid in Spin 1/2 J1-J2 Heisenberg Model on Square Lattice: Many-Variable Variational Monte Carlo Study Combined with Quantum-Number Projections

The nature of quantum spin liquids is studied for the spin-$1/2$ antiferromagnetic Heisenberg model on a square lattice containing exchange interactions between nearest-neighbor sites, $J_1$, and those between next-nearest-neighbor sites, $J_2$. We perform variational Monte Carlo simulations together with the quantum-number-projection technique and clarify the phase diagram in the ground state together with its excitation spectra. We obtain the nonmagnetic phase in the region $0.4< J_2/J_1\le 0.6$ sandwiched by the staggered and stripe antiferromagnetic (AF) phases. Our direct calculations of the spin gap support the notion that the triplet excitation from the singlet ground state is gapless in the range of $0.4 < J_2/J_1 \le 0.5$, while the gapped valence-bond-crystal (VBC) phase is stabilized for $0.5 < J_2/J_1 \le 0.6$. The VBC order is likely to have the columnar symmetry with a spontaneous symmetry breaking of the $C_{4v}$ symmetry. The power-law behaviors of the spin-spin and dimer-dimer correlation functions in the gapless region are consistent with the emergence of the algebraic quantum-spin-liquid phase (critical phase). The exponent of the spin correlation $\langle S(0)S(r)\rangle \propto 1/r^{z+η}$ at a long distance $r$ appears to increase from $z+η\sim 1$ at $J_2/J_1\sim0.4$ toward the continuous transition to the VBC phase at $J_1/J_1\sim0.5$. Our results, however, do not fully exclude the possibility of a direct quantum transition between the staggered AF and VBC phases with a wide critical region and deconfined criticality.

preprint2015arXiv

Variational Monte Carlo Method in the Presence of Spin-Orbit Interaction and Its Application to the Kitaev and Kitaev-Heisenberg Model

We propose an accurate variational Monte Carlo method applicable in the presence of the strong spin-orbit interaction. Our variational wave functions consist of generalized Pfaffian-Slater wave functions that involve mixtures of singlet and triplet Cooper pairs, Jastrow-Gutzwiller-type pro- jections, and quantum number projections. The generalized wave functions allow any symmetry breaking states, ranging from magnetic and/or charge ordered states to superconducting states and their fluctuations, on equal footing without any ad hoc ansatz for variational wave functions. We de- tail our optimization scheme for the generalized Pfaffian-Slater wave functions with complex-number variational parameters. Generalized quantum number projections are also introduced, which im- poses the conservation of not only spin and momentum quantum numbers but also Wilson loops. As a demonstration of the capability of the present variational Monte Carlo method, the accuracy and efficiency of the present method is tested for the Kitaev and Kitaev-Heisenberg models. The Kitaev model serves as a critical benchmark of the present method: The exact ground state of the model is a typical gapless quantum spin liquid far beyond the applicability range of simple mean-field wave functions. The newly introduced quantum number projections precisely reproduce the ground state degeneracy of the Kitaev spin liquids, in addition to their ground state energy. Our framework offers accurate solutions for the systems where strong electron correlation and spin-orbit interaction coexist.

preprint2014arXiv

Snapshot spectrum and critical phenomenon for two-dimensional classical spin systems

We investigate the eigenvalue distribution of the snapshot density matrix (SDM) generated by Monte Carlo simulation for two-dimensional classical spin systems. We find that the distribution in the high-temperature limit is well explained by the random-matrix theory, while that in the low-temperature limit can be characterized by the zero-eigenvalue condensation. At the critical point, we obtain the power-law distribution with a nontrivial exponent $α\equiv(2-η)/(1-η)$ and the asymptotic form of the snapshot entropy, on the basis of the relationship of the SDM with the correlation function matrix. The aspect-ratio dependence of the SDM spectrum is also mentioned.

preprint2013arXiv

Improved Multi-Variable Variational Monte Carlo Method Examined by High-Precision Calculations of One-Dimensional Hubbard Model

We revisit the accuracy of the variational Monte Carlo (VMC) method by taking an example of ground state properties for the one-dimensional Hubbard model. We start from the variational wave functions with the Gutzwiller and long-range Jastrow factor introduced by Capello et al. [Phys. Rev. B 72, 085121 (2005)] and further improve it by considering several quantum-number projections and a generalized one-body wave function. We find that the quantum spin projection and total momentum projection greatly improve the accuracy of the ground state energy within 0.5% error, for both small and large systems at half filling. Besides, the momentum distribution function n(k) at quarter filling calculated up to 196 sites allows us direct estimate of the critical exponents of the charge correlations from the power-law behavior of n(k) near the Fermi wave vector. Estimated critical exponents well reproduce those predicted by the Tomonaga-Luttinger theory.

preprint2011arXiv

First Simultaneous Observation of H-alpha Moreton Wave, EUV Wave, and Filament/Prominence Oscillations

We report on the first simultaneous observation of an H-alpha Moreton wave, the corresponding EUV fast coronal waves, and a slow and bright EUV wave (typical EIT wave). Associated with an X6.9 flare that occurred on 2011 August 9 at the active region NOAA 11263, we observed a Moreton wave in the H-alpha images taken by the Solar Magnetic Activity Research Telescope (SMART) at Hida Observatory of Kyoto University. In the EUV images obtained by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) we found not only the corresponding EUV fast "bright" coronal wave, but also the EUV fast "faint" wave that is not associated with the H-alpha Moreton wave. We also found a slow EUV wave, which corresponds to a typical EIT wave. Furthermore, we observed, for the first time, the oscillations of a prominence and a filament, simultaneously, both in the H-alpha and EUV images. To trigger the oscillations by the flare-associated coronal disturbance, we expect a coronal wave as fast as the fast-mode MHD wave with the velocity of about 570 - 800 km/s. These velocities are consistent with those of the observed Moreton wave and the EUV fast coronal wave.

preprint2010arXiv

Observations of Chromospheric Anemone Jets with Hinode SOT and Hida Ca II Spectroheliogram

We present the first simultaneous observations of chromospheric "anemone" jets in solar active regions with Hinode SOT Ca II H broadband filetergram and Ca II K spetroheliogram on the Domeless Solar Telescope (DST) at Hida Observatory. During the coordinated observation, 9 chromospheric anemone jets were simultaneously observed with the two instruments. These observations revealed three important features, i.e.: (1) the jets are generated in the lower chromosphere, (2) the length and lifetime of the jets are 0.4-5 Mm and 40-320 sec, (3) the apparent velocity of the jets with Hinode SOT are 3-24 km/s, while Ca II K3 component at the jets show blueshifts (in 5 events) in the range of 2- 6 km/s. The chromospheric anemone jets are associated with mixed polarity regions which are either small emerging flux regions or moving magnetic features. It is found that the Ca II K line often show red or blue asymmetry in K2/K1 component: the footpoint of the jets associated with emerging flux regions often show redshift (2-16 km/s), while the one with moving magnetic features show blueshift (around 5 km/s). Detailed analysis of magnetic evolution of the jet foaming regions revealed that the reconnection rate (or canceling rate) of the total magnetic flux at the footpoint of the jets are of order of 10^{16} Mx/s, and the resulting magnetic energy release rate (1.1-10) x 10^{24} erg/s, with the total energy release (1-13) x 10^{26} erg for the duration of the magnetic cancellations, 130s. These are comparable to the estimated total energy, 10^{26} erg, in a single chromospheric anemone jet. An observation-based physical model of the jet is presented. The relation between chromospheric anemone jets and Ellerman bombs is discussed.

preprint2007arXiv

Faster annealing schedules for quantum annealing

New annealing schedules for quantum annealing are proposed based on the adiabatic theorem. These schedules exhibit faster decrease of the excitation probability than a linear schedule. To derive this conclusion, the asymptotic form of the excitation probability for quantum annealing is explicitly obtained in the limit of long annealing time. Its first-order term, which is inversely proportional to the square of the annealing time, is shown to be determined only by the information at the initial and final times. Our annealing schedules make it possible to drop this term, thus leading to a higher order (smaller) excitation probability. We verify these results by solving numerically the time-dependent Schrodinger equation for small size systems