Bacterial spores respond to humidity similarly to hydrogels.

Bacterial spores have outstanding properties from the materials science perspective, which allow them to survive extreme environmental conditions. Recent work by [S. G. Harrellson et al., Nature 619, 500-505 (2023)] studied the mechanical properties of Bacillus subtilis spores and the evolution of these properties with the change of humidity. The experimental measurements were interpreted assuming that the spores behave as water-filled porous solids, subjected to hydration forces. Here, we revisit their experimental data using literature data on vapor sorption on spores and ideas from polymer physics. We demonstrate that upon the change of humidity, the spores behave like rubber with respect to their swelling, elasticity, and relaxation times. This picture is consistent with the knowledge of the materials comprising the bacterial cell walls-cross-linked peptidoglycan. Our results provide an interpretation of the mechanics of bacterial spores and can help in developing synthetic materials mimicking the mechanical properties of the spores.

Adsorption-Induced Deformation of Zeolites 4A and 13X: Experimental and Molecular Simulation Study.

Gas adsorption in zeolites leads to adsorption-induced deformation, which can significantly affect the adsorption and diffusive properties of the system. In this study, we conducted both experimental investigations and molecular simulations to understand the deformation of zeolites 13X and 4A during carbon dioxide adsorption at 273 K. To measure the sample's adsorption isotherm and strain simultaneously, we used a commercial sorption instrument with a custom-made sample holder equipped with a dilatometer. Our experimental data showed that while the zeolites 13X and 4A exhibited similar adsorption isotherms, their strain isotherms differed significantly. To gain more insight into the adsorption process and adsorption-induced deformation of these zeolites, we employed coupled Monte Carlo and molecular dynamics simulations with atomistically detailed models of the frameworks. Our modeling results were consistent with the experimental data and helped us identify the reasons behind the different deformation behaviors of the considered structures. Our study also revealed the sensitivity of the strain isotherm of zeolites to pore size and other structural and energetic features, suggesting that measuring adsorption-induced deformation could serve as a complementary method for material characterization and provide guidelines for related technical applications.

Surface Tension of Organophosphorus Compounds: Sarin and its Surrogates.

While the production and stockpiling of organophosphorus chemical warfare agents (CWAs), such as sarin, was banned three decades ago, CWAs have remained a threat. New approaches for decontamination and destruction of CWAs require detailed knowledge of their various physicochemical properties. In particular, surface tension is needed to describe the formation and evolution of hazardous aerosols when CWA liquids are dispersed in the air. Due to the extreme toxicity of sarin, most experimental studies are carried out using its surrogates─organophosphorus compounds which, while having similar structures, are much less toxic, e.g., dimethyl methylphosphonate (DMMP) and diisopropyl methylphosphonate (DIMP). However, not only for sarin, but also for its surrogates, literature data on the surface tension are scarce. Here we present experimental measurements and computational predictions of the surface tension of DMMP and DIMP. Classical molecular dynamics simulations using the Transferable Potentials for Phase Equilibria (TraPPE) force field produced an excellent agreement with the experimental results in the temperature range from 3 to 60 °C, validating the predictive capability of TraPPE. Consequently, we applied the TraPPE force field to sarin. Our modeled values for the sarin surface tension cover the range of temperatures from 0 to 85 °C, and the four experimental data points from the literature measured between 20 and 35 °C agree perfectly with our predictions. The temperature-dependent surface tension values for sarin and its surrogates obtained in our study can be used in models predicting the formation and evolution of aerosols made of these chemicals. Furthermore, our results justify the use of the TraPPE force field to derive the thermodynamic properties of other organophosphorus compounds with structures similar to the ones studied here.

The effect of interconnections on gas adsorption in materials with spherical mesopores: A Monte Carlo simulation study.

Gas adsorption is a standard method for measuring pore-size distributions of nanoporous materials. This method is often based on assuming the pores as separate entities of a certain simple shape: slit-like, cylindrical, or spherical. Here, we study the effect of interconnections on gas adsorption in materials with spherical pores, such as three-dimensionally ordered mesoporous (3DOm) carbons. We consider interconnected systems with two, four, and six windows of various sizes. We propose a simple method based on the integration of solid-fluid interactions to take into account these windows. We used Monte Carlo simulations to model argon adsorption at the normal boiling point and obtained adsorption isotherms for the range of systems. For a system with two windows, we obtained a remarkably smooth transition from the spherical to cylindrical isotherm. Depending on the size and number of windows, our system resembles both spherical and cylindrical pores. These windows can drastically shift the point of capillary condensation and result in pore-size distributions that are very different from the ones based on a spherical pore model. Our results can be further used for modeling fluids in a system of interconnected pores using Monte Carlo and density functional theory methods.

Pore-Size Distribution of Silica Colloidal Crystals from Nitrogen Adsorption Isotherms.

Silica colloidal crystals are face-centered cubic structures comprised of silica spheres with the diameters ranging between tens and hundreds of nanometers. The voids between the spheres form pores, which can be probed by nitrogen adsorption porosimetry. Here, we prepared two mesoporous samples and a macroporous reference sample and then measured nitrogen adsorption and desorption isotherms for further characterization. We proposed a straightforward procedure for calculation of the pore-size distribution of silica colloidal crystals from nitrogen adsorption isotherms. The procedure is based on the adsorption integral equation solution with a kernel of theoretical isotherms, consistent with the procedure used for many other porous materials. The solution is carried out using the non-negative least squares (NNLS) regression with Tikhonov regularization. The kernel of mesoporous isotherms is built on the basis of the macroscopic Derjaguin-Broekhoff-de Boer (DBdB) theory of capillary condensation considering the voids as a network of spheres. Application of our procedure for the analysis of the adsorption branches of experimental isotherms resulted in bimodal distributions, where the modes matched well with the sizes of the voids in the colloidal crystals face centered cubic structure: the main mode corresponds to the octahedral voids and the second mode to the tetrahedral voids. Furthermore, we modified the surface of the samples with organics and repeated the characterization procedure for the modified samples. The resulting pore-size distribution for the samples with the modified surface matched the original one quite closely. It demonstrates the procedure as a simple and efficient technique to estimate the pore-size distribution and justifies the spherical shape approximation for the voids in the silica colloidal crystals.

Porous Structure of Silica Colloidal Crystals.

We prepared silica colloidal crystals with different pore sizes using isothermal heating evaporation-induced self-assembly in quantities suitable for nitrogen porosimetry and studied their porous structure. We observed pores of two types in agreement with the description of silica colloidal crystals as face-centered cubic packed structures containing octahedral and tetrahedral voids. We calculated the sizes of these pores using the Derjaguin-Broekhoff-de Boer theory of capillary condensation for spherical pores. We also described the pore geometry mathematically and showed that the octahedral pore radii measured experimentally matches closely the radii of the spheres of the same volume. In the case of the tetrahedral pores, the proposed approach underestimated the pore radius by ca. 40%. Overall, this simple geometrical description provides a good representation of the porous system in silica colloidal crystals.

Mechanical Characterization of Hierarchical Structured Porous Silica by in Situ Dilatometry Measurements during Gas Adsorption.

Mechanical properties of hierarchically structured nanoporous materials are determined by the solid phase stiffness and the pore network morphology. We analyze the mechanical stiffness of hierarchically structured silica monoliths synthesized via a sol-gel process, which possess a macroporous scaffold built of interconnected struts with hexagonally ordered cylindrical mesopores. We consider samples with and without microporosity within the mesopore walls and analyze them on the macroscopic level as well as on the microscopic level of the mesopores. Untreated as-prepared samples still containing some organic components and the respective calcined and sintered counterparts of varying microporosity are investigated. To determine Young's moduli on the level of the macroscopic monoliths, we apply ultrasonic run time measurements, while Young's moduli of the mesopore walls are obtained by analysis of the in situ strain isotherms during N2 adsorption at 77 K. For the latter, we extended our previously reported theoretical approach for this type of materials by incorporating the micropore effects, which are clearly not negligible in the calcined and most of the sintered samples. The comparison of the macro- and microscopic Young's moduli reveals that both properties follow essentially the same trends, that is, calcination and sintering increase the mechanical stiffness on both levels. Consequently, stiffening of the monolithic samples can be primarily attributed to stiffening of the backbone material which is consistent with the fact that the morphology on the mesopore level is mainly preserved with the post-treatments applied.

In Situ Small-Angle Neutron Scattering Investigation of Adsorption-Induced Deformation in Silica with Hierarchical Porosity.

Adsorption-induced deformation of a series of silica samples with hierarchical porosity has been studied by in situ small-angle neutron scattering (SANS) and in situ dilatometry. Monolithic samples consisted of a disordered macroporous network of struts formed by a 2D lattice of hexagonally ordered cylindrical mesopores and disordered micropores within the mesopore walls. Strain isotherms were obtained at the mesopore level by analyzing the shift of the Bragg reflections from the ordered mesopore lattice in SANS data. Thus, SANS essentially measured the radial strain of the cylindrical mesopores including the volume changes of the mesopore walls due to micropore deformation. A H2O/D2O adsorbate with net zero coherent neutron scattering length density was employed in order to avoid apparent strain effects due to intensity changes during pore filling. In contrast to SANS, the strain isotherms obtained from in situ dilatometry result from a combination of axial and radial mesopore deformation together with micropore deformation. Strain data were quantitatively analyzed with a theoretical model for micro-/mesopore deformation by combining information from nitrogen and water adsorption isotherms to estimate the water-silica interaction. It was shown that in situ SANS provides complementary information to dilatometry and allows for a quantitative estimate of the elastic properties of the mesopore walls from water adsorption.

Molecular Simulations Shed Light on Potential Uses of Ultrasound in Nitrogen Adsorption Experiments.

Nitrogen adsorption is one of the main characterization techniques for nanoporous materials. The experimental adsorption isotherm provides information about the surface area and pore size distribution (PSD) for a sample. In this work we show that additional insight into PSD can be gained when the speed of sound propagation through a sample is measured during nitrogen adsorption experiment. We analyzed published experimental data on ultrasound propagation through a nanoporous Vycor glass sample during nitrogen adsorption experiment. Next, we calculated the change of the longitudinal and shear moduli of the sample as a function of relative vapor pressure. From this, we show that the shear modulus of the sample does not change upon filling the pores, evidencing that adsorbed nitrogen at 77 K has zero shear modulus, similarly to a bulk liquid. The longitudinal modulus of the sample behaves differently: it changes abruptly at the capillary condensation and keeps gradually increasing thereafter. We performed Monte Carlo molecular simulations to predict the compressibility of adsorbed nitrogen and then calculated the longitudinal modulus of the nitrogen-saturated Vycor using the Gassmann equation. Our theoretical predictions nicely matched the longitudinal modulus derived from the experimental data. Additionally, we performed molecular simulations to model nitrogen adsorbed in silica pores of sizes ranging from 2 to 8 nm. We found that the isothermal elastic modulus of adsorbed nitrogen depends linearly on the inverse pore size. This dependence, along with the proposed recipe for probing the modulus of adsorbed nitrogen, sets up the grounds for extracting additional information about the porous samples, when the nitrogen adsorption is combined with ultrasonic experiments.

Single Parameter for Predicting the Morphology of Atmospheric Black Carbon.

Black carbon (BC) from fuel combustion is an effective light absorber that contributes significantly to direct climate forcing. The forcing is altered when BC combines with other substances, which modify its mixing state and morphology, making the evaluation of its atmospheric lifetime and climate impact a challenge. To elucidate the associated mechanisms, we exposed BC aerosol to supersaturated vapors of different chemicals to form thin coatings and measured the coating mass required to induce the restructuring of BC aggregates. We found that studied chemicals fall into two distinct groups based on a single dimensionless parameter, χ, which depends on the diameter of BC monomer spheres and the coating material properties, including vapor supersaturation, molar volume, and surface tension. We show that when χ is small (low-volatility chemicals), the highly supersaturated vapor condenses uniformly over aggregates, including convex monomers and concave junctions in between monomers, but when χ is large (intermediate-volatility chemicals), junctions become preferred. The aggregates undergo prompt restructuring when condensation in the junctions dominates over condensation on monomer spheres. For a given monomer diameter, the coating distribution is mostly controlled by vapor supersaturation. The χ factor can be incorporated straightforwardly into atmospheric models to improve simulations of BC aging.

Effect of pore geometry on the compressibility of a confined simple fluid.

Fluids confined in nanopores exhibit properties different from the properties of the same fluids in bulk; among these properties is the isothermal compressibility or elastic modulus. The modulus of a fluid in nanopores can be extracted from ultrasonic experiments or calculated from molecular simulations. Using Monte Carlo simulations in the grand canonical ensemble, we calculated the modulus for liquid argon at its normal boiling point (87.3 K) adsorbed in model silica pores of two different morphologies and various sizes. For spherical pores, for all the pore sizes (diameters) exceeding 2 nm, we obtained a logarithmic dependence of fluid modulus on the vapor pressure. Calculation of the modulus at saturation showed that the modulus of the fluid in spherical pores is a linear function of the reciprocal pore size. The calculation of the modulus of the fluid in cylindrical pores appeared too scattered to make quantitative conclusions. We performed additional simulations at higher temperature (119.6 K), at which Monte Carlo insertions and removals become more efficient. The results of the simulations at higher temperature confirmed both regularities for cylindrical pores and showed quantitative difference between the fluid moduli in pores of different geometries. Both of the observed regularities for the modulus stem from the Tait-Murnaghan equation applied to the confined fluid. Our results, along with the development of the effective medium theories for nanoporous media, set the groundwork for analysis of the experimentally measured elastic properties of fluid-saturated nanoporous materials.

Adsorption-Induced Deformation of Hierarchically Structured Mesoporous Silica-Effect of Pore-Level Anisotropy.

The goal of this work is to understand adsorption-induced deformation of hierarchically structured porous silica exhibiting well-defined cylindrical mesopores. For this purpose, we performed an in situ dilatometry measurement on a calcined and sintered monolithic silica sample during the adsorption of N2 at 77 K. To analyze the experimental data, we extended the adsorption stress model to account for the anisotropy of cylindrical mesopores, i.e., we explicitly derived the adsorption stress tensor components in the axial and radial direction of the pore. For quantitative predictions of stresses and strains, we applied the theoretical framework of Derjaguin, Broekhoff, and de Boer for adsorption in mesopores and two mechanical models of silica rods with axially aligned pore channels: an idealized cylindrical tube model, which can be described analytically, and an ordered hexagonal array of cylindrical mesopores, whose mechanical response to adsorption stress was evaluated by 3D finite element calculations. The adsorption-induced strains predicted by both mechanical models are in good quantitative agreement making the cylindrical tube the preferable model for adsorption-induced strains due to its simple analytical nature. The theoretical results are compared with the in situ dilatometry data on a hierarchically structured silica monolith composed by a network of mesoporous struts of MCM-41 type morphology. Analyzing the experimental adsorption and strain data with the proposed theoretical framework, we find the adsorption-induced deformation of the monolithic sample being reasonably described by a superposition of axial and radial strains calculated on the mesopore level. The structural and mechanical parameters obtained from the model are in good agreement with expectations from independent measurements and literature, respectively.

Adsorption-Induced Surface Stresses of the Water/Quartz Interface: Ab Initio Molecular Dynamics Study.

Adsorption-induced deformation is expansion or contraction of a solid due to adsorption on its surface. This phenomenon is important for a wide range of applications, from chemomechanical sensors to methane recovery from geological formations. The strain of the solid is driven by the change of the surface stress due to adsorption. Using ab initio molecular dynamics, we calculate the surface stresses for the dry α-quartz surfaces, and investigate how these stresses change when the surfaces are exposed to water. We find that the nonhydroxylated surface shows small and approximately isotropic changes in stress, while the hydroxylated surface, which interacts more strongly with the polar water molecules, shows larger and qualitatively anisotropic (opposite sign in xx and yy) surface stress changes. All of these changes are several times larger than the surface tension of water itself. The anisotropy and possibility of positive surface stress change can explain experimentally observed surface area contraction due to adsorption.

Deformation of Microporous Carbons during N2, Ar, and CO2 Adsorption: Insight from the Density Functional Theory.

Using the nonlocal density functional theory, we investigate adsorption of N2 (77 K), Ar (77 K), and CO2 (273 K) and respective adsorption-induced deformation of microporous carbons. We show that the smallest micropores comparable in size and even smaller than the nominal molecular diameter of the adsorbate contribute significantly to the development of the adsorption stress. While pores of approximately the nominal adsorbate diameter exhibit no adsorption stress regardless of their filling level, the smaller pores cause expansive adsorption stresses up to almost 4 GPa. Accounting for this effect, we determined the pore-size distribution of a synthetic microporous carbon by simultaneously fitting its experimental CO2 adsorption isotherm (273 K) and corresponding adsorption-induced strain measured by in situ dilatometry. Based on the pore-size distribution and the elastic modulus fitted from CO2 data, we predicted the sample's strain isotherms during N2 and Ar adsorption (77 K), which were found to be in reasonable agreement with respective experimental data. The comparison of calculations and experimental results suggests that adsorption-induced deformation caused by micropores is not limited to the low relative pressures typically associated with the micropore filling, but is effective over the whole relative pressure range up to saturation pressure.

Revisiting Bangham's law of adsorption-induced deformation: changes of surface energy and surface stress.

When fluids are adsorbed on a solid surface they induce noticeable stresses, which cause the deformation of the solid. D. H. Bangham and co-authors performed a series of experimental measurements of adsorption-induced strains, and concluded that physisorption causes expansion, which is proportional to the lowering of the surface energy Δγ. This statement is referred to as the Bangham effect or Bangham's law. However, it is known that the quantity that controls the deformation is actually the change in surface stress Δf rather than surface energy Δγ, but this difference has not been considered in the context of adsorption-induced deformation of mesoporous materials. We use the Brunauer-Emmett-Teller (BET) theory to derive both values and show the difference between them. We find the condition when the difference between the two vanishes, and Bangham's law is applicable; it is likely that this condition is satisfied in most cases, and prediction of strain based on Δγ is a good approximation. We show that this is the case for adsorption of argon and water on Vycor glass. Finally, we show that the difference between Δγ and Δf can explain some of the experimental data that contradicts Bangham's law.

Modulus-pressure equation for confined fluids.

Ultrasonic experiments allow one to measure the elastic modulus of bulk solid or fluid samples. Recently such experiments have been carried out on fluid-saturated nanoporous glass to probe the modulus of a confined fluid. In our previous work [G. Y. Gor et al., J. Chem. Phys., 143, 194506 (2015)], using Monte Carlo simulations we showed that the elastic modulus K of a fluid confined in a mesopore is a function of the pore size. Here we focus on the modulus-pressure dependence K(P), which is linear for bulk materials, a relation known as the Tait-Murnaghan equation. Using transition-matrix Monte Carlo simulations we calculated the elastic modulus of bulk argon as a function of pressure and argon confined in silica mesopores as a function of Laplace pressure. Our calculations show that while the elastic modulus is strongly affected by confinement and temperature, the slope of the modulus versus pressure is not. Moreover, the calculated slope is in a good agreement with the reference data for bulk argon and experimental data for confined argon derived from ultrasonic experiments. We propose to use the value of the slope of K(P) to estimate the elastic moduli of an unknown porous medium.

Relation between pore size and the compressibility of a confined fluid.

When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here, we report a simple relation between the pore size and isothermal compressibility of argon confined in such pores. Compressibility is calculated from the fluctuations of the number of particles in the grand canonical ensemble using two different simulation techniques: conventional grand-canonical Monte Carlo and grand-canonical ensemble transition-matrix Monte Carlo. Our results provide a theoretical framework for extracting the information on the pore sizes of fluid-saturated samples by measuring the compressibility from ultrasonic experiments.

Adsorption stress changes the elasticity of liquid argon confined in a nanopore.

Knowledge of the elastic properties of a fluid is crucial for predicting its flow under high pressure, particularly in porous media. However, when a fluid is confined to a nanopore, many of its thermodynamic properties change as compared to bulk. Here we study the effect of confinement on the bulk modulus of liquid argon adsorbed in mesopores using classical density functional theory. We show that, at pressures lower than the saturation pressure, high adsorption stress in the pore causes the lowering of the fluid bulk modulus, a phenomenon which was recently observed experimentally [ Schappert, K.; Pelster, R. Europhys. Lett. 2014 , 105 , 5600 ]. Furthermore, we find that the pore size has a strong effect on the fluid bulk modulus, so that even at saturation, the elastic properties of nanoconfined fluid differ from the bulk values. We show that this difference is also due to the adsorption stress. Our results provide a basis for a new method for characterization of porous materials and have implications for modeling fluids in nanoporous geological formations, such as coal or shale.

Non-covalent interactions of nitrous oxide with aromatic compounds: spectroscopic and computational evidence for the formation of 1:1 complexes.

We present the first study of intermolecular interactions between nitrous oxide (N2O) and three representative aromatic compounds (ACs): phenol, cresol, and toluene. The infrared spectroscopic experiments were performed in a Ne matrix and were supported by high-level quantum chemical calculations. Comparisons of the calculated and experimental vibrational spectra provide direct identification and characterization of the 1:1 N2O-AC complexes. Our results show that N2O is capable of forming non-covalently bonded complexes with ACs. Complex formation is dominated by dispersion forces, and the interaction energies are relatively low (about -3 kcal mol(-1)); however, the complexes are clearly detected by frequency shifts of the characteristic bands. These results suggest that N2O can be bound to the amino-acid residues tyrosine or phenylalanine in the form of π complexes.