Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oligomerization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Metal-organic framework supported single-site nickel catalysts for butene dimerization

Homotopic sites in a well-controlled environment are not only ideal systems for mechanistic studies, but also allow optimal control of catalytic transformations. Sites having only a single metal cation and sites consisting of metal oxo complexes with few nickel (Ni) cations supported on the nodes of UiO-66 metal-organic framework (Ni-UiO-66) are studied for 1-butene dimerization. Monomeric Ni sites, which bind to the Zr 6 node via two Zr-OH(µ3) linkages, are active and selective for the dimerization of 1-butene to linear and mono-branched C 8 isomers. Ni oxo complexes with few Ni cations show lower activity and promote the oligomerization of transiently formed C 8 isomers. In conclusion, Kohn-Sham density function theory calculations combined with spectroscopic measurements and kinetic analyses indicate that dimerization follows a Cossee-Arlman reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intrinsic activation energies for ring contraction of allylic cations in zeolites

Cyclic carbocations are important intermediates in zeolite-catalyzed chemistries such as methanol-to-hydrocarbons conversion, naphthenes ring opening, or coke formation. While information about their thermodynamic stability exists, little is known about the kinetics of formation and transformation of cyclic carbocations in zeolites. To fill this knowledge gap, ring contraction of the 1,3,5,5-tetramethylcyclohexenyl cation (C 10 H 17 + ), a representative of 6-membered ring allylic cations, was investigated by in situ UV–vis and IR spectroscopy. Protonic forms of zeolites served as catalysts, at temperatures from 80 °C to 135 °C. Significant oligomerization and hydride transfer in BEA and FAU hampered kinetics analysis, whereas ring contraction dominated in the channels of MOR. The reactant cation, characterized by an electronic absorption at 314 nm and an allylic stretch at 1549 cm −1 , contracted to both a 1,3-alkyl-substituted cyclopentenyl cation (287 nm and 1506 cm −1 ) and a 1,2,3-alkyl-substituted cyclopentenyl cation (297 nm and 1489 cm −1 ). Collection of time-resolved IR spectra and fitting of the intensities with various kinetic models revealed a third transformation, which is expected from thermodynamics: the 1,3-alkyl-substituted cyclopentenyl cation isomerizes to the 1,2,3-alkyl-substituted cyclopentenyl cation. Series of IR spectra recorded at different temperatures delivered intrinsic activation enthalpies (entropies) in MOR of 67 ± 2 kJ mol −1 (−130 ± 6 J mol −1 K −1 ) and 90 ± 3 kJ mol −1 (−70 ± 8 J mol −1 K −1 ) for the contraction to 1,3- and 1,2,3-substituted species, and of 88 ± 5 kJ mol −1 (−84 ± 12 J mol −1 K −1 ) for the isomerization. The findings characterize one path – via contraction of larger rings – to different cyclopentenyl species in zeolites; and the associated, moderate activation energies suggest such transformations contribute to many complex hydrocarbon reaction networks.

Acid catalysis↗

Properties and Autoignition Reactivity of Diesel Boiling Range Ethers Produced from Guerbet Alcohols

We examine the properties of diesel boiling range ethers made from coupling of alcohols produced by oligomerization of ethanol (Guerbet alcohols) for their utility as low-carbon liquid fuel blendstocks. Basic properties of boiling point, flash point, freezing point, density and viscosity are well suited for blending into diesel fuels. For the mixture of ethers the lightest component, di-n-butyl ether, can be present at up to 20 vol% while still having adequately high flashpoint for safe handling. Soot formation tendency (as yield sooting index) is well below that of conventional diesel. The ethers have similar compatibility with elastomers as conventional diesel, based on Hansen solubility parameter analysis. Oxidation stability was assessed for 30 vol% blends of individual ethers in a conventional diesel fuel using a long-term storage test. Over 6 weeks we observed no formation of peroxides or degradation. n-alkyl ethers with carbon number of 8 or higher have cetane number over 100, which is outside the defined range of cetane number, while branched ethers are over 70. The ethers also blend antagonistically into conventional diesel for cetane number, meaning that the blend cetane value is lower than predicted based on a linear by volume, mass, or mole model. We show that aromatics and naphthenes likely act as radical scavengers to slow or shut down autoignition of the highly reactive ethers at low to medium blend levels. Overall, diesel boiling range ethers show significant promise as high quality low-net carbon diesel blendstocks.

09 BIOMASS FUELS↗

Domain crossover in the reductase subunit of NADPH-dependent assimilatory sulfite reductase

NADPH-dependent assimilatory sulfite reductase (SiR) from Escherichia coli performs a six-electron reduction of sulfite to the bioavailable sulfide. SiR is composed of a flavoprotein (SiRFP) reductase subunit and a hemoprotein (SiRHP) oxidase subunit. There is no known high-resolution structure of SiR or SiRFP, thus we do not yet fully understand how the subunits interact to perform their chemistry. Here, we used small-angle neutron scattering to understand the impact of conformationally restricting the highly mobile SiRFP octamer into an electron accepting (closed) or electron donating (open) conformation, showing that SiR remains active, flexible, and asymmetric even with these conformational restrictions. From these scattering data, we then model the first solution structure of SiRFP. Further, computational modeling of the N-terminal 52 amino acids that are responsible for SiRFP oligomerization suggests an eight-helical bundle tethers together the SiRFP subunits to form the SiR core. Finally, mass spectrometry analysis of the closed SiRFP variant show that SiRFP is capable of inter-molecular domain crossover, in which the electron donating domain from one polypeptide is able to interact directly with the electron accepting domain of another polypeptide. This structural characterization suggests that SiR performs its high-volume electron transfer through both inter- and intramolecular pathways between SiRFP domains and, thus, cis or trans transfer from reductase to oxidase subunits. Such highly redundant potential for electron transfer makes this system a potential target for designing synthetic enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Higher-order SPOP assembly reveals a basis for cancer mutant dysregulation

The speckle-type POZ protein (SPOP) functions in the Cullin3-RING ubiquitin ligase (CRL3) as a receptor for the recognition of substrates involved in cell growth, survival, and signaling. SPOP mutations have been attributed to the development of many types of cancers, including prostate and endometrial cancers. Prostate cancer mutations localize in the substrate-binding site of the substrate recognition (MATH) domain and reduce or prevent binding. However, most endometrial cancer mutations are dispersed in seemingly inconspicuous solvent-exposed regions of SPOP, offering no clear basis for their cancer-causing and peculiar gain-of-function properties. Herein, we present the first structure of SPOP in its oligomeric form, uncovering several new interfaces important for SPOP self-assembly and normal function. Given that many previously unaccounted-for cancer mutations are localized in these newly identified interfaces, we uncover molecular mechanisms underlying dysregulation of SPOP function, with effects ranging from gross structural changes to enhanced self-association, and heightened stability and activity.

59 BASIC BIOLOGICAL SCIENCES↗

Direct observation of key aluminum hydroxide prenucleation oligomers for gibbsite nucleation and crystallization in sodium aluminate solution by liquid ToF-SIMS

The mechanism of gibbsite (aluminum hydroxide) crystallization from highly alkaline solutions such as Bayer liquors remains poorly understood, where aluminum (Al) transforms from largely tetrahedrally coordinated aluminate monomers in sodium aluminate solutions into a network of octahedra in gibbsite crystals. A variety of traditional analytical approaches applied to this system do not readily reveal the presence of higher-order oligomeric intermediates. To overcome this limitation, we employed in-situ liquid Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) to examine the Al species present in concentrated sodium aluminate solutions favoring crystallization of gibbsite or sodium aluminate. A complex mixture of Al oligomers was found. By comparing the change in the relative concentration of Al oligomers with +1 and -1 charge, we were able to identify three major Al oligomer candidates, iso-tetramers, iso-pentamers, and cyclic-hexamers, for the nucleation and crystallization of gibbsite. The concentrations of iso-tetramers and iso-pentamers significantly surpass those of cyclic-hexamers. Time-dependent in-situ Raman spectroscopy analysis indicated that the appearance of gibbsite coincided with the peak concentration of these oligomers. The Density-functional theory (DFT) calculation suggests that the formation of iso-oligomers is more favorable than that of cyclic-hexamers. The combined results suggest that iso-tetramers and iso-pentamers play the most substantial role in the nucleation and growth of gibbsite in the sodium aluminate solutions. Our findings also suggest that the oligomers that promote gibbsite crystallization are more stable in dilute sodium aluminate solutions, making these solutions particularly suitable for efficient gibbsite crystallization. In conclusion, our study fills a major knowledge gap in understanding Al speciation that leads to the nucleation and crystallization of gibbsite in concentrated sodium aluminate solutions.

Bayer liquor↗

Expression, purification and characterization of human proton-coupled oligopeptide transporter 1 hPEPT1

The human peptide transporter hPEPT1 (SLC15A1) is responsible for uptake of dietary di- and tripeptides and a number of drugs from the small intestine by utilizing the proton electrochemical gradient, and hence an important target for peptide-like drug design and drug delivery. hPEPT1 belongs to the ubiquitous major facilitator superfamily that all contain a 12TM core structure, with global conformational changes occurring during the transport cycle. Several bacterial homologues of these transporters have been characterized, providing valuable insight into the transport mechanism of this family. Here we report the overexpression and purification of recombinant hPEPT1 in a detergent-solubilized state. Thermostability profiling of hPEPT1 at different pH values revealed that hPEPT1 is more stable at pH 6 as compared to pH 7 and 8. Micro-scale thermophoresis (MST) confirmed that the purified hPEPT1 was able to bind di- and tripeptides respectively. To assess the in-solution oligomeric state of hPEPT1, negative stain electron microscopy was performed, demonstrating a predominantly monomeric state.

59 BASIC BIOLOGICAL SCIENCES↗

Native Chemical Ligation of Peptoid Oligomers

Bioorganic chemists are inspired by natural biopolymers to design peptidomimetic oligomers that can exhibit sequence-structure-function relationships. Biomimetic polymers can be synthesized to incorporate a specific sequence of nonbiological monomer units using a variety of iterative solution-phase or solid-phase reaction schemes. These protocols generally provide access to a vast diversity of oligomeric compounds but are limited with respect to their ability to attain protein-like chain lengths. This constraint can preclude access to sequence-defined synthetic macromolecules with sufficient sizes required to exhibit tertiary structure and other protein-mimetic attributes. In contrast, peptide chemists have overcome this limitation by developing convergent synthetic methods, such as native chemical ligation, to join individual, smaller peptide chains together to make larger peptides or full proteins. A similar convergent approach is needed to establish efficient synthetic routes to non-natural sequence-defined macromolecules. Herein, we adapt the peptide native chemical ligation method to peptoid oligomers, demonstrating how short chains can be conjoined to create sequence-defined peptoid macromolecules. Nanosheet-forming peptoid polymers with distinct surface loop display domains were generated by sequential ligation of several discrete fragments. This method provides a reliable convergent ligation route for sequence-defined polypeptoids that results in a native amide bond joining the fragments. We envision that this strategy will be useful in synthesizing peptoid-based proteomimetics that incorporate diverse chemical features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling the Interface between Salts, Solvates, Co-crystals, and Ionic Liquids with Non-stoichiometric Protic Azolium Azolates [Plus Supporting Information]

In this work, a non-stoichiometric approach to control the solid-state behavior of protic ionic liquids (PILs) was demonstrated by direct mixing of 4,5-dicyanoimidazole (HDCNim) with either 1-ethylimidazole (C 2 im) or 1-butylimidazole (C 4 im) in different mole fractions. Isolation and characterization of three crystalline materials (all having melting points < 100 °C, thus fitting the PIL definition) revealed three different but closely related systems. The structure of [HC 2 im][DCNim] consists of a confused proton system, which is so named for the fast proton exchange between basic and acidic fragments and the formation of hydrogen-bonded anionic oligomers. The compound [HC 2 im][DCNim]·HDCNim consists of an oligomeric delocalized anion (a confused proton located between two azole fragments, in which one acts as a neutral moiety while the other as an ionic fragment, DCNim-H-DCNim) and [HC 4 im][DCNim]·HDCNim can be classified as a salt co-crystal. We believe that these observations will allow a deeper understanding of the behavior of “confused protons” and provide additional strategies for controlling the properties of important classes of materials such as active pharmaceutical ingredients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Elucidating Biomass-Derived Pyrolytic Lignin Structures from Demethylation Reactions through Density Functional Theory Calculations

Pyrolytic lignin is a fraction of pyrolysis oil that contains a wide range of phenolic compounds that can be used as intermediates to produce fuels and chemicals. However, the characteristics of the raw lignin structure make it difficult to establish a pyrolysis mechanism and determine pyrolytic lignin structures. Herein this study proposes dimer, trimer, and tetramer structures based on their relative thermodynamic stability for a hardwood lignin model in pyrolysis. Different configurations of oligomers were evaluated by varying the positions of the guaiacyl (G) and syringyl (S) units and the bonds βO4 and β5 in the hardwood model lignin through electronic structure calculations. The homolytic cleavage of βO4 bonds is assumed to occur and generate two free radical fragments. These can stabilize by taking hydrogen radicals that may be in solution during the intermediate liquid (pathway 1) formation before the thermal ejection. An alternative pathway (pathway 2) could occur when the radicals use intramolecular hydrogen, turning themselves into stable products. Subsequently, a demethylation reaction can take place, thus generating a methane molecule and new oligomeric lignin-derived molecules. The most probable resulting structures were studied. We used FTIR and NMR spectra of selected model compounds to evaluate our calculation approach. Thermophysical properties were calculated using group contribution methods. The results give insights into the lignin oligomer structures and how these molecules are formed. They also provide helpful information for the design of pyrolysis oil separation and upgrading equipment.

09 BIOMASS FUELS↗

Theory-Guided Inelastic Neutron Scattering of Crystalline Alkaline Aluminate Salts Bearing Principal Motifs of Solution-State Species

Aluminate salts precipitated from caustic alkaline solutions exhibit a correlation between the anionic speciation and the identity of the alkali cation in the precipitate, with the aluminate ions occurring either in monomeric (Al(OH) 4 – ) or dimeric (Al 2 O(OH) 6 2– ) forms. The origin of this correlation is poorly understood as are the roles that oligomeric aluminate species play in determining the solution structure, prenucleation clusters, and precipitation pathways. Characterization of aluminate solution speciation with vibrational spectroscopy results in spectra that are difficult to interpret because the ions access a diverse and dynamic configurational space. To investigate the Al(OH) 4 – and Al 2 O(OH) 6 2– anions within a well-defined crystal lattice, inelastic neutron scattering (INS) and Raman spectroscopic data were collected and simulated by density functional theory for K 2 [Al 2 O(OH) 6 ], Rb 2 [Al 2 O(OH) 6 ], and Cs[Al(OH) 4 ]·2H 2 O. These structures capture archetypal solution aluminate species: the first two salts contain dimeric Al 2 O(OH) 6 2– anions, while the third contains the monomeric Al(OH) 4 – anion. Here, comparisons were made to the INS and Raman spectra of sodium aluminate solutions frozen in a glassy state. In contrast to solution systems, the crystal lattice of the salts results in well-defined vibrations and associated resolved bands in the INS spectra. The use of a theory-guided analysis of the INS of this solid alkaline aluminate series revealed that differences were related to the nature of the hydrogen-bonding network and showed that INS is a sensitive probe of the degree of completeness and strength of the bond network in hydrogen-bonded materials. Results suggest that the ionic size may explain cation-specific differences in crystallization pathways in alkaline aluminate salts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of the Organic Solvent Polarity in Isolating Uranyl Peroxide Capsule Fragments

Uranyl peroxide capsules are a fascinating class of polyoxometalates (POMs), discovered only in the 21st century. Understanding the reactivity between peroxide, alkali cations, and uranyl in alkaline solutions is important in nuclear science disciplines including mineralogy, nuclear energy, and legacy nuclear wastes. In this work, we have developed a general procedure to isolate different fragments of the uranyl-peroxide POM capsules, using organic solvents to partially remove K + salts from crude solids of the monomer building block UO 2 (O 2 ) 3 4– (K–U 1 ), leading to stabilization of these reactive fragments. Higher polarity organic solvents remove more K + salts from the crude solid, owed to higher solubility, resulting in more extensive linking of uranyl peroxide building units. By this strategy we have isolated and structurally characterized a dimer K 6 [(UO 2 ) 2 (O 2 ) 4 (OH) 2 ]·7H 2 O (K–U 2 ) and a hexamer face frequently observed in the capsules, K 12 [(UO 2 ) 6 (O 2 ) 9 (OH) 6 ]$·x$H 2 O(K–U 6 ). Comparing experimental and computed Raman spectra shows that these intermediates crystallize by a solid-to-solid transformation, via polymerization of the monomer building block. By small-angle X-ray scattering (SAXS), we track the conversion of the fragments to POM capsules; the reaction rate increases from K–U 1 (days) < K–U 2 (hours) < K–U 6 (instantaneous). This study provides a general synthetic procedure to isolate metastable uranyl peroxide oligomers and control the oligomerization, which will be later applied to systems with the heavier alkalis that are even less stable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kinetics, Products, and Brown Carbon Formation by Aqueous-Phase Reactions of Glycolaldehyde with Atmospheric Amines and Ammonium Sulfate

Glycolaldehyde (GAld) is a C 2 water-soluble aldehyde produced during the atmospheric oxidation of isoprene and many other species and is commonly found in cloudwater. Previous work has established that glycolaldehyde evaporates more readily from drying aerosol droplets containing ammonium sulfate (AS) than does glyoxal, methylglyoxal, or hydroxyacetone, which implies that it does not oligomerize as quickly as these other species. Here, we report NMR measurements of glycolaldehyde’s aqueous-phase reactions with AS, methylamine, and glycine. Reaction rate constants are smaller than those of respective glyoxal and methylglyoxal reactions in the pH range of 3–6. In follow-up cloud chamber experiments, deliquesced glycine and AS seed particles were found to take up glycolaldehyde and methylamine and form brown carbon. At very high relative humidity, these changes were more than 2 orders of magnitude faster than predicted by our bulk liquid NMR kinetics measurements, suggesting that reactions involving surface-active species at crowded air–water interfaces may play an important role. The high-resolution liquid chromatography–electrospray ionization–mass spectrometric analysis of filter extracts of unprocessed AS + GAld seed particles identified sugar-like C 6 and C 12 GAld oligomers, including proposed product 3-deoxyglucosone, with and without modification by reactions with ammonia to diimine and imidazole forms. Chamber exposure to methylamine gas, cloud processing, and simulated sunlight increased the incorporation of both ammonia and methylamine into oligomers. Many C 4 –C 16 imidazole derivatives were detected in an extract of chamber-exposed aerosol along with a predominance of N-derivatized C 6 and C 12 glycolaldehyde oligomers, suggesting that GAld is capable of forming brown carbon SOA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hotspot Coevolution Is a Key Identifier of Near-Native Protein Complexes

Protein–protein interactions play a key role in mediating numerous biological functions, with more than half the proteins in living organisms existing as either homo- or hetero-oligomeric assemblies. Protein subunits that form oligomers minimize the free energy of the complex, but exhaustive computational search-based docking methods have not comprehensively addressed the challenge of distinguishing a natively bound complex from non-native forms. Current protein docking approaches address this problem by sampling multiple binding modes in proteins and scoring each mode, with the lowest-energy (or highest scoring) binding mode being regarded as a near-native complex. However, high-scoring modes often match poorly with the true bound form, suggesting a need for improvement of the scoring function. Here in this study, we propose a scoring function, KFC-E, that accounts for both conservation and coevolution of putative binding hotspot residues at protein–protein interfaces. We tested KFC-E on four benchmark sets of unbound examples and two benchmark sets of bound examples, with the results demonstrating a clear improvement over scores that examine conservation and coevolution across the entire interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Deciphering the Cofilin Oligomers via Intermolecular Disulfide Bond Formation: A Coarse-Grained Molecular Dynamics Approach to Understanding Cofilin’s Regulation on Actin Filaments

Cofilin, a key actin-binding protein, orchestrates the dynamics of the actomyosin network through its actin-severing activity and by promoting the recycling of actin monomers. Recent experimental work suggests that cofilin also forms functionally distinct oligomers through thiol post-translational modification (PTM) that encourages actin nucleation and assembly. Despite these advances, the structural conformations of cofilin oligomers that modulate actin activity remain elusive because there are combinatorial ways to oxidize thiols in cysteines to form disulfide bonds rapidly. This study employs molecular dynamics simulations to investigate human cofilin 1 as a case study for exploring cofilin dimers via disulfide bond formation. Using the free energy profiling, our simulations unveil a range of probable cofilin dimer structures not represented in current Protein Data Bank entries. These candidate dimers are characterized by their distinct population distributions and relative free energies. Of particular note is a dimer featuring an interface between cysteines 139 and 147 residues, which demonstrates stable free energy characteristics and intriguingly symmetrical geometry. In contrast, the experimentally proposed dimer structure exhibits a less stable free energy profile. Here, we also evaluate frustration quantification based on the energy landscape theory in the protein-protein interactions at the dimer interfaces. Notably, the 39-39 dimer configuration emerges as a promising candidate for forming cofilin tetramers, as substantiated by frustration analysis. Additionally, docking simulations with actin filaments further evaluate the stability of these cofilin dimer-actin complexes. Our findings thus offer a computational framework for understanding the role of thiol post-translational modification cofilin proteins in regulating oligomerization, and the subsequent cofilin-mediated actin dynamics in the actomyosin network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger–Weber potentials. In this study, the model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. In conclusion, the model enables simulation of early stage Al(OH) 3 nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗