Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Chemical Calculations”

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 19 records

Quantum chemically calculated Abraham parameters for quantifying and predicting polymer hydrophobicity

The leakage and accumulation of plastic in the environment is a significant and growing problem with numerous detrimental impacts and has led to a push toward the design and development of more environmentally benign materials. To this end, we have developed a quantum chemistry-based model for predicting the mobility of polymer materials from molecular structure. Hydrophobicity is used as a surrogate for mobility given that hydrophobic interactions drive much of the partitioning of contaminants in and out of various environmentally relevant compartments. To model polymer hydrophobicity, we adjusted a previously developed Quantum Chemically Calculated Abraham Parameter model to calculate Abraham parameters of small molecules from molecular structure information. The resulting model predicted the octanol-water partition coefficient (K OW ) of polymer repeating units with a root mean square error (RMSE) of 0.48 (log scale). Additionally, the hydrophobicity of high molecular weight polymer materials was captured through solubility parameters and Nile red staining experiments from the literature and predicted with RMSEs of 1.21 (J/cc) 0.5 and 3.42 nm, respectively. Finally, to test the environmental applicability of the model, the relative adsorption capacity of three polymers was predicted and used to unify sorption isotherms across multiple sorbates and polymer sorbents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

SF6 Negative Ion Formation in Charge Transfer Experiments

In the present work, we report an update and extension of the previous ion-pair formation study of Hubers, M.M.; Los, J. Chem. Phys. 1975, 10, 235–259, noting new fragment anions from time-of-flight mass spectrometry. The branching ratios obtained from the negative ions formed in K + SF6 collisions, in a wide energy range from 10.7 up to 213.1 eV in the centre-of-mass frame, show that the main anion is assigned to SF5− and contributing to more than 70% of the total ion yield, followed by the non-dissociated parent anion SF6− and F−. Other less intense anions amounting to <20% are assigned to SF3− and F2−, while a trace contribution at 32u is tentatively assigned to S− formation, although the rather complex intramolecular energy redistribution within the temporary negative ion is formed during the collision. An energy loss spectrum of potassium cation post-collision is recorded showing features that have been assigned with the help of theoretical calculations. Quantum chemical calculations for the lowest-lying unoccupied molecular orbitals in the presence of a potassium atom are performed to support the experimental findings. Apart from the role of the different resonances participating in the formation of different anions, the role of higher-lying electronic-excited states of Rydberg character are noted.

Biochemistry & Molecular Biology

Quantum Monte Carlo Calculations of Chemical Binding and Reactions

The auxiliary field quantum Monte Carlo method developed by the PIs has been shown to provide the most accurate description of strongly correlated electronic systems, from molecules to solids. Unlike other explicitly many‐body approaches, the quantum Monte Carlo method scales as a low order polynomial of systems size, similar to mean‐field methods such as density functional theory. However, the auxiliary field quantum Monte Carlo algorithm is significantly more expensive than traditional density functional calculations. This creates a bottleneck for applications to extended systems, such as large molecules and solids. One principal objective of this proposal was to develop new auxiliary field quantum Monte Carlo computational strategies to achieve improved scaling with system size, using downfolding and localization schemes, without sacrificing the predictive power of the calculations. A second goal is to extend the reach of auxiliary field quantum Monte Carlo to calculate excited states. This final report summarizes what has been achieved during the course the project toward these goals.

97 MATHEMATICS AND COMPUTING

Activation of methane by U + studied by guided ion beam tandem mass spectrometry and quantum chemistry

Reaction pathways of all products formed in the U + + CH 4 (CD 4 ) reaction were explored as a function of kinetic energy using guided ion beam tandem mass spectrometry and quantum chemical calculations. UH + , UC + , UCH + , UCH 2 + , and UCH 3 + (and their perdeuterated analogues) are formed in endothermic reactions. In both systems, the UCH 2 + (UCD 2 + ) dehydrogenated product was the dominant product in the low-energy region, whereas the UH + (UD + ) hydride product became predominant at high energies. The kinetic energy behavior of the various products is consistent with a common intermediate of H–U + –CH 3 (D–U + –CD 3 ). Here, the kinetic energy dependence of all product cross sections was modeled to obtain experimental bond dissociation energies at 0 K (in eV): D 0 (U + –H) = 2.42 ± 0.10, D 0 (U + –C) = 3.95 ± 0.12, D 0 (U + –CH) = 4.91 ± 0.09, D 0 (U + –CH 2 ) = 4.11 ± 0.04, and D 0 (U + –CH 3 ) = 2.41 ± 0.09. Quantum chemical calculations using the UCCSD(T) and UB3LYP approaches with the cc-pwCVXZ-PP basis set with MDF-60 pseudopotential for U + and the aug-cc-pCVXZ and aug-cc-pVXZ (X = T, Q) basis set for carbon and hydrogen, respectively, validate the experimental bond dissociation energies and outline the potential energy surface for all reactions observed. In addition, spin–orbit corrections of the bond energies for all products were calculated at a CASSCF-CASPT2-RASSI level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Advanced Method Optimization for Sampling and Analysis Instrumentation

This work presents a generalized approach for analytical method optimization that branches the gap between techniques historically employed and accurate modern optimization techniques suitable for various applications. The novelty of the described strategy is the utilization of multivariate, multiobjective optimization with Karush-Kuhn-Tucker conditions to bound the optimization space to solutions within the physical limitations of instrumentation. Briefly, the basic steps outlined in this paper are to (1) determine the objective(s) that should be maximized or minimized based on the goals of the analytical application, (2) conduct a screening experiment, (3) perform ANOVA to determine the parameters which have a statistically significant effect on the objective, (4) conduct an experiment (e.g., Box-Behnken design) to collect data for fitting the objective equation, and (5) determine the physical constraints of the parameters and solve the Lagrangian to determine the optimal method parameters. A broad approach to optimization target selection allows for robust method tuning to develop improved data sets amenable for chemometrics and machine learning algorithm development. Gas chromatography-mass spectrometry was selected as a use case due to its broad use across scientific fields and time-consuming method development involving numerous parameters. In conclusion, this strategy can reduce the cost of research, improve data quality, and enable the rapid development of new analytical technique.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Ground and excited state properties of ThBe and AcBe

In this work, the ground and excited states of ThBe and AcBe were investigated by performing high-level multireference and single-reference coupled-cluster quantum chemical calculations with large correlation consistent basis sets. Full potential energy curves (PECs), chemical bonding patterns, energetics, spectroscopic parameters (T e , r e , ω e , and ω e x e ), and spin–orbit effects of 13 and 8 electronic states of ThBe and AcBe, respectively, are reported. The ground electronic states of ThBe and AcBe are single-reference 1 3 Σ − (1σ 2 2σ 2 1π 2 ) and 1 2 Π (1σ 2 2σ 2 1π 1 ), respectively, and originate from their corresponding ground state fragments. The chemical bonding of ThBe (1 3 Σ − ) and AcBe (1 2 Π) are π-dative in character and are formed by d-electron transfers from Th/Ac to the empty 2p x and 2p y of the Be atom. The electron populations of the f-orbitals of both ThBe (1 3 Σ − ) and AcBe (1 2 Π) are minor which exhibit their “transition-metal-like” nature. The estimated bond energies of the spin–orbit ground states of ThBe (1 3 Σ − 0+ ) and AcBe (1 2 Π 1/2 ) are 12.79 and 11.02 kcal mol −1 , respectively. Finally, the bond energy of ThBe was used to estimate its heat of formation ΔH 0 f (298 K) of 869.61 ± 6 kJ mol −1 .

74 ATOMIC AND MOLECULAR PHYSICS

Advancing molecular machine learning representations with stereoelectronics-infused molecular graphs

Molecular representation is a critical element in our understanding of the physical world and the foundation for modern molecular machine learning. Previous molecular machine learning models have used strings, fingerprints, global features and simple molecular graphs that are inherently information-sparse representations. However, as the complexity of prediction tasks increases, the molecular representation needs to encode higher fidelity information. This work introduces a new approach to infusing quantum-chemical-rich information into molecular graphs via stereoelectronic effects, enhancing expressivity and interpretability. Learning to predict the stereoelectronics-infused representation with a tailored double graph neural network workflow enables its application to any downstream molecular machine learning task without expensive quantum-chemical calculations. We show that the explicit addition of stereoelectronic information substantially improves the performance of message-passing two-dimensional machine learning models for molecular property prediction. We show that the learned representations trained on small molecules can accurately extrapolate to much larger molecular structures, yielding chemical insight into orbital interactions for previously intractable systems, such as entire proteins, opening new avenues of molecular design. Finally, we have developed a web application (simg.cheme.cmu.edu) where users can rapidly explore stereoelectronic information for their own molecular systems.

Boiko, Daniil A

Activation of Semicore Electrons in Alkali Metals and Their Role in the B1–B2 Phase Transition under Pressure

Alkali fluorides are often thought of as archetypical ionic compounds whose structures can be understood in terms of the packing of rigid spheres. At ambient pressure, they assume the rocksalt (B1) structure, while under only a few GPa of pressure, the cesium chloride (B2) structure, with higher coordination numbers, is assumed for KF, RbF, and CsF. NaF requires almost 10 times more pressure to undergo this same phase transition, which has not been observed for LiF. Herein, we provide a detailed analysis, based upon quantum chemical calculations, explaining this behavior. We show that for the heavier alkali metals, the semicore p orbitals engage in metal-metal bonding in the B2 phase, facilitating the pressure-induced B1 → B2 structural transition. Furthermore, these findings suggest that the semicore orbitals of heavy alkali metals can be activated without the need of strong oxidants at very mild levels of compression, resulting in the formation of chemical bonds, challenging both traditional and modern core-valence distinctions. In addition, we argue that Cs 5p-5p bonding in CsCl occurs already at ambient pressure, stabilizing the B2 phase, and suggest experiments that may be able to detect signatures of such bonding.

Anions

Cryogenic Vibrational Spectroscopy of the Deprotonated Dimer of Phosphoric Acid

Phosphate-containing molecules are ubiquitous in nature, where they play crucial roles in biochemical processes. Further, they are of technical importance, for example, in certain batteries and in fuel cells, where a unique property of phosphoric acid is exploited its exceptionally high proton conductivity. Proton transport in phosphoric acid is known to involve proton shuttling; however, the elementary steps involved are not clear. To elucidate the hydrogen bonding preferences of phosphoric acid, we investigate the dihydrogen phosphate anion as well as the deprotonated dimer of phosphoric acid (H 3 PO 4 ·H 2 PO 4 − ) in the gas phase using infrared action spectroscopy in helium nanodroplets and infrared D 2 -tagging photodissociation spectroscopy, and the experimental spectra are compared to theoretical ones. Theory finds for H 3 PO 4 ·H 2 PO 4 − two different structures that are predicted to be nearly isoenergetic. The comparison to the experimental spectra, however, allows for a clear assignment and structure identification. The resulting structure has an interesting binding motif, which might be of relevance to interactions of phosphoric acid in the condensed phase and which can serve as a benchmark for quantum chemical calculations.

Chemical structure

Ultrafast solvent-to-solute proton transfer mediated by intermolecular coherent vibrations

Ultrafast photoinduced excited-state proton transfer (ESPT) plays a crucial role in protecting biomolecules and functional materials from photodamage. However, the influence of solute-solvent interactions on ESPT dynamics remains under active investigation. Here, we present an ultrafast spectroscopic study of ESPT in the photobase 2-(2´-pyridyl)benzimidazole (PBI) in methanol. Ultrafast absorption spectroscopy, supported by quantum chemical calculations, reveals three distinct kinetic steps: (1) a 2.2 ps solvent-to-solute proton transfer, (2) subsequent nonradiative relaxation to the ground state within 31 ps, producing a vibrationally hot ensemble with substantial excess kinetic energy, and (3) equilibration as this energy dissipates into the surrounding solvent bath over 186 ps. Femtosecond-resolved dynamics exhibit oscillatory signals indicative of coherent wavepacket motion on the S 1 potential energy surface. A phase flip in the excited-state absorption maximum confirms this assignment. Fourier analysis resolves two dominant periods (∼117 fs and ∼340 fs), corresponding to in-plane and out-of-plane vibrational modes coupled between PBI and the hydrogen-bonded methanol molecule. The rapid dephasing ( < 300 fs) suggests that the nuclear wavefunction evolves on an anharmonic potential energy surface while traversing the ESPT reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

A Four‐Coordinate Pr 4+ Imidophosphorane Complex

Abstract The imidophosphorane ligand, [NP t Bu 3 ] − ( t Bu= tert ‐butyl), enables isolation of a pseudo‐tetrahedral, tetravalent praseodymium complex, [Pr 4+ (NP t Bu 3 ) 4 ] ( 1‐Pr ), which is characterized by a suite of physical characterization methods including single‐crystal X‐ray diffraction, electron paramagnetic resonance, and L 3 ‐edge X‐ray near‐edge spectroscopies. Variable‐temperature direct‐current magnetic susceptibility data, supported by multiconfigurational quantum chemical calculations, demonstrate that the electronic structure diverges from the isoelectronic Ce 3+ analogue, driven by increased crystal field. The four‐coordinate environment around Pr 4+ in 1‐Pr , which is unparalleled in reported extended solid systems, provides a unique opportunity to study the interplay between crystal field splitting and spin‐orbit coupling in a molecular tetravalent lanthanide within a pseudo‐tetrahedral coordination geometry.

Boggiano, Andrew C.

A Four‐Coordinate Pr 4+ Imidophosphorane Complex

Abstract The imidophosphorane ligand, [NP t Bu 3 ] − ( t Bu= tert ‐butyl), enables isolation of a pseudo‐tetrahedral, tetravalent praseodymium complex, [Pr 4+ (NP t Bu 3 ) 4 ] ( 1‐Pr ), which is characterized by a suite of physical characterization methods including single‐crystal X‐ray diffraction, electron paramagnetic resonance, and L 3 ‐edge X‐ray near‐edge spectroscopies. Variable‐temperature direct‐current magnetic susceptibility data, supported by multiconfigurational quantum chemical calculations, demonstrate that the electronic structure diverges from the isoelectronic Ce 3+ analogue, driven by increased crystal field. The four‐coordinate environment around Pr 4+ in 1‐Pr , which is unparalleled in reported extended solid systems, provides a unique opportunity to study the interplay between crystal field splitting and spin‐orbit coupling in a molecular tetravalent lanthanide within a pseudo‐tetrahedral coordination geometry.

Boggiano, Andrew C.

Au147(SPh)30(PPh3)12: A Geometrically Closed, but Electronically Open Triple‐Shell Icosahedral Gold Cluster and its Geometrically Open Counterpart

The aesthetic platonic solids have been known since ancient times, and the structure of all five platonic solids is also found in chemical compounds. While gold sub-nanometer clusters and gold nanoparticles with an icosahedral structure have been known for a long time to exist, a multi-shell icosahedral gold cluster at the intermediate size between 13 and thousands of atoms has been elusive. Here we present the synthesis and crystallographic characterization of the first triple-shell icosahedral metal cluster, Au147(SPh)30(PPh3)12 1. The gold core in 1 is stabilized by phosphines and thiolates, but surprisingly no staple motifs are formed. A second cluster, Au146(SPh)30(PPh3)12 2, cocrystallizes and is identified as having a closed electronic shell but can be considered as a geometrically open pendant of 1. The unique clusters are characterized experimentally by EDX, UV/vis, DLS, and EPR and theoretically by quantum chemical calculations.

Strienz, Markus

A Thorough Characterization of the Tellurocyanate Anion

Tellurocyanate, [TeCN] − , is the heaviest group 16 congener of the cyanate anion, [OCN] − . Due to the relative instability of the C─Te bond, tellurocyanate chemistry has seen only scarce attention. Here, we present the facile synthesis and thorough characterization of [K@crypt-222][TeCN]. The anion is essentially linear with interatomic distances C─N = 1.150(6)Å and C─Te = 2.051(4)Å, thus approximating a C≡N triple bond and for C─Te a bond order between 1 and 2. Fully 13 C and 15 N labeled [Te 13 C 15 N] − allowed for the extraction of chemical shifts and all possible coupling constants ( 13 C = 77.8 ppm, 15 N = 285.7 ppm, 125 Te = −566 ppm, 1 J 13C-15N = 8 Hz, 1 J 13C-125Te = 748 Hz, 2 J1 5N-125Te = 55 Hz), which were also determined independently by quantum chemical calculations. In the series [ChCN] − (Ch = O─Te), [TeCN] − shows the strongest spin-orbit coupling (SOC) induced heavy-atom effect on the light-atom shielding (SO-HALA-effect). In contrast, 15 N shifts are also well described without considering relativistic effects and/or SOC. Negative-ion photoelectron spectroscopy was used to extract the electron affinity (EA = 3.034 eV) and spin-orbit splitting (3807 cm −1 ) of [TeCN] • . These values continue the trends of falling EA and rising SOC in the series [ChCN] • .

Bonding analysis

Redox‐Dependent Luminescence in a Pentavalent Uranyl Complex

The chemistry of pentavalent uranyl is of high interest due to the ion's central role as an intermediate in nuclear separations. However, due to the instability of [UO 2 ] + , little is known about its electronic structure. The synthesis and characterization of a pentavalent uranyl complex supported by a pyridine dipyrrolide pincer‐type ligand are reported. The pentavalent species is isolated following reduction of the hexavalent uranyl congener, allowing for direct comparisons to be made between the two species. The [UO 2 ] 2+ analog is nonemissive, but upon reduction, luminescence is “switched on” in the pentavalent uranyl complex, affording emission from a high‐energy state. Quantum chemical calculations support the experimental findings and help provide insight into the electronic structure of the excited states of the reduced complex. This work provides a unique example of ligand‐based non‐Kasha emission from a molecular pentavalent uranyl complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH