Population Balance Equations for Reactive Separation in Polymer Upcycling

Many polymer upcycling efforts aim to convert plastic waste into high-value liquid hydrocarbons. However, the subsequent cleavage of middle distillates to light gases can be problematic. The reactor often contains a vapor phase (light gases and middle distillates) and a liquid phase (molten polymers and waxes with a suspended or dissolved catalyst). Because the catalyst resides in the liquid phase, middle distillates that partition into the vapor phase are protected against further cleavage into light gases. In this paper, we consider a simple reactive separation strategy, in which a gas outflow removes the volatile products as they form. We combine vapor-liquid equilibrium models and population balance equations (PBEs) to describe polymer upcycling in a two-phase semibatch reactor. The results suggest that the temperature, headspace volume, and flow rate of the reactor can be used to tune selectivity toward the middle distillates, in addition to the molecular mechanism of catalysis. We anticipate that two-phase reactor models will be important in many polymer upcycling processes and that reactive separation strategies will provide ways to boost the yield of the desired products in these cases.

Engulfment Avalanches and Thermal Hysteresis for Antifreeze Proteins on Supercooled Ice

Antifreeze proteins (AFPs) bind to the ice-water surface and prevent ice growth at temperatures below 0 °C through a Gibbs-Thomson effect. Each adsorbed AFP creates a metastable depression on the surface that locally resists ice growth, until ice engulfs the AFP. We recently predicted the susceptibility to engulfment as a function of AFP size, distance between AFPs, and supercooling [ J. Chem. Phys. 2023, 158, 094501]. For an ensemble of AFPs adsorbed on the ice surface, the most isolated AFPs are the most susceptible, and when an isolated AFP gets engulfed, its former neighbors become more isolated and more susceptible to engulfment. Thus, an initial engulfment event can trigger an avalanche of subsequent engulfment events, leading to a sudden surge of unrestrained ice growth. This work develops a model to predict the supercooling at which the first engulfment event will occur for an ensemble of randomly distributed AFP pinning sites on an ice surface. Specifically, we formulate an inhomogeneous survival probability that accounts for the AFP coverage, the distribution of AFP neighbor distances, the resulting ensemble of engulfment rates, the ice surface area, and the cooling rate. We use the model to predict thermal hysteresis trends and compare with experimental data.

Site-Averaged Ab Initio Kinetics: Importance Learning for Multistep Reactions on Amorphous Supports

Single-atom centers on amorphous supports include catalysts for polymerization, partial oxidation, metathesis, hydrogenolysis, and more. The disordered environment makes each site different, and the kinetics exponentially magnifies these differences to make ab initio site-averaged kinetics calculations extremely difficult. This work extends the importance learning algorithm for efficient and precise site-averaged kinetics estimates to ab initio calculations and multistep reaction mechanisms. Specifically, we calculate site-averaged proton transfer relaxation rates on an ensemble of cluster models representing Brønsted acid sites on silica-alumina. We include direct and water-assisted proton transfer pathways and simultaneously estimate the water adsorption and activation enthalpies for forward and backward proton transfers. We use density functional theory (DFT) to obtain a site-averaged rate, somewhat like a turnover frequency, for the proton transfer relaxation rate. Finally, we show that importance learning can provide orders-of-magnitude acceleration over standard sampling methods for site-averaged rate calculations in cases where the rate is dominated by a few highly active sites.

Broken bond models, magic-sized clusters, and nucleation theory in nanoparticle synthesis

Magic clusters are metastable faceted nanoparticles that are thought to be important and, sometimes, observable intermediates in the nucleation of certain faceted crystallites. This work develops a broken bond model for spheres with a face-centered-cubic packing that form tetrahedral magic clusters. With just one bond strength parameter, statistical thermodynamics yield a chemical potential driving force, an interfacial free energy, and free energy vs magic cluster size. These properties exactly correspond to those from a previous model by Mule et al. [J. Am. Chem. Soc. 143, 2037 (2021)]. Interestingly, a Tolman length emerges (for both models) when the interfacial area, density, and volume are treated consistently. To describe the kinetic barriers between magic cluster sizes, Mule et al. invoked an energy parameter to penalize the two-dimensional nucleation and growth of new layers in each facet of the tetrahedra. According to the broken bond model, barriers between magic clusters are insignificant without the additional edge energy penalty. We estimate the overall nucleation rate without predicting the rates of formation for intermediate magic clusters by using the Becker-Döring equations. Our results provide a blueprint for constructing free energy models and rate theories for nucleation via magic clusters starting from only atomic-scale interactions and geometric considerations.

Free energy barriers for anti-freeze protein engulfment in ice: Effects of supercooling, footprint size, and spatial separation

Anti-freeze proteins (AFPs) protect organisms at freezing conditions by attaching to the ice surface and arresting its growth. Each adsorbed AFP locally pins the ice surface, resulting in a metastable dimple for which the interfacial forces counteract the driving force for growth. As supercooling increases, these metastable dimples become deeper, until metastability is lost in an engulfment event where the ice irreversibly swallows the AFP. Engulfment resembles nucleation in some respects, and this paper develops a model for the "critical profile" and free energy barrier for the engulfment process. Specifically, we variationally optimize the ice-water interface and estimate the free energy barrier as a function of the supercooling, the AFP footprint size, and the distance to neighboring AFPs on the ice surface. Finally, we use symbolic regression to derive a simple closed-form expression for the free energy barrier as a function of two physically interpretable, dimensionless parameters.

Zirconium-Catalyzed C-H Alumination of Polyolefins, Paraffins, and Methane

C-H/Et-Al exchange in zirconium-catalyzed reactions of saturated hydrocarbons and AlEt₃ affords versatile organoaluminum compounds and ethane. The grafting of commercially available Zr(OtBu)₄ on silica/alumina gives monopodal ≡SiO-Zr(OtBu)₃ surface pre-catalyst sites that are activated in situ by ligand exchange with AlEt₃. The catalytic C-H alumination of dodecane at 150 °C followed by quenching in air affords n-dodecanol as the major product, revealing selectivity for methyl group activation. Shorter hydrocarbon or alcohol products were not detected under these conditions. Catalytic reactions of cyclooctane and AlEt₃, however, afford ring-opened products, indicating that C-C bond cleavage occurs readily in methyl group-free reactants. This selectivity for methyl group alumination enables the C-H alumination of polyethylenes, polypropylene, polystyrene, and poly-α-olefin oils without significant chain deconstruction. In addition, the smallest hydrocarbon, methane, undergoes selective mono-alumination under solvent-free catalytic conditions, providing a direct route to Al-Me species.

Oriented attachment kinetics for rod-like particles at a flat surface: Buffon's needle at the nanoscale

The adsorption of large rod-like molecules or crystallites on a flat crystal face, similar to Buffon's needle, requires the rods to "land," with their binding sites in precise orientational alignment with matching sites on the surface. An example is provided by long, helical antifreeze proteins (AFPs), which bind at specific facets and orientations on the ice surface. The alignment constraint for adsorption, in combination with the loss in orientational freedom as the molecule diffuses toward the surface, results in an entropic barrier that hinders the adsorption. Prior kinetic models do not factor in the complete geometry of the molecule, nor explicitly enforce orientational constraints for adsorption. Here, we develop a diffusion-controlled adsorption theory for AFP molecules binding at specific orientations to flat ice surfaces. We formulate the diffusion equation with relevant boundary conditions and present analytical solutions to the attachment rate constant. The resulting rate constant is a function of the length and aspect ratio of the AFP, the distance threshold associated with binding, and solvent conditions such as temperature and viscosity. These results and methods of calculation may also be useful for predicting the kinetics of crystal growth through oriented attachment.

Kinetic coefficient for ice-water interface from simulated non-equilibrium relaxation at coexistence

In the theory of solidification, the kinetic coefficient multiplies the local supercooling to give the solid-liquid interface velocity. The same coefficient should drive interface migration at the coexistence temperature in proportion to a curvature force. This work computes the ice-water kinetic coefficient from molecular simulations starting from a sinusoidal ice-water interface at the coexistence temperature. We apply this method to the basal and prismatic ice planes and compare results to previous estimates from equilibriumcorrelation functions and simulations at controlled supercooling.

Persistent nucleation and size dependent attachment kinetics produce monodisperse PbS nanocrystals

Modern syntheses of colloidal nanocrystals yield extraordinarily narrow size distributions that are believed to result from a rapid "burst of nucleation" (La Mer, JACS, 1950, 72(11), 4847-4854) followed by diffusion limited growth and size distribution focusing (Reiss, J. Chem. Phys., 1951, 19, 482). Using a combination of in situ X-ray scattering, optical absorption, and 13C nuclear magnetic resonance (NMR) spectroscopy, we monitor the kinetics of PbS solute generation, nucleation, and crystal growth from three thiourea precursors whose conversion reactivity spans a 2-fold range. In all three cases, nucleation is found to be slow and continues during >50% of the precipitation. A population balance model based on a size dependent growth law (1/r) fits the data with a single growth rate constant (k G) across all three precursors. However, the magnitude of the k G and the lack of solvent viscosity dependence indicates that the rate limiting step is not diffusion from solution to the nanoparticle surface. Several surface reaction limited mechanisms and a ligand penetration model that fits data our experiments using a single fit parameter are proposed to explain the results.

Size-Controlled Nanoparticles Embedded in a Mesoporous Architecture Leading to Efficient and Selective Hydrogenolysis of Polyolefins

A catalytic architecture, comprising a mesoporous silica shell surrounding platinum nanoparticles (NPs) supported on a solid silica sphere (mSiO₂/Pt-X/SiO₂; X is the mean NP diameter), catalyzes hydrogenolysis of melt-phase polyethylene (PE) into a narrow C₂₃-centered distribution of hydrocarbons in high yield using very low Pt loadings (∼10^-5 g Pt/g PE). During catalysis, a polymer chain enters a pore and contacts a Pt NP where the C-C bond cleavage occurs and then the smaller fragment exits the pore. mSiO₂/Pt/SiO₂ resists sintering or leaching of Pt and provides high yields of liquids; however, many structural and chemical effects on catalysis are not yet resolved. Here, we report the effects of Pt NP size on activity and selectivity in PE hydrogenolysis. Time-dependent conversion and yields and a lumped kinetics model based on the competitive adsorption of long vs short chains reveal that the activity of catalytic material is highest with the smallest NPs, consistent with a structure-sensitive reaction. Remarkably, the three mSiO₂/Pt-X/SiO₂ catalysts give equivalent selectivity. We propose that mesoscale pores in the catalytic architecture template the C₂₃-centered distribution, whereas the active Pt sites influence the carbon-carbon bond cleavage rate. This conclusion provides a framework for catalyst design by separating the C-C bond cleavage activity at catalytic sites from selectivity for chain lengths of the products influenced by the structure of the catalytic architecture. The increased activity, selectivity, efficiency, and lifetime obtained using this architecture highlight the benefits of localized and confined environments for isolated catalytic particles under condensed-phase reaction conditions.

Diffusion Attachment Model for Long Helical Antifreeze Proteins to Ice

Some of the most potent antifreeze proteins (AFPs) are approximately rigid helical structures that bind with one side in contact with the ice surface at specific orientations. These AFPs take random orientations in solution; however, most orientations become sterically inaccessible as the AFP approaches the ice surface. The effect of these inaccessible orientations on the rate of adsorption of AFP to ice has never been explored. Here, we present a diffusion-controlled theory of adsorption kinetics that accounts for these orientational restrictions to predict a rate constant for adsorption (k_on, in m/s) as a function of the length and width of the AFP molecules. We find that k_on decreases with length and diameter of the AFP and is almost proportional to the inverse of the area of the binding surface. We demonstrate that the restricted orientations create an entropic barrier to AFP adsorption, which we compute to be approximately 7 k_BT for most AFPs and up to 9 k_BT for Maxi, the largest known AFP. We compare the entropic resistance 1/k_on to resistances for diffusion through boundary layers and across typical distances in the extracellular matrix and find that these entropic and diffusion resistances could become comparable in the small confined spaces of biological environments.

Free energies of crystals computed using Einstein crystal with fixed center of mass and differing spring constants

Free energies of crystals computed using a center of mass constraint require a finite-size correction, as shown in previous work by Polson et al. [J. Chem. Phys. 112, 5339-5342 (2000)]. Their reference system is an Einstein crystal with equal spring constants. In this paper, we extend the work of Polson et al. [J. Chem. Phys. 112, 5339-5342 (2000)] to the case of differing spring constants. The generalization is convenient for constraining the center of mass in crystals with atoms of differing masses, and it helps to optimize the free energy calculations. To test the theory, we compare the free energies of LiI and NaCl crystals from calculations with differing spring constants to those computed using equal spring constants. Using these center of mass finite size corrections, we compute the true free energies of these crystals for different system sizes to eliminate the intrinsic finite-size effects. These calculations help demonstrate the size of these finite-size corrections relative to other contributions to the absolute free energy of the crystals.

Analysis of a large cohort of African/Afro-Caribbean patients by cardiac magnetic resonance

Funding Acknowledgements

Type of funding sources: None.

Background

Limited data is available in the literature on Cardiac Magnetic Resonance (CMR) features of African/Afro-Caribbeans, in particular exploring differences between hypertensive heart disease (HHD) and hypertrophic cardiomyopathy (HCM). 

Purpose

To describe the cardiac morphology and pattern of late gadolinium enhancement (LGE) in patients of African/Afro-Caribbean origin. 

Methods

We retrospectively analysed African/Afro-Caribbean patients who underwent clinical CMR at a tertiary centre. Three groups were specifically investigated: HHD, HCM and combination of HHD and HCM or ambiguous (HHD&HCM). 

Results

Overall, 166 consecutive patients (58% male, mean age 55 ± 14yo) were analysed. One-hundred fifty-four (93%) had history of arterial hypertension (HTN), including 17 with uncontrolled/malignant HTN. Overall, 28(17%) had normal scans, 70(42%) HHD, 15(9%) apical HCM, 10(6%) classical septal HCM, 19(11%) dilated cardiomyopathy, 6(4%) cardiac amyloidosis, 3(2%) ischemic heart disease, 4(2%) myocarditis, 2 sarcoidosis and 1 valvular disease. In 7(4%) the diagnosis was ambiguous between HHD and HCM and 1 uncertain. Forty-four (27%) had dual pathology, most frequently HHD, bystander myocardial infarction (MI) and embolic MI. LGE was detected in 95(57%), 26 with ischemic pattern, 13 diffuse mid-wall, 29 focal non-ischemic and 28 diffuse/multifocal non-ischemic. CMR features and correlations between subgroups are reported in the Table. Patients with HHD had significantly higher left ventricular (LV) end-diastolic (ED) volume indexed (Vi) and LV end-systolic (ES) Vi, but lower LV ejection fraction (EF) and LV maximum wall thickness (MWT) compared to HCM patients. HHD&HCM had higher LVEF and MWT compared to HHD and higher LVEDVi compared to HCM. LGE was more frequently seen in HCM and HHD&HCM as focal non-ischemic (6vs5vs10,p = 0.049) and diffuse multifocal(5vs6vs9,p = 0.009). A history of uncontrolled/malignant HTN was more frequent in HHD and HHD&HCM (11vs1vs5,p = 0.025) and was associated with diffuse LGE with lateral wall involvement (p < 0.0001) (Figure). Conclusions: CMR findings in African/Afro-Caribbeans may overlap between aetiologies. A specific pattern of diffuse non-ischaemic LGE involving the lateral wall appears though to be more often associated with severe uncontrolled HTN.

Table HHD(n = 70) HCM(n = 19) HHD&HCM or ambiguous (n = 13) p LVEDVi, ml/m2 88 ± 33 62 ± 11 77 ± 6 0.003 LVESVi, ml/m2 36 ± 26 16 ± 10 26 ± 19 0.002 LVEF, % 62 ± 16 79 ± 6 74 ± 12 0.567 LVMWT, mm 13 ± 2 15 ± 4 14 ± 3 <0.0001 CMR features and correlations between subgroups Abstract Figure

Short- and Long-term Outcome After Interventional VSD Closure: A Single-Center Experience in Pediatric and Adult Patients

Interventional closure of congenital ventricular septal defects (VSD) is recording a continuous rise in acceptance. Complete atrioventricular block (cAVB) and residual shunting are major concerns during follow-up, but long-term data for both are still limited. We retrospectively evaluated the outcome of patients with interventional VSD closure and focused on long-term results (> 1 year follow-up). Transcatheter VSD closures were performed between 1993 and 2015, in 149 patients requiring 155 procedures (104 perimembranous, 29 muscular, 19 residual post-surgical VSDs, and 3 with multiple defects). The following devices were used: 65 × Amplatzer™ Membranous VSD Occluder, 33 × Duct Occluder II, 27 × Muscular VSD Occluder, 3 × Duct Occluder I, 24 × PFM-Nit-Occlud®, and 3 × Rashkind-Occluder. The median age at time of implantation was 6.2 (0.01-66.1) years, median height 117 (49-188) cm, and median weight 20.9 (3.2-117) kg. Median follow-up time was 6.2 (1.1-21.3) years and closure rate was 86.2% at last follow-up. Complications resulting in device explantation include one case of cAVB with a Membranous VSD occluder 7 days after implantation and four cases due to residual shunt/malposition. Six (4%) deaths occurred during follow-up with only one procedural related death from a hybrid VSD closure. Overall, our reported results of interventional VSD closure show favorable outcomes with only one (0.7%) episode of cAVB. Interventional closure offers a good alternative to surgical closure and shows improved performance by using softer devices. However, prospective long-term data in the current era with different devices are still mandatory to assess the effectiveness and safety of this procedure.

Diabat method for polymorph free energies: Extension to molecular crystals

Lattice-switch Monte Carlo and the related diabat methods have emerged as efficient and accurate ways to compute free energy differences between polymorphs. In this work, we introduce a one-to-one mapping from the reference positions and displacements in one molecular crystal to the positions and displacements in another. Two features of the mapping facilitate lattice-switch Monte Carlo and related diabat methods for computing polymorph free energy differences. First, the mapping is unitary so that its Jacobian does not complicate the free energy calculations. Second, the mapping is easily implemented for molecular crystals of arbitrary complexity. We demonstrate the mapping by computing free energy differences between polymorphs of benzene and carbamazepine. Free energy calculations for thermodynamic cycles, each involving three independently computed polymorph free energy differences, all return to the starting free energy with a high degree of precision. The calculations thus provide a force field independent validation of the method and allow us to estimate the precision of the individual free energy differences.

Importance learning estimator for the site-averaged turnover frequency of a disordered solid catalyst

For disordered catalysts such as atomically dispersed "single-atom" metals on amorphous silica, the active sites inherit different properties from their quenched-disordered local environments. The observed kinetics are site-averages, typically dominated by a small fraction of highly active sites. Standard sampling methods require expensive ab initio calculations at an intractable number of sites to converge on the site-averaged kinetics. We present a new method that efficiently estimates the site-averaged turnover frequency (TOF). The new estimator uses the same importance learning algorithm [Vandervelden et al., React. Chem. Eng. 5, 77 (2020)] that we previously used to compute the site-averaged activation energy. We demonstrate the method by computing the site-averaged TOF for a simple disordered lattice model of an amorphous catalyst. The results show that with the importance learning algorithm, the site-averaged TOF and activation energy can now be obtained concurrently with orders of magnitude reduction in required ab initio calculations.

Absolute chemical potentials for complex molecules in fluid phases: A centroid reference for predicting phase equilibria

Solid-fluid phase equilibria are difficult to predict in simulations because bound degrees of freedom in the crystal phase must be converted to free translations and rotations in the fluid phase. Here, we avoid the solid-to-fluid transformation step by starting with chemical potentials for two reference systems, one for the fluid phase and one for the solid phase. For the solid, we start from the Einstein crystal and transform to the fully interacting molecular crystal. For the fluid phase, we introduce a new reference system, the "centroid," and then transform to gas phase molecules. We illustrate the new calculations by predicting the sublimation vapor pressure of succinic acid in the temperature range of 300 K-350 K.

Is mapping the BCG vaccine-induced immune responses the key to improving the efficacy against tuberculosis?

In recent years, the century-old Mycobacterium bovis Bacillus Calmette-Guérin (BCG) vaccine against tuberculosis (TB) has been re-evaluated for its capacity to stem the global tide of TB. There is increasing evidence that the efficacy of BCG can be improved by the modified administration methods and schedules. Here, we first discuss recent approaches of vaccine administration, revaccination or boosting that have been used to try to improve the efficacy of BCG against TB. We then dive deeper into studies investigating the immune correlates of protection and describe studies that have investigated BCG-specific T-cell responses and the influence of environmental exposures. These studies all highlight that there is still a lot to learn about the immune response induced by BCG, both in terms of phenotype and specificity, which has been surprisingly understudied. We argue that several critical gaps in knowledge exist and must be addressed by future research to rationally improve the efficacy of BCG, including comprehensive, proteome-wide understanding of the epitopes derived from BCG recognized by BCG-vaccinated individuals, the phenotype of responding antigen-specific T cells and how previous exposure to environmental mycobacteria affect these parameters and thus influence vaccine efficacy. The development of modern techniques allows us to answer some of these questions to better understand how BCG works in terms of both protection against TB and the immune response that it triggers.

Polyethylene upcycling to long-chain alkylaromatics by tandem hydrogenolysis/aromatization

The current scale of plastics production and the accompanying waste disposal problems represent a largely untapped opportunity for chemical upcycling. Tandem catalytic conversion by platinum supported on γ-alumina converts various polyethylene grades in high yields (up to 80 weight percent) to low-molecular-weight liquid/wax products, in the absence of added solvent or molecular hydrogen, with little production of light gases. The major components are valuable long-chain alkylaromatics and alkylnaphthenes (average ~C30, dispersity Ð = 1.1). Coupling exothermic hydrogenolysis with endothermic aromatization renders the overall transformation thermodynamically accessible despite the moderate reaction temperature of 280°C. This approach demonstrates how waste polyolefins can be a viable feedstock for the generation of molecular hydrocarbon products.

[A boy with white sponge naevus of the buccal mucosa]

White sponge naevus (WSN) is a rare, autosomal dominant disorder that causes various complaints WSN is most commonly found on the buccal mucosa. Clinically, the white, slightly elevated lesions of WSN may be confused with other disorders on oral mucosa. We report a case of WSN in a 14-year-old boy who had complaints for a considerable period of time. WSN is caused by mutations in KRT4 and KRT13.

Secondary effectiveness factors for catalytic reactions in series: extension to slab, cylindrical, and spherical geometries

‘Secondary’ effectiveness factors are developed for reactions in series that are useful for controlling the yield of the intermediate product (B).

Grafting metal complexes onto amorphous supports: from elementary steps to catalyst site populationsviakernel regression

We use transition state theory, kernel regression, and population balance modeling techniques to model the grafting of metal complexes onto amorphous catalyst supports.

Site-averaged kinetics for catalysts on amorphous supports: an importance learning algorithm

We combine importance sampling and kernel regression techniques to efficiently predict site-averaged kinetics for isolated catalyst sites on amorphous supports.

Female patients are less satisfied with biological treatment for psoriasis and experience more side-effects than male patients: results from the prospective BioCAPTURE registry

BACKGROUND

Female sex has been reported as a predictor for treatment discontinuation with biological therapies for psoriasis, although reasons remain unclear. It can be hypothesized that lower satisfaction with biological treatment in women might add to the lower drug survival rates.

OBJECTIVES

To identify possible differences in satisfaction with biological treatment between female and male patients using the Treatment Satisfaction Questionnaire for Medication (TSQM).

METHODS

Data of psoriasis patients treated with biologics were obtained from the prospective, multicentre, daily-practice BioCAPTURE registry. Longitudinal TSQM data were analysed by linear mixed models. Relevant patient characteristics were incorporated as possible confounding factors. Post hoc analysis of adverse events was performed in order to investigate differences between sexes.

RESULTS

We included 315 patients with 396 corresponding treatment episodes (137 adalimumab, 90 etanercept, 137 ustekinumab, 24 secukinumab and 8 infliximab). Almost forty per cent of the patients were female. Women had significantly lower baseline PASI scores (P = 0.01). Longitudinal analyses demonstrated lower TSQM scores for 'side-effects' (P = 0.05) and 'global satisfaction' (P = 0.01) in female patients compared with male patients over 1 year of treatment. Women reported more relevant adverse events in the context of biologic treatment compared to men (rate ratio 1.79; P < 0.001), with more fungal (rate ratio 2.20; P = 0.001) and herpes simplex infections (rate ratio 3.25; P = 0.005).

CONCLUSIONS

This study provides a prospective, longitudinal analysis of treatment satisfaction with biologics in female and male patients with psoriasis. Women were slightly less satisfied with treatment regarding side-effects and global satisfaction. Differences in treatment satisfaction and side-effects might add to the fact that women discontinue biological treatments more often.

Gibbs free-energy differences between polymorphs via a diabat approach

Polymorph free-energy differences are critical to several applications. A recently proposed diabat interpolation framework estimated free-energy differences between polymorphs by quadratic interpolation of diabats. This work extends the Zwanzig-Bennett relation to the NPT ensemble so that the diabats directly give Gibbs free-energy differences. We also demonstrate how the approach can be used in cases where the diabats are not parabolic. We illustrate the diabat method for Gibbs free-energy difference of zirconium (BCC and HCP phases) and compare it with the conventional lattice switch Monte Carlo approach.

The Clathrate-Water Interface Is Oleophilic

The slow nucleation of clathrate hydrates is a central challenge for their use in the storage and transportation of natural gas. Molecules that strongly adsorb to the clathrate-water interface decrease the crystal-water surface tension, lowering the barrier for clathrate nucleation. Surfactants are widely used to promote the nucleation and growth of clathrate hydrates. It has been proposed that these amphiphilic molecules bind to the clathrate surface via hydrogen bonding. However, recent studies reveal that PVCap, an amphiphilic polymer, binds to clathrates through hydrophobic moieties. Here we use molecular dynamic simulations and theory to investigate the mode and strength of binding of surfactants to the clathrate-water interface and their effect on the nucleation rate. We find that the surfactants bind to the clathrate-water interface exclusively through their hydrophobic tails. The binding is strong, driven by the entropy of dehydration of the alkyl chain, as it penetrates empty cavities at the hydrate surface. The hydrophobic attraction of alkyl groups to the clathrate surface also results in strong adsorption of alkanes. We identify two regimes for the binding of surfactants as a function of their density at the hydrate surface, which we interpret to correspond to the two steps of the Langmuir adsorption isotherm observed in experiments. Our results indicate that hydrophobic attraction to the clathrate-water interface is key for the design of soluble additives that promote the nucleation of hydrates. We use the calculated adsorption coefficients to estimate the concentration of sodium dodecyl sulfate (SDS) required to reach nucleation rates for methane hydrate consistent with those measured in experiments. To our knowledge, this study is the first to quantify the effect of surfactant concentration in the nucleation rate of clathrate hydrates.

HIV: ageing, cognition and neuroimaging at 4-year follow-up

OBJECTIVES

The aim of the study was to investigate the hypothesis of accelerated cognitive ageing in HIV-positive individuals using longitudinal assessment of cognitive performance and quantitative magnetic resonance imaging (MRI).

METHODS

We assessed a broad cognitive battery and quantitative MRI metrics [voxel-based morphometry (VBM) and diffusion tensor imaging (DTI)] in asymptomatic HIV-positive men who have sex with men (15 aged 20-40 years and 15 aged ≥ 50 years), and HIV-seronegative matched controls (nine aged 20-40 years and 16 aged ≥ 50 years).

RESULTS

Being HIV positive was associated with greater decreases in executive function and global cognition. Additionally, using DTI, we found that the HIV-positive group had a greater increase in mean diffusivity, but we did not find group differences in volume change using VBM. With respect to the HIV status by age group interaction, this was statistically significant for change in global cognition, with older HIV-positive individuals showing greater global cognitive decline, but there were no significant interaction effects on other measures. Lastly, change in cognitive performance was correlated with change in the DTI measures, and this effect was stronger for the HIV-positive participants.

CONCLUSIONS

In the present study, we found some evidence for accelerated ageing in HIV-positive individuals, with a statistically significant HIV status by age group interaction in global cognition, although this interaction could not be explained by the imaging findings. Moreover, we also found that change in cognitive performance was correlated with change in the DTI measures, and this effect was stronger for the HIV-positive participants. This will need replication in larger studies using a similarly lengthy follow-up period.