Abstract: In recent decades new recycling technologies for mixed plastic waste have emerged. In pyrolysis, the polymer chains are thermally broken (pyrolyzed) to obtain hydrocarbon materials of different molecular weights such as naphtha, oil or waxes, whose yields can be controlled by varying the reaction parameters. Naphtha represents a closed-loop recycling process as it is a feedstock for (poly)olefins; while the co-production of waxes, having several applications in e.g. the construction industry, exemplifies an open-loop recycling process. This paper compares the economic performance of the pyrolysis of mixed polyolefin waste in a closed-loop and open-loop scheme, including a probabilistic approach to the most important variables. From an economic perspective, open-loop pyrolysis as presented outperforms closed-loop recycling, due to the high prices of wax. However, the results present a high dispersion caused by the volatility of the prices of crude oil and its derivates. Considering the current oil price projections, our case study analysis showed that for open-loop recycling there is a future probability of almost a 98 % of observing positive results and around 57 % of probability in the case of closed-loop recycling, under the assumptions made. Yet, in a future scenario where decarbonized electricity would decrease oil prices, the probability of a positive outcome reduces to 57 % for the open-loop case and to less than 8 % in the case of closed-loop recycling. To make these pathways attractive to investors, the nameplate capacity should be at least 70 kt/year for open-loop recycling and 115 kt/year for closed-loop recycling. A 120 kt/year plant should operate minimally at 80 % of its capacity for open-loop recycling, while closed-loop recycling would demand running close to maximum capacity. Security of feedstock supply therefore is required.

Abstract: The valuable properties of carbon steel slag are currently underexploited. To date, research mainly focusses on valorising a single property of the slag. In this study an ex-ante life cycle assessment (LCA) was applied to evaluate the environmental profile of a novel technological pathway aimed at the extraction of chromium from carbon steel slag in combination with high quality valorisation of the residual matrix material. A comparison with current practice was made, not only by calculating the environmental impact of the lab scale observations, but more importantly by estimating the impact on an industrial scale. Practical guidance on ex-ante LCA is limited, so this study contributes by incorporating simulations on thermodynamic behaviour, complemented with empirical calculation rules and including information derived from similar technologies to perform the upscaling. These principles of ex-ante LCA were applied to the lab results of two consecutive research iterations. Substantial improvements of the environmental profile were observed: ex-ante results turned out to be a factor 20 lower compared to the results from the lab observations after the first iteration and had decreased by a factor 2 compared to the small pilot scale of the second iteration. All upscaled results are better than those from the worst case reference scenario (landfill). Based on the experience gained after this iterative research cycle, a practical recommendation is that at a low technology readiness level using more simple calculation rules in combination with a flowsheet based on elementary design principles for processes at an industrial scale is a more efficient way of modelling compared to a fully-fledged process design from the start.

Abstract: Uniform acceptance force biased Monte Carlo (UFMC) simulations have previously been shown to be a powerful tool to simulate atomic scale processes, enabling one to follow the dynamical path during the simulation. In this contribution, we present a simple proof to demonstrate that this uniform acceptance still complies with the condition of detailed balance, on the condition that the characteristic parameter lambda = 1/2 and that the maximum allowed step size is chosen to be sufficiently small. Furthermore, the relation to Metropolis Monte Carlo (MMC) is also established, and it is shown that UFMC reduces to MMC by choosing the characteristic parameter lambda = 0 [Rao, M. et al. Mol. Phys. 1979, 37, 1773]. Finally, a simple example compares the UFMC and MMC methods.

Abstract: We introduce a simple enhanced sampling approach for the calculation of free energy differences and barriers along a one-dimensional reaction coordinate. First, a small number of short nonequilibrium simulations are carried out along the reaction coordinate, and the Jarzynski equality is used to learn an approximate free energy surface from the nonequilibrium work distribution. This free energy estimate is represented in a compact form as an artificial neural network and used as an external bias potential to accelerate rare events in a subsequent molecular dynamics simulation. The final free energy estimate is then obtained by reweighting the equilibrium probability distribution of the reaction coordinate sampled under the influence of the external bias. We apply our reweighted Jarzynski sampling recipe to four processes of varying scales and complexities.spanning chemical reaction in the gas phase, pair association in solution, and droplet nucleation in supersaturated vapor. In all cases, we find reweighted Jarzynski sampling to be a very efficient strategy, resulting in rapid convergence of the free energy to high precision.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are “hot spots” for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

Abstract: MFI zeolite nanosheets with thickness of 2 and 8 nm were synthesized, transformed into bifunctional catalysts by loading with platinum and tested in n-decane isomerization and hydrocracking. Detailed analysis of skeletal isomers and hydrocracked products revealed that the MFI nanosheets display transition-state shape-selectivity similar to bulk MFI zeolite crystals. The suppressed formation of bulky skeletal isomers and C5 cracking products are observed both in the nanosheets and the bulk crystals grown in three dimensions. This is typical for restricted transition-state shape-selectivity, characteristic for the MFI type pores. It is a first clear example of transition-state shape-selectivity inside a zeolitic nanosheet. Owing to the short diffusion path across the sheets, expression of diffusion-based discrimination of reaction products in the MFI nanosheets was limited. The 2-methylnonane formation among monobranched C10 isomers and 2,7-dimethyloctane among dibranched C10 isomers, which in MFI zeolite are favored by product diffusion, was much less favored on the nanosheets compared to the reference bulk ZSM-5 material.

Abstract: Different zeolites were tested for selective removal of methyl elaidate (trans isomer) from an equimolar mixture with methyl oleate (cis isomer). Sorption experiments of the geometric isomers show that only ZSM-5 samples with reduced Al content in the framework are able to discriminate among the bent cis and the linear trans fatty acid methyl esters. Hydrogenation experiments of equimolar methyl oleate and elaidate mixtures at low temperature (65 °C) and high hydrogen pressure (6.0 MPa), using Pt catalysts, confirm this result. Only with a Pt/NaZSM-5 catalyst outspoken selectivity for the hydrogenation of the trans isomer is obtained. In order to prepare a selective Pt/ZSM-5 catalyst, the influence of Pt addition (impregnation, ion-exchange and competitive ion-exchange) and Pt activation (different calcination and reduction temperatures) on the Pt-distribution and Pt particle size was investigated using SEM, bright-field and HR TEM, EDX, electron tomography, CO-chemisorption, XPS, XRD, and UVvis measurements. The best result in terms of hydrogenation activity and selectivity is obtained with a Pt/ZSM-5 catalyst, which is prepared via competitive ion-exchange, followed by slow calcination up to 350 °C under high O2 flow and a reduction up to 500 °C under H2. This preparation method leads to a Pt/ZSM-5 catalyst with the best Pt distribution and the smallest Pt clusters occluded in the zeolite structure. Finally, the influence of zeolite crystal size, morphology, and elemental composition of ZSM-5 on hydrogenation activity and selectivity was investigated in detail.