Hybrid evolutionary oligopolies and the dynamics of corporate social responsibility

The diffusion of corporate social responsibility is investigated by employing a hybrid evolutionary game where a firm chooses between being either socially responsible, which implies devoting a fraction of its profit to social projects, or non-socially responsible. Consumers prize socially responsible companies by paying a higher reservation price for their products. The hybrid evolutionary framework is characterized by a quantity dynamics that describes the oligopolistic competition given firms’ belief about the composition of the industry. At regular intervals of time, this belief is endogenously updated by a retrospective comparison on the profits obtained and on the basis of an evolutionary mechanism. Assuming that firms are Nash players, that is at each instant of time they produce the Nash equilibrium-in-belief quantity, the investigation of the model reveals that an industry homogeneously populated by socially responsible firms is a stable equilibrium when the fraction of profits earmarked for socially responsible activities is sufficiently limited. However, the extra marginal profits of a socially responsible firm are reduced when the number of competitors increases, impeding the diffusion of socially responsible companies. In particular, the trade-off between a higher net margin on sales obtained by socially responsible firms and a lower level of production that reduces the profit gap between a socially responsible firm and the rest of the market shows that an increased size of the industry favors mixed oligopolies. Moreover, imposing the hypothesis of neutrality of CSR activities, the model reveals that being socially responsible is an evolutionarily stable strategy for firms and is convenient for customers. Relaxing the hypothesis of Nash players by introducing boundedly rational firms that decide their level of production according to a partial adjustment toward the best reply, the robustness of these results is confirmed

Quasi-equilibrium binary black hole sequences for puncture data derived from helical Killing vector conditions

We construct a sequence of binary black hole puncture data derived under the assumptions (i) that the ADM mass of each puncture as measured in the asymptotically flat space at the puncture stays constant along the sequence, and (ii) that the orbits along the sequence are quasi-circular in the sense that several necessary conditions for the existence of a helical Killing vector are satisfied. These conditions are equality of ADM and Komar mass at infinity and equality of the ADM and a rescaled Komar mass at each puncture. In this paper we explicitly give results for the case of an equal mass black hole binary without spin, but our approach can also be applied in the general case. We find that up to numerical accuracy the apparent horizon mass also remains constant along the sequence and that the prediction for the innermost stable circular orbit is similar to what has been found with the effective potential method.Comment: 6 pages, 3 figures, 1 tabl

A new numerical method to construct binary neutron star initial data

We present a new numerical method for the generation of binary neutron star initial data using a method along the lines of the the Wilson-Mathews or the closely related conformal thin sandwich approach. Our method uses six different computational domains, which include spatial infinity. Each domain has its own coordinates which are chosen such that the star surfaces always coincide with domain boundaries. These properties facilitate the imposition of boundary conditions. Since all our fields are smooth inside each domain, we are able to use an efficient pseudospectral method to solve the elliptic equations associated with the conformal thin sandwich approach. Currently we have implemented corotating configurations with arbitrary mass ratios, but an extension to arbitrary spins is possible. The main purpose of this paper is to introduce our new method and to test our code for several different configurations.Comment: 18 pages, 8 figures, 1 tabl

Black hole puncture initial data with realistic gravitational wave content

We present improved post-Newtonian-inspired initial data for non-spinning black-hole binaries, suitable for numerical evolution with punctures. We revisit the work of Tichy et al. [W. Tichy, B. Bruegmann, M. Campanelli, and P. Diener, Phys. Rev. D 67, 064008 (2003)], explicitly calculating the remaining integral terms. These terms improve accuracy in the far zone and, for the first time, include realistic gravitational waves in the initial data. We investigate the behavior of these data both at the center of mass and in the far zone, demonstrating agreement of the transverse-traceless parts of the new metric with quadrupole-approximation waveforms. These data can be used for numerical evolutions, enabling a direct connection between the merger waveforms and the post-Newtonian inspiral waveforms.Comment: 13 pages, 7 figures; replaced with published versio

A minimal no-radiation approximation to Einstein's field equations

An approximation to Einstein's field equations in Arnowitt-Deser-Misner (ADM) canonical formalism is presented which corresponds to the magneto-hydrodynamics (MHD) approximation in electrodynamics. It results in coupled elliptic equations which represent the maximum of elliptic-type structure of Einstein's theory and naturally generalizes previous conformal-flat truncations of the theory. The Hamiltonian, in this approximation, is identical with the non-dissipative part of the Einsteinian one through the third post-Newtonian order. The proposed scheme, where stationary spacetimes are exactly reproduced, should be useful to construct {\em realistic} initial data for general relativistic simulations as well as to model astrophysical scenarios, where gravitational radiation reaction can be neglected.Comment: 9 page

High-accuracy high-mass ratio simulations for binary neutron stars and their comparison to existing waveform models

The subsequent observing runs of the advanced gravitational-wave detector network will likely provide us with various gravitational-wave observations of binary neutron star systems. For an accurate interpretation of these detections, we need reliable gravitational-wave models. To test and to point out how existing models could be improved, we perform a set of high-resolution numerical-relativity simulations for four different physical setups with mass ratios $q$ = $1.25$, $1.50$, $1.75$, $2.00$, and total gravitational mass $M = 2.7M_\odot$ . Each configuration is simulated with five different resolutions to allow a proper error assessment. Overall, we find approximately 2nd order converging results for the dominant $(2,2)$, but also subdominant $(2,1)$, $(3,3)$, $(4,4)$ modes, while, generally, the convergence order reduces slightly for an increasing mass ratio. Our simulations allow us to validate waveform models, where we find generally good agreement between state-of-the-art models and our data, and to prove that scaling relations for higher modes currently employed for binary black hole waveform modeling also apply for the tidal contribution. Finally, we also test if the current NRTidal model to describe tidal effects is a valid description for high-mass ratio systems. We hope that our simulation results can be used to further improve and test waveform models in preparation for the next observing runs

Constructing binary neutron star initial data with high spins, high compactnesses, and high mass ratios

The construction of accurate and consistent initial data for various binary parameters is a critical ingredient for numerical relativity simulations of the compact binary coalescence. In this article, we present an upgrade of the pseudospectral SGRID code, which enables us to access even larger regions of the binary neutron star parameter space. As a proof of principle, we present a selected set of first simulations based on initial configurations computed with the new code version. In particular, we simulate two millisecond pulsars close to their breakup spin, highly compact neutron stars with masses at about $98\%$ of the maximum supported mass of the employed equation of state, and an unequal mass systems with mass ratios even outside the range predicted by population synthesis models ($q = 2.03$). The discussed code extension will help us to simulate previously unexplored binary configurations. This is a necessary step to construct and test new gravitational wave approximants and to interpret upcoming binary neutron star merger observations. When we construct initial data, one has to specify various parameters, such as a rotation parameter for each star. Some of these parameters do not have direct physical meaning, which makes comparisons with other methods or models difficult. To facilitate this, we introduce simple estimates for the initial spin, momentum, mass, and center of mass of each individual star

Developing a Predictive Metric to Assess School Viability

This article examines a wide range of parish school indicators that can be used to predict long-term viability

Feasibility of approximating spatial and local entanglement in long-range interacting systems using the extended Hubbard model

We investigate the extended Hubbard model as an approximation to the local and spatial entanglement of a one-dimensional chain of nanostructures where the particles interact via a long range interaction represented by a `soft' Coulomb potential. In the process we design a protocol to calculate the particle-particle spatial entanglement for the Hubbard model and show that, in striking contrast with the loss of spatial degrees of freedom, the predictions are reasonably accurate. We also compare results for the local entanglement with previous results found using a contact interaction (PRA, 81 (2010) 052321) and show that while the extended Hubbard model recovers a better agreement with the entanglement of a long-range interacting system, there remain realistic parameter regions where it fails to predict the quantitative and qualitative behaviour of the entanglement in the nanostructure system.Comment: 6 pages, 5 figures and 1 table; added results with correlated hopping term; accepted by EP

Binary black hole initial data from matched asymptotic expansions

We present an approximate metric for a binary black hole spacetime to construct initial data for numerical relativity. This metric is obtained by asymptotically matching a post-Newtonian metric for a binary system to a perturbed Schwarzschild metric for each hole. In the inner zone near each hole, the metric is given by the Schwarzschild solution plus a quadrupolar perturbation corresponding to an external tidal gravitational field. In the near zone, well outside each black hole but less than a reduced wavelength from the center of mass of the binary, the metric is given by a post-Newtonian expansion including the lowest-order deviations from flat spacetime. When the near zone overlaps each inner zone in a buffer zone, the post-Newtonian and perturbed Schwarzschild metrics can be asymptotically matched to each other. By demanding matching (over a 4-volume in the buffer zone) rather than patching (choosing a particular 2-surface in the buffer zone), we guarantee that the errors are small in all zones. The resulting piecewise metric is made formally $C^\infty$ with smooth transition functions so as to obtain the finite extrinsic curvature of a 3-slice. In addition to the metric and extrinsic curvature, we present explicit results for the lapse and the shift, which can be used as initial data for numerical simulations. This initial data is not accurate all the way to the asymptotically flat ends inside each hole, and therefore must be used with evolution codes which employ black hole excision rather than puncture methods. This paper lays the foundations of a method that can be sraightforwardly iterated to obtain initial data to higher perturbative order.Comment: 24 pages, 15 figures. Replaced with published version. Major editing of text, no major change to the physic