Search NASA⌕ Search

Engineering topics

Colón, Yamil J.

Publications and source records attributed to Colón, Yamil J..

Combining High-Throughput Experiments and Active Learning to Characterize Deep Eutectic Solvents

The high tunability of deep eutectic solvents (DESs) stems from the ease of changing their precursors and relative compositions. However, measuring the physicochemical properties across large composition and temperature ranges, necessary to properly design target-specific DESs, is tedious and error-prone and represents a bottleneck in the advancement and scalability of DES-based applications. As such, active learning (AL) methodologies based on Gaussian processes (GPs) were developed in this work to minimize the experimental effort necessary to characterize DESs. Owing to its importance for large-scale applications, the reduction of DES viscosity through the addition of a low-molecular-weight solvent was explored as a case study. A high-throughput experimental screening was initially performed on nine different ternary DESs. Then, GPs were successfully trained to predict DES viscosity from its composition and temperature, showcasing the ability of these stochastic, nonparametric models to accurately describe the physicochemical properties of complex mixtures. Finally, the ability of GPs to provide estimates of their own uncertainty was leveraged through an AL framework to minimize the number of data points necessary to obtain accurate viscosity modes. This led to a significant reduction in data requirements, with many systems requiring only five independent viscosity data points to be properly described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Expanding the Design Space of Polymer–Metal Organic Framework (MOF) Gels by Understanding Polymer–MOF Interactions

The fabrication of polymer-MOF composite gels holds great potential to provide emergent properties for drug delivery, environmental remediation, and catalysis. To leverage the full potential of these composites, we investigated how the presence and chemistry of polymers impact MOF formation within the composites and, in turn, how MOFs impact polymer gelation. We show that polymers with a high density of strongly metal-binding carboxylic acids inhibit MOF formation; however, reducing the density of carboxylic acids or substituting them with weaker metal-binding hydroxyl groups permits both MOF formation and gelation within composites. Preparing composites with poly(ethylene glycol) (PEG), which does not bind MOF zirconium (Zr)-oxo clusters, and observing gelation suggests that MOFs can entrap polymer chains to create cross-links in addition to cross-linking them through polymer-Zr-oxo interactions. Both simulations and experiments show composite hydrogels formed with poly(vinyl alcohol) (PVA) to be more stable than those made with PEG, which can reptate through MOF pores upon heating. We demonstrate the generalizability of this composite formation process across different Zr-based MOFs (UiO-66, NU-901, UiO-67, and MOF-525) and by spin-coating gels into conformable films. PVA-UiO-66 composite hydrogels demonstrated high sorption and sustained release of methylene blue relative to the polymer alone (3× loading, 28× slower release), and PVA-MOF-525 composite hydrogels capably sorb the therapeutic peptide Angiotensin 1–7. By understanding the influence of polymer-MOF interactions on the structure and properties of composite gels, this work informs and expands the design space of this emerging class of materials.

36 MATERIALS SCIENCE↗

Stochastic machine learning via sigma profiles to build a digital chemical space

This work establishes a different paradigm on digital molecular spaces and their efficient navigation by exploiting sigma profiles. To do so, the remarkable capability of Gaussian processes (GPs), a type of stochastic machine learning model, to correlate and predict physicochemical properties from sigma profiles is demonstrated, outperforming state-of-the-art neural networks previously published. The amount of chemical information encoded in sigma profiles eases the learning burden of machine learning models, permitting the training of GPs on small datasets which, due to their negligible computational cost and ease of implementation, are ideal models to be combined with optimization tools such as gradient search or Bayesian optimization (BO). Gradient search is used to efficiently navigate the sigma profile digital space, quickly converging to local extrema of target physicochemical properties. While this requires the availability of pretrained GP models on existing datasets, such limitations are eliminated with the implementation of BO, which can find global extrema with a limited number of iterations. A remarkable example of this is that of BO toward boiling temperature optimization. Holding no knowledge of chemistry except for the sigma profile and boiling temperature of carbon monoxide (the worst possible initial guess), BO finds the global maximum of the available boiling temperature dataset (over 1,000 molecules encompassing more than 40 families of organic and inorganic compounds) in just 15 iterations (i.e., 15 property measurements), cementing sigma profiles as a powerful digital chemical space for molecular optimization and discovery, particularly when little to no experimental data is initially available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Activity coefficient acquisition with thermodynamics–informed active learning for phase diagram construction

This work explores the use of thermodynamics-informed Gaussian processes (GPs) and active learning (AL) to model activity coefficients and construct phase diagrams. Relying on synthetic data generated from an excess Gibbs energy model, GPs were found to accurately describe the activity coefficients of several binary mixtures across large composition and temperature ranges. Moreover, GPs could estimate their own uncertainty and identify composition/temperature regions where activity coefficient data provide the most information to the models. This was leveraged to build AL algorithms targeted at modeling phase equilibria. In many cases, a single active-learning-acquired data point was sufficient to describe the phase diagrams studied. Lastly, the ability of AL to greatly reduce the amount of data needed to obtain accurate models was further verified on experimental case studies, namely individual ion activity coefficients, the solid–liquid and vapor–liquid equilibrium of deep eutectic solvents, and phase equilibria in ternary mixtures.

25 ENERGY STORAGE↗

How Reproducible are Surface Areas Calculated from the BET Equation?

Abstract Porosity and surface area analysis play a prominent role in modern materials science. At the heart of this sits the Brunauer–Emmett–Teller (BET) theory, which has been a remarkably successful contribution to the field of materials science. The BET method was developed in the 1930s for open surfaces but is now the most widely used metric for the estimation of surface areas of micro‐ and mesoporous materials. Despite its widespread use, the calculation of BET surface areas causes a spread in reported areas, resulting in reproducibility problems in both academia and industry. To prove this, for this analysis, 18 already‐measured raw adsorption isotherms were provided to sixty‐one labs, who were asked to calculate the corresponding BET areas. This round‐robin exercise resulted in a wide range of values. Here, the reproducibility of BET area determination from identical isotherms is demonstrated to be a largely ignored issue, raising critical concerns over the reliability of reported BET areas. To solve this major issue, a new computational approach to accurately and systematically determine the BET area of nanoporous materials is developed. The software, called “BET surface identification” (BETSI), expands on the well‐known Rouquerol criteria and makes an unambiguous BET area assignment possible.

36 MATERIALS SCIENCE↗

The effect of co-adsorbed solvent molecules on H 2 binding to metal alkoxides

The introduction of metal alkoxides has been proposed as an attractive option to enhance hydrogen binding energies in porous materials such as metal–organic frameworks (MOFs) for room-temperature hydrogen storage applications. The presence of residual solvent molecules from MOF synthesis can, however, affect the performance of these functional groups. We perform quantum chemical calculations to predict solvent binding energies onto the metal-alkoxides and the temperatures required to drive off the solvent molecules and successfully activate porous materials with these moieties. Calculations are performed for Li, Mg, Zn, Cu, and Ni alkoxides and chloroform (CHCl 3 ), dimethylformamide (DMF), ethanol, methanol, and water solvent molecules. We identify CHCl 3 as a promising solvent that can be removed from these alkoxides at mild temperatures, whereas DMF binds strongly to the metal alkoxides and removal would require temperatures above the present upper bound of thermal stability in MOFs. As a second objective, we calculated the binding energies of hydrogen to metal alkoxide–solvent complexes to explore the effect of any solvent molecules that cannot be removed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Free energy of metal-organic framework self-assembly

Metal-organic frameworks (MOFs) represent an important class of materials. Careful selection of building blocks allows for tailoring of the properties of the resulting framework. The self-assembly process, however, is not understood, and without detailed knowledge of the underlying molecular mechanism, it is difficult to anticipate whether a particular design can be realized, or whether the material adopts a metastable, kinetically arrested state. We present a detailed examination of early-stage self-assembly pathways of the MOF-5. Enhanced sampling techniques are used to model a self-assembly in an explicit solvent (dimethylformamide, DMF). We identify several free energy barriers encountered during the assembly of the final MOF, which arise from structural rearrangements preceding MOF formation and from disrupted MOF-solvent interactions as formation proceeds. In all cases considered here, MOFs exhibit favorable entropic gains during the assembly. More generally, the strategy presented provides a step toward the experimental design characterizing the formation of ordered frameworks and possible sources of polymorphism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗