Fast Algorithm for Finding the Eigenvalue Distribution of Very Large Matrices

A theoretical analysis is given of the equation of motion method, due to Alben et al., to compute the eigenvalue distribution (density of states) of very large matrices. The salient feature of this method is that for matrices of the kind encountered in quantum physics the memory and CPU requirements of this method scale linearly with the dimension of the matrix. We derive a rigorous estimate of the statistical error, supporting earlier observations that the computational efficiency of this approach increases with matrix size. We use this method and an imaginary-time version of it to compute the energy and the specific heat of three different, exactly solvable, spin-1/2 models and compare with the exact results to study the dependence of the statistical errors on sample and matrix size.Comment: 24 pages, 24 figure

Computer simulation of Wheeler's delayed choice experiment with photons

We present a computer simulation model of Wheeler's delayed choice experiment that is a one-to-one copy of an experiment reported recently (V. Jacques {\sl et al.}, Science 315, 966 (2007)). The model is solely based on experimental facts, satisfies Einstein's criterion of local causality and does not rely on any concept of quantum theory. Nevertheless, the simulation model reproduces the averages as obtained from the quantum theoretical description of Wheeler's delayed choice experiment. Our results prove that it is possible to give a particle-only description of Wheeler's delayed choice experiment which reproduces the averages calculated from quantum theory and which does not defy common sense.Comment: Europhysics Letters (in press

Beyond the poor man's implementation of unconditionally stable algorithms to solve the time-dependent Maxwell Equations

For the recently introduced algorithms to solve the time-dependent Maxwell equations (see Phys.Rev.E Vol.64 p.066705 (2001)), we construct a variable grid implementation and an improved spatial discretization implementation that preserve the property of the algorithms to be unconditionally stable by construction. We find that the performance and accuracy of the corresponding algorithms are significant and illustrate their practical relevance by simulating various physical model systems.Comment: 18 pages, 16 figure

Photon and spin dependence of the resonance lines shape in the strong coupling regime

We study the quantum dynamics of a spin ensemble coupled to cavity photons. Recently, related experimental results have been reported, showing the existence of the strong coupling regime in such systems. We study the eigenenergy distribution of the multi-spin system (following the Tavis-Cummings model) which shows a peculiar structure as a function of the number of cavity photons and of spins. We study how this structure causes changes in the spectrum of the admittance in the linear response theory, and also the frequency dependence of the excited quantities in the stationary state under a probing field. In particular, we investigate how the structure of the higher excited energy levels changes the spectrum from a double-peak structure (the so-called vacuum field Rabi splitting) to a single peak structure. We also point out that the spin dynamics in the region of the double-peak structure corresponds to recent experiments using cavity ringing while in region of the single peak structure, it corresponds to the coherent Rabi oscillation in a driving electromagnetic filed. Using a standard Lindblad type mechanism, we study the effect of dissipations on the line width and separation in the computed spectra. In particular, we study the relaxation of the total spin in the general case of a spin ensemble in which the total spin of the system is not specified. The theoretical results are correlated with experimental evidence of the strong coupling regime, achieved with a spin 1/2 ensemble

Lattice dynamics of random and quasiperiodic heterostructures

We report on the quantum-mechanical displacement form factor in quasiperiodic and random heterostructures. A one-dimensional treatment is adopted to describe the longitudinal displacement along the growth axis. Elastic properties are assumed to be homogeneous, while the inhomogeneous mass density characterizes the heterostructure. In the low-frequency limit, the peak structure can be attributed to acoustic phonons, whereas for higher frequencies the quasiperiodic and random cases differ markedly. In the quasiperiodic case and constant momentum transfer, resonances separated by gaps occur and their number depends on the resolution in the frequency domain. The random case is dominated by an acoustic resonance becoming broader with increasing frequency

Logical Hidden Markov Models

Logical hidden Markov models (LOHMMs) upgrade traditional hidden Markov models to deal with sequences of structured symbols in the form of logical atoms, rather than flat characters. This note formally introduces LOHMMs and presents solutions to the three central inference problems for LOHMMs: evaluation, most likely hidden state sequence and parameter estimation. The resulting representation and algorithms are experimentally evaluated on problems from the domain of bioinformatics

A constrained stochastic state selection method applied to quantum spin systems

We describe a further development of the stochastic state selection method, which is a kind of Monte Carlo method we have proposed in order to numerically study large quantum spin systems. In the stochastic state selection method we make a sampling which is simultaneous for many states. This feature enables us to modify the method so that a number of given constraints are satisfied in each sampling. In this paper we discuss this modified stochastic state selection method that will be called the constrained stochastic state selection method in distinction from the previously proposed one (the conventional stochastic state selection method) in this paper. We argue that in virtue of the constrained sampling some quantities obtained in each sampling become more reliable, i.e. their statistical fluctuations are less than those from the conventional stochastic state selection method. In numerical calculations of the spin-1/2 quantum Heisenberg antiferromagnet on a 36-site triangular lattice we explicitly show that data errors in our estimation of the ground state energy are reduced. Then we successfully evaluate several low-lying energy eigenvalues of the model on a 48-site lattice. Our results support that this system can be described by the theory based on the spontaneous symmetry breaking in the semiclassical Neel ordered antiferromagnet.Comment: 15 pgaes, 5 figure

The Stochastic State Selection Method Combined with the Lanczos Approach to Eigenvalues in Quantum Spin Systems

We describe a further development of the stochastic state selection method, a new Monte Carlo method we have proposed recently to make numerical calculations in large quantum spin systems. Making recursive use of the stochastic state selection technique in the Lanczos approach, we estimate the ground state energy of the spin-1/2 quantum Heisenberg antiferromagnet on a 48-site triangular lattice. Our result for the upper bound of the ground state energy is -0.1833 +/- 0.0003 per bond. This value, being compatible with values from other work, indicates that our method is efficient in calculating energy eigenvalues of frustrated quantum spin systems on large lattices.Comment: 11 page

Possible Experience: from Boole to Bell

Mainstream interpretations of quantum theory maintain that violations of the Bell inequalities deny at least either realism or Einstein locality. Here we investigate the premises of the Bell-type inequalities by returning to earlier inequalities presented by Boole and the findings of Vorob'ev as related to these inequalities. These findings together with a space-time generalization of Boole's elements of logic lead us to a completely transparent Einstein local counterexample from everyday life that violates certain variations of the Bell inequalities. We show that the counterexample suggests an interpretation of the Born rule as a pre-measure of probability that can be transformed into a Kolmogorov probability measure by certain Einstein local space-time characterizations of the involved random variables.Comment: Published in: EPL, 87 (2009) 6000