2 months ago
Adaptive Allocation of Monte Carlo Samples for Efficient, Multifidelity Computational Screening of Metal–Organic Frameworks
For applications in gas sensing, purification, and capture, we often wish to search a large set of metal–organic frameworks (MOFs) for the top-K in terms of their Henry coefficients for an adsorbate. A molecular simulation to predict the Henry coefficient of a MOF constitutes a Monte Carlo integration where each sample consists of inserting an adsorbate in the MOF at a random position, orientation, and configuration, then calculating the MOF–adsorbate interaction energy. Our idea is to leverage top-K arm identification algorithms, developed for the multi-armed bandit problem in reinforcement learning, to sequentially and adaptively allocate adsorbate insertions among the MOFs, in a data-driven manner, to obtain the most accurate top-K subset under a fixed insertion budget. By analogy, each MOF is a slot machine in a casino that, upon pulling its arm (inserting an adsorbate), offers a stochastic reward (a noisy estimate of its Henry coefficient) sampled from a static, unknown probability distribution. Each adaptive allocation algorithm (1) proceeds in a feedback loop of (i) allocate adsorbate insertions to MOF(s), (ii) update the running estimates of the Henry coefficients of the MOF(s), then (iii) judiciously allocate adsorbate insertions to the next MOF(s); (2) sequentially dials-up the fidelities of ongoing molecular simulations in the MOFs, giving a multifidelity computational screening; and (3) circumvents the need to hand-craft structural or chemical features of the MOFs for decision making. As a case study, we implement, benchmark, and analyze the sequential halving, successive accepts and rejects, and narrowing exploration (our proposed heuristic) algorithms to adaptively allocate xenon insertions to screen a set of ca. 300 MOFs for the top-K Xe Henry coefficient subset over differing insertion budgets. Provided with a sufficient budget, we find that these adaptive insertion algorithms can significantly reduce (by a factor of 2–3) the simple regret (sum of true minus empirical top-K true Henry coefficients) and error in the top-K subset of MOFs output by a computational screening. By another metric, adaptive insertion allocation provided a ca. 60% discount on the computational cost to identify the top-K MOFs with less than 5% error. We thereby demonstrate that top-K arm identification algorithms may generally be useful for more efficiently screening materials for various properties via Monte Carlo molecular simulations. This efficiency improvement is especially important when adopting more computationally expensive, sophisticated force fields or even ab initio calculations for the potential energy of configurations to lend higher-fidelity screenings.
check out our new paper on adaptively allocating Monte Carlo samples of MOF-adsorbate configurations for efficient, multi-fidelity computational screening of MOFs for an adsorption property using molecular simulations.
pubs.acs.org/doi/full/10....
6
1
0
0
2 months ago
I'll see you at the AIChE conference in 2075!
2
0
1
0
3 months ago
the singular value decomposition is my favorite matrix factorization by far.
if I were to get a tattoo, it would be “A = UΣVᵀ".
cliché for a professor teaching SVD, but in my grad-level “math for chemical engineers” class, I compressed a photo of my dog using the SVD in Julia. 🐶
2
0
0
0
3 months ago
interesting point! beauty/simplicity/convenience => finds more applications. thinking of where I've encountered symmetric matrices: kernels (Gram matrix), adjacency matrix for an undirected graph, Hessian matrix, description of ellipse... there, the symmetry seems natural. SVD/PCA, less.
0
0
0
0
3 months ago
"it is no exaggeration to say that symmetric matrices are the most important matrices the world will ever see."
"if symmetry makes a matrix important, [the] extra property [of having all positive eigenvalues] makes it truly special."
- Gilbert Strang
4
0
1
0
3 months ago
thank you for the spotlight! 😀
1
0
0
0
3 months ago
Conductive Covalent Organic Frameworks as Chemiresistive Sensor Arrays for the Detection and Differentiation of Gasotransmitters
This paper describes a chemiresistive sensor array using four structurally analogous, but chemically distinct, conductive covalent organic frameworks (COFs) (M-COF-DC-8, M = Fe, Co, Ni, and Cu) capable of detecting and differentiating four important gaseous analytes: nitric oxide (NO), carbon monoxide (CO), hydrogen sulfide (H2S), and ammonia (NH3). The COFs were synthesized from the condensation of 2,3,9,10,16,17,23,24-octaamino-metallophthalocyanine precursors with pyrenetetraone linkers resulting in chemically robust and electrically conductive materials. Chemiresistive sensing experiments, together with machine learning to parse the response pattern of the sensor array, show that the M-COF-DC-8 (M = Fe, Co, Ni, Cu) materials can detect and differentiate this suite of oxidizing and reducing gases at parts-per-million concentrations, with theoretical limits of detection (LOD) in the parts-per-billion range in dry N2. Importantly, the COF array containing M-COF-DC-8 (M = Co, Ni, Cu) retains its ability to detect and differentiate these analytes in air and humidity under low power consumption. Spectroscopic investigations reveal that the synthetic control over the identity of the metallophthalocyanine core efficiently tunes material–analyte interactions and, therefore, emergent device performance. The use of highly tunable COFs as the active material in sensor arrays enables low-power, sensitive, and real-time gas detection with future applications in healthcare and personal protection.
a sensor array of conductive COFs, made by Prof. Kat Mirica's group at Dartmouth, can distinguish between NO, CO, NH₃, and H₂S. cool for us to contribute with PCA and k-NN. 😀
pubs.acs.org/doi/10.1021/...
1
0
0
0
3 months ago
it'd be a special X-mas seminar! jello molds served! j.k.
0
0
0
0
3 months ago
cool! (if I remember correctly, you are from Oregon, right? if so, please reach out next time you’re back home, to visit Oregon State and give a seminar!)
1
0
1
0
3 months ago
CC @rociomer.bsky.social @bessvlai.bsky.social
0
0
1
0
3 months ago
after eight years as a ChemE prof., I had a fantastic day when my PhD advisor Prof. Berend Smit visited Oregon State University! 😁
8
0
2
0
4 months ago
my PhD student G. Fabusola trained and tested machine learning algorithms to parse the response pattern of a conductive-MOF sensor array from K. Mirica's group!
👃 the electronic nose could detect and differentiate toxic gases and H₂S/SO₂ mixtures at ppm-levels.
pubs.acs.org/doi/10.1021/...
4
2
0
0
5 months ago
We interrupt our regular programming to announce…
80
23
1
3
5 months ago
pretty cool!
0
0
0
0
5 months ago
in a “it’s a small world” moment, I ran into @bessvlai.bsky.social at Case Western Reserve University in Cleveland, OH. she was giving a seminar in the chemistry department; me, in chemical engineering. great to see you, Bess!
3
0
1
1
6 months ago
new preprint,
"adaptive allocation of Monte Carlo samples for efficient, multi-fidelity computational screening of metal-organic frameworks"
feedback welcome!
chemrxiv.org/engage/chemr...
5
1
0
0
7 months ago
"guidelines for multi-fidelity Bayesian optimization of molecules and materials"
our News & Views article in Nature Computational Science.
rdcu.be/ext6h
4
0
0
0
8 months ago
beavers are cool. glad our mascot is a beaver.
> The fur trade transformed North America but it nearly destroyed the population of several fur-bearers like muskrats and beavers who are critically important to their ecosystem.
🥺
1
0
0
0
9 months ago
🍷solving a linear program for optimal wine blending in Julia
simonensemble.github.io/pluto_nbs/wi...
2
1
0
0
10 months ago
😅 yeah, I think he wanted to go on a walk!
1
0
0
0
10 months ago
a post-fermentation blend of *nine* white wines from Oregon! and a linear program for wine blending.
5
0
1
0
10 months ago
🚰 "Optimizing mixtures of metal–organic frameworks for robust and bespoke passive atmospheric water harvesting" by C. Harriman, Q. Ke, T. Vlugt, A. Howarth, C. Simon.
feedback welcome on our ChemRxiv preprint:
chemrxiv.org/engage/chemr...
6
0
0
0
10 months ago
fascinating: atmospheric water harvesting by indigenous populations on the Canary Islands long ago.
Kennedy & Boreyko. “Bio‐inspired fog harvesting meshes: a review”. Advanced Functional Materials. 2023.
3
1
2
0
10 months ago
finally got to meet Mark Allendorf from Sandia National Lab! currently co-director of the DOE Hydrogen Materials – Advanced Research Consortium (HyMARC). been following his work since grad school.
20
1
0
0
11 months ago
with @cgbischak.bsky.social and @shijingsun.bsky.social at the Automating Chemical Labs Scialog in Tucson! 🌵
4
0
0
0