Researcher profile

Creighton K. Thomas

Creighton K. Thomas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
12works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

12 published item(s)

preprint2016arXiv

Charge avalanches and depinning in the Coulomb glass: The role of long-range interactions

We explore the stability of far-from-equilibrium metastable states of a three-dimensional Coulomb glass at zero temperature by studying charge avalanches triggered by a slowly varying external electric field. Surprisingly, we identify a sharply defined dynamical ("depinning") phase transition from stationary to nonstationary charge displacement at a critical value of the external electric field. Using particle-conserving dynamics, scale-free system-spanning avalanches are observed only at the critical field. We show that the qualitative features of this depinning transition are completely different for an equivalent short-range model, highlighting the key importance of long-range interactions for nonequilibrium dynamics of Coulomb glasses.

preprint2015arXiv

Coulomb energy of uniformly-charged spheroidal shell systems

We provide exact expressions for the electrostatic energy of uniformly-charged prolate and oblate spheroidal shells. We find that uniformly-charged prolate spheroids of eccentricity greater than 0.9 have lower Coulomb energy than a sphere of the same area. For the volume-constrained case, we find that a sphere has the highest Coulomb energy among all spheroidal shells. Further, we derive the change in the Coulomb energy of a uniformly-charged shell due to small, area-conserving perturbations on the spherical shape. Our perturbation calculations show that buckling-type deformations on a sphere can lower the Coulomb energy. Finally, we consider the possibility of counterion condensation on the spheroidal shell surface. We employ a Manning-Oosawa two-state model approximation to evaluate the renormalized charge and analyze the behavior of the equilibrium free energy as a function of the shell's aspect ratio for both area-constrained and volume-constrained cases. Counterion condensation is seen to favor the formation of spheroidal structures over a sphere of equal area for high values of shell volume fractions.

preprint2014arXiv

Electrostatics-driven shape transitions in soft shells

Manipulating the shape of nanoscale objects in a controllable fashion is at the heart of designing materials that act as building blocks for self-assembly or serve as targeted drug delivery carriers. Inducing shape deformations by controlling external parameters is also an important way of designing biomimetic membranes. In this paper, we demonstrate that electrostatics can be used as a tool to manipulate the shape of soft, closed membranes by tuning environmental conditions such as the electrolyte concentration in the medium. Using a molecular dynamics-based simulated annealing procedure, we investigate charged elastic shells that do not exchange material with their environment, such as elastic membranes formed in emulsions or synthetic nanocontainers. We find that by decreasing the salt concentration or increasing the total charge on the shell's surface, the spherical symmetry is broken, leading to the formation of ellipsoids, discs, and bowls. Shape changes are accompanied by a significant lowering of the electrostatic energy and a rise in the surface area of the shell. To substantiate our simulation findings, we show analytically that a uniformly charged disc has a lower Coulomb energy than a sphere of the same volume. Further, we test the robustness of our results by including the effects of charge renormalization in the analysis of the shape transitions and find the latter to be feasible for a wide range of shell volume fractions.

preprint2013arXiv

Novel disordering mechanism in ferromagnetic systems with competing interactions

Ferromagnetic Ising systems with competing interactions are considered in the presence of a random field. We find that in three space dimensions the ferromagnetic phase is disordered by a random field which is considerably smaller than the typical interaction strength between the spins. This is the result of a novel disordering mechanism triggered by an underlying spin-glass phase. Calculations for the specific case of the long-range dipolar LiHo_xY_{1-x}F_4 compound suggest that the above mechanism is responsible for the peculiar dependence of the critical temperature on the strength of the random field and the broadening of the susceptibility peaks as temperature is decreased, as found in recent experiments by Silevitch et al. [Nature (London) 448, 567 (2007)]. Our results thus emphasize the need to go beyond the standard Imry-Ma argument when studying general random-field systems.

preprint2013arXiv

Numerically exact correlations and sampling in the two-dimensional Ising spin glass

A powerful existing technique for evaluating statistical mechanical quantities in two-dimensional Ising models is based on constructing a matrix representing the nearest neighbor spin couplings and then evaluating the Pfaffian of the matrix. Utilizing this technique and other more recent developments in evaluating elements of inverse matrices and exact sampling, a method and computer code for studying two-dimensional Ising models is developed. The formulation of this method is convenient and fast for computing the partition function and spin correlations. It is also useful for exact sampling, where configurations are directly generated with probability given by the Boltzmann distribution. These methods apply to Ising model samples with arbitrary nearest-neighbor couplings and can also be applied to general dimer models. Example results of computations are described, including comparisons with analytic results for the ferromagnetic Ising model, and timing information is provided.

preprint2013arXiv

Shapes of pored membranes

We study the shapes of pored membranes within the framework of the Helfrich theory under the constraints of fixed area and pore size. We show that the mean curvature term leads to a budding- like structure, while the Gaussian curvature term tends to flatten the membrane near the pore; this is corroborated by simulation. We propose a scheme to deduce the ratio of the Gaussian rigidity to the bending rigidity simply by observing the shape of the pored membrane. This ratio is usually difficult to measure experimentally. In addition, we briefly discuss the stability of a pore by relaxing the constraint of a fixed pore size and adding the line tension. Finally, the flattening effect due to the Gaussian curvature as found in studying pored membranes is extended to two-component membranes. We find that sufficiently high contrast between the components' Gaussian rigidities leads to budding which is distinct from that due to the line tension.

preprint2012arXiv

Boolean decision problems with competing interactions on scale-free networks: Critical thermodynamics

We study the critical behavior of Boolean variables on scale-free networks with competing interactions (Ising spin glasses). Our analytical results for the disorder-network-decay-exponent phase diagram are verified using Monte Carlo simulations. When the probability of positive (ferromagnetic) and negative (antiferromagnetic) interactions is the same, the system undergoes a finite-temperature spin-glass transition if the exponent that describes the decay of the interaction degree in the scale-free graph is strictly larger than 3. However, when the exponent is equal to or less than 3, a spin-glass phase is stable for all temperatures. The robustness of both the ferromagnetic and spin-glass phases suggests that Boolean decision problems on scale-free networks are quite stable to local perturbations. Finally, we show that for a given decay exponent spin glasses on scale-free networks seem to obey universality. Furthermore, when the decay exponent of the interaction degree is larger than 4 in the spin-glass sector, the universality class is the same as for the mean-field Sherrington-Kirkpatrick Ising spin glass.

preprint2011arXiv

Optimizing glassy p-spin models

Computing the ground state of Ising spin-glass models with p-spin interactions is, in general, an NP-hard problem. In this work we show that unlike in the case of the standard Ising spin glass with two-spin interactions, computing ground states with p=3 is an NP-hard problem even in two space dimensions. Furthermore, we present generic exact and heuristic algorithms for finding ground states of p-spin models with high confidence for systems of up to several thousand spins.

preprint2011arXiv

Sampling the ground-state magnetization of d-dimensional p-body Ising models

We demonstrate that a recently introduced heuristic optimization algorithm [Phys. Rev. E 83, 046709 (2011)] that combines a local search with triadic crossover genetic updates is capable of sampling nearly uniformly among ground-state configurations in spin-glass-like Hamiltonians with p-spin interactions in d space dimensions that have highly degenerate ground states. Using this algorithm we probe the zero-temperature ferromagnet to spin-glass transition point q_c of two example models, the disordered version of the two-dimensional three-spin Baxter-Wu model [q_c = 0.1072(1)] and the three-dimensional Edwards-Anderson model [q_c = 0.2253(7)], by computing the Binder ratio of the ground-state magnetization.

preprint2011arXiv

Simplest model to study reentrance in physical systems

We numerically investigate the necessary ingredients for reentrant behavior in the phase diagram of physical systems. Studies on the possibly simplest model that exhibits reentrance, the two-dimensional random bond Ising model, show that reentrant behavior is generic whenever frustration is present in the model. For both discrete and continuous disorder distributions, the phase diagram in the disorder-temperature plane is found to be reentrant, where for some disorder strengths a paramagnetic phase exists at both high and low temperatures, but an ordered ferromagnetic phase exists for intermediate temperatures.

preprint2011arXiv

Universality in phase boundary slopes for spin glasses on self dual lattices

We study the effects of disorder on the slope of the disorder--temperature phase boundary near the Onsager point (Tc = 2.269...) in spin-glass models. So far, studies have focused on marginal or irrelevant cases of disorder. Using duality arguments, as well as exact Pfaffian techniques we reproduce these analytical estimates. In addition, we obtain different estimates for spin-glass models on hierarchical lattices where the effects of disorder are relevant. We show that the phase-boundary slope near the Onsager point can be used to probe for the relevance of disorder effects.

preprint2010arXiv

Chaos and universality in two-dimensional Ising spin glasses

Recently extended precise numerical methods and droplet scaling arguments allow for a coherent picture of the glassy states of two-dimensional Ising spin glasses to be assembled. The length scale at which entropy becomes important and produces "chaos", the extreme sensitivity of the state to temperature, is found to depend on the type of randomness. For the $\pm J$ model this length scale dominates the low-temperature specific heat. Although there is a type of universality, some critical exponents do depend on the distribution of disorder.