Source author record

P. A. Rikvold

P. A. Rikvold 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

16works
11topics
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

16 published item(s)

preprint2016arXiv

Equilibrium, Metastability, and Hysteresis in a Model Spin-crossover Material with Nearest-neighbor Antiferromagnetic-like and Long-range Ferromagnetic-like Interactions

Phase diagrams and hysteresis loops were obtained by Monte Carlo simulations and a mean-field method for a simplified model of a spin-crossover material with a two-step transition between the high-spin and low-spin states. This model is a mapping onto a square-lattice $S=1/2$ Ising model with antiferromagnetic nearest-neighbor and ferromagnetic Husimi-Temperley (equivalent-neighbor) long-range interactions. Phase diagrams obtained by the two methods for weak and strong long-range interactions are found to be similar. However, for intermediate-strength long-range interactions, the Monte Carlo simulations show that tricritical points decompose into pairs of critical endpoints and mean-field critical points surrounded by horn-shaped regions of metastability. Hysteresis loops along paths traversing the horn regions are strongly reminiscent of thermal two-step transition loops with hysteresis, recently observed experimentally in several spin-crossover materials. We believe analogous phenomena should be observable in experiments and simulations for many systems that exhibit competition between local antiferromagnetic-like interactions and long-range ferromagnetic-like interactions caused by elastic distortions.

preprint2015arXiv

Monte Carlo Simulations of the Critical Properties of a Ziff-Gulari-Barshad model of Catalytic CO Oxidation with Long-range Reactivity

The Ziff-Gulari-Barshad (ZGB) model, a simplified description of the oxidation of carbon monoxide (CO) on a catalyst surface, is widely used to study properties of nonequilibrium phase transitions. In particular, it exhibits a nonequilibrium, discontinuous transition between a reactive and a CO poisoned phase. If one allows a nonzero rate of CO desorption ($k$), the line of phase transitions terminates at a critical point ($k_{c}$). In this work, instead of restricting the CO and atomic oxygen (O) to react to form carbon dioxide (CO$_{2}$) only when they are adsorbed in close proximity, we consider a modified model that includes an adjustable probability for adsorbed CO and O atoms located far apart on the lattice to react. We employ large-scale Monte Carlo simulations for system sizes up to 240$\times$240 lattice sites, using the crossing of fourth-order cumulants to study the critical properties of this system. We find that the nonequilibrium critical point changes from the two-dimensional Ising universality class to the mean-field universality class upon introducing even a weak long-range reactivity mechanism. This conclusion is supported by measurements of cumulant fixed-point values, cluster percolation probabilities, correlation-length finite-size scaling properties, and the critical exponent ratio $β/ ν$. The observed behavior is consistent with that of the \emph{equilibrium} Ising ferromagnet with additional weak long-range interactions [T. Nakada, P. A. Rikvold, T. Mori, M. Nishino, and S. Miyashita, Phys. Rev. B 84, 054433 (2011)]. The large system sizes and the use of fourth-order cumulants also enable determination with improved accuracy of the critical point of the original ZGB model with CO desorption.

preprint2013arXiv

Model for the catalytic oxidation of CO, including gas-phase impurities and CO desorption

We present results of kinetic Monte Carlo simulations of a modified Ziff-Gulari-Barshad model for the reaction CO+O --> CO_2 on a catalytic surface. Our model includes impurities in the gas phase, CO desorption, and a modification known to eliminate the unphysical O poisoned phase. The impurities can adsorb and desorb on the surface, but otherwise remain inert. In a previous work that did not include CO desorption [G. M. Buendia and P. A. Rikvold, Phys. Rev. E, 85 031143 (2012)], we found that the impurities have very distinctive effects on the phase diagram and greatly diminish the reactivity of the system. If the impurities do not desorb, once the system reaches a stationary state, the CO_2 production disappears. When the impurities are allowed to desorb, there are regions where the CO_2 reaction window reappears, although greatly reduced. Following experimental evidence that indicates that temperature effects are crucial in many catalytic processes, here we further analyze these effects by including a CO desorption rate. We find that the CO desorption has the effect to smooth the transition between the reactive and the CO rich phase, and most importantly it can counteract the negative effects of the presence of impurities by widening the reactive window such that now the system remains catalytically active in the whole range of CO pressures.

preprint2012arXiv

Effects of inert species in the gas phase in a model for the catalytic oxidation of CO

We study by kinetic Monte Carlo simulations the catalytic oxidation of carbon monoxide on a surface in the presence of contaminants in the gas phase. The process is simulated by a Ziff-Gulari-Barshad (ZGB) model that has been modified to include the effect of the contaminants and to eliminate the unphysical oxygen-poisoned phase. The impurities can adsorb and desorb on the surface, but otherwise remain inert. We find that, if the impurities can not desorb, no matter how small their proportion in the gas mixture, the first order transition and the reactive window that characterize the ZGB model disappear. The coverages become continuous, and once the surface has reached a steady state there is no production of CO$_2$. This is quite different from the behavior of a system in which the surface presents a fixed percentage of impurities. When the contaminants are allowed to desorb, the reactive window appears again, and disappears at a value that depends on the proportion of contaminants in the gas and on their desorption rate.

preprint2012arXiv

Optimization of Survivability Analysis for Large-Scale Engineering Networks

Engineering networks fall into the category of large-scale networks with heterogeneous nodes such as sources and sinks. The survivability analysis of such networks requires the analysis of the connectivity of the network components for every possible combination of faults to determine a network response to each combination of faults. From the computational complexity point of view, the problem belongs to the class of exponential time problems at least. Partially, the problem complexity can be reduced by mapping the initial topology of a complex large-scale network with multiple sources and multiple sinks onto a set of smaller sub-topologies with multiple sources and a single sink connected to the network of sources by a single link. In this paper, the mapping procedure is applied to the Florida power grid.

preprint2011arXiv

Kinetic Monte Carlo Simulations of a Model for Heat-assisted Magnetization Reversal in Ultrathin Films

To develop practically useful systems for ultra-high-density information recording with densities above terabits/cm$^2$, it is necessary to simultaneously achieve high thermal stability at room temperature and high recording rates. One method that has been proposed to reach this goal is heat-assisted magnetization reversal (HAMR). In this method, one applies a high-coercivity material, whose coercivity is temporarily lowered during the writing process through localized heating. Here we present kinetic Monte Carlo simulations of a model of HAMR for ultrathin films, in which the temperature in the central part of the film is momentarily increased above the critical temperature, for example by a laser pulse. We observe that the speed-up achieved by this method, relative to the switching time at a constant, subcritical temperature, is optimal for an intermediate strength of the writing field. This effect is explained using the theory of nucleation-induced magnetization switching in finite systems. Our results should be particularly relevant to recording media with strong perpendicular anisotropy, such as ultrathin Co/Pt or Co/Pd multilayers.

preprint2010arXiv

A new battery-charging method suggested by molecular dynamics simulations

Based on large-scale molecular dynamics simulations, we propose a new charging method that should be capable of charging a Lithium-ion battery in a fraction of the time needed when using traditional methods. This charging method uses an additional applied oscillatory electric field. Our simulation results show that this charging method offers a great reduction in the average intercalation time for Li+ ions, which dominates the charging time. The oscillating field not only increases the diffusion rate of Li+ ions in the electrolyte but, more importantly, also enhances intercalation by lowering the corresponding overall energy barrier.

preprint2010arXiv

A New Charging Method for Li-ion Batteries: Dependence of the charging time on the Direction of an Additional Oscillating Field

We have recently proposed a new method for charging Li-ion batteries based on large-scale molecular dynamics studies (I. Abou Hamad et al, Phys. Chem. Chem. Phys., 12, 2740 (2010)). Applying an additional oscillating electric field in the direction perpendicular to the graphite sheets of the anode showed an exponential decrease in charging time with increasing amplitude of the applied oscillating field. Here we present new results exploring the effect on the charging time of changing the orientation of the oscillating field. Results for oscillating fields in three orthogonal directions are compared.

preprint2007arXiv

Two Modes of Magnetization Switching in a Simulated Iron Nanopillar in an Obliquely Oriented Field

Finite-temperature micromagnetics simulations are employed to study the magnetization-switching dynamics driven by a field applied at an angle to the long axis of an iron nanopillar. A bi-modal distribution in the switching times is observed, and evidence for two competing modes of magnetization-switching dynamics is presented. For the conditions studied here, temperature $T = 20$ K and the reversal field 3160 Oe at an angle of 75$^\circ$ to the long axis, approximately 70% of the switches involve unstable decay (no free-energy barrier) and 30% involve metastable decay (a free-energy barrier is crossed). The latter are indistinguishable from switches which are constrained to start at a metastable free-energy minimum. Competition between unstable and metastable decay could greatly complicate applications involving magnetization switches near the coercive field.

preprint2001arXiv

Halide Adsorption on Single-crystal Silver Substrates: Dynamic Simulations and ab-initio Density-functional Theory

We investigate the static and dynamic behaviors of a Br adlayer electrochemically deposited onto single-crystal Ag(100) using an off-lattice model of the adlayer. Unlike previous studies using a lattice-gas model, the off-lattice model allows adparticles to be located at any position within a two-dimensional approximation to the substrate. Interactions with the substrate are approximated by a corrugation potential. Using Density Functional Theory (DFT) to calculate surface binding energies, a sinusoidal approximation to the corrugation potential is constructed. A variety of techniques, including Monte Carlo and Langevin simulations, are used to study the behavior of the adlayer. The lateral root-mean-square (rms) deviation of the adparticles from the binding sites is presented along with equilibrium coverage isotherms, and the thermally activated Arrhenius barrier-hopping model used in previous dynamic Monte Carlo simulations is tested.

preprint1998arXiv

A projection method for statics and dynamics of lattice spin systems

A method based on Monte Carlo sampling of the probability flows projected onto the subspace of one or more slow variables is proposed for investigation of dynamic and static properties of lattice spin systems. We illustrate the method by applying it, with projection onto the order-parameter subspace, to the three-dimensional 3-state Potts model in equilibrium and to metastable decay in a three-dimensional 3-state kinetic Potts model.

preprint1995arXiv

A Method to Study Relaxation of Metastable Phases: Macroscopic Mean-Field Dynamics

We propose two different macroscopic dynamics to describe the decay of metastable phases in many-particle systems with local interactions. These dynamics depend on the macroscopic order parameter $m$ through the restricted free energy $F(m)$ and are designed to give the correct equilibrium distribution for $m$. The connection between macroscopic dynamics and the underlying microscopic dynamic are considered in the context of a projection- operator formalism. Application to the square-lattice nearest-neighbor Ising ferromagnet gives good agreement with droplet theory and Monte Carlo simulations of the underlying microscopic dynamic. This includes quantitative agreement for the exponential dependence of the lifetime on the inverse of the applied field $H$, and the observation of distinct field regions in which the derivative of the lifetime with respect to $1/H$ depends differently on $H$. In addition, at very low temperatures we observe oscillatory behavior of this derivative with respect to $H$, due to the discreteness of the lattice and in agreement with rigorous results. Similarities and differences between this work and earlier works on finite Ising models in the fixed-magnetization ensemble are discussed.

preprint1995arXiv

Coadsorption of Copper and Sulfate on Au(111) Electrodes: Monte Carlo Simulation of a Lattice-Gas Model

We report ground-state calculations and Monte Carlo simulations for a lattice-gas model of the underpotential deposition of copper on Au(111) in sulfate-containing electrolytes. In a potential range of approximately 100$\sim$150 mV, this system exhibits a $(\sqrt3\!\times\!\sqrt3)$ mixed phase with 2/3 monolayer (ML) copper and 1/3 ML sulfate. Our simulation results agree well with experimental results and with other theoretical work.

preprint1994arXiv

Application of a constrained-transfer-matrix method to metastability in the $d$$=$$2$ Ising ferromagnet

Applying a numerical transfer-matrix formalism, we obtain complex-valued constrained free energies for the two-dimensional square-lattice nearest-neighbor Ising ferromagnet below its critical temperature and in an external magnetic field. In particular, we study the imaginary part of the constrained free-energy branch that corresponds to the metastable phase. Although droplets are not introduced explicitly, the metastable free energy is obtained in excellent agreement with field-theoretical droplet-model predictions. The finite-size scaling properties are different in the weak-field and intermediate-field regimes, and we identify the corresponding different critical-droplet shapes. For intermediate fields, we show that the surface free energy of the critical droplet is given by a Wulff construction with the equilibrium surface tension. We also find a prefactor exponent in complete agreement with the field-theoretical droplet model. Our results extend the region of validity for known results of this field-theoretical droplet model, and they indicate that this transfer-matrix approach provides a nonperturbative numerical continuation of the equilibrium free energy into the metastable phase.

preprint1994arXiv

Recent Results on the Decay of Metastable Phases

We review some aspects of current knowledge regarding the decay of metastable phases in many-particle systems. In particular we emphasize recent theoretical and computational developments and numerical results regarding homogeneous nucleation and growth in kinetic Ising and lattice-gas models. An introductory discussion of the droplet theory of homogeneous nucleation is followed by a discussion of Monte Carlo and transfer-matrix methods commonly used for numerical study of metastable decay, including some new algorithms. Next we discuss specific classes of systems. These include a brief discussion of recent progress for fluids, and more exhaustive considerations of ferromagnetic Ising models ({\it i.e.}, attractive lattice-gas models) with weak long-range interactions and with short-range interactions. Whereas weak-long-range-force (WLRF) models have infinitely long-lived metastable phases in the infinite-range limit, metastable phases in short-range-force (SRF) models eventually decay, albeit extremely slowly. Recent results on the finite-size scaling of metastable lifetimes in SRF models are reviewed, and it is pointed out that such effects may be experimentally observable.