MPI parallel programming of mixed integer optimization problems using CPLEX with COIN-OR

The aim of this technical report is to present some detailed explanations in order to help to understand and use the Message Passing Interface (MPI) parallel programming for solving several mixed integer optimization problems. We have developed a C++ experimental code that uses the IBM ILOG CPLEX optimizer within the COmputational INfrastructure for Operations Research (COIN-OR) and MPI parallel computing for solving the optimization models under UNIX-like systems. The computational experience illustrates how can we solve 44 optimization problems which are asymmetric with respect to the number of integer and continuous variables and the number of constraints. We also report a comparative with the speedup and efficiency of several strategies implemented for some available number of threads.

Scenario Cluster Lagrangian Decomposition in two stage stochastic mixed 0-1 optimization

In this paper we introduce four scenario Cluster based Lagrangian Decomposition (CLD) procedures for obtaining strong lower bounds to the (optimal) solution value of two-stage stochastic mixed 0-1 problems. At each iteration of the Lagrangian based procedures, the traditional aim consists of obtaining the solution value of the corresponding Lagrangian dual via solving scenario submodels once the nonanticipativity constraints have been dualized. Instead of considering a splitting variable representation over the set of scenarios, we propose to decompose the model into a set of scenario clusters. We compare the computational performance of the four Lagrange multiplier updating procedures, namely the Subgradient Method, the Volume Algorithm, the Progressive Hedging Algorithm and the Dynamic Constrained Cutting Plane scheme for different numbers of scenario clusters and different dimensions of the original problem. Our computational experience shows that the CLD bound and its computational effort depend on the number of scenario clusters to consider. In any case, our results show that the CLD procedures outperform the traditional LD scheme for single scenarios both in the quality of the bounds and computational effort. All the procedures have been implemented in a C++ experimental code. A broad computational experience is reported on a test of randomly generated instances by using the MIP solvers COIN-OR and CPLEX for the auxiliary mixed 0-1 cluster submodels, this last solver within the open source engine COIN-OR. We also give computational evidence of the model tightening effect that the preprocessing techniques, cut generation and appending and parallel computing tools have in stochastic integer optimization. Finally, we have observed that the plain use of both solvers does not provide the optimal solution of the instances included in the testbed with which we have experimented but for two toy instances in affordable elapsed time. On the other hand the proposed procedures provide strong lower bounds (or the same solution value) in a considerably shorter elapsed time for the quasi-optimal solution obtained by other means for the original stochastic problem.

Reference Points Based on Dynamic Optimisation: A Versatil Algorithm for Mixed Fishery Management with Bio-economic Agestructured Models

Single-species management objectives may not be consistent within mixed fisheries. They may lead species to unsafe situations, promote discarding of over-quota and/or misreporting of catches. We provide an algorithm for characterising bio-economic reference points for a mixed fishery as the steady-state solution of a dynamic optimal management problem. The optimisation problem takes into account: i) that species are fishing simultaneously in unselective fishing operations and ii)intertemporal discounting and fleet costs to relate reference points to discounted economic profits along optimal trajectories. We illustrate how the algorithm can be implemented by applying it to the European Northern Stock of Hake (Merluccius merluccius), where fleets also capture Northern megrim (Lepidorhombus whiffiagonis) and Northern anglerfish (Lophius piscatorius and Lophius budegassa). We find that optimal mixed management leads to a target reference point that is quite similar to the 2/3 of the Fmsy single-species (hake) target. Mixed management is superior to singlespecies management because it leads the fishery to higher discounted profits with higher long-term SSB for all species. We calculate that the losses due to the use of the Fmsy single-species (hake) target in this mixed fishery account for 11.4% of total discounted profits.

Generating cluster submodels from a multistage stochastic mixed integer optimization model using break stage

We present a scheme to generate clusters submodels with stage ordering from a (symmetric or a nonsymmetric one) multistage stochastic mixed integer optimization model using break stage. We consider a stochastic model in compact representation and MPS format with a known scenario tree. The cluster submodels are built by storing first the 0-1 the variables, stage by stage, and then the continuous ones, also stage by stage. A C++ experimental code has been implemented for reordering the stochastic model as well as the cluster decomposition after the relaxation of the non-anticipativiy constraints until the so-called breakstage. The computational experience shows better performance of the stage ordering in terms of elapsed time in a randomly generated testbed of multistage stochastic mixed integer problems.

Long-lived oscillatory incoherent electron dynamics in molecules: trans-polyacetylene oligomers

We identify an intriguing feature of the electron-vibrational dynamics of molecular systems via a computational examination of trans-polyacetylene oligomers. Here, via the vibronic interactions, the decay of an electron in the conduction band resonantly excites an electron in the valence band, and vice versa, leading to oscillatory exchange of electronic population between two distinct electronic states that lives for up to tens of picoseconds. The oscillatory structure is reminiscent of beating patterns between quantum states and is strongly suggestive of the presence of long-lived molecular electronic coherence. Significantly, however, a detailed analysis of the electronic coherence properties shows that the oscillatory structure arises from a purely incoherent process. These results were obtained by propagating the coupled dynamics of electronic and vibrational degrees of freedom in a mixed quantum-classical study of the Su-Schrieffer-Heeger Hamiltonian for polyacetylene. The incoherent process is shown to occur between degenerate electronic states with distinct electronic configurations that are indirectly coupled via a third auxiliary state by vibronic interactions. A discussion of how to construct electronic superposition states in molecules that are truly robust to decoherence is also presented

Some experiments on solving multistage stochastic mixed 0-1 programs with time stochastic dominance constraints

In this work we extend to the multistage case two recent risk averse measures for two-stage stochastic programs based on first- and second-order stochastic dominance constraints induced by mixed-integer linear recourse. Additionally, we consider Time Stochastic Dominance (TSD) along a given horizon. Given the dimensions of medium-sized problems augmented by the new variables and constraints required by those risk measures, it is unrealistic to solve the problem up to optimality by plain use of MIP solvers in a reasonable computing time, at least. Instead of it, decomposition algorithms of some type should be used. We present an extension of our Branch-and-Fix Coordination algorithm, so named BFC-TSD, where a special treatment is given to cross scenario group constraints that link variables from different scenario groups. A broad computational experience is presented by comparing the risk neutral approach and the tested risk averse strategies. The performance of the new version of the BFC algorithm versus the plain use of a state-of-the-artMIP solver is also reported.

Efficiency vs. Stability in a Mixed Network Formation Model

Documentos de Trabajo

Metallic thin fi lms on stepped surfaces: lateral scattering of quantum well states

Quantum well states of Ag films grown on stepped Au(111) surfaces are shown to undergo lateral scattering, in analogy with surface states of vicinal Ag(111). Applying angle resolved photoemission spectroscopy we observe quantum well bands with zone-folding and gap openings driven by surface/interface step lattice scattering. Experiments performed on a curved Au(111) substrate allow us to determine a subtle terrace-size effect, i.e., a fine step-density-dependent upward shift of quantum well bands. This energy shift is explained as mainly due to the periodically stepped crystal potential offset at the interface side of the film. Finally, the surface state of the stepped Ag film is analyzed with both photoemission and scanning tunneling microscopy. We observe that the stepped film interface also affects the surface state energy, which exhibits a larger terrace-size effect compared to surface states of bulk vicinal Ag(111) crystals

Spin-helical Dirac states in graphene induced by polar-substrate surfaces with giant spin-orbit interaction: a new platform for spintronics

Spintronics, or spin electronics, is aimed at efficient control and manipulation of spin degrees of freedom in electron systems. To comply with demands of nowaday spintronics, the studies of electron systems hosting giant spin-orbit-split electron states have become one of the most important problems providing us with a basis for desirable spintronics devices. In construction of such devices, it is also tempting to involve graphene, which has attracted great attention because of its unique and remarkable electronic properties and was recognized as a viable replacement for silicon in electronics. In this case, a challenging goal is to lift spin degeneracy of graphene Dirac states. Here, we propose a novel pathway to achieve this goal by means of coupling of graphene and polar-substrate surface states with giant Rashba-type spin-splitting. We theoretically demonstrate it by constructing the graphene@BiTeCl system, which appears to possess spin-helical graphene Dirac states caused by the strong interaction of Dirac and Rashba electrons. We anticipate that our findings will stimulate rapid growth in theoretical and experimental investigations of graphene Dirac states with real spin-momentum locking, which can revolutionize the graphene spintronics and become a reliable base for prospective spintronics applications.

Coincidence point theorems in quasi-metric spaces without assuming the mixed monotone property and consequences in G-metric spaces

In this paper, we present some coincidence point theorems in the setting of quasi-metric spaces that can be applied to operators which not necessarily have the mixed monotone property. As a consequence, we particularize our results to the field of metric spaces, partially ordered metric spaces and G-metric spaces, obtaining some very recent results. Finally, we show how to use our main theorems to obtain coupled, tripled, quadrupled and multidimensional coincidence point results.

Microwave absorption properties of permalloy nanodots in the vortex and quasi-uniform magnetization states

When the in-plane bias magnetic field acting on a flat circular magnetic dot is smaller than the saturation field, there are two stable competing magnetization configurations of the dot: the vortex and the quasi-uniform (C-state). We measured microwave absorption properties in an array of non-interacting permalloy dots in the frequency range 1-8 GHz when the in-plane bias magnetic field was varied in the region of the dot magnetization state bi-stability. We found that the microwave absorption properties in the vortex and quasi-uniform stable states are substantially different, so that switching between these states in a fixed bias field can be used for the development of reconfigurable microwave magnetic materials.

Analysis of the differences between the rates of inflation associated with two aggregate price level: The case of the euro area and the United States

Idioma: Inglés Abstract: This project focuses on two indicators of prices, the GDP deflator and the consumer price index (CPI), and analyzes the differences and similarities they present. These price indexes have been chosen taking into account its great representativeness and importance to economic and social level, and its direct relationship to the overall functioning of the economy and, in particular, inflation. It should be also mentioned that this study was conducted for cases of the euro area and the United States, as the impact of these economies in the economic and social situation at international level is very significant.

One-dimensional potential for image-potential states on graphene

In the framework of dielectric theory, the static non-local self-energy of an electron near an ultra-thin polarizable layer has been calculated and applied to study binding energies of image-potential states near free-standing graphene. The corresponding series of eigenvalues and eigenfunctions have been obtained by numerically solving the one-dimensional Schrodinger equation. The imagepotential state wave functions accumulate most of their probability outside the slab. We find that the random phase approximation (RPA) for the nonlocal dielectric function yields a superior description for the potential inside the slab, but a simple Fermi-Thomas theory can be used to get a reasonable quasi-analytical approximation to the full RPA result that can be computed very economically. Binding energies of the image-potential states follow a pattern close to the Rydberg series for a perfect metal with the addition of intermediate states due to the added symmetry of the potential. The formalism only requires a minimal set of free parameters: the slab width and the electronic density. The theoretical calculations are compared with experimental results for the work function and image-potential states obtained by two-photon photoemission.

Efficacy of a cognitive and behavioural psychotherapy applied by primary care psychologists in patients with mixed anxiety-depressive disorder: a research protocol

Background: In contrast with the recommendations of clinical practice guidelines, the most common treatment for anxiety and depressive disorders in primary care is pharmacological. The aim of this study is to assess the efficacy of a cognitive-behavioural psychological intervention, delivered by primary care psychologists in patients with mixed anxiety-depressive disorder compared to usual care. Methods/Design: This is an open-label, multicentre, randomized, and controlled study with two parallel groups. A random sample of 246 patients will be recruited with mild-to-moderate mixed anxiety-depressive disorder, from the target population on the lists of 41 primary care doctors. Patients will be randomly assigned to the intervention group, who will receive standardised cognitive-behavioural therapy delivered by psychologists together with usual care, or to a control group, who will receive usual care alone. The cognitive-behavioural therapy intervention is composed of eight individual 60-minute face-to face sessions conducted in eight consecutive weeks. A follow-up session will be conducted over the telephone, for reinforcement or referral as appropriate, 6 months after the intervention, as required. The primary outcome variable will be the change in scores on the Short Form-36 General Health Survey. We will also measure the change in the frequency and intensity of anxiety symptoms (State-Trait Anxiety Inventory) and depression (Beck Depression Inventory) at baseline, and 3, 6 and 12 months later. Additionally, we will collect information on the use of drugs and health care services. Discussion: The aim of this study is to assess the efficacy of a primary care-based cognitive-behavioural psychological intervention in patients with mixed anxiety-depressive disorder. The international scientific evidence has demonstrated the need for psychologists in primary care. However, given the differences between health policies and health services, it is important to test the effect of these psychological interventions in our geographical setting.