Researcher profile

Owen Bradley

Owen Bradley contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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 graph slice

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

Published work

4 published item(s)

preprint2022arXiv

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. We perform a variety of tests, and in all cases find rapid convergence of the chemical potential -- inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

preprint2022arXiv

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

preprint2022arXiv

Instabilities of spin-1 Kitaev spin liquid phase in presence of single-ion anisotropies

We study the spin-one Kitaev model on the honeycomb lattice in the presence of single-ion anisotropies. We consider two types of single ion anisotropies: A $D_{111}$ anisotropy which preserves the symmetry between $X$, $Y$, and $Z$ bonds but violates flux conservation and a $D_{100}$ anisotropy that breaks the symmetry between $X$, $Y$, and $Z$ bonds but preserves flux conservation. We use series expansion methods, degenerate perturbation theory, and exact diagonalization to study these systems. Large positive $D_{111}$ anisotropy leads to a simple product ground state with conventional magnon-like excitations, while large negative $D_{111}$ leads to a broken symmetry and degenerate ground states. For both signs there is a phase transition at a small $|D_{111}| \approx 0.12$ separating the more conventional phases from the Kitaev spin liquid phase. With large $D_{100}$ anisotropy, the ground state is a simple product state, but the model lacks conventional dispersive excitations due to the large number of conservation laws. Large negative $D_{100}$ leads to decoupled one-dimensional systems and many degenerate ground states. No evidence of a phase transition is seen in our numerical studies at any finite $D_{100}$. Convergence of the series expansion extrapolations all the way to $D_{100}=0$ suggests that the nontrivial Kitaev spin liquid is a singular limit of this type of single-ion anisotropy going to zero, which also restores symmetry between the $X$, $Y$, and $Z$ bonds.

preprint2021arXiv

Thermodynamic behavior of modified integer-spin Kitaev models on the honeycomb lattice

We study the thermodynamic properties of modified spin-$S$ Kitaev models introduced by Baskaran, Sen and Shankar (Phys. Rev. B 78, 115116 (2008)). These models have the property that for half-odd-integer spins their eigenstates map on to those of spin-1/2 Kitaev models, with well-known highly entangled quantum spin-liquid states and Majorana fermions. For integer spins, the Hamiltonian is made out of commuting local operators. Thus, the eigenstates can be chosen to be completely unentangled between different sites, though with a significant degeneracy for each eigenstate. For half-odd-integer spins, the thermodynamic properties can be related to the spin-1/2 Kitaev models apart from an additional degeneracy. Hence we focus here on the case of integer spins. We use transfer matrix methods, high temperature expansions and Monte Carlo simulations to study the thermodynamic properties of ferromagnetic and antiferromagnetic models with spin $S=1$ and $S=2$. Apart from large residual entropies, which all the models have, we find that they can have a variety of different behaviors. Transfer matrix calculations show that for the different models, the correlation lengths can be finite as $T\to 0$, become critical as $T\to 0$ or diverge exponentially as $T\to 0$. There is a conserved $Z_2$ flux variable associated with each hexagonal plaquette which saturates at the value $+1$ as $T\rightarrow0$ in all models except the $S=1$ antiferromagnet where the mean flux remains zero as $T\to 0$. We provide qualitative explanations for these results.