Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive quantum chemistry”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Extreme Nanoconfinement Dramatically Enhances Small Molecule Solubility in Nonpolar Polymers

Elucidating gas solubility in confined polymer systems addresses a fundamental gap in polymer physics and has important implications for gas barrier and separation technologies and polymer upcycling reactions. In this study, we examine the solubility of methanol and n-hexane in polystyrene and low-density polyethylene confined within the interstitial pores of disordered silica nanoparticle packings. Using capillary rise infiltration, these polymers are infiltrated into the nanoparticle packings. Gas solubility in these confined polymers is measured using a quartz crystal microbalance. Remarkably, confinement leads to a 10- to ∼100-fold increase in gas solubility. Systematic experiments reveal that pore size plays a dominant role in increasing solubility. In contrast, the molecular weight of the polymer and the surface wetting characteristics of the nanoparticles, achieved by hydrophobically modifying the nanoparticles, have minimal effects. By systematically varying penetrant polarity, polymer crystallinity, pore size, and nanopore surface chemistry, this study isolates geometric nanoconfinement as the dominant factor governing solubility enhancement in confined polymers. Atomistic simulations revealed that confinement and surface−polymer interactions both contribute to enhanced solubility, with polymer packing playing a significant role in modulating gas uptake. Here, these results suggest that the changes in the molecular arrangements of polymer segments underlie the observed trends. This study highlights the potential of confined polymers in engineering the separation performance of membranes and heterogeneously catalyzed polymer upcycling reactions.

alcohols↗
DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2021

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Time-resolved photoelectron imaging of complex resonances in molecular nitrogen

We have used the FERMI free-electron laser to perform time-resolved photoelectron imaging experiments on a complex group of resonances near 15.38 eV in the absorption spectrum of molecular nitrogen, N-2, under jet-cooled conditions. The new data complement and extend the earlier work of Fushitani et al. [Opt. Express 27, 19702-19711 (2019)], who recorded time-resolved photoelectron spectra for this same group of resonances. Time-dependent oscillations are observed in both the photoelectron yields and the photoelectron angular distributions, providing insight into the interactions among the resonant intermediate states. In addition, for most states, we observe an exponential decay of the photoelectron yield that depends on the ionic final state. This observation can be rationalized by the different lifetimes for the intermediate states contributing to a particular ionization channel. Finally, although there are nine resonances within the group, we show that by detecting individual photoelectron final states and their angular dependence, we can identify and differentiate quantum pathways within this complex system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Preparation and Characterization of Inclusion Complexes of β-Cyclodextrin and Phenolics from Wheat Bran by Combination of Experimental and Computational Techniques

Bitterness often associated with whole wheat products may be related to phenolics in the bran. Cyclodextrins (CDs) are known to form inclusion complexes. The objective was to form inclusion complexes between β-CD and wheat phenolics. Pure phenolic acids (trans-ferulic acid (FA), caffeic acid (CA), and p-coumaric acid (CO)) and phenolic acids from wheat bran were used to investigate complex formation potential. Complexes were characterized by spectroscopy techniques, and a computational and molecular modeling study was carried out. The relative amount of complex formation between β-CD and wheat bran extract was CA > CO > FA. The phenolic compounds formed inclusion complexes with β-CDs by non-covalent bonds. The quantum-mechanical calculations supported the experimental results. The most stable complex was CO/β-CD complex. The ΔH value for CO/β-CD complex was -11.72 kcal/mol and was about 3 kcal/mol more stable than the other complexes. The QSPR model showed good correlation between binding energy and 1H NMR shift for the H5 signal. This research shows that phenolics and β-CD inclusion complexes could be utilized to improve the perception of whole meal food products since inclusion complexes have the potential to mask the bitter flavor and enhance the stability of the phenolics in wheat bran.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A family of ionic supersalts with covalent-like directionality and unconventional multiferroicity

Ionic crystals composed of elemental ions such as NaCl are non-polar due to directionless ionic bonding interactions. Here, we show that these can develop polarity by changing their building blocks from elemental ions to superalkali and superhalogen cluster-ions, which mimic the chemistry of alkali and halogen atoms, respectively. Due to the non-spherical geometries of these cluster ions, corresponding supersalts form anisotropic polar structures with ionic bonding, yet covalent-like directionality, akin to sp 3 hybridized systems. Using density functional theory and extensive structure searches, we predict a series of stable ferroelectric/ferroelastic supersalts, PnH 4 MX 4 (Pn = N, P; M = B, Al, Fe; X = Cl, Br) composed of superalkali PnH 4 and superhalogen MX 4 ions. Unlike traditional ferroelectric/ferroelastic materials, the cluster-ion based supersalts possess ultra-low switching barrier and can endure large ion displacements and reversible strain. In particular, PH 4 FeBr 4 exhibits triferroic coupling of ferroelectricity, ferroelasticity, and antiferromagnetism with controllable spin directions via either ferroelastic or 90-degree ferroelectric switching.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Optical properties of a diamond NV color center from capped embedded multiconfigurational correlated wavefunction theory

Diamond defects are among the most promising qubits. Modeling their properties through accurate quantum mechanical simulations can further their development into robust units of information. We use the recently developed capped density functional embedding theory (capped-DFET) with the multiconfigurational n-electron valence second-order perturbation theory to characterize the electronic excitation energies for different spin manifolds of the well-characterized negatively charged substitutional N defect adjacent to a vacancy (V C ) in diamond (N C V C − ). We successfully reproduce vertical excitation energies for both triplet and singlet states of N C V C − with errors < 0.1 eV. Unlike other embedding methods, capped-DFET exhibits robust predictions that are approximately independent of the embedded cluster size: it only requires a cluster to contain the defect atoms and their nearest neighbors (as small as a 40-atom capped cluster). Furthermore, our method is free from slowly converging Coulomb interactions between charged defects, and thus also only weakly dependent on supercell size.

Chemistry↗
DOE OSTI2025

TaTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-titanium (Ta-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Size-Dependent Janus-Ligand Shell Formation on PbS Quantum Dots

We studied the size-dependent Janus ligand shell formation on PbS QDs employing an X-type ligand exchange reaction between native oleate ligands and two substituted cinnamic acid ligands, trifluoromethyl- and dimethyl amino-cinnamic acid, representing electron donating and electron withdrawing ligands. The exchange reactions become significantly more favorable for both electron donating and withdrawing ligands (..delta..G becomes more negative) in the smaller QDs compared to the larger QDs likely because the ligand density is smaller on the larger QDs reducing the strength of the ligand-ligand interactions. We found that Janus-ligand shells form more readily on smaller QDs than on bigger QDs with electron donating ligands. We also observed a dependence on the QD concentration that should be considered when forming Janus-ligand shells. Two-dimensional solution nuclear magnetic resonance spectroscopy (2D-NMR) shows evidence of pronounced phase segregation between oleate and electron donating ligands on the smaller QDs consistent with the enhanced ligand-ligand interactions. This study broadens our understanding of how to construct Janus and patchy ligand shell morphologies on small QDs.

14 SOLAR ENERGY↗
DOE OSTI2026

Hole Transfer to Carbazole Derivatives: Untold Story of “Self-Assembled Monolayers” Employed in “Halide Perovskite Solar Cells”

Carbazole derivatives, such as MeO-2PACz and 2PACz are known to improve the performance of halide perovskite solar cells by facilitating hole transfer. To assess their interaction with halide perovskites, this work probes the hole transfer from excited CsPbBr3 quantum dots to MeO-2PACz and 2PACz using emission spectroscopic and transient absorption techniques. The different oxidation potentials of these two carbazoles result in divergent interactions with CsPbBr 3 QDs. Whereas MeO-2PACz quenches the emission of CsPbBr 3 QDs, 2PACz enhances the emission by remediating the surface traps. Transient absorption studies confirm the formation of MeO-2PACz +• cation radical with characteristic absorption in the near IR region. No such oxidation process was observed with 2PACz. The mechanistic insights into the interaction of the two carbazole derivatives with excited perovskite nanocrystals will add another piece to the untold story behind the improved performance of perovskite photovoltaic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Frustrated Magnetism in FeGe 3 O 4 with a Chiral Trillium Network

The discovery of new magnetic ground states in geometrically frustrated lattices remains a central challenge in materials science. Here, we report the synthesis, structural characterization, and frustrated magnetic properties of FeGe 3 O 4 , a newly identified compound that crystallizes in the noncentrosymmetric cubic space group P 2 1 3. In this structure, Fe atoms form an intricate double-trillium lattice with nearest-neighbor Fe−Fe distances of ∼4.2 Å, while Ge 2+ ions mediate magnetic interactions through Fe− Ge−Fe pathways. Field-dependent magnetization at 2 K shows a pronounced nonlinearity, reaching a maximum moment of 2.55(3) μ B /Fe 2+ at 70 kOe without evidence of saturation. Magnetic susceptibility, heat capacity, and neutron scattering collectively reveal the onset of short-range magnetic interactions near 5 K, with no longrange ordering detected down to 0.06 K. Specific heat measurements demonstrate strong frustration: only ∼34% of the expected magnetic entropy is recovered at 2.4 K. Taken together, these results establish FeGe 3 O 4 as a rare example of a geometrically frustrated trillium lattice magnet, offering a promising platform for exploring exotic quantum magnetic phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Linear versus Bent Uranium(II) Metallocenes—A Local Vibrational Mode Study

Uranium metallocenes have recently attracted attention driven by their use as catalysts in organometallic synthesis. In addition to bent U(IV) and U(III), an U(II) metallocene [(η 5 -C 5 i Pr 5 ) 2 U] was synthesized with an unusual linear Cp–U–Cp angle. Here, in this work, we investigated 22 U(II) metallocenes, (i) assessing the intrinsic strength of the U-ring interactions in these complexes with a novel bond strength measure based on our local vibrational mode analysis and (ii) systematically exploring what makes these U(II) metallocenes bent. We included relativistic effects through the NESCau Hamiltonian and complemented the local mode analysis with natural bonding orbital (NBO) and quantum theory of atoms in molecules (QTAIM) data. Our study led to the following results: (i) reduction of bulky U-ring ligand substituents does not lead to bent complexes for alkyl substituents ( i Pr and i Bu) in contrast to SiMe 3 ring substituents, which are all bent. (ii) The most bent complexes are [(η 5 -C 5 H 4 SiMe 3 ) 2 U] (130°) and [η 5 -P 5 H 5 ) 2 U] (143°). (iii) Linear complexes showed one hybridized NBO with s/d character, while bent structures were characterized by s/d/f mixing. (iv) We did not observe a correlation between the strength of the U-ring interaction and the amount of the ring-U-ring bend; the strongest interaction was found for [η 5 -Cp) 2 U] and the weakest for [η 5 -P 5 H 5 ) 2 U]. In conclusion, our results provide a foundation for the design of U(II) metallocenes with specific physicochemical properties and increased reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Atoms and interatomic bonding synergism inherent in molecular electronic wave functions

The chemical model of matter consists of atoms held together by bonds. The success of this model implies that the physical interactions of the electrons and nuclei in molecules combine into compound interactions that create the bonding. In the quantum mechanical description, the modified atoms in molecules and the bonding synergism are contained in the molecular electronic wave function. So far, only part of this information has been recovered from the wave function. Notably, the atoms have remained unidentified in the wave function. One reason is that conventional energy decomposition analyses formulate separate model wave functions, independent of the actual wave function, to represent “prepared atoms” and preconceived interactions and, then, intuitively catenate the parts. In the present work, the embedded modified atoms and the inherent physical synergisms between them are determined by a unified derivation entirely from the actual molecular valence space wave function. By means of a series of intrinsic orbital and configurational transformations of the wave function, the energy of formation of a molecule is additively resolved in terms of intra-atomic energy changes, interference energies, and quasi-classical, non-classical, and charge-transfer Coulombic interactions. Here, the analysis furnishes an algorithm for the quantitative resolution of the energy of formation, which enables analyses elucidating reaction energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Magnetic‐Field‐Switchable Laser via Optical Pumping of Rubrene

Abstract Volumetric optical imaging of magnetic fields is challenging with existing magneto‐optical materials, motivating the search for dyes with strong magnetic field interactions, distinct emission spectra, and an ability to withstand high photon flux and incorporation within samples. Here, the magnetic field effect on singlet‐exciton fission is exploited to demonstrate spatial imaging of magnetic fields in a thin film of rubrene. Doping rubrene with the high‐quantum yield dye dibenzotetraphenylperiflanthene (DBP) is shown to enable optically pumped, slab waveguide lasing. This laser is magnetic‐field‐switchable: when operated just below the lasing threshold, application of a 0.4 T magnetic field switches the device between nonlasing and lasing modes, accompanied by an intensity modulation of +360%. This is thought to be the first demonstration of a magnetically switchable laser, as well as the largest magnetically induced change in emission brightness in a singlet‐fission material to date. These results demonstrate that singlet‐fission materials are promising materials for magnetic sensing applications and could inspire a new class of magneto‐optical modulators.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Mechanisms of Degradation of Toxic Nerve Agents: Quantum-chemical Insight into Interactions of Sarin and Soman with Molybdenum Dioxide

With the urgent need to provide an efficient and reliable protection for people from fatal chemical and biological weapons, our fundamental understanding of how toxins interact with filters is far from complete. The situation is further complicated by natural difficulties of performing experimental measurements with lethal toxins. Unlike experiments, computational modeling offers an attractive and safe yet reliable way of studying behavior of toxic agents on a variety of substrates at a great level of detail. Here, we report DFT-based quantum-chemical calculations of adsorption and decomposition of DMMP, sarin and soman on MoO 2 (011) surface. Our calculations show that MoO 2 strongly adsorbs toxic nerve agents and quickly decomposes them. Decomposition of DMMP on the MoO 2 (011) surface proceeds via the PO-CH 3 bond breaking and a formation of a surface methoxy group. The calculated activation barrier for this reaction is 131.5 kJ mol -1 . Unlike DMMP, decomposition of sarin and soman proceeds via the dealkylation reaction yielding propene and 3,3-dimethyl-1-butene, respectively. Decomposition of sarin requires a remarkably low energy (53.7 kJ mol- 1 ), whereas the similar reaction in soman requires 50 kJ mol -1 more energy. We also make specific predictions to guide Ambient-Pressure X-ray Photoelectron Spectroscopy (APXPS) experiments on sarin interaction with MoO 2 samples. We conclude that MoO 2 serves as an efficient substrate able of degrading toxins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗