Monte Carlo simulations for two-dimensional quantum systems

The Trotter-Suzuki transformation has been used to obtain the classical representation ford-dimensional lattice systems with boson and fermion degrees of freedom. A Monte Carlo algorithm for the equivalent (d+1)-dimensional classical system is presented. Numerical results are shown for the Heisenberg-spin-glass, the XY model and the spinless fermion lattice gas in two dimensions.     

Monte Carlo study of hot electrons in quantum wells

We have performed an ensemble Monte Carlo calculation of the high field response at 77 K of two different two-dimensional systems, one a modulation doped (MD) single square well of GaAs/AlGaAs and the other a δ-doping layer formed by a sheet of donor impurities embedded in bulk GaAs. Results of this calculation show a reduced low field mobility and a decrease in the transient overshoot velocity in the δ-doping versus MD system due to the effect of impurity scattering. At high fields, however, the steady state saturated drift velocities appear to be similar in the two systems. Because of the higher carrier density, the performance of short-channel devices based on the δ-layer will be advantageous where high output currents are necessary.     

Variational quantum Monte Carlo simulations with tensor-network states

We show that the formalism of tensor-network states, such as the matrix-product states (MPS), can be used as a basis for variational quantum Monte Carlo simulations. Using a stochastic optimization method, we demonstrate the potential of this approach by explicit MPS calculations for the transverse Ising chain with up to N=256 spins at criticality, using periodic boundary conditions and D x D matrices with D up to 48. The computational cost of our scheme formally scales as ND3, whereas standard MPS approaches and the related density matrix renormalization group method scale as ND5 and ND6, respectively, for periodic systems.     

Pfaffian pairing wave functions in electronic-structure quantum Monte Carlo simulations

We investigate the accuracy of trial wave functions for quantum Monte Carlo based on Pfaffian functional form with singlet and triplet pairing. Using a set of first row atoms and molecules we find that these wave functions provide very consistent and systematic behavior in recovering the correlation energies on the level of 95%. In order to get beyond this limit we explore the possibilities of multi-Pfaffian pairing wave functions. We show that a small number of Pfaffians recovers another large fraction of the missing correlation energy comparable to the larger-scale configuration interaction wave functions. We also find that Pfaffians lead to substantial improvements in fermion nodes when compared to Hartree-Fock wave functions.     

Kinetic Monte Carlo simulations of three-dimensional self-assembled quantum dot islands

By three-dimensional kinetic Monte Carlo simulations, the effects of the temperature, the flux rate, the total coverage and the interruption time on the distribution and the number of self-assembled InAs/GaAs(001) quantum dot(QD) islands are studied, which shows that a higher temperature, a lower flux rate and a longer growth time correspond to a better island distribution. The relations between the number of islands and the temperature and the flux rate are also successfully simulated. It is observed that for the total coverage lower than 0.5 ML, the number of islands decreases with the temperature increasing and other growth parameters fixed and the number of islands increases with the flux rate increasing when the deposition is lower than 0.6 ML and the other parameters are fixed.     

Loop updates for quantum Monte Carlo simulations in the canonical ensemble

We present a new nonlocal updating scheme for quantum Monte Carlo simulations, which conserves particle number and other symmetries. It allows exact symmetry projection and direct evaluation of the equal-time Green's function and other observables in the canonical ensemble. The method is applicable to a wide variety of systems. We show results for bosonic atoms in optical lattices, neutron pairs in atomic nuclei, and electron pairs in ultrasmall superconducting grains.     

Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations

Quantum Monte Carlo simulations, while being efficient for bosons, suffer from the "negative sign problem" when applied to fermions--causing an exponential increase of the computing time with the number of particles. A polynomial time solution to the sign problem is highly desired since it would provide an unbiased and numerically exact method to simulate correlated quantum systems. Here we show that such a solution is almost certainly unattainable by proving that the sign problem is nondeterministic polynomial (NP) hard, implying that a generic solution of the sign problem would also solve all problems in the complexity class NP in polynomial time.     

Supercooled superfluids in Monte Carlo simulations

We perform path integral Monte Carlo simulations to study the imaginary time dynamics of metastable supercooled superfluid states and nearly superglassy states of a one component fluid of spinless bosons square wells. Our study shows that the identity of the particles and the exchange symmetry is crucial for the frustration necessary to obtain metastable states in the quantum regime. Whereas the simulation time has to be chosen to determine whether we are in a metastable state or not, the imaginary time dynamics tells us if we are or not close to an arrested glassy state.     

EPR and ferromagnetism in diluted magnetic semiconductor quantum wells

Motivated by recent measurements of electron paramagnetic resonance spectra in modulation-doped CdMnTe quantum wells [Phys. Rev. Lett. 91, 077201 (2003)]], we develop a theory of collective spin excitations in quasi-two-dimensional diluted magnetic semiconductors. Our theory explains the anomalously large Knight shift found in these experiments as a consequence of collective coupling between Mn-ion local moments and itinerant-electron spins. We use this theory to discuss the physics of ferromagnetism in (II,Mn)VI quantum wells and to speculate on the temperature at which it is likely to be observed in n-type modulation-doped systems.     

Convergence behaviour of Green's function quantum Monte Carlo simulations of pi electron systems

Green's function quantum Monte Carlo (GF QMC) simulations of fermionic ensembles require the definition of a so-called spectral parameter W to yield trustworthy estimates of the fully correlated ground state energy E. In the present work we discuss the influence of the shift parameter, W, on the convergence behaviour of GF QMC simulations. As model systems we have considered some π molecules which are studied in the framework of the Pariser–Parr–Pople (PPP) Hamiltonian. The influence of the W parameter on the convergence of the GF QMC simulations in many-electron systems with an odd number of electronic permutations within one spin direction exceeds the influence observed in systems without such odd permutations. They do not occur in polyenes and Huckel annulenes with an electron count of (4n + 2) (n = 0, 1, 2…). The GF QMC technique adopted as a computational tool has been developed to study π molecules with fermionic sign problems owing to odd electronic interchanges within one spin direction.     

Equation of state and Raman frequency of diamond from quantum Monte Carlo simulations

We report variational and diffusion quantum Monte Carlo (VMC and DMC) calculations of the equation of state and Raman frequency of diamond. The pressure derivative of the bulk modulus, which has not been accurately determined experimentally, is calculated to be 3.8 (VMC) and 3.7 (DMC). The values of the Raman frequency calculated at the experimental volume are 1373 cm-1(VMC) and 1359 cm-1 (DMC), in good agreement with the experimental value of 1333 cm-1.     

Monte Carlo simulations of sexual reproduction

Modifying the Redfield model of sexual reproduction and the Penna model of biological aging, we compare reproduction with and without recombination in age-structured populations. In constrast to Redfield and in agreement with Bernardes we find sexual reproduction to be preferred to asexual one. In particular, the presence of old but still reproducing males helps the survival of younger females beyond their reproductive age.     

Monte Carlo simulations of parapatric speciation

Parapatric speciation is studied using an individual-based model with sexual reproduction. We combine the theory of mutation accumulation for biological ageing with an environmental selection pressure that varies according to the individuals geographical positions and phenotypic traits. Fluctuations and genetic diversity of large populations are crucial ingredients to model the features of evolutionary branching and are intrinsic properties of the model. Its implementation on a spatial lattice gives interesting insights into the population dynamics of speciation on a geographical landscape and the disruptive selection that leads to the divergence of phenotypes. Our results suggest that assortative mating is not an obligatory ingredient to obtain speciation in large populations at low gene flow.