Some examples of modeling in chemical engineering: thermal treatment of polymers, phase equilibrium calculations and process design

Presentation submitted to PSE Seminar, Chemical Engineering Department, Center for Advanced Process Design-making (CAPD), Carnegie Mellon University, Pittsburgh (USA), October 2012.

Emergence of half-metallicity in suspended NiO chains: Ab initio electronic structure and quantum transport calculations

Contrary to the antiferromagnetic and insulating character of bulk NiO, one-dimensional chains of this material can become half metallic due to the lower coordination of their atoms. Here we present ab initio electronic structure and quantum transport calculations of ideal infinitely long NiO chains and of more realistic short ones suspended between Ni electrodes. While infinite chains are insulating, short suspended chains are half-metallic minority-spin conductors that displays very large magnetoresistance and a spin-valve behavior controlled by a single atom.

VO: a strongly correlated metal close to a Mott-Hubbard transition

Here, we present experimental and computational evidences to support that rocksalt cubic VO is a strongly correlated metal with non-Fermi-liquid thermodynamics and an unusually strong spin-lattice coupling. An unexpected change of sign of metallic thermopower with composition is tentatively ascribed to the presence of a pseudogap in the density of states. These properties are discussed as signatures of the proximity to a magnetic quantum phase transition. The results are summarized in an electronic phase diagram for the 3d monoxides, which resembles that of other strongly correlated systems. The structural and electronic simplicity of 3d monoxides makes them ideal candidates to progress in the understanding of highly correlated electron systems.

Psychometric behaviour of the strengths and difficulties questionnaire (SDQ) in the Spanish national health survey 2006

Background: The Strengths and Difficulties Questionnaire (SDQ) is a tool to measure the risk for mental disorders in children. The aim of this study is to describe the diagnostic efficiency and internal structure of the SDQ in the sample of children studied in the Spanish National Health Survey 2006. Methods: A representative sample of 6,773 children aged 4 to 15 years was studied. The data were obtained using the Minors Questionnaire in the Spanish National Health Survey 2006. The ROC curve was constructed and calculations made of the area under the curve, sensitivity, specificity and the Youden J indices. The factorial structure was studied using models of exploratory factorial analysis (EFA) and confirmatory factorial analysis (CFA). Results: The prevalence of behavioural disorders varied between 0.47% and 1.18% according to the requisites of the diagnostic definition. The area under the ROC curve varied from 0.84 to 0.91 according to the diagnosis. Factor models were cross-validated by means of two different random subsamples for EFA and CFA. An EFA suggested a three correlated factor model. CFA confirmed this model. A five-factor model according to EFA and the theoretical five-factor model described in the bibliography were also confirmed. The reliabilities of the factors of the different models were acceptable (>0.70, except for one factor with reliability 0.62). Conclusions: The diagnostic behaviour of the SDQ in the Spanish population is within the working limits described in other countries. According to the results obtained in this study, the diagnostic efficiency of the questionnaire is adequate to identify probable cases of psychiatric disorders in low prevalence populations. Regarding the factorial structure we found that both the five and the three factor models fit the data with acceptable goodness of fit indexes, the latter including an externalization and internalization dimension and perhaps a meaningful positive social dimension. Accordingly, we recommend studying whether these differences depend on sociocultural factors or are, in fact, due to methodological questions.

Handling of Uncertainty in Life Cycle Inventory by Correlated Multivariate Lognormal Distributions: Application to the Design of Supply Chain Networks

Poster presented in the 24th European Symposium on Computer Aided Process Engineering (ESCAPE 24), Budapest, Hungary, June 15-18, 2014.

Role of potassium orbitals in the metallic behavior of K3picene

Detailed electronic structure calculations of picene clusters doped by potassium modeling the crystalline K3picene structure show that while two electrons are completely transferred from potassium atoms to the lowest-energy unoccupied molecular orbital of pristine picene, the third one remains closely attached to both material components. Multiconfigurational analysis is necessary to show that many structures of almost degenerate total energies compete to define the cluster ground state. Our results prove that the 4s orbital of potassium should be included in any interaction model describing the material. We propose a quarter-filled two-orbital model as the most simple model capable of describing the electronic structure of K-intercalated picene. Precise solutions obtained by a development of the Lanczos method show low-energy electronic excitations involving orbitals located at different positions. Consequently, metallic transport is possible in spite of the clear dominance of interaction over hopping.

Handling of Uncertainty in Life Cycle Inventory by Correlated Multivariate Lognormal Distributions: Application to the Design of Supply Chain Networks

In this work, we analyze the effect of incorporating life cycle inventory (LCI) uncertainty on the multi-objective optimization of chemical supply chains (SC) considering simultaneously their economic and environmental performance. To this end, we present a stochastic multi-scenario mixed-integer linear programming (MILP) coupled with a two-step transformation scenario generation algorithm with the unique feature of providing scenarios where the LCI random variables are correlated and each one of them has the desired lognormal marginal distribution. The environmental performance is quantified following life cycle assessment (LCA) principles, which are represented in the model formulation through standard algebraic equations. The capabilities of our approach are illustrated through a case study of a petrochemical supply chain. We show that the stochastic solution improves the economic performance of the SC in comparison with the deterministic one at any level of the environmental impact, and moreover the correlation among environmental burdens provides more realistic scenarios for the decision making process.

Improvements in DFT Calculations of Spin–Spin Coupling Constants

Different types of spin–spin coupling constants (SSCCs) for several representative small molecules are evaluated and analyzed using a combination of 10 exchange functionals with 12 correlation functionals. For comparison, calculations performed using MCSCF, SOPPA, other common DFT methods, and also experimental data are considered. A detailed study of the percentage of Hartree–Fock exchange energy in SSCCs and in its four contributions is carried out. From the above analysis, a combined functional formed with local Slater (34%), Hartree–Fock exchange (66%), and P86 correlation functional (S66P86) is proposed in this paper. The accuracy of the values obtained with this hybrid functional (mean absolute deviation of 4.5 Hz) is similar to that of the SOPPA method (mean absolute deviation of 4.6 Hz).

Cleavage of the C–C Bond in the Ethanol Oxidation Reaction on Platinum. Insight from Experiments and Calculations

Using a combination of experimental and computational methods, mainly FTIR and DFT calculations, new insights are provided here in order to better understand the cleavage of the C–C bond taking place during the complete oxidation of ethanol on platinum stepped surfaces. First, new experimental results pointing out that platinum stepped surfaces having (111) terraces promote the C–C bond breaking are presented. Second, it is computationally shown that the special adsorption properties of the atoms in the step are able to promote the C–C scission, provided that no other adsorbed species are present on the step, which is in agreement with the experimental results. In comparison with the (111) terrace, the cleavage of the C–C bond on the step has a significantly lower activation energy, which would provide an explanation for the observed experimental results. Finally, reactivity differences under acidic and alkaline conditions are discussed using the new experimental and theoretical evidence.

Thermodynamic calculations for the salt crystallisation damage in porous built heritage using PHREEQC

This work considers the crystallisation mechanisms of the most common and aggressive salts that generate stress in porous building stones as a result of changing ambient conditions. These mechanisms include the salt crystallisation that result from decreasing relative humidity and changes in temperature and, in hydrated salts, the dissolution of the lower hydrated form and the subsequent precipitation of the hydrated salt. We propose a new methodology for thermodynamic calculations using PHREEQC that includes these crystallisation mechanisms. This approach permits the calculation of the equilibrium relative humidity and the parameterization of the critical relative humidity and crystallisation pressures for the dissolution–precipitation transitions. The influence of other salts on the effectives of salt crystallisation and chemical weathering is also assessed. We review the sodium and magnesium sulphate and sodium chloride systems, in both single and multicomponent solutions, and they are compared to the sodium carbonate and calcium carbonate systems. The variation of crystallisation pressure, the formation of new minerals and the chemical dissolution by the presence of other salts is also evaluated. Results for hydrated salt systems show that high crystallisation pressures are possible as lower hydrated salts dissolve and more hydrated salts precipitate. High stresses may be also produced by decreasing temperature, although it requires that porous materials are wet for long periods of time. The presence of other salts changes the temperature and relative humidity of salt transitions that generates stress rather than reducing the pressure of crystallisation, if any salt has previously precipitated. Several practical conclusions derive from proposed methodology and provide conservators and architects with information on the potential weathering activity of soluble salts. Furthermore, the model calculations might be coupled with projections of future climate to give as improved understanding of the likely changes in the frequency of phase transitions in salts within porous stone.