Source author record

Seiji Miyashita

Seiji Miyashita 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

44works
14topics
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

44 published item(s)

preprint2022arXiv

Finite-temperature properties of extended Nagaoka ferromagnetism

We study finite-temperature properties of a Hubbard model including sites of a particle bath which was proposed as a microscopic model to show itinerant ferromagnetism at finite electron density. We use direct numerical methods, such as exact diagonalization and random vector methods. The temperature dependence of quantities is surveyed in the full range of the temperature. We find that the specific heat has several peaks, which correspond to ordering processes in different energy scales. In particular, magnetic order appears at very low temperature. Depending on the chemical potential of the particle bath and the Coulomb interaction, the system exhibits an itinerant ferromagnetic state or an antiferromagnetic state of the Mott insulator. Microscopically the competition between these two types of orderings causes a peculiar ordering process of local spin correlations. Some local ferromagnetic correlations are found to be robust, which indicates that the ferromagnetic correlation originates from the motion of itinerant electrons in a short-range cluster.

preprint2020arXiv

Nontrivial temperature dependence of ferromagnetic resonance frequency for spin reorientation transitions

We find unusual temperature dependence of the ferromagnetic resonance (FMR) frequency $f_{\rm R}$ for the spin-reorientation (SR) transition, in which the easy axis changes depending on temperature, observed in the Nd permanent magnet, Nd$_2$Fe$_{14}$B: $f_{\rm R} \sim 0$ below the SR transition temperature ($T_{\rm R}$), drastic increase of $f_{\rm R}$ around $T_{\rm R}$, and decrease from a peak at higher temperatures. It is nontrivial that the SR transition causes the unusual behavior of the FMR frequency in a wide temperature region. We show the mechanism of the temperature dependence by theoretical and computational analyses. We derive a general relation between $f_{\rm R}$ and magnetizations to help the understanding of the mechanism, and clarify that the fluctuation of the transverse magnetization is a key ingredient for the resonance in all temperature regions.

preprint2020arXiv

The effect of the surface magnetic anisotropy of the neodymium atoms on the coercivity in the neodymium permanent magnet

The Nd permanent magnet (Nd$_{2}$Fe$_{14}$B) is an indispensable material used in modern energy conversion devices. The realization of high coercivity at finite temperatures is a burning issue. One of the important ingredients for controlling the coercive force is the surface property of magnetic grains. It has been reported by first-principles studies that the Nd atoms in the first (001) surface layer facing the vacuum have in-plane anisotropy perpendicular to the $c$ axis, which may decrease the coercivity. Focusing on the surface anisotropy effect on the coercivity, we examine the coercivity at zero and finite temperatures by using an atomistic model reflecting the lattice structure of the Nd magnet with a stochastic Landau-Lifshitz-Gilbert equation method. We study general three cases, in which the Nd atoms in surface layers have (1) no anisotropy, (2) in-plane anisotropy, and (3) reinforced anisotropy for two types of surfaces, (001) and (100) surfaces. We find that in contrast to the zero-temperature case, due to the thermal fluctuation effect, the modification of only the first surface layer has little effect on the coercivity at finite temperatures. However, the modification of a few layers results in significant effects. We discuss the details of the dependence of the coercivity on temperature, type of surface, and modified layer depth, and also the features of domain growth in magnetization reversal.

preprint2019arXiv

Dynamical phase transition in Floquet optical bistable systems: An approach from finite-size quantum systems

We study a dynamical phase transition in optical bistable systems subject to a time-periodic driving field. The phase transition occurs in the structure of limit cycle as a function of the frequency of the driving field. In the thermodynamic limit, a single limit cycle is divided into two separated limit cycles at the transition point. In finite-size systems, however, there is always a single limit cycle due to the quantum tunneling effect. We use a Floquet dissipative map, which is a time-evolution operator over one period in a dynamics given by a quantum master equation, and discuss the decay rate of relaxation dynamics into the limit cycle based on the dominant eigenvalue of the map. We found that the decay rate exhibits qualitatively different system-size dependence before and after the phase transition, and it shows a finite-size scaling of spinodal phenomena around the transition point. The present work provides a systematic way of studying dynamical phase transition observed in time-periodically driven open systems in terms of the Floquet dissipative map.

preprint2016arXiv

Effective Floquet-Gibbs states for dissipative quantum systems

A periodically driven quantum system, when coupled to a heat bath, relaxes to a non-equilibrium asymptotic state. In the general situation, the retrieval of this asymptotic state presents a rather non-trivial task. It was recently shown that in the limit of an infinitesimal coupling, using so-called rotating wave approximation (RWA), and under strict conditions imposed on the time-dependent system Hamiltonian, the asymptotic state can attain the Gibbs form. A Floquet-Gibbs state is characterized by a density matrix which is diagonal in the Floquet basis of the system Hamiltonian with the diagonal elements obeying a Gibbs distribution, being parametrized by the corresponding Floquet quasi-energies. Addressing the non-adiabatic driving regime, upon using the Magnus expansion, we employ the concept of a corresponding effective Floquet Hamiltonian. In doing so we go beyond the conventionally used RWA and demonstrate that the idea of Floquet-Gibbs states can be extended to the realistic case of a weak, although finite system-bath coupling, herein termed effective Floquet-Gibbs states.

preprint2016arXiv

Monte Carlo analysis for finite temperature magnetism of Nd$_2$Fe$_{14}$B permanent magnet

We investigate the effects of magnetic inhomogeneities and thermal fluctuations on the magnetic properties of a rare earth intermetallic compound, Nd$_2$Fe$_{14}$B. The constrained Monte Carlo method is applied to a Nd$_2$Fe$_{14}$B bulk system to realize the experimentally observed spin reorientation and magnetic anisotropy constants $K^{\rm A}_m (m=1, 2, 4)$ at finite temperatures. % Subsequently, it is found that the temperature dependence of $K^{\rm A}_1$ deviates from the Callen--Callen law, $K^{\rm A}_1(T) \propto M(T)^3$, even above room temperature, $T_{\rm R}\sim 300\rm\, K$, when the Fe (Nd) anisotropy terms are removed to leave only the Nd (Fe) anisotropy terms. This is because the exchange couplings between Nd moments and Fe spins are much smaller than those between Fe spins. % It is also found that the exponent $n$ in the external magnetic field $H_{\rm ext}$ response of barrier height $\mathcal{F}_{\rm B}=\mathcal{F}_{\rm B}^0(1-H_{\rm ext}/H_0)^n$ is less than $2$ in the low-temperature region below $T_{\rm R}$, whereas $n$ approaches $2$ when $T>T_{\rm R}$, indicating the presence of Stoner--Wohlfarth-type magnetization rotation. This reflects the fact that the magnetic anisotropy is mainly governed by the $K^{\rm A}_1$ term in the $T>T_{\rm R}$ region.

preprint2016arXiv

Two-step asymptotics of scaled Dunkl processes

Dunkl processes are generalizations of Brownian motion obtained by using the differential-difference operators known as Dunkl operators as a replacement of spatial partial derivatives in the heat equation. Special cases of these processes include Dyson's Brownian motion model and the Wishart-Laguerre eigenvalue processes, which are well-known in random matrix theory. It is known that the dynamics of Dunkl processes is obtained by transforming the heat kernel using Dunkl's intertwining operator. It is also known that, under an appropriate scaling, their distribution function converges to a steady-state distribution which depends only on the coupling parameter $β$ as the process time $t$ tends to infinity. We study scaled Dunkl processes starting from an arbitrary initial distribution, and we derive expressions for the intertwining operator in order to calculate the asymptotics of the distribution function in two limiting situations. In the first one, $β$ is fixed and $t$ tends to infinity (approach to the steady state), and in the second one, $t$ is fixed and $β$ tends to infinity (strong-coupling limit). We obtain the deviations from the limiting distributions in both of the above situations, and we find that they are caused by the two different mechanisms which drive the process, namely, the drift and exchange mechanisms. We find that the deviation due to the drift mechanism decays as $t^{-1}$, while the deviation due to the exchange mechanism decays as $t^{-1/2}$.

preprint2015arXiv

Condition for emergence of the Floquet-Gibbs state in periodically driven open systems

We study probability distribution of a steady state of a periodically driven system coupled to a thermal bath by using a quantum master equation in the weak coupling limit. It is proved that, even when the external field is strong, the probability distribution is independent of the detailed nature of the thermal bath under the following conditions: (i) the Hamiltonian of the relevant system is bounded and the period of the driving field is short, (ii) the Hamiltonians for the driving field at different times commute, and (iii) the Hamiltonians of the driving field and of the interaction between the relevant system and the thermal bath commute. It is shown that the steady state is described by the Gibbs distribution of the Floquet states of the relevant system at the temperature of the thermal bath.

preprint2015arXiv

Doping control of realization of an extended Nagaoka ferromagnetic state from the Mott state

Inspired by the Nagaoka ferromagnetism, we propose an itinerant model to study the transition between the Mott singlet state and a ferromagnetic state by emulating a doping process in finite lattices. In the Nagaoka ferromagnetism, the total spin of the system takes the maximum value when an electron is removed from the half-filled system. To incorporate a procedure of the electron removal, our model contains extra sites as a reservoir of electrons, and the chemical potential of the reservoir controls the distribution of electrons. As a function of the chemical potential, the system exhibits ground-state phase transitions among various values of the total spin, including a saturated ferromagnetic state due to the Nagaoka mechanism at finite hole density. We discuss the nature of the ferromagnetism by measuring various physical quantities, such as the distribution of electrons, the spin correlation functions, the magnetization process in the magnetic field, and also the entanglement entropy.

preprint2015arXiv

Effect of the short-range interaction on critical phenomena in elastic interaction systems

The elastic interaction, induced by the lattice distortion due to the difference of the molecular size, causes an effective long-range interaction. In spin-crossover (SC) compounds, local bistable states, i.e., high-spin and low-spin states have different molecular sizes, and the elastic interaction is important. In bipartite lattices, e.g., the square lattice, the ground state can be two types of phases: ferromagnetic-like (FM) and antiferromagnetic-like (AFM) phases. In systems like SC compounds, the former phase consists of all small or large molecules, and the latter phase has the configuration of alternating small and large molecules. In fact, both cases are observed in SC systems. In this paper, we have studied the effect of the short-range interaction in the elastic system on the properties of those order-disorder phase transitions. We have obtained a phase diagram in the coordinates of the temperature and the strength of the short-range interaction, including the metastable structures. We show that effects of the short-range interaction are essentially different for FM and AFM phase transitions. In the FM transition, the long-range interaction of elasticity is relevant, and the system exhibits a phase transition in the mean-filed universality class. In this case, the long-range interaction strongly enhances the FM order, and it works cooperatively with the short-range interaction. In contrast, in the AFM transition, the elastic interaction slightly enhances the AFM order, but essentially it does not contribute to the ordering, and the system shows a transition in the Ising universality class. We have found that in the border region between FM and AFM phases, the AFM phase has an advantage at finite temperatures. We discuss the critical properties of two-step SC transitions with comparison between the elastic interaction model and conventional SC models (Ising-like models).

preprint2015arXiv

Numerical study on ESR by making use of Wiener-Khinchin relation in time domain

To evaluate ESR spectrum at finite temperatures for specified spatial configurations of spins is very important issue to study quantum spin systems. Although a direct numerical estimation of the Kubo formula provides exact data, the application is limited to small size of the system because of the restriction of the computer capacity. The method of the Fourier transform of the autocorrelation function improved the restriction. As an extension of the method, we propose a new method for numerical calculation of the ESR spectrum from the time evolution of the magnetization by making use of the Wiener-Khinchin theorem.

preprint2015arXiv

Quantum Stoner-Wohlfarth model

The quantum mechanical counterpart of the famous Stoner-Wohlfarth model -- an easy-axis magnet in a tilted magnetic field -- is studied theoretically and through simulations, as a function of the spin-size $S$ in a sweeping longitudinal field. Beyond the classical Stoner-Wohlfarth transition, the sweeping field-induced adiabatic change of states slows down as $S$ increases, leading to a dynamical quantum phase transition. This result is described as a critical phenomenon associated with Landau-Zener tunneling gaps at metastable quasi-avoided crossings. Furthermore, a beating of the magnetization is discovered after the Stoner-Wohlfarth transition. The period of the beating, obtained analytically, arises from a new type of quantum phase factor.

preprint2015arXiv

Realization of the thermal equilibrium in inhomogeneous magnetic systems by the Landau-Lifshitz-Gilbert equation with stochastic noise, and its dynamical aspects

It is crucially important to investigate effects of temperature on magnetic properties such as critical phenomena, nucleation, pinning, domain wall motion, coercivity, etc. The Landau-Lifshitz-Gilbert (LLG) equation has been applied extensively to study dynamics of magnetic properties. Approaches of Langevin noises have been developed to introduce the temperature effect into the LLG equation. To have the thermal equilibrium state (canonical distribution) as the steady state, the system parameters must satisfy some condition known as the fluctuation-dissipation relation. In inhomogeneous magnetic systems in which spin magnitudes are different at sites, the condition requires that the ratio between the amplitude of the random noise and the damping parameter depends on the magnitude of the magnetic moment at each site. Focused on inhomogeneous magnetic systems, we systematically showed agreement between the stationary state of the stochastic LLG equation and the corresponding equilibrium state obtained by Monte Carlo simulations in various magnetic systems including dipole-dipole interactions. We demonstrated how violations of the condition result in deviations from the true equilibrium state. We also studied the characteristic features of the dynamics depending on the choice of the parameter set. All the parameter sets satisfying the condition realize the same stationary state (equilibrium state). In contrast, different choices of parameter set cause seriously different relaxation processes. We show two relaxation types, i.e., magnetization reversals with uniform rotation and with nucleation.

preprint2015arXiv

Termination of the Berezinskii-Kosterlitz-Thouless phase with a new critical universality in spin-crossover systems

Two dimensional systems with U(1) symmetry exhibit a peculiar phase, i.e., the Berezinskii-Kosterlitz-Thouless (BKT) phase. In particular situations, the BKT phase exists as an intermediate temperature phase. There have been scenarios for the phase transitions at the two endpoints of the intermediate BKT phase, i.e., the phase transition at the low-temperature endpoint is a BKT transition and that at the high-temperature endpoint is either a BKT transition or a first-order transition. The present study gives a novel scenario, i.e., a second-order transition with a new critical universality and a BKT transition. We found that this new phase transition is realized in spin-crossover systems on a triangular lattice with an antiferromagnetic short-range interaction. At the low-temperature transition the elastic interaction plays as a ferromagnetic infinite-range interaction and encourages the breaking of $Z_2$ symmetry between high-spin rich and low-spin rich states.

preprint2014arXiv

Monte Carlo Studies of the Ising Antiferromagnet with a Ferromagnetic Mean-field Term

The unusual thermodynamic properties of the Ising antiferromagnet supplemented with a ferromagnetic, mean-field term are outlined. This simple model is inspired by more realistic models of spin-crossover materials. The phase diagram is estimated using Metropolis Monte Carlo methods, and differences with preliminary Wang-Landau Monte Carlo results for small systems are noted.

preprint2014arXiv

Novel symmetry-broken phase in a driven cavity system in the thermodynamic limit

We study non-equilibrium stationary states of a cavity system consisting of many atoms interacting with a quantized cavity field mode, under a driving field in a dissipative environment. We derive a quantum master equation which is suitable for treating systems with a strong driving field and a strong atom-photon interaction. We do this by making use of the fact that the mean-field dynamics are exact in the thermodynamic limit thanks to a uniform coupling between atoms and photons. We find ordered states with symmetry-broken components of the photon field and atomic excitation driven by the external field. The mechanism by which these ordered states arise is discussed from the viewpoint of the quantum interference effect.

preprint2014arXiv

Quantum Annealing for Variational Bayes Inference

This paper presents studies on a deterministic annealing algorithm based on quantum annealing for variational Bayes (QAVB) inference, which can be seen as an extension of the simulated annealing for variational Bayes (SAVB) inference. QAVB is as easy as SAVB to implement. Experiments revealed QAVB finds a better local optimum than SAVB in terms of the variational free energy in latent Dirichlet allocation (LDA).

preprint2014arXiv

Two limiting regimes of interacting Bessel processes

We consider the interacting Bessel processes, a family of multiple-particle systems in one dimension where particles evolve as individual Bessel processes and repel each other via a log-potential. We consider two limiting regimes for this family on its two main parameters: the inverse temperature beta and the Bessel index nu. We obtain the time-scaled steady-state distributions of the processes for the cases where beta or nu are large but finite. In particular, for large beta we show that the steady-state distribution of the system corresponds to the eigenvalue distribution of the beta-Laguerre ensembles of random matrices. We also estimate the relaxation time to the steady state in both cases. We find that in the freezing regime beta->infinity, the scaled final positions of the particles are locked at the square root of the zeroes of the Laguerre polynomial of parameter nu-1/2 for any initial configuration, while in the regime nu->infinity, we prove that the scaled final positions of the particles converge to a single point. In order to obtain our results, we use the theory of Dunkl operators, in particular the intertwining operator of type B. We derive a previously unknown expression for this operator and study its behaviour in both limiting regimes. By using these limiting forms of the intertwining operator, we derive the steady-state distributions, the estimations of the relaxation times and the limiting behaviour of the processes.

preprint2013arXiv

Asymptotic forms and scaling properties of the relaxation time near threshold points in spinodal-type dynamical phase transitions

We study critical properties of the relaxation time at a threshold point in switching processes between bistable states under change of external fields. In particular, we investigate the relaxation processes near the spinodal point of the infinitely long-range interaction model (the Husimi-Temperley model) by analyzing the scaling properties of the corresponding Fokker-Planck equation. We also confirm the obtained scaling relations by direct numerical solution of the original master equation, and by kinetic Monte Carlo simulation of the stochastic decay process. In particular, we study the asymptotic forms of the divergence of the relaxation time near the spinodal point, and reexamine its scaling properties. We further extend the analysis to transient critical phenomena such as a threshold behavior with diverging switching time under a general external driving perturbation. This models photo-excitation processes in spin-crossover materials. In the ongoing development of nano-size fabrication, such size-dependence of switching processes should be an important issue, and the properties obtained here will be applicable to a wide range of physical processes.

preprint2013arXiv

Complex Magnetic Phase Diagram of a Geometrically Frustrated Sm Lattice: SmPd2Al3 case

Magnetism in SmPd2Al3 was investigated on a single crystal by magnetometry and neutron diffraction. SmPd2Al3 represents a distinctive example of the Sm magnetism exhibiting complex magnetic behavior at low temperatures with four consecutive magnetic phase transitions at 3.4, 3.9, 4.3 and 12.5 K. The rich magnetic phase diagram of this compound reflects the specific features of the Sm3+ ion, namely the energy nearness of the ground-state multiplet J = 5/2 and the first excited multiplet J = 7/2 in conjunction with strong crystal field influence. Consequently, a significantly reduced Sm magnetic moment in comparison with the theoretical Sm3+ free-ion value is observed. Despite the strong neutron absorption by natural samarium and the small Sm magnetic moment (~ 0.2 μB) we have successfully determined the magnetic k-vector (1/3, 1/3, 0) of the phase existing in the temperature interval 12.5 - 4.3 K. This observation classifies the SmPd2Al3 compound as a magnetically frustrated system. The complex magnetic behavior of this material is further illustrated by kinetic effects of the magnetization inducing rather complicated magnetic structure with various metastable states.

preprint2013arXiv

Quantum Annealing for Dirichlet Process Mixture Models with Applications to Network Clustering

We developed a new quantum annealing (QA) algorithm for Dirichlet process mixture (DPM) models based on the Chinese restaurant process (CRP). QA is a parallelized extension of simulated annealing (SA), i.e., it is a parallel stochastic optimization technique. Existing approaches [Kurihara et al. UAI2009, Sato et al. UAI2009] and cannot be applied to the CRP because their QA framework is formulated using a fixed number of mixture components. The proposed QA algorithm can handle an unfixed number of classes in mixture models. We applied QA to a DPM model for clustering vertices in a network where a CRP seating arrangement indicates a network partition. A multi core processor was used for running QA in experiments, the results of which show that QA is better than SA, Markov chain Monte Carlo inference, and beam search at finding a maximum a posteriori estimation of a seating arrangement in the CRP. Since our QA algorithm is as easy as to implement the SA algorithm, it is suitable for a wide range of applications.

preprint2013arXiv

Quantum Decoherence Scaling with Bath Size: Importance of Dynamics, Connectivity, and Randomness

The decoherence of a quantum system $S$ coupled to a quantum environment $E$ is considered. For states chosen uniformly at random from the unit hypersphere in the Hilbert space of the closed system $S+E$ we derive a scaling relationship for the sum of the off-diagonal elements of the reduced density matrix of $S$ as a function of the size $D_E$ of the Hilbert space of $E$. This sum decreases as $1/\sqrt{D_E}$ as long as $D_E\gg 1$. This scaling prediction is tested by performing large-scale simulations which solve the time-dependent Schr{ö}dinger equation for a ring of spin-1/2 particles, four of them belonging to $S$ and the others to $E$. Provided that the time evolution drives the whole system from the initial state toward a state which has similar properties as states belonging to the class of quantum states for which we derived the scaling relationship, the scaling prediction holds. For systems which do not exhibit this feature, it is shown that increasing the complexity (in terms of connections) of the environment or introducing a small amount of randomness in the interactions in the environment suffices to observe the predicted scaling behavior.

preprint2012arXiv

Calogero-Moser Systems as a Diffusion-Scaling Transform of Dunkl Processes on the Line

The Calogero-Moser systems are a series of interacting particle systems on one dimension that are both classically and quantum-mechanically integrable. Their integrability has been established through the use of Dunkl operators (a series of differential-difference operators that depend on the choice of an abstract set of vectors, or root system). At the same time, Dunkl operators are used to define a family of stochastic processes called Dunkl processes. We showed in a previous paper that when the coupling constant of interaction of the symmetric Dunkl process on the root system A(N-1) goes to infinity (the freezing regime), its final configuration is proportional to the roots of the Hermite polynomials. It is also known that the positions of the particles of the Calogero-Moser system with particle exchange become fixed at the roots of the Hermite polynomials in the freezing regime. Although both systems present a freezing behaviour that depends on the roots of the Hermite polynomials, the reason for this similarity has been an open problem until now. In the present work, we introduce a new type of similarity transformation called the diffusion-scaling transformation, in which a new space variable is given by a diffusion-scaling variable constructed using the original space and time variables. We prove that the abstract Calogero-Moser system on an arbitrary root system is a diffusion-scaling transform of the Dunkl process on the same root system. With this, we prove that the similar freezing behaviour of the two systems on A(N-1) stems from their similar mathematical structure.

preprint2012arXiv

Critical temperature and correlation length of an elastic interaction model for spin-crossover materials

It has previously been pointed out that the coexistence of infinite-range and short-range interactions causes a system to have a phase transition of the mean-field universality class, in which the cluster size is finite even at the critical point. In the present paper, we study this property in a model of bistable molecules, whose size changes depending on the bistable states. The molecules can move in space, interacting via an elastic interaction. It is known that due to the different sizes, an effective long-range interaction between the spins appears, and thus this model has a mean-field type of phase transition. It is found that the scaling properties of the shift of the critical temperature from the pure short-range limit in the model with infinite-range and short-range interactions hold also in the present model, regarding the ratio of the size of the two states as a control parameter for the strength of the long-range interaction. By studying the structure factor, it is shown that the dependence of the cluster size at the critical temperature also shows the same scaling properties as a previously studied model with both infinite-range and short-range interactions. We therefore conclude that these scaling relations hold universally in hybrid models with both short-range and weak long-range interactions.

preprint2012arXiv

Equilibration and Thermalization of Classical Systems

It is demonstrated that the canonical distribution for a subsystem of a closed system follows directly from the solution of the time-reversible Newtonian equation of motion in which the total energy is strictly conserved. It is shown that this conclusion holds for both integrable or nonintegrable systems even though the whole system may contain as little as a few thousand particles. In other words, we demonstrate that the canonical distribution holds for subsystems of experimentally relevant sizes and observation times.

preprint2012arXiv

Interacting Particles on the Line and Dunkl Intertwining Operator of Type A: Application to the Freezing Regime

We consider a one-dimensional system of Brownian particles that repel each other through a logarithmic potential. We study two formulations for the system and the relation between them. The first, Dyson's Brownian motion model, has an interaction coupling constant determined by the parameter beta > 0. When beta = 1,2 and 4, this model can be regarded as a stochastic realization of the eigenvalue statistics of Gaussian random matrices. The second system comes from Dunkl processes, which are defined using differential-difference operators (Dunkl operators) associated with finite abstract vector sets called root systems. When the type-A root system is specified, Dunkl processes constitute a one-parameter system similar to Dyson's model, with the difference that its particles interchange positions spontaneously. We prove that the type-A Dunkl processes with parameter k > 0 starting from any symmetric initial configuration are equivalent to Dyson's model with the parameter beta = 2k. We focus on the intertwining operators, since they play a central role in the mathematical theory of Dunkl operators, but their general closed form is not yet known. Using the equivalence between symmetric Dunkl processes and Dyson's model, we extract the effect of the intertwining operator of type A on symmetric polynomials from these processes' transition probability densities. In the strong coupling limit, the intertwining operator maps all symmetric polynomials onto a function of the sum of their variables. In this limit, Dyson's model freezes, and it becomes a deterministic process with a final configuration proportional to the roots of the Hermite polynomials multiplied by the square root of the process time, while being independent of the initial configuration.

preprint2012arXiv

Macroscopic nucleation phenomena in continuum media with long-range interactions

Nucleation, commonly associated with discontinuous transformations between metastable and stable phases, is crucial in fields as diverse as atmospheric science and nanoscale electronics. Traditionally, it is considered a microscopic process (at most nano-meter), implying the formation of a microscopic nucleus of the stable phase. Here we show for the first time, that considering long-range interactions mediated by elastic distortions, nucleation can be a macroscopic process, with the size of the critical nucleus proportional to the total system size. This provides a new concept of "macroscopic barrier-crossing nucleation". We demonstrate the effect in molecular dynamics simulations of a model spin-crossover system with two molecular states of different sizes, causing elastic distortions.

preprint2011arXiv

Crossover between a Short-range and a Long-range Ising model

Recently, it has been found that an effective long-range interaction is realized among local bistable variables (spins) in systems where the elastic interaction causes ordering of the spins. In such systems, generally we expect both long-range and short-range interactions to exist. In the short-range Ising model, the correlation length diverges at the critical point. In contrast, in the long-range interacting model the spin configuration is always uniform and the correlation length is zero. As long as a system has non-zero long-range interactions, it shows criticality in the mean-field universality class, and the spin configuration is uniform beyond a certain scale. Here we study the crossover from the pure short-range interacting model to the long-range interacting model. We investigate the infinite-range model (Husimi-Temperley model) as a prototype of this competition, and we study how the critical temperature changes as a function of the strength of the long-range interaction. This model can also be interpreted as an approximation for the Ising model on a small-world network. We derive a formula for the critical temperature as a function of the strength of the long-range interaction. We also propose a scaling form for the spin correlation length at the critical point, which is finite as long as the long-range interaction is included, though it diverges in the limit of the pure short-range model. These properties are confirmed by extensive Monte Carlo simulations.

preprint2010arXiv

Approach to equilbrium in nano-scale systems at finite temperatur

We study the time evolution of the reduced density matrix of a system of spin-1/2 particles interacting with an environment of spin-1/2 particles. The initial state of the composite system is taken to be a product state of a pure state of the system and a pure state of the environment. The latter pure state is prepared such that it represents the environment at a given finite temperature in the canonical ensemble. The state of the composite system evolves according to the time-dependent Schr{ö}dinger equation, the interaction creating entanglement between the system and the environment. It is shown that independent of the strength of the interaction and the initial temperature of the environment, all the eigenvalues of the reduced density matrix converge to their stationary values, implying that also the entropy of the system relaxes to a stationary value. We demonstrate that the difference between the canonical density matrix and the reduced density matrix in the stationary state increases as the initial temperature of the environment decreases. As our numerical simulations are necessarily restricted to a modest number of spin-1/2 particles ($<36$), but do not rely on time-averaging of observables nor on the assumption that the coupling between system and environment is weak, they suggest that the stationary state of the system directly follows from the time evolution of a pure state of the composite system, even if the size of the latter cannot be regarded as being close to the thermodynamic limit.

preprint2010arXiv

Corpuscular model of two-beam interference and double-slit experiments with single photons

We introduce an event-based corpuscular simulation model that reproduces the wave mechanical results of single-photon double slit and two-beam interference experiments and (of a one-to-one copy of an experimental realization) of a single-photon interference experiment with a Fresnel biprism. The simulation comprises models that capture the essential features of the apparatuses used in the experiment, including the single-photon detectors recording individual detector clicks. We demonstrate that incorporating in the detector model, simple and minimalistic processes mimicking the memory and threshold behavior of single-photon detectors is sufficient to produce multipath interference patterns. These multipath interference patterns are built up by individual particles taking one single path to the detector where they arrive one-by-one. The particles in our model are not corpuscular in the standard, classical physics sense in that they are information carriers that exchange information with the apparatuses of the experimental set-up. The interference pattern is the final, collective outcome of the information exchanges of many particles with these apparatuses. The interference patterns are produced without making reference to the solution of a wave equation and without introducing signalling or non-local interactions between the particles or between different detection points on the detector screen.

preprint2010arXiv

Non-monotonic Dynamics in Frustrated Ising Model with Time-Dependent Transverse Field

We study how the degree of ordering depends on the strength of the thermal and quantum fluctuations in frustrated systems by investigating the correlation function of the order parameter. Concretely, we compare the equilibrium spin correlation function in a frustrated lattice which exhibits a non-monotonic temperature dependence (reentrant type dependence) with that in the ground state as a function of the transverse field that causes the quantum fluctuation. We find the correlation function in the ground state also shows a non-monotonic dependence on the strength of the transverse field. We also study the real-time dynamics of the spin correlation function under a time-dependent field. After sudden decrease of the temperature, we found non-monotonic changes of the correlation function reflecting the static temperature dependence, which indicates that an effective temperature of the system changes gradually. For the quantum system, we study the dependence of changes of the correlation function on the sweeping speed of the transverse field. Contrary to the classical case, the correlation function varies little in a rapid change of the field, though it shows a non-monotonic change when we sweep the field slowly.

preprint2010arXiv

Quantum Field Induced Orderings in Fully Frustrated Ising Spin Systems

We study ordering mechanisms which are induced by the quantum fluctuation in fully frustrated Ising spin systems. Since there are many degenerated states in frustrated systems, "order by thermal disorder" often takes place due to a kind of entropy effect. To consider "order by quantum disorder" in fully frustrated Ising spin systems, we apply transverse field as quantum fluctuation. There exists a ferromagnetic correlation in each sublattice. The sublattice correlation at zero temperature is enlarged due to transverse field. The quantum fluctuation enhances the solid order at zero temperatures. This is an example of quantum field induced ordering in fully frustrated systems. We also study a case in which the transverse field induces a reentrant behavior as another type of order by quantum disorder, and compare correspondent cases in the classical systems.

preprint2009arXiv

Comparison among Various Expressions of Complex Admittance for Quantum System in Contact with Heat Reservoir

Relation among various expressions of the complex admittance for quantum systems in contact with heat reservoir is studied. Exact expressions of the complex admittance are derived in various types of formulations of equations of motion under contact with heat reservoir. Namely, the complex admittance is studied in the relaxation method and the external-field method. In the former method, the admittance is calculated using the Kubo formula for quantum systems in contact with heat reservoir in no external driving fields, while in the latter method the admittance is directly calculated from equations of motion with external driving terms. In each method, two types of equation of motions are considered, i.e., the time-convolution (TC) equation and time-convolutionless (TCL) equation. That is, the full of the four cases are studied. It is turned out that the expression of the complex admittance obtained by using the relaxation method with the TC equation exactly coincides with that obtained by using the external-field method with the TC equation, while other two methods give different forms. It is also explicitly demonstrated that all the expressions of the complex admittance coincide with each other in the lowest Born approximation for the systemreservoir interaction. The formulae necessary for the higher order expansions in powers of the system-reservoir interaction are derived, and also the expressions of the admittance in the n-th order approximation are given. To characterize the TC and TCL methods, we study the expressions of the admittances of two exactly solvable models. Each exact form of admittance is compared with the results of the two methods in the lowest Born approximation. It is found that depending on the model, either of TC and TCL would be the better method.

preprint2009arXiv

Intrinsic effects of the boundary condition on switching processes in spin-crossover solids

We investigated domain growth in switching processes between the low-spin and high-spin phases in thermally induced hysteresis loops of spin-crossover (SC) solids. Elastic interactions among the molecules induce effective long-range interactions, and thus the boundary condition plays a significant role in the dynamics. In contrast to SC systems with periodic boundary conditions, where uniform configurations are maintained during the switching process, we found that domain structures appear with open boundary conditions. Unlike Ising-like models with short-range interactions, domains always grow from the corners of the system. The present clustering mechanism provides an insight into the switching dynamics of SC solids, in particular, in nano-scale systems.

preprint2009arXiv

Master equation approach to line shape in dissipative systems

We propose a formulation to obtain the line shape of a magnetic response with dissipative effects that directly reflects the nature of the environment. Making use of the fact that the time evolution of a response function is described by the same equation as the reduced density operator, we formulate a full description of the complex susceptibility. We describe the dynamics using the equation of motion for the reduced density operator, including the term for the initial correlation between the system and a thermal bath. In this formalism, we treat the full description of non-Markovian dynamics, including the initial correlation. We present an explicit and compact formula up to the second order of cumulants, which can be applied in a straightforward way to multiple spin systems. We also take into account the frequency shift by the system-bath interaction. We study the dependence of the line shape on the type of interaction between the system and the thermal bath. We demonstrate that the present formalism is a powerful tool for investigating various kinds of systems, and we show how it is applied to spin systems, including those with up to three spins. We distinguish the contributions of the initial correlation and the frequency shift, and make clear the role of each contribution in the Ohmic coupling spectral function. As examples of applications to multispin systems, we obtain the dependence of the line shape on the spatial orientation in relation to the direction of the static field (Nagata-Tazuke effect), including the effects of the thermal environment, in a two-spin system, along with the dependence on the arrangement of a triangle in a three-spin system.

preprint2009arXiv

Nagaoka ferromagnetism in large-spin systems -Fermion and Boson systems--

We study magnetic properties of itinerant quantum magnetic particles described by a generalized Hubbard model with large spin ($S>1/2$) which may be realized in optical lattices of laser-cooled atom systems. In fermion systems (half-integer spins), an extended form of Nagaoka ferromagnetism may be realized. However, as novel aspects of the large spin cases, we found that the condition on the lattice connectivity is more stringent than in the case of $S=1/2$ particles, and that the system shows a peculiar degenerate structure of the ground statein which the ferromagnetic state is included. In contrast, it turns out that the ground state of itinerant bosonic systems (integer spins) has a degenerate structure similar to that of fermion system with $S>1/2$ regardless of the shape, connectivity or filling of the lattice, and that the state with the maximum total spin is always one of the ground states.

preprint2009arXiv

Threshold Phenomena under Photo Excitation of Spin-crossover Materials with Cooperativity due to Elastic Interactions

Photo-induced switching from the low-spin state to the high-spin state is studied in a model of spin-crossover materials, in which long-range interactions are induced by elastic distortions due to different molecular sizes the two spin states. At a threshold value of the light intensity we observe nonequilibrium critical behavior corresponding to a mean-field spinodal point. Finite-size scaling of the divergence of the relaxation time is revealed by analysis of kinetic Monte Carlo simulations.

preprint2005arXiv

Dynamical properties of temperature chaos and memory effect

In the frustrated interaction systems, the nature of ordered configuration can be intrinsically temperature dependent. There, the idea of effective coupling of decorated and frustrated bond plays an important role. The idea of effective coupling is verified in the study of equilibrium properties. We will study how the effective coupling is realized in the dynamical processes, e.g. the relaxation time after a change of temperature. We also study the memory phenomena which have been found in the spin glass.

preprint2005arXiv

Temperature dependence of ESR intensity for the nanoscale molecular magnet V15

The electron spin resonance (ESR) of nanoscale molecular magnet ${\rm V}_{15}$ is studied. Since the Hamiltonian of ${\rm V}_{15}$ has a large Hilbert space and numerical calculations of the ESR signal evaluating the Kubo formula with exact diagonalization method is difficult, we implement the formula with the help of the random vector technique and the Chebyshev polynominal expansion, which we name the double Chebyshev expansion method. We calculate the temperature dependence of the ESR intensity of ${\rm V}_{15}$ and compare it with the data obtained in experiment. As another complementary approach, we also implement the Kubo formula with the subspace iteration method taking only important low-lying states into account. We study the ESR absorption curve below $100{\rm K}$ by means of both methods. We find that side peaks appear due to the Dzyaloshinsky-Moriya interaction and these peaks grows as temperature decreases.

preprint2001arXiv

Noise effect on the nonadiabatic transition and correction to the tunneling energy gap estimated by the Landau-Zener-Stuckelberg formula

We study effect of noise on the tunneling process of the magnetization in a sweeping field from the viewpoint of the nonadiabatic transition. The tunneling probability depends on the properties of the noise, e.g. the direction, the amplitude, and the relaxation time of the autocorrelation function. We investigate the tunneling process in the presence of noise by solving the time-dependent Schr$\rm{\ddot{o}}$dinger equation numerically, and analyze its parameter dependence. We estimate the effective tunneling gap making use of the Landau-Zener-St$\rm{\ddot{u}}$ckelberg (LZS) formula. The dependence of the gap on the sweeping velocity reproduces the corresponding experimental data qualitatively. For fast sweeping field the estimated gap has no velocity dependence and this gap seems to give the true gap. We, however, show that a correction for the estimated energy gap is necessary and we propose an correction formula.

preprint2000arXiv

Dynamical process for switching dynamics between metastable ordered magnetic state and nonmagnetic ground state in photoinduced phase transition

We propose a dynamical mechanism of the two-way switching between the metastable state and the stable state, which has been found in experiments of photoinduced reversible magnetization and photoinduced structural phase transition. We find that the two-way switching with a non-symmetry breaking perturbation such as illumination is possible only in systems with appropriate parameters. We make it clear that the existence of two time scales in the dynamical process is important for the two-way switching.

preprint1998arXiv

Non-Adiabatic Transition in Spin-Boson Model and Generalization of the Landau-Zener Formula

Non-adiabatic transitions are studied in a spin-boson model with multiple scattering points. In order to generalize the Landau-Zener formula, which describes the case of a single scattering point, we define an ``effective gap'' for a set of scattering points. The generalized formula agrees very well with numerical results of the non-adiabatic dynamics, which we obtained by a direct numerical method. This will make the Landau-Zener formula yet more useful in analyzing experimental data of magnetic-moment inversion.

preprint1997arXiv

Ground-State and Thermodynamic Properties of the Quantum Mixed Spin-1/2-1/2-1-1 Chain

We investigate both analytically and numerically the ground-state and thermodynamic properties of the quantum mixed spin-1/2-1/2-1-1 chain described by the Hamiltonian $H=\sum_{\ell=1}^{N/4} (J_1\vecs_{4\ell-3}\cdot \vecs_{4\ell-2}+J_2\vecs_{4\ell-2}\cdot\vecS_{4\ell-1}+J_3\vecS_{4\ell-1}\cdot \vecS_{4\ell}+J_2\vecS_{4\ell}\cdot\vecs_{4\ell+1})$, where two $S=1/2$ spins ($\vecs_{4\ell-3}$ and $\vecs_{4\ell-2}$) and two S=1 spins ($\vecS_{4\ell-1}$ and $\vecS_{4\ell}$) are arranged alternatively. Inseveral limiting cases of $J_1$, $J_2$, and $J_3$ we apply the Wigner-Eckart theorem and carry out a perturbation calculation to examine the behavior of the massless lines where the energy gap vanishes. Performing a quantum Monte Carlo calculation without global flips at a sufficiently low temperature for the case where $J_1=J_3=1.0$ and $J_2>0$, we find that the ground state of the present system in this case undergoes a second-order phase transition accompanying the vanishing of the energy gap at $J_2=J_{2c$ with $J_{2 c} = 0.77 \pm 0.01$. We also find that the ground states for both $J_2 < J_{2 c$ and $J_2 > J_{2c$ can be understood by means of the valence-bond-solid picture. A quantum Monte Carlo calculation which takes the global flips along the Trotter direction into account is carried out to elucidate the temperature dependences of the specific heat and the magnetic susceptibility. In particular, it is found that the susceptibility per unit cell for $J_2=0.77$ with $J_1=J_3=1.0$ takes a finite value at absolute zero temperature and that the specific heat per unit cell versus temperature curve for $J_2=5.0$ with $J_1=J_3=1.0$ has a double peak.