Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “General Physics”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Molecular Dynamics Study of Combustion Reactions in a Supercritical Environment. Part 2: Boxed MD Study of CO + OH → CO 2 + H Reaction Kinetics

Oxy-fuel combustion technology holds a great promise in both increasing the efficiency of the energy conversion and reducing environmental impact. Yet, effects of the higher pressures and replacement of the nitrogen with carbon dioxide diluent are not well understood at present. The title reaction is one of the most important processes in combustion. Despite numerous studies, the effects of supercritical carbon dioxide environment did not receive much attention in the past. In this work, we report the results of boxed molecular dynamics simulations of these effects at QM/MM theory level with periodical boundary conditions. The free energy barriers for HOCO intermediate formation and decomposition were tabulated in a wide range of pressures (1–1000 atm) and temperatures (400–1600 K). Pressure dependence of calculated rate constants for these reaction steps and overall reaction were analyzed. We discovered that the CO 2 environment may increase these rate constants up to a factor of 25, at near critical conditions. At higher temperatures, this effect weakens significantly. Numerical values for parameters of extended Arrhenius equation, suitable for combustion kinetic modeling are reported.

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

Molecular Dynamics Study of Combustion Reactions in Supercritical Environment. Part 3: Boxed MD Study of CH 3 + HO 2 → CH 3 O + OH Reaction Kinetics

The kinetics of reaction CH 3 + HO 2 → CH 3 O + OH in supercritical carbon dioxide media at pressures from 0.3 to 1000 atm in the temperature range (600–1600) K was studied using boxed molecular dynamics simulations at QM/MM theory level with periodical boundary conditions. The mechanism of this process includes two consecutive steps: formation and decomposition of CH 3 OOH intermediate. We calculated the activation free energies and rate constants of each step, then used Bodenstein’s quasistationary concentrations approximation to estimate the rate constants of the reaction. On the basis of the temperature dependence of the rate constants, parameters in the extended Arrhenius equation were determined. Furthermore, we found that reaction rate of each step, as well as overall reaction, increases with increasing CO 2 pressure in the system. The most capable zone for the process is T = 1000–1200 K, and the CO 2 pressure is about 100 atm.

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

Effect of Coulomb Correlation on the Magnetic Properties of Mn Clusters

In spite of decades of research, a fundamental understanding of the unusual magnetic behavior of small Mn clusters remains a challenge. Experiments show that Mn 2 is antiferromagnetic while small clusters containing up to five Mn atoms are ferromagnetic with magnetic moments of 5 μB/atom and become ferrimagnetic as they grow further. Theoretical studies based on density functional theory (DFT), however, find Mn2 to be ferromagnetic, with ferrimagnetic order setting in at different sizes that depend upon the computational methods used. While quantum chemical techniques correctly account for the antiferromagnetic ground state of Mn 2 , they are computationally too demanding to treat larger clusters, making it difficult to understand the evolution of magnetism. These studies clearly point to the importance of correlation and the need to find ways to treat it effectively for larger clusters and nanostructures. In this report we show that the DFT+U method can be used to account for strong correlation. We determine the on-site Coulomb correlation, Hubbard U self-consistently by using the linear response theory and study its effect on the magnetic coupling of Mn clusters containing up to five atoms. With a calculated U value of 4.8 eV, we show that the ground state of Mn 2 is antiferromagnetic with a Mn–Mn distance of 3.34 Å, which agrees well with the electron spin resonance experiment. Equally important, we show that on-site Coulomb correlation also plays an important role in the evolution of magnetic coupling in larger clusters, as the results differ significantly from standard DFT calculations. We conclude that for a proper understanding of magnetism of Mn nanostructures (clusters, chains, and layers) one must take into account the effect of strong correlation.

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

Strong Effect of Organic Ligands on the Electronic Structure of Metal-Chalcogenide Clusters

It is shown that multiple ionization energies of metal-chalcogenide clusters can be substantially reduced by adding ligands that form charge transfer complexes. In this report we demonstrate this intriguing phenomenon by considering metal-chalcogenide clusters including cases where a cluster has a filled electronic shell with a large gap between the occupied and unoccupied states reminiscent of stable species. The studies include a Co 6 Se 8 core ligated with tri-ethylphosphine (PEt 3 ) ligands forming a stable Co 6 Se 8 (PEt 3 ) 6 species. All of the ligated clusters have a first ionization energy in the range for alkali atoms and multiple ionization energies that are considerably lower than those for the non-ligated clusters. The change in electronic behavior upon ligation can be associated with a shift in the electronic spectrum via a crystal field like effect due to attaching ligands that form charge transfer complexes. We also show that metal-chalcogenide species can be programmed by proper ligand replacement to promote dimerization by first forming the Co 6 Se 8 (PEt 3 ) n (CO) 6-n (n = 0–6) clusters where the CO ligands could be replaced by diisocyanide (CNC 6 H 4 NC) ligands. The diisocyanide ligand acts as a rigid linker between the metallic cores, enabling the formation of a Co 6 Se 8 (PEt 3 ) 5 (CNC 6 H 4 NC)Co 6 Se 8 (PEt 3 ) 5 superatomic molecule (SM), and we examine the electronic and magnetic properties of the recently synthesized SM via studies on an analogous SM with smaller ligands.

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

Simple Flux-Side Formulation of State-Resolved Thermal Reaction Rates for Ring-Polymer Surface Hopping

Employing the recently developed isomorphic Hamiltonian framework for including nuclear quantum effects in mixed quantum–classical nonadiabatic dynamics, we present a flux-side formulation of state-resolved thermal reaction rates for ring-polymer surface hopping (iso-RPSH). An appealing aspect of the new approach is that calculation of multiple state-resolved nonadiabatic thermal reaction rates is enabled with only a single free-energy surface calculation, whereas previous nonadiabatic flux-side formulations for surface hopping involve multiple free-energy surface calculations. In conclusion, the method is shown to be robust and straightforwardly implemented, and numerical results reveal that RPSH in the isomorphic Hamiltonian framework leads to better dividing surface independence than alternative RPSH methods due to improved preservation of the path-integral statistics.

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

Neural Network Water Model Based on the MB-Pol Many-Body Potential

The MB-pol many-body potential accurately predicts many properties of water, including cluster, liquid phase, and vapor–liquid equilibrium properties, but its high computational cost can make applying it in large-scale simulations quite challenging. Here, in order to address this limitation, we developed a “deep potential” neural network (DPMD) model based on the MB-pol potential for water. We find that a DPMD model trained on mostly liquid configurations yields a good description of the bulk liquid phase but severely underpredicts vapor–liquid coexistence densities. By contrast, adding cluster configurations to the neural network training set leads to a good agreement for the vapor coexistence densities. Liquid phase densities under supercooled conditions are also represented well, even though they were not included in the training set. These results confirm that neural network models can combine accuracy and transferability if sufficient attention is given to the construction of a representative training set for the target system.

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

Characterizing the Mechanisms of Ca and Mg Carbonate Ion-Pair Formation with Multi-Level Molecular Dynamics/Quantum Mechanics Simulations

The carbonate minerals of Ca and Mg are abundant throughout the lithosphere and have recently garnered significant research interest as possible long-term carbon sinks in the sequestration of atmospheric carbon dioxide. Nonetheless, an understanding of the atomic-level processes comprising their mineralization remains limited. Furthermore, we characterize and contrast the mechanisms of contact ion-pair formation in aqueous Ca and Mg carbonate systems, which represents the most fundamental step leading to the formation of their mineral solids. Utilizing multilevel embedded correlated wavefunction-based ab initio molecular dynamics/quantum mechanics simulations, we characterize not only the dynamics of these processes but also factors arising from the electronic structure of the involved species, revealing further details of the fundamentally different mechanisms for the interconversion between the contact ion-pairs and solvent-shared ion-pairs of Ca versus Mg carbonate.

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

Effect of K + Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol

Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion–ion and ion–solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K + ion. Diffusion-based Nernst–Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.

25 ENERGY STORAGE↗
DOE OSTI2024

Molecular Dynamical and Quantum Mechanical Exploration of the Site-Specific Dynamics of Cy3 Dimers Internally Linked to dsDNA

Performing spectroscopic measurements on biomolecules labeled with fluorescent probes is a powerful approach to locating the molecular behavior and dynamics of large systems at specific sites within their local environments. The indocarbocyanine dye Cy3 has emerged as one of the most commonly used chromophores. The incorporation of Cy3 dimers into DNA enhances experimental resolution owing to the spectral characteristics influenced by the geometric orientation of excitonically coupled monomeric units. Various theoretical models and simulations have been utilized to aid in the interpretation of the experimental spectra. In this study, we employ all-atom molecular dynamics simulations to study the structural dynamics of Cy3 dimers internally linked to the dsDNA backbone. We used quantum mechanical calculations to derive insights from both the linear absorption spectra and the circular dichroism data. Furthermore, we explore potential limitations within a commonly used force field for cyanine dyes. The molecular dynamics simulations suggest the presence of four possible Cy3 dimeric populations. The spectral simulations on the four populations show one of them to agree better with the experimental signatures, suggesting it to be the dominant population. Furthermore, the relative orientation of Cy3 in this population compares very well with previous predictions from the Holstein–Frenkel Hamiltonian model.

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

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

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

Charged nanoparticle attraction in multivalent salt solution: A classical-fluids density functional theory and molecular dynamics study

Here, negatively charged nanoparticles (NPs) in 1:1, 1:2, and 1:3 electrolyte solutions are studied in a primitive ion model using molecular dynamics (MD) simulations and classical density functional theory (DFT). We determine the conditions for attractive interactions between the like-charged NPs. Ion density profiles and NP–NP interaction free energies are compared between the two methods and are found to be in qualitative agreement. The NP interaction free energy is purely repulsive for monovalent counterions, but can be attractive for divalent and trivalent counterions. Using DFT, the NP interaction free energy for different NP diameters and charges is calculated. The depth and location of the minimum in the interaction depend strongly on the NPs’ charge. For certain parameters, the depth of the attractive well can reach 8–10 k B T, indicating that kinetic arrest and aggregation of the NPs due to electrostatic interactions is possible. Rich behavior arises from the geometric constraints of counterion packing at the NP surface. Layering of counterions around the NPs is observed and, as secondary counterion layers form the minimum of the NP–NP interaction free energy shifts to larger separation, and the depth of the free energy minimum varies dramatically. We find that attractive interactions occur with and without NP overcharging.

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

Molecular Theory of Hydration at Different Temperatures

Solvation plays an important role in diverse chemical processes ranging from reaction kinetics to molecular recognition, solubility, and phase separations. Despite a long-history of theoretical exploration, quantitative prediction of solvation remains a theoretical challenge without relying on the macroscopic properties of the solvent as an input. Here we present a molecular density functional theory that provides a self-consistent description of the solvation structure and thermodynamic properties of small organic molecules in liquid water at different temperatures. Based on the solute configuration and force-field parameters generated from first-principles calculations, the theoretical predictions are found in good agreement with experimental data for the hydration free energies of 197 organic molecules in a temperature range from 0 to 40 °C. In addition to calibration with experimental results, the theoretical predictions are compared with recent molecular dynamics simulations for the hydration of five highly explosive nitrotoluenes. Lastly, this work demonstrates the potential of the classical density functional theory for high-throughput prediction of solvation properties over a broad range of temperatures.

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

Restrained Proton Indicator in Combined Quantum-Mechanics/Molecular-Mechanics Dynamics Simulations of Proton Transfer through a Carbon Nanotube

Recently, a collective variable “proton indicator” was purposed for tracking an excess proton solvated in bulk water in molecular dynamics simulations. In this work, we demonstrate the feasibility of utilizing the position of this proton indicator as a reaction coordinate to model an excess proton migrating through a hydrophobic carbon nanotube in combined quantum-mechanics/molecular-mechanics simulations. Our results indicate that applying a harmonic restraint to the proton indicator in the bulk solvent near the nanotube pore entrance leads to the recruitment of water molecules into the pore. This is consistent with an earlier study that employed a multistate empirical valence bond potential and a different representation (center of excess charge) of the proton. We attribute this water recruitment to the delocalized nature of the solvated proton, which prefers to be in high-dielectric bulk solvent. While water recruitment into the pore is considered an artifact in the present simulations (because of the artificially imposed restraint on the proton), if the proton were naturally restrained, it could assist in building water wires prior to proton transfer through the pore. The potential of mean force for a proton translocation through the water-filled pore was computed by umbrella sampling, where the bias potentials were applied to the proton indicator. The free energy curve and barrier heights agree reasonably with those in the literature. Finally, the results suggest that the proton indicator can be used as a reaction coordinate in simulations of proton transport in confined environments.

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

Water Lone Pair Delocalization in Classical and Quantum Descriptions of the Hydration of Model Ions

Understanding the nature of ionic hydration at a fundamental level has eluded scientists despite intense interest for nearly a century. In particular, the microscopic origins of the asymmetry of ion solvation thermodynamics with respect to the sign of the ionic charge remains a mystery. Here, we determine the response of accurate quantum mechanical water models to strong nanoscale solvation forces arising from excluded volumes and ionic electrostatic fields. This is compared to the predictions of two important limiting classes of classical models of water with fixed point changes, differing in their treatment of “lone pair” electrons. Using the quantum water model as our standard of accuracy, we find that a single fixed classical treatment of lone pair electrons cannot accurately describe solvation of both apolar and cationic solutes, emphasizing the need for a more flexible description of local electronic effects in solvation processes. However, we explicitly show that all water models studied respond to weak long-ranged electrostatic perturbations in a manner that follows macroscopic dielectric continuum models, as would be expected. Here, we emphasize the importance of these findings in the context of realistic ion models, using density functional theory and empirical models, and discuss the implications of our results for quantitatively accurate reduced descriptions of solvation in dielectric media.

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

A Hybrid Quantum–Classical Study of Ion Adsorption at the Copper Electrode

Here, we present a set of computational studies that surveys the molecular structure and thermodynamics of simple ions at the electrode–electrolyte interface. A hybrid quantum–classical approach is used to obtain free energy profiles and solvent structure as a function of the ion’s position relative to the electrochemical liquid/solid interface and system potential. Model systems consist of a copper electrode and single halide or alkali metal ion in contact with an implicit solvent and reveal new insights into the interplay of forces that drive or resist ion absorption at the electrode/ electrolyte interface. We also expand on subtleties of absorption site preference and describe the complex relationship between the energetic barrier to absorption and the ion–electrode binding energy as a function of electrode potential.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface Engineering of Copper Catalyst through CO* Adsorbate

The electrochemical reduction of CO 2 with Cu-based catalysts depends intimately on the instantaneous local chemical environment of the catalyst-electrolyte interface. This microenvironment fluctuates according to the concentration of surface-adsorbed competing reaction intermediates and the applied electrode potential. In practice, disentangling these factors is exceedingly challenging, yet they critically determine the electrocatalyst efficiency and selectivity. Using grand canonical quantum-classical hybrid calculations, we quantify the complex interdependence between electrode potential, CO* coverage, and the interfacial field strength. Here we show that the often overlooked CO* coverage effect in fact strongly influences the field strength, with a magnitude change exceeding 1V/Å at certain potentials; among other effects, this change should lower the CO* dimerization barrier that dictates selectivity toward multi-carbon products. Beyond showcasing the importance of surface coverage for CO2 reduction, our results highlight the power of surface additives to modulate interfacial fields toward tailored electrochemical pathways.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2015

Carbyne fiber synthesis within evaporating metallic liquid carbon

Carbyne (e.g., linear chains of sp-bonded carbon) has been the subject of intense research focus due to its presence in astrophysical bodies, as well as its potential for use as a nanoelectronic device and superhard material. In this work, we discuss the formation of carbyne fiber bundles over a nanosecond time scale in laser pulse melting studies, using a previously determined density functional tight binding model for carbon coupled with a new correction for the dispersion energy. We determine our dispersion energy model by optimizing a modified Lennard-Jones potential to an experimentally determined equation of state for graphite, yielding excellent results for the bulk modulus and density under ambient conditions. We then simulate previous experiments by heating graphite to high temperature, followed by expanding the ensuing liquid phase to low density. Our results indicate that the initial, hot liquid phase mainly consists of sp2-bonded carbon atoms, which form a system of sp-bonded strands bound together via dispersion interactions upon achieving low density and temperature. Lastly, the high computational efficiency of our approach allows for direct comparison with experiments that span a wide range of thermodynamic conditions and can help determine parameters for synthesis of carbon-based materials with potentially exotic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗