Source author record

Makoto Yamashita

Makoto Yamashita 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

28works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

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

Published work

28 published item(s)

preprint2026arXiv

An inexact infeasible arc-search interior-point method for linear optimization problems

We propose an inexact infeasible arc-search interior-point method for solving linear optimization problems. The method combines an arc-search strategy with inexact solutions to Newton systems and admits a polynomial iteration complexity bound. In existing inexact infeasible interior-point methods, both the linearization error of the central path and the inexactness of the Newton system accumulate along the search direction, which forces the algorithm to take very small steps. The proposed method mitigates this effect by using an arc-search strategy: the curved search path provides a more accurate approximation of the central path, so the step size can remain larger even when the Newton system is solved inexactly. As a result, the proposed method achieves a provably tighter worst-case iteration bound than existing inexact infeasible line-search methods. Numerical experiments on NETLIB benchmark problems demonstrate that the proposed method reduces both the number of iterations and the computation time.

preprint2022arXiv

Exact SDP relaxations for quadratic programs with bipartite graph structures

For nonconvex quadratically constrained quadratic programs (QCQPs), we first show that, under certain feasibility conditions, the standard semidefinite (SDP) relaxation is exact for QCQPs with bipartite graph structures. The exact optimal solutions are obtained by examining the dual SDP relaxation and the rank of the optimal solution of this dual SDP relaxation under strong duality. Our results on the QCQPs generalize the results on QCQP with sign-definite bipartite graph structures, QCQPs with forest structures, and QCQPs with nonpositive off-diagonal data elements. Second, we propose a conversion method from QCQPs with no particular structure to the ones with bipartite graph structures. As a result, we demonstrate that a wider class of QCQPs can be exactly solved by the SDP relaxation. Numerical instances are presented for illustration.

preprint2021arXiv

Homology and K-theory of dynamical systems. I. torsion-free ample groupoids

Given an ample groupoid, we construct a spectral sequence with groupoid homology with integer coefficients on the second sheet, converging to the K-groups of the (reduced) groupoid C*-algebra, provided the groupoid has torsion-free stabilizers and satisfies a strong form of the Baum-Connes conjecture. The construction is based on the triangulated category approach to the Baum-Connes conjecture developed by Meyer and Nest. We also present a few applications to topological dynamics and discuss the HK conjecture of Matui.

preprint2020arXiv

Exact SDP relaxations of quadratically constrained quadratic programs with forest structures

We study the exactness of the semidefinite programming (SDP) relaxation of quadratically constrained quadratic programs (QCQPs). With the aggregate sparsity matrix from the data matrices of a QCQP with $n$ variables, the rank and positive semidefiniteness of the matrix are examined. We prove that if the rank of the aggregate sparsity matrix is not less than $n-1$ and the matrix remains positive semidefinite after replacing some off-diagonal nonzero elements with zeros, then the standard SDP relaxation provides an exact optimal solution for the QCQP under feasibility assumptions. In particular, we demonstrate that QCQPs with forest-structured aggregate sparsity matrix, such as the tridiagonal or arrow-type matrix, satisfy the exactness condition on the rank. The exactness is attained by considering the feasibility of the dual SDP relaxation, the strong duality of SDPs, and a sequence of QCQPs with perturbed objective functions, under the assumption that the feasible region is compact. We generalize our result for a wider class of QCQPs by applying simultaneous tridiagonalization on the data matrices. Moreover, simultaneous tridiagonalization is applied to a matrix pencil so that QCQPs with two constraints can be solved exactly by the SDP relaxation.

preprint2016arXiv

A Successive LP Approach with C-VaR Type Constraints for IMRT Optimization

Radiation therapy is considered to be one of important treatment protocols for cancers. Radiation therapy employs several beams of ionizing radiation to kill cancer tumors, but such irradiation also causes damage to normal tissues. Therefore, a treatment plan should satisfy dose-volume constraints (DVCs). Intensity-modulated radiotherapy treatment (IMRT) enables to control the beam intensities and gives more flexibility for the treatment plan to satisfy the DVCs. Romeijn et al. (2003) replaced the DVCs in an IMRT optimization with C-VaR (Conditional Value-at-Risk) type constraints, and proposed a numerical method based on linear programming (LP). Their approach reduced the computation cost of the original DVCs, but the feasible region of their LP problems was much narrow compared to the DVCs, therefore, their approach often failed to find a feasible plan even when the DVCs were not so tight. In this paper, we propose a successive LP approach with the C-VaR type constraints. We detect outliers form the solution of LP problems, and remove them from the domain of the C-VaR type constraints. This eases the sensitivity of C-VaR type constraints to outliers and we can search feasible plans from wider regions. Furthermore, we can give a mathematical proof that if the optimal value of the LP problem in the proposed approach is non-positive, the corresponding optimal solution satisfies all the DVCs. From numerical experiments on test data sets, we observed that our approach found feasible solutions more appropriately than existing LP approaches. In addition, our approach required fewer LP problems, and this led to a short computation time.

preprint2016arXiv

Dynamical instability in the S=1 Bose-Hubbard model

We study the dynamical instabilities of superfluid flows in the S=1 Bose-Hubbard model. The time evolution of each spin component in a condensate is calculated based on the dynamical Gutzwiller approximation for a wide range of interactions, from a weakly correlated regime to a strongly correlated regime near the Mott-insulator transition. Owing to the spin-dependent interactions, the superfluid flow of the spin-1 condensate decays at a different critical momentum from a spinless case when the interaction strength is the same. We furthermore calculate the dynamical phase diagram of this model and clarify that the obtained phase boundary has very different features depending on whether the average number of particles per site is even or odd. Finally, we analyze the density and spin modulations that appear in association with the dynamical instability. We find that spin modulations are highly sensitive to the presence of a uniform magnetic field.

preprint2016arXiv

Spin-1 Quantum Walks

We study the quantum walks of two interacting spin-1 bosons. We derive an exact solution for the time-dependent wave function that describes the two-particle dynamics governed by the one-dimensional spin-1 Bose-Hubbard model. We show that propagation dynamics in real space and mixing dynamics in spin space are correlated via the spin-dependent interaction in this system. The spin-mixing dynamics has two characteristic frequencies in the limit of large spin-dependent interactions. One of the characteristic frequencies is determined by the energy difference between two bound states, and the other frequency relates to the cotunneling process of a pair of spin-1 bosons. Furthermore, we numerically analyze the growth of the spin correlations in quantum walks. We find that long-range spin correlations emerge showing a clear dependence on the sign of the spin-dependent interaction and the initial state.

preprint2015arXiv

An efficient second-order cone programming approach for optimal selection in tree breeding

An important problem in tree breeding is optimal selection from candidate pedigree members to produce the highest performance in seed orchards, while conserving essential genetic diversity. The most beneficial members should contribute as much as possible, but such selection of orchard parents would reduce performance of the orchard progeny due to serious inbreeding. To avoid inbreeding, we should include a constraint on the numerator relationship matrix to keep a group coancestry under an appropriate threshold. Though an SDP (semidefinite programming) approach proposed by Pong-Wong and Woolliams gave an accurate optimal value, it required rather long computation time. In this paper, we propose an SOCP (second-order cone programming) approach to reduce this computation time. We demonstrate that the same solution is attained by the SOCP formulation, but requires much less time. Since a simple SOCP formulation is not much more efficient compared to the SDP approach, we exploit a sparsity structure of the numerator relationship matrix, and formulate the SOCP constraint using Henderson's algorithm. Numerical results show that the proposed SOCP approach reduced computation time in a case study from 39,200 seconds under the SDP approach to less than 2 seconds.

preprint2015arXiv

An iterative method using boundary distance for box-constrained nonlinear semidefinite programs

We propose an iterative method for nonlinear semidefinite programs with box constraints. The search direction in the proposed method utilizes the distance from the current point to the boundary of a feasible set. The computation of the search direction exploits the second derivative of the objective function only in a quadratic form, and this property saves the computation cost compared to an evaluation of the whole entries of the second derivative. We compute a step length in an interval determined by a radius and we update the radius using a quadratic approximation function. In this paper, we also discuss convergence properties of the proposed method based on structures of the search direction. Numerical tests show that the proposed method solves problems in which the size of a variable matrix is larger than 5,000 and that it is faster than a feasible direction method for objective functions with strong nonlinearity.

preprint2015arXiv

BCS superconducting transitions in lattice fermions

We develop a general description of the superconductivity of lattice fermions based on the BCS theory. We propose a modeling of the density of states (DOS) of lattice fermions, where divergent and semi-metallic structures are described by asymptotic expansions around the Fermi energy. This modeling leads to a unified representation of the transition temperature $T_c$ at half filling, which reproduces asymptotic forms of $T_c$ derived in several lattices, such as the square, honeycomb, and Lieb lattices, for the weakly interacting limit. The derived asymptotic forms of $T_c$ are categorized into four types, which is attributed to the different responses of the degenerate fermions depending on the DOS structures. The DOS with a delta-functional singularity induces the highest $T_c$ in the weakly interacting region, where $T_c$ is linearly proportional to the pairing interaction $U$. Three kinds of universal ratios defined in the BCS theory no longer reduce to constants independent of the system parameters but can be parameterized with a certain variable that characterizes the singular structures of the DOS. We find universal relationship among thermodynamic quantities that holds for all parameter regions. Further, we numerically demonstrate that in multi-band systems the correlation effects can induce an effective delta-functional singularity. This phenomenon generally appears in multi-energy (or multi-gap) systems and may provide a plausible guideline for material designs of high-$T_c$ superconductor.

preprint2015arXiv

Magnetism in the three-dimensional layered Lieb lattice: Enhanced transition temperature via flat-band and Van Hove singularities

We describe the enhanced magnetic transition temperatures $T_c$ of two-component fermions in three-dimensional layered Lieb lattices, which are created in cold atom experiments. We determine the phase diagram at half-filling using the dynamical mean-field theory. The dominant mechanism of enhanced $T_c$ gradually changes from the (delta-functional) flat-band to the (logarithmic) Van Hove singularity as the interlayer hopping increases. We elucidate that the interaction induces an effective flat-band singularity from a dispersive flat (or narrow) band. We offer a general analytical framework for investigating the singularity effects, where a singularity is treated as one parameter in the density of states. This framework provides a unified description of the singularity-induced phase transitions, such as magnetism and superconductivity, where the weight of the singularity characterizes physical quantities. This treatment of the flat-band provides the transition temperature and magnetization as a universal form (i.e., including the Lambert function). We also elucidate a specific feature of the magnetic crossover in magnetization at finite temperatures.

preprint2015arXiv

Theoretical Analysis on Spectroscopy of Atomic Bose-Hubbard Systems

We provide a numerical method to calculate comprehensively the microwave and the laser spectra of ultracold bosonic atoms in optical lattices at finite temperatures. Our formulation is built up with the sum rules, up to the second order, derived from the general principle of spectroscopy. The sum rule approach allows us to discuss the physical origins of a spectral peak shift and also a peak broadening. We find that a spectral broadening of superfluid atoms can be determined from number fluctuations of atoms, while that of normal-state atoms is mainly attributed to quantum fluctuations resulting from hopping of atoms. To calculate spectra at finite temperatures, based on the sum rule approach, we provide a two-mode approximation assuming that spectra of the superfluid and normal state atoms can be calculated separately. Our method can properly deal with multi-peak structures of spectra resulting from thermal fluctuations and also coexisting of the superfluid and the normal states. By combining the two-mode approximation with a finite temperature Gutzwiller approximation, we calculate spectra at finite temperatures by considering realistic systems, and the calculated spectra show nice agreements with those in experiments.

preprint2014arXiv

Density Modulations Associated with the Dynamical Instability in the Bose-Hubbard Model

A superfluid flow beyond a critical momentum in an optical lattice decays drastically by the interplay between nonlinearity due to the interparticle interactions in Bose-Einstein condensate and periodicity of the lattice; this instability is called dynamical instability. The complex density modulational profiles after the condensate becomes unstable observed experimentally is not completely understood, while the dynamical instability has been studied theoretically and experimentally. In this paper, we analyze the density modulation of condensates as a precursor of the dynamical instability in the two-dimensional Bose-Hubbard model. Our analysis has clarified the unexplored properties of the density modulations associated with the dynamical instability at low filling and in a wide range of interactions, while the previous works have analyzed the density modulation on the basis of Gross-Pitaevskii equation under the specific condition that one-dimensional optical lattice is very shallow and the filling is very large. The numerical simulations based on the dynamical Gutzwiller approximation elucidate that the principal mode of density modulation highly depends on interaction strength U and the momentum acceleration rate. We briefly discuss these features with the stability phase diagram calculated on the basis of the Bogoliubov theory.

preprint2014arXiv

Fast implementation for semidefinite programs with positive matrix completion

Solving semidefinite programs (SDP) in a short time is the key to managing various mathematical optimization problems. The matrix-completion primal-dual interior-point method (MC-PDIPM) extracts a sparse structure of input SDP by factorizing the variable matrices. In this paper, we propose a new factorization based on the inverse of the variable matrix to enhance the performance of MC-PDIPM. We also use multithreaded parallel computing to deal with the major bottlenecks in MC-PDIPM. Numerical results show that the new factorization and multithreaded computing reduce the computation time for SDPs that have structural sparsity.

preprint2014arXiv

Flat-band ferromagnetism in the multilayer Lieb optical lattice

We theoretically study magnetic properties of two-component cold fermions in half-filled multilayer Lieb optical lattices, i.e., two, three, and several layers, using the dynamical mean-field theory. We clarify that the magnetic properties of this system become quite different depending on whether the number of layers is odd or even. In odd-number-th layers in an odd-number-layer system, finite magnetization emerges even with an infinitesimal interaction. This is a striking feature of the flatband ferromagnetic state in multilayer systems as a consequence of the Lieb theorem. In contrast, in even-number layers, magnetization develops from zero on a finite interaction. These different magnetic behaviours are triggered by the flat bands in the local density of states and become identical in the limit of the infinite-layer (i.e., three-dimensional) system. We also address how interlayer hopping affects the magnetization process. Further, we point out that layer magnetization, which is a population imbalance between up and down atoms on a layer, can be employed to detect the emergence of the flat-band ferromagnetic state without addressing sublattice magnetization.

preprint2014arXiv

High-fidelity cluster state generation for ultracold atoms in an optical lattice

We propose a method for generating high-fidelity multipartite spin-entanglement of ultracold atoms in an optical lattice in a short operation time with a scalable manner, which is suitable for measurement-based quantum computation. To perform the desired operations based on the perturbative spin-spin interactions, we propose to actively utilize the extra degrees of freedom (DOFs) usually neglected in the perturbative treatment but included in the Hubbard Hamiltonian of atoms, such as, (pseudo-)charge and orbital DOFs. Our method simultaneously achieves high fidelity, short operation time, and scalability by overcoming the following fundamental problem: enhancing the interaction strength for shortening operation time breaks the perturbative condition of the interaction and inevitably induces unwanted correlations among the spin and extra DOFs.

preprint2013arXiv

Equivariant comparison of quantum homogeneous spaces

We prove the deformation invariance of the quantum homogeneous spaces of the q-deformation of simply connected simple compact Lie groups over the Poisson-Lie quantum subgroups, in the equivariant KK-theory with respect to the translation action by maximal tori. This extends a result of Neshveyev-Tuset to the equivariant setting. As applications, we prove the ring isomorphism of the K-group of Gq with respect to the coproduct of C(Gq), and an analogue of the Borsuk-Ulam theorem for quantum spheres.

preprint2013arXiv

Operator algebra of foliations with projectively invariant transverse measure

We study the structure of operator algebras associated with the foliations which have projectively invariant measures. When a certain ergodicity condition on the measure preserving holonomies holds, the lack of holonomy invariant transverse measure can be established in terms of a cyclic cohomology class associated with the transverse fundamental cocycle and the modular automorphism group.

preprint2013arXiv

Twisted Cyclic Cohomology and Modular Fredholm Modules

Connes and Cuntz showed in [Comm. Math. Phys. 114 (1988), 515-526] that suitable cyclic cocycles can be represented as Chern characters of finitely summable semifinite Fredholm modules. We show an analogous result in twisted cyclic cohomology using Chern characters of modular Fredholm modules. We present examples of modular Fredholm modules arising from Podleś spheres and from $SU_q(2)$.

preprint2013arXiv

Twisting the q-deformations of compact semisimple Lie groups

Given a compact semisimple Lie group $G$ of rank $r$, and a parameter $q>0$, we can define new associativity morphisms in Rep(Gq) using a 3-cocycle $Φ$ on the dual of the center of G, thus getting a new tensor category Rep(Gq)$^Φ$. For a class of cocycles $Φ$ we construct compact quantum groups $G^τ_q$ with representation categories Rep(Gq)$^Φ$. The construction depends on the choice of an r-tuple $τ$ of elements in the center of G. In the simplest case of G=SU(2) and $τ=-1$, our construction produces Woronowicz's quantum group SU_{-q}(2) out of SUq(2). More generally, for G=SU(n), we get quantum group realizations of the Kazhdan-Wenzl categories.

preprint2012arXiv

Mott Transition and Spin Structures of Spin-1 Bosons in Two-Dimensional Optical Lattice at Unit Filling

We study the ground state properties of spin-1 bosons in a two-dimensional optical lattice, by applying a variational Monte Carlo method to the S=1 Bose-Hubbard model on a square lattice at unit filling. A doublon-holon binding factor introduced in the trial state provides a noticeable improvement in the variational energy over the conventional Gutzwiller wave function and allows us to deal effectively with the inter-site correlations of particle densities and spins. We systematically show how spin-dependent interactions modify the superfluid-Mott insulator transitions in the S=1 Bose-Hubbard model due to the interplay between the density and spin fluctuations of bosons. Furthermore, regarding the magnetic phases in the Mott region, the calculated spin structure factor elucidates the emergence of nematic and ferromagnetic spin orders for antiferromagnetic ($U_2>0$) and ferromagnetic ($U_2<0$) couplings, respectively.

preprint2012arXiv

The second-order reduced density matrix method and the two-dimensional Hubbard model

The second-order reduced density matrix method (the RDM method) has performed well in determining energies and properties of atomic and molecular systems, achieving coupled-cluster singles and doubles with perturbative triples (CC SD(T)) accuracy without using the wave-function. One question that arises is how well does the RDM method perform with the same conditions that result in CCSD(T) accuracy in the strong correlation limit. The simplest and a theoretically important model for strongly correlated electronic systems is the Hubbard model. In this paper, we establish the utility of the RDM method when employing the $P$, $Q$, $G$, $T1$ and $T2^\prime$ conditions in the two-dimension al Hubbard model case and we conduct a thorough study applying the $4\times 4$ Hubbard model employing a coefficients. Within the Hubbard Hamilt onian we found that even in the intermediate setting, where $U/t$ is between 4 and 10, the $P$, $Q$, $G$, $T1$ and $T2^\prime$ conditions re produced good ground state energies.

preprint2011arXiv

On subfactors arising from asymptotic representations of symmetric groups

We consider the infinite symmetric group and its infinite index subgroup given as the stabilizer subgroup of one element under the natural action on a countable set. This inclusion of discrete groups induces a hyperfinite subfactor for each finite factorial representation of the larger group. We compute subfactor invariants of this construction in terms of the Thoma parameter.

preprint2010arXiv

Thermodynamic properties of two-component fermionic atoms trapped in a two-dimensional optical lattice

We study the finite temperature properties of two-component fermionic atoms trapped in a two-dimensional optical lattice. We apply the self-energy functional approach to the two-dimensional Hubbard model with a harmonic trapping potential, and systematically investigate the thermodynamic properties of this system. We find that entropy and grand potential provide evidence of a crossover between the Mott insulating and metallic phases at certain temperatures. In addition, we find that entropy exhibits a cusp-like anomaly at lower temperatures, suggesting a second or higher order antiferromagnetic transition. We estimate the antiferromagnetic transition temperatures, and clarify how the trapping potential affects this magnetic transition.

preprint2010arXiv

Time-of-flight imaging method to observe signatures of antiferromagnetically ordered states of fermionic atoms in an optical lattice

We propose a simple method to detect the antiferromagnetic (AF) state of fermionic atoms in an optical lattice by combining a time-of-flight (TOF) imaging method and a Feshbach resonance. In this scheme, the nontrivial dynamics of fermionic atoms during the imaging process works as a probe with respect to the breaking of the translational symmetry in the AF state. Precise numerical simulations demonstrate that the characteristic oscillatory dynamics induced by the scattering process that transfers an AF ordering vector appears in TOF images, which can be easily observed experimentally.

preprint2006arXiv

Fermionic atoms trapped in one-dimensional optical superlattice with harmonic confinement

We study the ground-state properties of spin-1/2 fermionic atoms confined in a one-dimensional optical superlattice with harmonic confinement by using the density matrix renormalization group method. For this purpose, we consider an ionic Hubbard model that has superlattice potentials with 2-site periodicity. We find that several different types of insulating regimes coexist even if the number of atoms at each site is not an integer, but its average within the unit cell is an integer or half integer. This is contrasted to the coexisting phase of the metallic and Mott-insulating regimes known for the ordinary Hubbard model in an optical lattice. The phase characteristics are elucidated by investigating the profiles of the atom density, the local density/spin fluctuations, the double occupation probability and the spin correlations in detail.