The magmatic to hydrothermal evolution of the intrusive Mont Saint-Hilaire Complex: Insights into the late-stage evolution of peralkaline rocks

The Cretaceous Mont Saint-Hilaire complex (Quebec, Canada) comprises three major rock units that were emplaced in the following sequence: (I) gabbros; (II) diorites; (III) diverse partly agpaitic foid syenites. The major element compositions of the rock-forming minerals, age-corrected Nd and oxygen isotope data for mineral separates and trace element data of Fe-Mg silicates from the various lithologies imply a common source for all units. The distribution of the rare earth elements in clinopyroxene from the gabbros indicates an ocean island basalt type composition for the parental magma. Gabbros record temperatures of 1200 to 800 degrees C, variable silica activities between 0 center dot 7 and 0 center dot 3, and f(O2) values between -0 center dot 5 and +0 center dot 7 (log delta FMQ, where FMQ is fayalite-magnetite-quartz). The diorites crystallized under uniform a(SiO2) (a(SiO2) = 0 center dot 4-0 center dot 5) and more reduced f(O2) conditions (log delta FMQ similar to-1) between similar to 1100 and similar to 800 degrees C. Phase equilibria in various foid syenites indicate that silica activities decrease from 0 center dot 6-0 center dot 3 at similar to 1000 degrees C to < 0 center dot 3 at similar to 550 degrees C. Release of an aqueous fluid during the transition to the hydrothermal stage caused a(SiO2) to drop to very low values, which results from reduced SiO(2) solubilities in aqueous fluids compared with silicate melts. During the hydrothermal stage, high water activities stabilized zeolite-group minerals. Fluid inclusions record a complex post-magmatic history, which includes trapping of an aqueous fluid that unmixed from the restitic foid syenitic magma. Cogenetic aqueous and carbonic fluid inclusions reflect heterogeneous trapping of coexisting immiscible external fluids in the latest evolutionary stage. The O and C isotope characteristics of fluid-inclusion hosted CO(2) and late-stage carbonates imply that the surrounding limestones were the source of the external fluids. The mineral-rich syenitic rocks at Mont Saint-Hilaire evolved as follows: first, alkalis, high field strength and large ion lithophile elements were pre-enriched in the (late) magmatic and subsequent hydrothermal stages; second, percolation of external fluids in equilibrium with the carbonate host-rocks and mixing processes with internal fluids as well as fluid-rock interaction governed dissolution of pre-existing minerals, element transport and precipitation of mineral assemblages determined by locally variable parameters. It is this hydrothermal interplay between internal and external fluids that is responsible for the mineral wealth found at Mont Saint-Hilaire.

The crew-scheduling module in the GIST system

The public transportation is gaining importance every year basically duethe population growth, environmental policies and, route and streetcongestion. Too able an efficient management of all the resources relatedto public transportation, several techniques from different areas are beingapplied and several projects in Transportation Planning Systems, indifferent countries, are being developed. In this work, we present theGIST Planning Transportation Systems, a Portuguese project involving twouniversities and six public transportation companies. We describe indetail one of the most relevant modules of this project, the crew-scheduling module. The crew-scheduling module is based on the application of meta-heuristics, in particular GRASP, tabu search and geneticalgorithm to solve the bus-driver-scheduling problem. The metaheuristicshave been successfully incorporated in the GIST Planning TransportationSystems and are actually used by several companies.

On certain greedoid polyhedra, partially indexable scheduling problems and extended restless bandit allocation indices

We present a polyhedral framework for establishing general structural properties on optimal solutions of stochastic scheduling problems, where multiple job classes vie for service resources: the existence of an optimal priority policy in a given family, characterized by a greedoid(whose feasible class subsets may receive higher priority), where optimal priorities are determined by class-ranking indices, under restricted linear performance objectives (partial indexability). This framework extends that of Bertsimas and Niño-Mora (1996), which explained the optimality of priority-index policies under all linear objectives (general indexability). We show that, if performance measures satisfy partial conservation laws (with respect to the greedoid), which extend previous generalized conservation laws, then theproblem admits a strong LP relaxation over a so-called extended greedoid polytope, which has strong structural and algorithmic properties. We present an adaptive-greedy algorithm (which extends Klimov's) taking as input the linear objective coefficients, which (1) determines whether the optimal LP solution is achievable by a policy in the given family; and (2) if so, computes a set of class-ranking indices that characterize optimal priority policies in the family. In the special case of project scheduling, we show that, under additional conditions, the optimal indices can be computed separately for each project (index decomposition). We further apply the framework to the important restless bandit model (two-action Markov decision chains), obtaining new index policies, that extend Whittle's (1988), and simple sufficient conditions for their validity. These results highlight the power of polyhedral methods (the so-called achievable region approach) in dynamic and stochastic optimization.

Driver scheduling problem modelling

The Drivers Scheduling Problem (DSP) consists of selecting a set of duties for vehicle drivers, for example buses, trains, plane or boat drivers or pilots, for the transportation of passengers or goods. This is a complex problem because it involves several constraints related to labour and company rules and can also present different evaluation criteria and objectives. Being able to develop an adequate model for this problem that can represent the real problem as close as possible is an important research area.The main objective of this research work is to present new mathematical models to the DSP problem that represent all the complexity of the drivers scheduling problem, and also demonstrate that the solutions of these models can be easily implemented in real situations. This issue has been recognized by several authors and as important problem in Public Transportation. The most well-known and general formulation for the DSP is a Set Partition/Set Covering Model (SPP/SCP). However, to a large extend these models simplify some of the specific business aspects and issues of real problems. This makes it difficult to use these models as automatic planning systems because the schedules obtained must be modified manually to be implemented in real situations. Based on extensive passenger transportation experience in bus companies in Portugal, we propose new alternative models to formulate the DSP problem. These models are also based on Set Partitioning/Covering Models; however, they take into account the bus operator issues and the perspective opinions and environment of the user.We follow the steps of the Operations Research Methodology which consist of: Identify the Problem; Understand the System; Formulate a Mathematical Model; Verify the Model; Select the Best Alternative; Present the Results of theAnalysis and Implement and Evaluate. All the processes are done with close participation and involvement of the final users from different transportation companies. The planner s opinion and main criticisms are used to improve the proposed model in a continuous enrichment process. The final objective is to have a model that can be incorporated into an information system to be used as an automatic tool to produce driver schedules. Therefore, the criteria for evaluating the models is the capacity to generate real and useful schedules that can be implemented without many manual adjustments or modifications. We have considered the following as measures of the quality of the model: simplicity, solution quality and applicability. We tested the alternative models with a set of real data obtained from several different transportation companies and analyzed the optimal schedules obtained with respect to the applicability of the solution to the real situation. To do this, the schedules were analyzed by the planners to determine their quality and applicability. The main result of this work is the proposition of new mathematical models for the DSP that better represent the realities of the passenger transportation operators and lead to better schedules that can be implemented directly in real situations.

Composition of electricity generation portfolios, pivotal dynamics and market prices

We use a simulation model to study how the diversification of electricity generation portfoliosinfluences wholesale prices. We find that technological diversification generally leads to lower market prices but that the relationship is mediated by the supply to demand ratio. In each demand case there is a threshold where pivotal dynamics change. Pivotal dynamics pre- and post-threshold are the cause of non-linearities in the influence of diversification on market prices. The findings are robust to our choice of behavioural parameters and match close-form solutions where those are available.

Generation of the primary hair follicle pattern.

Hair follicles are spaced apart from one another at regular intervals through the skin. Although follicles are predominantly epidermal structures, classical tissue recombination experiments indicated that the underlying dermis defines their location during development. Although many molecules involved in hair follicle formation have been identified, the molecular interactions that determine the emergent property of pattern formation have remained elusive. We have used embryonic skin cultures to dissect signaling responses and patterning outcomes as the skin spatially organizes itself. We find that ectodysplasin receptor (Edar)-bone morphogenetic protein (BMP) signaling and transcriptional interactions are central to generation of the primary hair follicle pattern, with restriction of responsiveness, rather than localization of an inducing ligand, being the key driver in this process. The crux of this patterning mechanism is rapid Edar-positive feedback in the epidermis coupled with induction of dermal BMP4/7. The BMPs in turn repress epidermal Edar and hence follicle fate. Edar activation also induces connective tissue growth factor, an inhibitor of BMP signaling, allowing BMP action only at a distance from their site of synthesis. Consistent with this model, transgenic hyperactivation of Edar signaling leads to widespread overproduction of hair follicles. This Edar-BMP activation-inhibition mechanism appears to operate alongside a labile prepattern, suggesting that Edar-mediated stabilization of beta-catenin active foci is a key event in determining definitive follicle locations.

SwissParam: a fast force field generation tool for small organic molecules.

The drug discovery process has been deeply transformed recently by the use of computational ligand-based or structure-based methods, helping the lead compounds identification and optimization, and finally the delivery of new drug candidates more quickly and at lower cost. Structure-based computational methods for drug discovery mainly involve ligand-protein docking and rapid binding free energy estimation, both of which require force field parameterization for many drug candidates. Here, we present a fast force field generation tool, called SwissParam, able to generate, for arbitrary small organic molecule, topologies, and parameters based on the Merck molecular force field, but in a functional form that is compatible with the CHARMM force field. Output files can be used with CHARMM or GROMACS. The topologies and parameters generated by SwissParam are used by the docking software EADock2 and EADock DSS to describe the small molecules to be docked, whereas the protein is described by the CHARMM force field, and allow them to reach success rates ranging from 56 to 78%. We have also developed a rapid binding free energy estimation approach, using SwissParam for ligands and CHARMM22/27 for proteins, which requires only a short minimization to reproduce the experimental binding free energy of 214 ligand-protein complexes involving 62 different proteins, with a standard error of 2.0 kcal mol(-1), and a correlation coefficient of 0.74. Together, these results demonstrate the relevance of using SwissParam topologies and parameters to describe small organic molecules in computer-aided drug design applications, together with a CHARMM22/27 description of the target protein. SwissParam is available free of charge for academic users at www.swissparam.ch.

A simple optimised search heuristic for the job-shop scheduling problem

This paper presents a simple Optimised Search Heuristic for the Job Shop Scheduling problem that combines a GRASP heuristic with a branch-and-bound algorithm. The proposed method is compared with similar approaches and leads to better results in terms of solution quality and computing times.

Metaheuristics for the bus-driver scheduling problem

We present new metaheuristics for solving real crew scheduling problemsin a public transportation bus company. Since the crews of thesecompanies are drivers, we will designate the problem by the bus-driverscheduling problem. Crew scheduling problems are well known and severalmathematical programming based techniques have been proposed to solvethem, in particular using the set-covering formulation. However, inpractice, there exists the need for improvement in terms of computationalefficiency and capacity of solving large-scale instances. Moreover, thereal bus-driver scheduling problems that we consider can present variantaspects of the set covering, as for example a different objectivefunction, implying that alternative solutions methods have to bedeveloped. We propose metaheuristics based on the following approaches:GRASP (greedy randomized adaptive search procedure), tabu search andgenetic algorithms. These metaheuristics also present some innovationfeatures based on and genetic algorithms. These metaheuristics alsopresent some innovation features based on the structure of the crewscheduling problem, that guide the search efficiently and able them tofind good solutions. Some of these new features can also be applied inthe development of heuristics to other combinatorial optimizationproblems. A summary of computational results with real-data problems ispresented.

Solving a concrete sleepers production scheduling by genetic algorithms

PRECON S.A is a manufacturing company dedicated to produce prefabricatedconcrete parts to several industries as rail transportation andagricultural industries.Recently, PRECON signed a contract with RENFE,the Spanish Nnational Rail Transportation Company to manufacturepre-stressed concrete sleepers for siding of the new railways of the highspeed train AVE. The scheduling problem associated with the manufacturingprocess of the sleepers is very complex since it involves severalconstraints and objectives. The constraints are related with productioncapacity, the quantity of available moulds, satisfying demand and otheroperational constraints. The two main objectives are related withmaximizing the usage of the manufacturing resources and minimizing themoulds movements. We developed a deterministic crowding genetic algorithmfor this multiobjective problem. The algorithm has proved to be a powerfuland flexible tool to solve the large-scale instance of this complex realscheduling problem.

Carbon and oxygen isotope study of hydrothermal carbonates in the zinc-lead deposits of the San Vicente district, central Peru: A quantitative modeling on mixing processes and CO2 degassing

Carbon and oxygen isotope studies of the host and gangue carbonates of Mississippi Valley-type zinc-lead deposits in the San Vicente District hosted in the Upper Triassic to Lower Jurassic dolostones of the Pucara basin (central Peru) were used to constrain models of the ore formation. A mixing model between an incoming hot saline slightly acidic radiogenic (Pb, Sr) fluid and the native formation water explains the overall isotopic variation (delta(13)C = - 11.5 to + 2.5 parts per thousand relative to PDB and delta(18)O = + 18.0 to + 24.3 parts per thousand relative to SMOW) of the carbonate generations. The dolomites formed during the main ore stage show a narrower range (delta(13)C = - 0.1 to + 1.7 parts per thousand and delta(18)O = + 18.7 to + 23.4 parts per thousand) which is explained by exchange between the mineralizing fluids and the host carbonates combined with changes in temperature and pressure. This model of fluid-rock interaction explains the pervasive alteration of the host dolomite I and precipitation of sphalerite I. The open-space filling hydrothermal white sparry dolomite and the coexisting sphalerite II formed by prolonged fluid-host dolomite interaction and limited CO2 degassing. Late void-filling dolomite III (or calcite) and the associated sphalerite III formed as the consequence of CO2 degassing and concomitant pH increase of a slightly acidic ore fluid. Widespread brecciation is associated to CO2 outgassing. Consequently, pressure variability plays a major role in the ore precipitation during the late hydrothermal events in San Vicente. The presence of native sulfur associated with extremely carbon-light calcites replacing evaporitic sulfates (e.g., delta(13)C = - 11.5 parts per thousand), altered native organic matter and heavier hydrothermal bitumen (from - 27.0 to - 23.0 parts per thousand delta(13)C) points to thermochemical reduction of sulfate and/or thiosulfate. The delta(13)C- and delta(18)O-values of the altered host dolostone and hydrothermal carbonates, and the carbon isotope composition of the associated organic matter show a strong regional homogeneity. These results coupled with the strong mineralogical and petrographic similarities of the different MVT occurrences perhaps reflects the fact that the mineralizing processes were similar in the whole San Vicente belt, suggesting the existence of a common regional mineralizing hydrothermal system with interconnected plumbing.

A glutathione precursor, N-acetyl-cysteine, modulates mismatch negativity generation in schizophrenia patients

Schizophrenia patients exhibit deficits in low-level processing, including pitch discrimination. This deficiency manifests in auditory evoked potentials (AEPs) as an impaired mismatch negativity (MMN), an electrophysiological response to infrequent target stimuli interspersed among frequent standard stimuli that typically peaks ~100ms post-stimulus onset. NMDA receptor antagonists have been shown to block MMN generation in both animals and humans, and NMDA dysfunction has been linked to the underlying pathophysiology of schizophrenia. A parallel line of evidence indicates that glutathione (GSH) regulation is perturbed in schizophrenia patients at the gene, protein and functional levels (Tosic et al., 2006). This GSH dysregulation leads to NMDA receptors' hypofunction through interaction with their redox site (Steullet et al., 2006). The present study aimed to modulate GSH levels in schizophrenia patients and assessed the effects of such a modulation on MMN generation mechanisms. N-acetyl-cysteine (NAC), a GSH precursor, was administered to schizophrenia patients, using a double-blind cross-over protocol. One group received NAC (2g/day) for 60 days and then placebo for another 60 days, and vice-versa for the second group. AEPs from patients were recorded at the onset of the protocol, at the point of cross-over, and at the end of the study. Participants were instructed to manually respond to target stimuli (2kHz pure tones occurring 20% of the time among 1kHz pure tones). Analyses of AEPs recorded at protocol onset indicated that patients (n=11) were significantly impaired in generating the MMN relative to age-matched controls (n=11). Specifically, the global field power (GFP), an index of AEP magnitude, was measured over the 70- 155ms post-stimulus interval and submitted to an analysis of variance (ANOVA). There was a significant interaction between population and stimulus frequency, indicating impaired MMN generation in patients at protocol onset. Analyses of AEPs recorded during administration of NAC (n=7) versus placebo (n=7) revealed the efficacy of this GSH precursor in modulating MMN generation mechanisms. ANOVA of GFP over the 70- 155ms post-stimulus interval, using stimulus frequency and treatment as within-participants variables, revealed a significant interaction and indicated that NAC can ameliorate MMN generation. We discuss these results in terms of potential therapeutic strategies for schizophrenia.