Towards More Accurate Molecular Dynamics Calculation of Thermal Conductivity. Case Study: GaN Bulk Crystals

Significant differences exist among literature for thermal conductivity of various systems computed using molecular dynamics simulation. In some cases, unphysical results, for example, negative thermal conductivity, have been found. Using GaN as an example case and the direct non-equilibrium method, extensive molecular dynamics simulations and Monte Carlo analysis of the results have been carried out to quantify the uncertainty level of the molecular dynamics methods and to identify the conditions that can yield sufficiently accurate calculations of thermal conductivity. We found that the errors of the calculations are mainly due to the statistical thermal fluctuations. Extrapolating results to the limit of an infinite-size system tend to magnify the errors and occasionally lead to unphysical results. The error in bulk estimates can be reduced by performing longer time averages using properly selected systems over a range of sample lengths. If the errors in the conductivity estimates associated with each of the sample lengths are kept below a certain threshold, the likelihood of obtaining unphysical bulk values becomes insignificant. Using a Monte-Carlo approach developed here, we have determined the probability distributions for the bulk thermal conductivities obtained using the direct method. We also have observed a nonlinear effect that can become a source of significant errors. For the extremely accurate results presented here, we predict a [0001] GaN thermal conductivity of 185 $\rm{W/K \cdot m}$ at 300 K, 102 $\rm{W/K \cdot m}$ at 500 K, and 74 $\rm{W/K \cdot m}$ at 800 K. Using the insights obtained in the work, we have achieved a corresponding error level (standard deviation) for the bulk (infinite sample length) GaN thermal conductivity of less than 10 $\rm{W/K \cdot m}$, 5 $\rm{W/K \cdot m}$, and 15 $\rm{W/K \cdot m}$ at 300 K, 500 K, and 800 K respectively

Effective Free Energy for Individual Dynamics

Physics and economics are two disciplines that share the common challenge of linking microscopic and macroscopic behaviors. However, while physics is based on collective dynamics, economics is based on individual choices. This conceptual difference is one of the main obstacles one has to overcome in order to characterize analytically economic models. In this paper, we build both on statistical mechanics and the game theory notion of Potential Function to introduce a rigorous generalization of the physicist's free energy, which includes individual dynamics. Our approach paves the way to analytical treatments of a wide range of socio-economic models and might bring new insights into them. As first examples, we derive solutions for a congestion model and a residential segregation model.Comment: 8 pages, 2 figures, presented at the ECCS'10 conferenc

Dynamic scaling regimes of collective decision making

We investigate a social system of agents faced with a binary choice. We assume there is a correct, or beneficial, outcome of this choice. Furthermore, we assume agents are influenced by others in making their decision, and that the agents can obtain information that may guide them towards making a correct decision. The dynamic model we propose is of nonequilibrium type, converging to a final decision. We run it on random graphs and scale-free networks. On random graphs, we find two distinct regions in terms of the "finalizing time" -- the time until all agents have finalized their decisions. On scale-free networks on the other hand, there does not seem to be any such distinct scaling regions

On Spatial Consensus Formation: Is the Sznajd Model Different from a Voter Model?

In this paper, we investigate the so-called ``Sznajd Model'' (SM) in one dimension, which is a simple cellular automata approach to consensus formation among two opposite opinions (described by spin up or down). To elucidate the SM dynamics, we first provide results of computer simulations for the spatio-temporal evolution of the opinion distribution $L(t)$, the evolution of magnetization $m(t)$, the distribution of decision times $P(\tau)$ and relaxation times $P(\mu)$. In the main part of the paper, it is shown that the SM can be completely reformulated in terms of a linear VM, where the transition rates towards a given opinion are directly proportional to frequency of the respective opinion of the second-nearest neighbors (no matter what the nearest neighbors are). So, the SM dynamics can be reduced to one rule, ``Just follow your second-nearest neighbor''. The equivalence is demonstrated by extensive computer simulations that show the same behavior between SM and VM in terms of $L(t)$, $m(t)$, $P(\tau)$, $P(\mu)$, and the final attractor statistics. The reformulation of the SM in terms of a VM involves a new parameter $\sigma$, to bias between anti- and ferromagnetic decisions in the case of frustration. We show that $\sigma$ plays a crucial role in explaining the phase transition observed in SM. We further explore the role of synchronous versus asynchronous update rules on the intermediate dynamics and the final attractors. Compared to the original SM, we find three additional attractors, two of them related to an asymmetric coexistence between the opposite opinions.Comment: 22 pages, 20 figures. For related publications see http://www.ais.fraunhofer.de/~fran

Enhanced Anandamide Plasma Levels in Patients with Complex Regional Pain Syndrome following Traumatic Injury: A Preliminary Report

The complex regional pain syndrome (CRPS) is a disabling neuropathic pain condition that may develop following injuries of the extremities. The pathogenesis of this syndrome is not clear; however, it includes complex interactions between the nervous and the immune system resulting in chronic inflammation, pain and trophic changes. This interaction may be mediated by chronic stress which is thought to activate the endogenous cannabinoid (endocannabinoid) system (ECS). We conducted an open, prospective, comparative clinical study to determine plasma level of the endocannabinoid anandamide by high-performance liquid chromatography and a tandem mass spectrometry system in 10 patients with CRPS type I versus 10 age- and sex-matched healthy controls. As compared to healthy controls, CRPS patients showed significantly higher plasma concentrations of anandamide. These results indicate that the peripheral ECS is activated in CRPS. Further studies are warranted to evaluate the role of the ECS in the limitation of inflammation and pain. Copyright (C) 2009 S. Karger AG, Base

Reinforced communication and social navigation generate groups in model networks

To investigate the role of information flow in group formation, we introduce a model of communication and social navigation. We let agents gather information in an idealized network society, and demonstrate that heterogeneous groups can evolve without presuming that individuals have different interests. In our scenario, individuals' access to global information is constrained by local communication with the nearest neighbors on a dynamic network. The result is reinforced interests among like-minded agents in modular networks; the flow of information works as a glue that keeps individuals together. The model explains group formation in terms of limited information access and highlights global broadcasting of information as a way to counterbalance this fragmentation. To illustrate how the information constraints imposed by the communication structure affects future development of real-world systems, we extrapolate dynamics from the topology of four social networks.Comment: 7 pages, 3 figure

Statistical physics of the Schelling model of segregation

We investigate the static and dynamic properties of a celebrated model of social segregation, providing a complete explanation of the mechanisms leading to segregation both in one- and two-dimensional systems. Standard statistical physics methods shed light on the rich phenomenology of this simple model, exhibiting static phase transitions typical of kinetic constrained models, nontrivial coarsening like in driven-particle systems and percolation-related phenomena.Comment: 4 pages, 3 figure

Role of Activity in Human Dynamics

The human society is a very complex system; still, there are several non-trivial, general features. One type of them is the presence of power-law distributed quantities in temporal statistics. In this Letter, we focus on the origin of power-laws in rating of movies. We present a systematic empirical exploration of the time between two consecutive ratings of movies (the interevent time). At an aggregate level, we find a monotonous relation between the activity of individuals and the power-law exponent of the interevent-time distribution. At an individual level, we observe a heavy-tailed distribution for each user, as well as a negative correlation between the activity and the width of the distribution. We support these findings by a similar data set from mobile phone text-message communication. Our results demonstrate a significant role of the activity of individuals on the society-level patterns of human behavior. We believe this is a common character in the interest-driven human dynamics, corresponding to (but different from) the universality classes of task-driven dynamics.Comment: 5 pages, 6 figures. Accepted by EP

Quantum Games

In these lecture notes we investigate the implications of the identification of strategies with quantum operations in game theory beyond the results presented in [J. Eisert, M. Wilkens, and M. Lewenstein, Phys. Rev. Lett. 83, 3077 (1999)]. After introducing a general framework, we study quantum games with a classical analogue in order to flesh out the peculiarities of game theoretical settings in the quantum domain. Special emphasis is given to a detailed investigation of different sets of quantum strategies.Comment: 13 pages (LaTeX), 3 figure