[N]pT ensemble and finite-size-scaling study of the critical isostructural transition in the generalized exponential model of index 4.

First-order transitions of system where both lattice site occupancy and lattice spacing fluctuate, such as cluster crystals, cannot be efficiently studied by traditional simulation methods, which necessarily fix one of these two degrees of freedom. The difficulty, however, can be surmounted by the generalized [N]pT ensemble [J. Chem. Phys. 136, 214106 (2012)]. Here we show that histogram reweighting and the [N]pT ensemble can be used to study an isostructural transition between cluster crystals of different occupancy in the generalized exponential model of index 4 (GEM-4). Extending this scheme to finite-size scaling studies also allows us to accurately determine the critical point parameters and to verify that it belongs to the Ising universality class.

Mixtures of g-priors in Generalized Linear Models

Mixtures of Zellner's g-priors have been studied extensively in linear models and have been shown to have numerous desirable properties for Bayesian variable selection and model averaging. Several extensions of g-priors to Generalized Linear Models (GLMs) have been proposed in the literature; however, the choice of prior distribution of g and resulting properties for inference have received considerably less attention. In this paper, we extend mixtures of g-priors to GLMs by assigning the truncated Compound Confluent Hypergeometric (tCCH) distribution to 1/(1+g) and illustrate how this prior distribution encompasses several special cases of mixtures of g-priors in the literature, such as the Hyper-g, truncated Gamma, Beta-prime, and the Robust prior. Under an integrated Laplace approximation to the likelihood, the posterior distribution of 1/(1+g) is in turn a tCCH distribution, and approximate marginal likelihoods are thus available analytically. We discuss the local geometric properties of the g-prior in GLMs and show that specific choices of the hyper-parameters satisfy the various desiderata for model selection proposed by Bayarri et al, such as asymptotic model selection consistency, information consistency, intrinsic consistency, and measurement invariance. We also illustrate inference using these priors and contrast them to others in the literature via simulation and real examples.

Transferability of Regional Permafrost Disturbance Susceptibility Modelling Using Generalized Linear and Generalized Additive Models

To effectively assess and mitigate risk of permafrost disturbance, disturbance-p rone areas can be predicted through the application of susceptibility models. In this study we developed regional susceptibility models for permafrost disturbances using a field disturbance inventory to test the transferability of the model to a broader region in the Canadian High Arctic. Resulting maps of susceptibility were then used to explore the effect of terrain variables on the occurrence of disturbances within this region. To account for a large range of landscape charac- teristics, the model was calibrated using two locations: Sabine Peninsula, Melville Island, NU, and Fosheim Pen- insula, Ellesmere Island, NU. Spatial patterns of disturbance were predicted with a generalized linear model (GLM) and generalized additive model (GAM), each calibrated using disturbed and randomized undisturbed lo- cations from both locations and GIS-derived terrain predictor variables including slope, potential incoming solar radiation, wetness index, topographic position index, elevation, and distance to water. Each model was validated for the Sabine and Fosheim Peninsulas using independent data sets while the transferability of the model to an independent site was assessed at Cape Bounty, Melville Island, NU. The regional GLM and GAM validated well for both calibration sites (Sabine and Fosheim) with the area under the receiver operating curves (AUROC) N 0.79. Both models were applied directly to Cape Bounty without calibration and validated equally with AUROC's of 0.76; however, each model predicted disturbed and undisturbed samples differently. Addition- ally, the sensitivity of the transferred model was assessed using data sets with different sample sizes. Results in- dicated that models based on larger sample sizes transferred more consistently and captured the variability within the terrain attributes in the respective study areas. Terrain attributes associated with the initiation of dis- turbances were similar regardless of the location. Disturbances commonly occurred on slopes between 4 and 15°, below Holocene marine limit, and in areas with low potential incoming solar radiation

Parallel nonconvex generalized Benders decomposition for natural gas production network planning under uncertainty

A scenario-based two-stage stochastic programming model for gas production network planning under uncertainty is usually a large-scale nonconvex mixed-integer nonlinear programme (MINLP), which can be efficiently solved to global optimality with nonconvex generalized Benders decomposition (NGBD). This paper is concerned with the parallelization of NGBD to exploit multiple available computing resources. Three parallelization strategies are proposed, namely, naive scenario parallelization, adaptive scenario parallelization, and adaptive scenario and bounding parallelization. Case study of two industrial natural gas production network planning problems shows that, while the NGBD without parallelization is already faster than a state-of-the-art global optimization solver by an order of magnitude, the parallelization can improve the efficiency by several times on computers with multicore processors. The adaptive scenario and bounding parallelization achieves the best overall performance among the three proposed parallelization strategies.

A generalized model for population dynamics where interactions change over time

La dinámica demográfica ha sido modelada con ecuaciones diferenciales desde que Malthus comenzó sus estudios hace más de doscientos años atrás. Los modelos convencionales siempre tratan relaciones entre especies como estáticas, denotando sólo su dependencia durante un período fijo del tiempo, aunque sea conocido que las relaciones entre especies pueden cambiar con el tiempo. Aquí proponemos un modelo para la dinámica demográfica que incorpora la evolución con el tiempo de las interacciones entre especies. Este modelo incluye una amplia gama de interacciones, de depredador-presa a las relaciones mutualistas, ya sea obligada o facultativa. El mecanismo que describimos permite la transición de una clase de relación entre especies a algún otro, según algunos parámetros externos fijados por el contexto. Estas transiciones podrían evitar la extinción de una de las especies, si esto termina por depender demasiado del ambiente o su relación con las otras especies.

A Generalized Methodology to Characterize Composite Materials for Pyrolysis Models

The predictive capabilities of computational fire models have improved in recent years such that models have become an integral part of many research efforts. Models improve the understanding of the fire risk of materials and may decrease the number of expensive experiments required to assess the fire hazard of a specific material or designed space. A critical component of a predictive fire model is the pyrolysis sub-model that provides a mathematical representation of the rate of gaseous fuel production from condensed phase fuels given a heat flux incident to the material surface. The modern, comprehensive pyrolysis sub-models that are common today require the definition of many model parameters to accurately represent the physical description of materials that are ubiquitous in the built environment. Coupled with the increase in the number of parameters required to accurately represent the pyrolysis of materials is the increasing prevalence in the built environment of engineered composite materials that have never been measured or modeled. The motivation behind this project is to develop a systematic, generalized methodology to determine the requisite parameters to generate pyrolysis models with predictive capabilities for layered composite materials that are common in industrial and commercial applications. This methodology has been applied to four common composites in this work that exhibit a range of material structures and component materials. The methodology utilizes a multi-scale experimental approach in which each test is designed to isolate and determine a specific subset of the parameters required to define a material in the model. Data collected in simultaneous thermogravimetry and differential scanning calorimetry experiments were analyzed to determine the reaction kinetics, thermodynamic properties, and energetics of decomposition for each component of the composite. Data collected in microscale combustion calorimetry experiments were analyzed to determine the heats of complete combustion of the volatiles produced in each reaction. Inverse analyses were conducted on sample temperature data collected in bench-scale tests to determine the thermal transport parameters of each component through degradation. Simulations of quasi-one-dimensional bench-scale gasification tests generated from the resultant models using the ThermaKin modeling environment were compared to experimental data to independently validate the models.

The generating function of the generalized Fibonacci sequence

Using tools of the theory of orthogonal polynomials we obtain the generating function of the generalized Fibonacci sequence established by Petronilho for a sequence of real or complex numbers {Qn} defined by Q0 = 0, Q1 = 1, Qm = ajQm−1 + bjQm−2, m ≡ j (mod k), where k ≥ 3 is a fixed integer, and a0, a1, . . . , ak−1, b0, b1, . . . , bk−1 are 2k given real or complex numbers, with bj #0 for 0 ≤ j ≤ k−1. For this sequence some convergence proprieties are obtained.

Solvabilty of generalized Hammerstein integral equations on unbounded domains, with sign-changing kernels

In this work we study an Hammerstein generalized integral equation u(t)=∫_{-∞}^{+∞}k(t,s) f(s,u(s),u′(s),...,u^{(m)}(s))ds, where k:ℝ²→ℝ is a W^{m,∞}(ℝ²), m∈ℕ, kernel function and f:ℝ^{m+2}→ℝ is a L¹-Carathéodory function. To the best of our knowledge, this paper is the first one to consider discontinuous nonlinearities with derivatives dependence, without monotone or asymptotic assumptions, on the whole real line. Our method is applied to a fourth order nonlinear boundary value problem, which models moderately large deflections of infinite nonlinear beams resting on elastic foundations under localized external loads.

Beam Energy Dependence Of Moments Of The Net-charge Multiplicity Distributions In Au+au Collisions At Rhic.

We report the first measurements of the moments--mean (M), variance (σ(2)), skewness (S), and kurtosis (κ)--of the net-charge multiplicity distributions at midrapidity in Au+Au collisions at seven energies, ranging from sqrt[sNN]=7.7 to 200 GeV, as a part of the Beam Energy Scan program at RHIC. The moments are related to the thermodynamic susceptibilities of net charge, and are sensitive to the location of the QCD critical point. We compare the products of the moments, σ(2)/M, Sσ, and κσ(2), with the expectations from Poisson and negative binomial distributions (NBDs). The Sσ values deviate from the Poisson baseline and are close to the NBD baseline, while the κσ(2) values tend to lie between the two. Within the present uncertainties, our data do not show nonmonotonic behavior as a function of collision energy. These measurements provide a valuable tool to extract the freeze-out parameters in heavy-ion collisions by comparing with theoretical models.

Partitioning The Net Effect Of Host Diversity On An Emerging Amphibian Pathogen.

The 'dilution effect' (DE) hypothesis predicts that diverse host communities will show reduced disease. The underlying causes of pathogen dilution are complex, because they involve non-additive (driven by host interactions and differential habitat use) and additive (controlled by host species composition) mechanisms. Here, we used measures of complementarity and selection traditionally employed in the field of biodiversity-ecosystem function (BEF) to quantify the net effect of host diversity on disease dynamics of the amphibian-killing fungus Batrachochytrium dendrobatidis (Bd). Complementarity occurs when average infection load in diverse host assemblages departs from that of each component species in uniform populations. Selection measures the disproportionate impact of a particular species in diverse assemblages compared with its performance in uniform populations, and therefore has strong additive and non-additive properties. We experimentally infected tropical amphibian species of varying life histories, in single- and multi-host treatments, and measured individual Bd infection loads. Host diversity reduced Bd infection in amphibians through a mechanism analogous to complementarity (sensu BEF), potentially by reducing shared habitat use and transmission among hosts. Additionally, the selection component indicated that one particular terrestrial species showed reduced infection loads in diverse assemblages at the expense of neighbouring aquatic hosts becoming heavily infected. By partitioning components of diversity, our findings underscore the importance of additive and non-additive mechanisms underlying the DE.

Dependence of the crossover exponent with the diffusion rate in the generalized contact process model

We study how the crossover exponent, phi, between the directed percolation (DP) and compact directed percolation (CDP) behaves as a function of the diffusion rate in a model that generalizes the contact process. Our conclusions are based in results pointed by perturbative series expansions and numerical simulations, and are consistent with a value phi = 2 for finite diffusion rates and phi = 1 in the limit of infinite diffusion rate.

Fracionamento dos carboidratos pelas equações do Cornell Net Carbohydrate and Protein System de três cultivares de girassol na presença ou não de irrigação

Objetivou-se quantificar as frações de carboidratos pelas equações do Cornell Net Carbohydrate and Protein System (CNCPS) de três cultivares de girassol (Helianthus annuus L.) cultivados na presença ou não de irrigação. A utilização de uma preparação fibrosa, denominada parede celular (PC), nas equações da CNCPS, em substituição à fibra em detergente neutro (FDN) não promoveu diferenças nas frações de carboidratos B1 e C, mas influenciou as frações A e B2. Como os valores da fração B1, obtidos pelo modelo CNCPS foram menores que os teores de amido e pectina determinados em laboratório, supõe-se que a pectina e outros oligossacarídeos da parede celular, solubilizados pela solução de detergente neutro (fibra solúvel), nunca fizeram parte da fração B1, e sim da fração A. Apesar de os carboidratos da fibra solúvel apresentarem elevadas taxas de degradação, não parece adequada a caracterização da fibra solúvel na fração A. Parece mais adequado que a fibra solúvel (que inclui a pectina) seja alocada a uma fração exclusivamente sua, que pode ser a fração B2, e que seja criada uma nova fração, a B3, para os carboidratos digeríveis da parede celular. Assim, a fração B1 seria composta apenas de amido. A equação da fração C, que estima os carboidratos indigeríveis da parede celular, pode ser simplificada, relacionando a fração indigerível ao teor de lignina na matéria seca, e não à FDN isenta de cinzas e proteína, como atualmente utilizado. Esta proposta tem implicações práticas, uma vez que a fração indigerível da parede celular tem sido expressa em relação à FDN, e não na MS, com base no fato de que os efeitos inibitórios da lignina ocorrem sobre os componentes fibrosos da parede celular vegetal, e não sobre o conteúdo celular.