Quantum-number projection in the path-integral renormalization group method

We present a quantum-number projection technique which enables us to exactly treat spin, momentum and other symmetries embedded in the Hubbard model. By combining this projection technique, we extend the path-integral renormalization group method to improve the efficiency of numerical computations. By taking numerical calculations for the standard Hubbard model and the Hubbard model with next nearest neighbor transfer, we show that the present extended method can extremely enhance numerical accuracy and that it can handle excited states, in addition to the ground state.Comment: 11 pages, 7 figures, submitted to Phys. Rev.

Quantum Phase Transitions to Charge Order and Wigner Crystal Under Interplay of Lattice Commensurability and Long-Range Coulomb Interaction

Relationship among Wigner crystal, charge order and Mott insulator is studied by the path-integral renormalization group method for two-dimensional lattices with long-range Coulomb interaction. In contrast to Hartree-Fock results, the solid stability drastically increases with lattice commensurability. The transition to liquid occurs at the electron gas parameter $r_s \sim 2$ for the filling $n=1/2$ showing large reduction from $r_s \sim 35$ in the continuum limit. Correct account of quantum fluctuations are crucial to understand charge-order stability generally observed only at simple fractional fillings and nature of quantum liquids away from them.Comment: 4 pages including 7 figure

Drude Weight of the Two-Dimensional Hubbard Model -- Reexamination of Finite-Size Effect in Exact Diagonalization Study --

The Drude weight of the Hubbard model on the two-dimensional square lattice is studied by the exact diagonalizations applied to clusters up to 20 sites. We carefully examine finite-size effects by consideration of the appropriate shapes of clusters and the appropriate boundary condition beyond the imitation of employing only the simple periodic boundary condition. We successfully capture the behavior of the Drude weight that is proportional to the squared hole doping concentration. Our present result gives a consistent understanding of the transition between the Mott insulator and doped metals. We also find, in the frequency dependence of the optical conductivity, that the mid-gap incoherent part emerges more quickly than the coherent part and rather insensitive to the doping concentration in accordance with the scaling of the Drude weight.Comment: 9 pages with 10 figures and 1 table. accepted in J. Phys. Soc. Jp

Suppressed Coherence due to Orbital Correlations in the Ferromagnetically Ordered Metallic Phase of Mn Compounds

Small Drude weight $D$ together with small specific heat coefficient $\gamma$ observed in the ferromagnetic phase of R$_{1-x}$A$_x$MnO$_3$ (R=La, Pr, Nd, Sm; A=Ca, Sr, Ba) are analyzed in terms of a proximity effect of the Mott insulator. The scaling theory for the metal-insulator transition with the critical enhancement of orbital correlations toward the staggered ordering of two $e_g$ orbitals such as $3x^2-r^2$ and $3y^2-r^2$ symmetries may lead to the critical exponents of $D \propto \delta^{u}$ and $\gamma \propto \delta^v$ with $u=3/2$ and $v=0$. The result agrees with the experimental indications.Comment: 4 pages LaTeX using jpsj.sty. To appear in J. Phys. Soc. Jpn. 67(1998)No.

Improvement of solar cycle prediction: Plateau of solar axial dipole moment

Aims. We report the small temporal variation of the axial dipole moment near the solar minimum and its application to the solar cycle prediction by the surface flux transport (SFT) model. Methods. We measure the axial dipole moment using the photospheric synoptic magnetogram observed by the Wilcox Solar Observatory (WSO), the ESA/NASA Solar and Heliospheric Observatory Michelson Doppler Imager (MDI), and the NASA Solar Dynamics Observatory Helioseismic and Magnetic Imager (HMI). We also use the surface flux transport model for the interpretation and prediction of the observed axial dipole moment. Results. We find that the observed axial dipole moment becomes approximately constant during the period of several years before each cycle minimum, which we call the axial dipole moment plateau. The cross-equatorial magnetic flux transport is found to be small during the period, although the significant number of sunspots are still emerging. The results indicates that the newly emerged magnetic flux does not contributes to the build up of the axial dipole moment near the end of each cycle. This is confirmed by showing that the time variation of the observed axial dipole moment agrees well with that predicted by the SFT model without introducing new emergence of magnetic flux. These results allows us to predict the axial dipole moment in Cycle 24/25 minimum using the SFT model without introducing new flux emergence. The predicted axial dipole moment of Cycle 24/25 minimum is 60--80 percent of Cycle 23/24 minimum, which suggests the amplitude of Cycle 25 even weaker than the current Cycle 24. Conclusions. The plateau of the solar axial dipole moment is an important feature for the longer prediction of the solar cycle based on the SFT model.Comment: 5 pages, 3 figures, accepted for publication in A&A Lette

Quantification of propagation modes in an astronomical instrument from its radiation pattern

Context. Understanding complex phenomena and unsolved problems in modern astronomy requires wider-bandwidth observations. The current technique for designing and fabricating an astronomical instrument potentially provides such observations with higher efficiency and precision than in the past. Higher-order modes in an instrument associated with wider bandwidths have been reported, which may degrade observation precision. Aims. To reduce the unfavorable degradation, we need to quantify the higher-order propagation modes, though their power is too difficult to measure directly. Instead of the direct mode measurement, we aim at developing a method based on measurable radiation patterns from an instrument of interest. Method. 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 inside modes and the outside radiated field and calculate the pseudo-inverse matrix. To understand the estimation accuracy of the proposed method, we demonstrate two cases with numerical simulations, axially-corrugated horn case and offset Cassegrain antenna case, and investigate the effect of random errors on the accuracy. Results. Both cases showed the estimated 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. The demonstrated method works independently of the details of a system. Conclusions. The method can quantify the propagation modes inside an instrument and will be applicable to most of linear components and antennas. This method can be employed for a general purpose, such as diagnosis of feed alignment and higher-performance feed design.Comment: 7 pages, 11 figure

Fate of Quasiparticle at Mott Transition and Interplay with Lifshitz Transition Studied by Correlator Projection Method

Filling-control metal-insulator transition on the two-dimensional Hubbard model is investigated by using the correlator projection method, which takes into account momentum dependence of the free energy beyond the dynamical mean-field theory. The phase diagram of metals and Mott insulators is analyzed. Lifshitz transitions occur simultaneously with metal-insulator transitions at large Coulomb repulsion. On the other hand, they are separated each other for lower Coulomb repulsion, where the phase sandwiched by the Lifshitz and metal-insulator transitions appears to show violation of the Luttinger sum rule. Through the metal-insulator transition, quasiparticles retain nonzero renormalization factor and finite quasi-particle weight in the both sides of the transition. This supports that the metal-insulator transition is caused not by the vanishing renormalization factor but by the relative shift of the Fermi level into the Mott gap away from the quasiparticle band, in sharp contrast with the original dynamical mean-field theory. Charge compressibility diverges at the critical end point of the first-order Lifshitz transition at finite temperatures. The origin of the divergence is ascribed to singular momentum dependence of the quasiparticle dispersion.Comment: 24 pages including 10 figure

Spin-gap phase in nearly-half-filled one-dimensional conductors coupled with phonons

Asymptotic properties of nearly-half-filled one-dimensional conductors coupled with phonons are studied through a renormalization group method. Due to spin-charge coupling via electron-phonon interaction, the spin correlation varies with filling as well as the charge correlation. Depending on the relation between cut-off energy scales of the Umklapp process and of the electron-phonon interaction, various phases appear. We found a metallic phase with a spin gap and a dominant charge- density-wave correlation near half filling between a gapless density-wave phase (like in the doped repulsive Hubbard model) and a superconductor phase with a spin gap. The spin gap is produced by phonon-assisted backward scatterings which are interfered with the Umklapp process constructively or destructively depending on the character of electron-phonon coupling.Comment: 14 pages, revtex, replaced 5 ps figures, published in PR