Source author record

Per Arne Rikvold

Per Arne 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

29works
14topics
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

29 published item(s)

preprint2022arXiv

Competing gravitational and viscous efects in 3D two-phase flow investigated with a table-top optical scanner

We present experiments and theory describing the transition from viscosity-stabilized flow to gravitationally unstable fingering for two-phase flow in a 3D synthetic porous medium. Observation is made possible by the use of our newly developed table-top 3D-scanner based on optical index matching and laser-induced fluorescence, which is described in detail. In the experiment, a more dense, more viscous fluid injected at a fixed flow-rate from a point source at the top of the flow cell displaces a less viscous, less dense fluid. We observe a stable invasion zone near the inlet, which increases in size with increasing flow rates, and presents initially a close to hemispherical shape. At later times, the invasion front transits to an unstable mode and a fingering flow regime. The transition occurs at a predicted critical radius, Rc, corresponding to the zero of the combined viscous and gravitational pressure gradient.

preprint2020arXiv

Absorbing Random Walks Interpolating Between Centrality Measures on Complex Networks

Centrality, which quantifies the "importance" of individual nodes, is among the most essential concepts in modern network theory. As there are many ways in which a node can be important, many different centrality measures are in use. Here, we concentrate on versions of the common betweenness and it closeness centralities. The former measures the fraction of paths between pairs of nodes that go through a given node, while the latter measures an average inverse distance between a particular node and all other nodes. Both centralities only consider shortest paths (i.e., geodesics) between pairs of nodes. Here we develop a method, based on absorbing Markov chains, that enables us to continuously interpolate both of these centrality measures away from the geodesic limit and toward a limit where no restriction is placed on the length of the paths the walkers can explore. At this second limit, the interpolated betweenness and closeness centralities reduce, respectively, to the well-known it current betweenness and resistance closeness (information) centralities. The method is tested numerically on four real networks, revealing complex changes in node centrality rankings with respect to the value of the interpolation parameter. Non-monotonic betweenness behaviors are found to characterize nodes that lie close to inter-community boundaries in the studied networks.

preprint2020arXiv

Visualization by optical fluorescence of two-phase flow in a three-dimensional porous medium

Slow flow of a single fluid through a porous medium is well understood on a macroscopic level through Darcy's law, a linear relation between flow rate and a combination of pressure differences, viscosity, and gravitational forces. Two-phase flow is complicated by the interface separating the fluids, but understanding of two-dimensional, two-phase flow has been obtained from experiments using transparent cells. In most three-dimensional media, however, visual observation is difficult. Here, we present preliminary results of experiments on a model medium consisting of randomly packed glass spheres, in which one fluorescent liquid invades another. By refractive index matching and scanning with a sheet-shaped laser beam, we obtain slices of the flow patterns, which we combine into three-dimensional pictures. We observe a compact region of invading fluid, surrounded by finger-like protrusions. The compact region becomes more dominant with increasing invader flow rate. The patterns are theoretically analyzed in terms of the interplay between gravitational, viscous, and capillary forces.

preprint2015arXiv

A model for the catalytic oxidation of CO that includes CO desorption and diffusion, O repulsion, and impurities in the gas phase

We present kinetic Monte Carlo simulations exploring the nonequilibrium phase diagram of a modified Ziff-Gulari-Barshad (ZGB) dynamic lattice-gas model for the catalytic oxidation of carbon monoxide (CO) on a surface. The modified model includes the simultaneous presence of contaminants in the gas phase, CO desorption, CO diffusion, and strong repulsion between adsorbed oxygen (O) atoms; all of which have been observed in experimental systems. We find that the strong O-O repulsion produces higher reaction rates, albeit in a reduced reactive pressure window. In systems with impurities, the CO$_2$ production rate is greatly reduced, but this effect is mitigated by CO desorption and diffusion. CO desorption has the effect of widening the reactive pressure window, while CO diffusion has the effect of increasing the reaction rate. In some parameter regimes the presence of impurities destroys the discontinuous transition between the reactive and high CO coverage phases.

preprint2015arXiv

A Universal Lifetime Distribution for Multi-Species Systems

Lifetime distributions of social entities, such as enterprises, products, and media contents, are one of the fundamental statistics characterizing the social dynamics. To investigate the lifetime distribution of mutually interacting systems, simple models having a rule for additions and deletions of entities are investigated. We found a quite universal lifetime distribution for various kinds of inter-entity interactions, and it is well fitted by a stretched-exponential function with an exponent close to 1/2. We propose a "modified Red-Queen" hypothesis to explain this distribution. We also review empirical studies on the lifetime distribution of social entities, and discussed the applicability of the model.

preprint2015arXiv

Centrality Fingerprints for Power Grid Network Growth Models

In our previous work, we have shown that many of the properties of the Florida power grid are reproduced by deterministic network growth models based on the minimization of energy dissipation $E_\mathrm{diss}$. As there is no $a~ priori$ best $E_\mathrm{diss}$ minimizing growth model, we here present a tool, called the "centrality fingerprint," for probing the behavior of different growth models. The centrality fingerprints are comparisons of the current flow into/out of the network with the values of various centrality measures calculated at every step of the growth process. Finally, we discuss applications to the Maryland power grid.

preprint2014arXiv

Architecture of the Florida Power Grid as a Complex Network

We study the Florida high-voltage power grid as a technological network embedded in space. Measurements of geographical lengths of transmission lines, the mixing of generators and loads, the weighted clustering coefficient, as well as the organization of edge conductance weights show a complex architecture quite different from random-graph models usually considered. In particular, we introduce a parametrized mixing matrix to characterize the mixing pattern of generators and loads in the Florida Grid, which is intermediate between the random mixing case and the semi-bipartite case where generator-generator transmission lines are forbidden. Our observations motivate an investigation of optimization (design) principles leading to the structural organization of power grids. We thus propose two network optimization models for the Florida Grid as a case study. Our results show that the Florida Grid is optimized not only by reducing the construction cost (measured by the total length of power lines), but also through reducing the total pairwise edge resistance in the grid, which increases the robustness of power transmission between generators and loads against random line failures. We then embed our models in spatial areas of different aspect ratios and study how this geometric factor affects the network structure, as well as the box-counting fractal dimension of the grids generated by our models.

preprint2014arXiv

Monte Carlo Studies of the Ising Antiferromagnet with a Ferromagnetic Mean-field Term

The unusual thermodynamic properties of the Ising antiferromagnet supplemented with a ferromagnetic, mean-field term are outlined. This simple model is inspired by more realistic models of spin-crossover materials. The phase diagram is estimated using Metropolis Monte Carlo methods, and differences with preliminary Wang-Landau Monte Carlo results for small systems are noted.

preprint2013arXiv

Asymptotic forms and scaling properties of the relaxation time near threshold points in spinodal-type dynamical phase transitions

We study critical properties of the relaxation time at a threshold point in switching processes between bistable states under change of external fields. In particular, we investigate the relaxation processes near the spinodal point of the infinitely long-range interaction model (the Husimi-Temperley model) by analyzing the scaling properties of the corresponding Fokker-Planck equation. We also confirm the obtained scaling relations by direct numerical solution of the original master equation, and by kinetic Monte Carlo simulation of the stochastic decay process. In particular, we study the asymptotic forms of the divergence of the relaxation time near the spinodal point, and reexamine its scaling properties. We further extend the analysis to transient critical phenomena such as a threshold behavior with diverging switching time under a general external driving perturbation. This models photo-excitation processes in spin-crossover materials. In the ongoing development of nano-size fabrication, such size-dependence of switching processes should be an important issue, and the properties obtained here will be applicable to a wide range of physical processes.

preprint2013arXiv

Dynamics of Desorption with Lateral Diffusion

The dynamics of desorption from a submonolayer of adsorbed atoms or ions are significantly influenced by the absence or presence of lateral diffusion of the adsorbed particles. When diffusion is present, the adsorbate configuration is simultaneously changed by two distinct processes, proceeding in parallel: adsorption/desorption, which changes the total adsorbate coverage, and lateral diffusion, which is coverage conserving. Inspired by experimental results, we here study the effects of these competing processes by kinetic Monte Carlo simulations of a simple lattice-gas model. In order to untangle the various effects, we perform large-scale simulations, in which we monitor coverage, correlation length, and cluster-size distributions, as well as the behavior of representative individual clusters, during desorption. For each initial adsorbate configuration, we perform multiple, independent simulations, without and with diffusion, respectively. We find that, compared to desorption without diffusion, the coverage-conserving diffusion process produces two competing effects: a retardation of the desorption rate, which is associated with a coarsening of the adsorbate configuration, and an acceleration due to desorption of monomers "evaporated" from the cluster perimeters. The balance between these two effects is governed by the structure of the adsorbate layer at the beginning of the desorption process. Deceleration and coarsening are predominant for configurations dominated by monomers and small clusters, while acceleration is predominant for configurations dominated by large clusters.

preprint2013arXiv

Phase separation on a hyperbolic lattice

We report a preliminary numerical study by kinetic Monte Carlo simulation of the dynamics of phase separation following a quench from high to low temperature in a system with a single, conserved, scalar order parameter (a kinetic Ising ferromagnet) confined to a hyperbolic lattice. The results are compared with simulations of the same system on two different, Euclidean lattices, in which cases we observe power-law domain growth with an exponent near the theoretically known value of 1/3. For the hyperbolic lattice we observe much slower domain growth, consistent to within our current accuracy with power-law growth with a much smaller exponent near 0.13. The paper also includes a brief introduction to non-Euclidean lattices and their mapping to the Euclidean plane.

preprint2012arXiv

Critical temperature and correlation length of an elastic interaction model for spin-crossover materials

It has previously been pointed out that the coexistence of infinite-range and short-range interactions causes a system to have a phase transition of the mean-field universality class, in which the cluster size is finite even at the critical point. In the present paper, we study this property in a model of bistable molecules, whose size changes depending on the bistable states. The molecules can move in space, interacting via an elastic interaction. It is known that due to the different sizes, an effective long-range interaction between the spins appears, and thus this model has a mean-field type of phase transition. It is found that the scaling properties of the shift of the critical temperature from the pure short-range limit in the model with infinite-range and short-range interactions hold also in the present model, regarding the ratio of the size of the two states as a control parameter for the strength of the long-range interaction. By studying the structure factor, it is shown that the dependence of the cluster size at the critical temperature also shows the same scaling properties as a previously studied model with both infinite-range and short-range interactions. We therefore conclude that these scaling relations hold universally in hybrid models with both short-range and weak long-range interactions.

preprint2012arXiv

Effects of Lateral Diffusion on the Dynamics of Desorption

The adsorbate dynamics during simultaneous action of desorption and lateral adsorbate diffusion is studied in a simple lattice-gas model by kinetic Monte Carlo simulations. It is found that the action of the coverage-conserving diffusion process during the course of the desorption has two distinct, competing effects: a general acceleration of the desorption process, and a coarsening of the adsorbate configuration through Ostwald ripening. The balance between these two effects is governed by the structure of the adsorbate layer at the beginning of the desorption process.

preprint2012arXiv

Macroscopic nucleation phenomena in continuum media with long-range interactions

Nucleation, commonly associated with discontinuous transformations between metastable and stable phases, is crucial in fields as diverse as atmospheric science and nanoscale electronics. Traditionally, it is considered a microscopic process (at most nano-meter), implying the formation of a microscopic nucleus of the stable phase. Here we show for the first time, that considering long-range interactions mediated by elastic distortions, nucleation can be a macroscopic process, with the size of the critical nucleus proportional to the total system size. This provides a new concept of "macroscopic barrier-crossing nucleation". We demonstrate the effect in molecular dynamics simulations of a model spin-crossover system with two molecular states of different sizes, causing elastic distortions.

preprint2011arXiv

Crossover between a Short-range and a Long-range Ising model

Recently, it has been found that an effective long-range interaction is realized among local bistable variables (spins) in systems where the elastic interaction causes ordering of the spins. In such systems, generally we expect both long-range and short-range interactions to exist. In the short-range Ising model, the correlation length diverges at the critical point. In contrast, in the long-range interacting model the spin configuration is always uniform and the correlation length is zero. As long as a system has non-zero long-range interactions, it shows criticality in the mean-field universality class, and the spin configuration is uniform beyond a certain scale. Here we study the crossover from the pure short-range interacting model to the long-range interacting model. We investigate the infinite-range model (Husimi-Temperley model) as a prototype of this competition, and we study how the critical temperature changes as a function of the strength of the long-range interaction. This model can also be interpreted as an approximation for the Ising model on a small-world network. We derive a formula for the critical temperature as a function of the strength of the long-range interaction. We also propose a scaling form for the spin correlation length at the critical point, which is finite as long as the long-range interaction is included, though it diverges in the limit of the pure short-range model. These properties are confirmed by extensive Monte Carlo simulations.

preprint2011arXiv

Floridian high-voltage power-grid network partitioning and cluster optimization using simulated annealing

Many partitioning methods may be used to partition a network into smaller clusters while minimizing the number of cuts needed. However, other considerations must also be taken into account when a network represents a real system such as a power grid. In this paper we use a simulated annealing Monte Carlo (MC) method to optimize initial clusters on the Florida high-voltage power-grid network that were formed by associating each load with its "closest" generator. The clusters are optimized to maximize internal connectivity within the individual clusters and minimize the power deficiency or surplus that clusters may otherwise have.

preprint2011arXiv

Parameter Estimation by Density Functional Theory for a Lattice-gas Model of Br and Cl Chemisorption on Ag(100)

We study Bromine and Chlorine chemisorption on a Ag(100) surface, using a lattice-gas model and the quantum-mechanical Density Functional Theory (DFT) method. In this model the Br and Cl ions adsorb at the fourfold hollow sites of the Ag(100) surface, which can be represented by a square lattice of adsorption sites. Five different coverages were used for each kind of adsorbate. For each adsorbate and coverage, we obtained the minimum-energy configuration, its energy, and its charge distribution. From these data we calculated dipole moments, lateral interaction energies, and binding energies. Our results show that for Br the lateral interactions obtained by fitting to the adsorption energies obtained from the DFT calculation are consistent with long-range dipole-dipole lateral interactions obtained using the dipole moments calculated from the DFT charge distribution. For Cl we found that, while the long-range dipole-dipole lateral interactions are important, short-range attractive interactions are also present. Our results are overall consistent with parameter estimates previously obtained by fitting room-temperature Monte Carlo simulations to electrochemical adsorption isotherms [I.\ Abou Hamad et al., J.\ Electroanal. Chem.\ 554 (2003), 211; Electrochim.\ Acta 50 (2005), 5518].

preprint2010arXiv

Effects of demographic stochasticity on biological community assembly on evolutionary time scales

We study the effects of demographic stochasticity on the long-term dynamics of biological coevolution models of community assembly. The noise is induced in order to check the validity of deterministic population dynamics. While mutualistic communities show little dependence on the stochastic population fluctuations, predator-prey models show strong dependence on the stochasticity, indicating the relevance of the finiteness of the populations. For a predator-prey model, the noise causes drastic decreases in diversity and total population size. The communities that emerge under influence of the noise consist of species strongly coupled with each other and have stronger linear stability around the fixed-point populations than the corresponding noiseless model. The dynamics on evolutionary time scales for the predator-prey model are also altered by the noise. Approximate $1/f$ fluctuations are observed with noise, while $1/f^{2}$ fluctuations are found for the model without demographic noise.

preprint2010arXiv

Positive interactions and the emergence of community structure in metacommunities

The significant role of space in maintaining species coexistence and determining community structure and function is well established. However, community ecology studies have mainly focused on simple competition and predation systems, and the relative impact of positive interspecific interactions in shaping communities in a spatial context is not well understood. Here we employ a spatially explicit metacommunity model to investigate the effect of local dispersal on the structure and function of communities in which species are linked through an interaction web comprising mutualism, competition and exploitation. Our results show that function, diversity and interspecific interactions of locally linked communities undergo a phase transition with changes in the rate of species dispersal. We find that low spatial interconnectedness favors the spontaneous emergence of strongly mutualistic communities which are more stable but less productive and diverse. On the other hand, high spatial interconnectedness promotes local biodiversity at the expense of local stability and supports communities with a wide range of interspecific interactions. We argue that investigations of the relationship between spatial processes and the self-organization of complex interaction webs are critical to understanding the geographic structure of interactions in real landscapes.

preprint2010arXiv

Random walk in genome space : A key ingredient of intermittent dynamics of community assembly on evolutionary time scales

Community assembly is studied using individual-based multispecies models. The models have stochastic population dynamics with mutation, migration, and extinction of species. Mutants appear as a result of mutation of the resident species, while migrants have no correlation with the resident species. It is found that the dynamics of community assembly with mutations are quite different from the case with migrations. In contrast to mutation models, which show intermittent dynamics of quasi-steady states interrupted by sudden reorganizations of the community, migration models show smooth and gradual renewal of the community. As a consequence, instead of the $1/f$ diversity fluctuations found for the mutation models, $1/f^2$, random-walk like fluctuations are observed for the migration models. In addition, a characteristic species-lifetime distribution is found: a power law that is cut off by a "skewed" distribution in the long-lifetime regime. The latter has a longer tail than a simple exponential function, which indicates an age-dependent species-mortality function. Since this characteristic profile has been observed, both in fossil data and in several other mathematical models, we conclude that it is a universal feature of macroevolution.

preprint2010arXiv

Spectral matrix methods for partitioning power grids: Applications to the Italian and Floridian high-voltage networks

Intentional islanding is used to limit cascading power failures by isolating highly connected "islands" with local generating capacity. To efficiently isolate an island, one should break as few power lines as possible. This is a graph partitioning problem, and here we give preliminary results on islanding of the Italian and Floridian high-voltage grids by spectral matrix methods.

preprint2009arXiv

Degree Correlations in a Dynamically Generated Model Food Web

We explore aspects of the community structures generated by a simple predator-prey model of biological coevolution, using large-scale kinetic Monte Carlo simulations. The model accounts for interspecies and intraspecies competition for resources, as well as adaptive foraging behavior. It produces a metastable low-diversity phase and a stable high-diversity phase. The structures and joint indegree-outdegree distributions of the food webs generated in the latter phase are discussed.

preprint2009arXiv

Threshold Phenomena under Photo Excitation of Spin-crossover Materials with Cooperativity due to Elastic Interactions

Photo-induced switching from the low-spin state to the high-spin state is studied in a model of spin-crossover materials, in which long-range interactions are induced by elastic distortions due to different molecular sizes the two spin states. At a threshold value of the light intensity we observe nonequilibrium critical behavior corresponding to a mean-field spinodal point. Finite-size scaling of the divergence of the relaxation time is revealed by analysis of kinetic Monte Carlo simulations.

preprint2008arXiv

Are Genetically Robust Regulatory Networks Dynamically Different from Random Ones?

We study a genetic regulatory network model developed to demonstrate that genetic robustness can evolve through stabilizing selection for optimal phenotypes. We report preliminary results on whether such selection could result in a reorganization of the state space of the system. For the chosen parameters, the evolution moves the system slightly toward the more ordered part of the phase diagram. We also find that strong memory effects cause the Derrida annealed approximation to give erroneous predictions about the model's phase diagram.

preprint2008arXiv

Effects of stochastic population fluctuations in two models of biological macroevolution

Two mathematical models of macroevolution are studied. These models have population dynamics at the species level, and mutations and extinction of species are also included. The population dynamics are updated by difference equations with stochastic noise terms that characterize population fluctuations. The effects of the stochastic population fluctuations on diversity and total population sizes on evolutionary time scales are studied. In one model, species can make either predator-prey, mutualistic, or competitive interactions, while the other model allows only predator-prey interactions. When the noise in the population dynamics is strong enough, both models show intermittent behavior and their power spectral densities show approximate $1/f$ fluctuations. In the noiseless limit, the two models have different power spectral densities. For the predator-prey model, $1/f^2$ fluctuations appears, indicating random-walk like behavior, while the other model still shows $1/f$ noise. These results indicate that stochastic population fluctuations may significantly affect long-time evolutionary dynamics.

preprint2007arXiv

Self-optimization, community stability, and fluctuations in two individual-based models of biological coevolution

We compare and contrast the long-time dynamical properties of two individual-based models of biological coevolution. Selection occurs via multispecies, stochastic population dynamics with reproduction probabilities that depend nonlinearly on the population densities of all species resident in the community. New species are introduced through mutation. Both models are amenable to exact linear stability analysis, and we compare the analytic results with large-scale kinetic Monte Carlo simulations, obtaining the population size as a function of an average interspecies interaction strength. Over time, the models self-optimize through mutation and selection to approximately maximize a community fitness function, subject only to constraints internal to the particular model. If the interspecies interactions are randomly distributed on an interval including positive values, the system evolves toward self-sustaining, mutualistic communities. In contrast, for the predator-prey case the matrix of interactions is antisymmetric, and a nonzero population size must be sustained by an external resource. Time series of the diversity and population size for both models show approximate 1/f noise and power-law distributions for the lifetimes of communities and species. For the mutualistic model, these two lifetime distributions have the same exponent, while their exponents are different for the predator-prey model. The difference is probably due to greater resilience toward mass extinctions in the food-web like communities produced by the predator-prey model.

preprint2002arXiv

Ab initio calculations for bromine adlayers on the Ag(100) and Au(100) surfaces: the c(2x2) structure

Ab initio total-energy density-functional methods with supercell models have been employed to calculate the c(2x2) structure of the Br-adsorbed Ag(100) and Au(100) surfaces. The atomic geometries of the surfaces and the preferred bonding sites of the bromine have been determined. The bonding character of bromine with the substrates has also been studied by analyzing the electronic density of states and the charge transfer. The calculations show that while the four-fold hollow-site configuration is more stable than the two-fold bridge-site topology on the Ag(100) surface, bromine prefers the bridge site on the Au(100) surface. The one-fold on-top configuration is the least stable configuration on both surfaces. It is also observed that the second layer of the Ag substrate undergoes a small buckling as a consequence of the adsorption of Br. Our results provide a theoretical explanation for the experimental observations that the adsorption of bromine on the Ag(100) and Au(100) surfaces results in different bonding configurations.

preprint1998arXiv

Test of the Kolmogorov-Johnson-Mehl-Avrami picture of metastable decay in a model with microscopic dynamics

The Kolmogorov-Johnson-Mehl-Avrami (KJMA) theory for the time evolution of the order parameter in systems undergoing first-order phase transformations has been extended by Sekimoto to the level of two-point correlation functions. Here, this extended KJMA theory is applied to a kinetic Ising lattice-gas model, in which the elementary kinetic processes act on microscopic length and time scales. The theoretical framework is used to analyze data from extensive Monte Carlo simulations. The theory is inherently a mesoscopic continuum picture, and in principle it requires a large separation between the microscopic scales and the mesoscopic scales characteristic of the evolving two-phase structure. Nevertheless, we find excellent quantitative agreement with the simulations in a large parameter regime, extending remarkably far towards strong fields (large supersaturations) and correspondingly small nucleation barriers. The original KJMA theory permits direct measurement of the order parameter in the metastable phase, and using the extension to correlation functions one can also perform separate measurements of the nucleation rate and the average velocity of the convoluted interface between the metastable and stable phase regions. The values obtained for all three quantities are verified by other theoretical and computational methods. As these quantities are often difficult to measure directly during a process of phase transformation, data analysis using the extended KJMA theory may provide a useful experimental alternative.

preprint1996arXiv

Numerical simulations of scattering speckle from phase ordering systems

The scattering of coherent X-rays from dynamically evolving systems is currently becoming experimentally feasible. The scattered beam produces a pattern of bright and dark speckles, which fluctuate almost independently in time and can be used to study the dynamics of the system. Here we report large-scale computer simulations of the speckle dynamics for a phase ordering system, using a two-dimensional model quenched through an order--disorder transition into the two-phase regime. The intensity at each wave vector ${\bf k}$ is found to be an exponentially distributed random variable. The scaling hypothesis is extended to the two-time correlation function of the scattering intensity at a given wave vector, $Corr({\bf k};t_1,t_2)$. The characteristic decay time difference for the correlation function, $|t_1 - t_2|_c$, is found to scale as $k|t_1 + t_2|^{1/2}$.