Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular potentials”

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 145 records · Page 8

Converting Second‐Order Saddle Points to Transition States: New Principles for the Design of 4π Photoswitches

Abstract Molecular photoswitches have demonstrated potential for storing solar energy at the molecular level, with power densities comparable to commercial batteries and hydroelectric energy storage. However, development of efficient photoswitches is hindered by limitations in cyclability and optical properties of existing materials. We here demonstrate that certain limitations in photoswitches based on electrocyclizations stem from the issue of controlling competition between Woodward‐Hoffmann allowed and forbidden pathways. Our approach moves beyond the traditional view of activation barriers and reveals that second‐order saddle points are crucial in dictating the competition between disrotatory and conrotatory pathways. These insights suggest new opportunities to manipulate the competition between these pathways through geometric constraints, fundamentally altering the connectivity of the potential energy surface. Our study also emphasizes the necessity of multi‐reference methods and the need to conduct higher‐dimensional explorations for competing pathways beyond photoswitch design.

Chemistry↗

AladynPi – Adaptive Neural Network Molecular Dynamics Simulation Code with Physically Informed Potential: Computational Materials Mini-Application

This report provides an overview and description of commands used in the Computational Materials mini-application, AladynPi. AladynPi is an extension of a previously released mini-application, Aladyn (https://github.com/nasa/aladyn; Yamakov, V.I., and Glaessgen, E.H., NASA/TM-2018-220104). Aladyn and AladynPi are basic molecular dynamics codes written in FORTRAN 2003, which are designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory method. An input for the ANN is a set of structure coefficients, characterizing the local atomic environment of each atom, for which the atomic energy is obtained in the ANN inference process. In Aladyn, the ANN gives directly the energy of interatomic interactions. In AladynPi, the ANN gives optimized parameters for a predefined empirical function, known as bond-order-potential (BOP). The parameterized BOP function is then used to calculate the energy. AladynPi code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing manycore architectures based on graphic processing units (GPUs). The effort is supported by the High Performance Computing incubator (HPCi) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗

Ab initio investigation of a hypersonic double cone experiment

This article presents a direct molecular simulation (DMS) of a reactive Mach 8.2 oxygen flow over a double cone geometry. The free stream conditions and article configuration generate a flow with thermal and chemical nonequilibrium, which are common attributes of hypersonic flight. This scenario was first studied experimentally at Calspan University of Buffalo Research Center’s test facility. DMS is a particle method that uses quantum mechanically derived interaction potentials to simulate molecular collisions within a flow field. Since these interaction potentials are the only modeling inputs used in the simulation, all flow features can solely be attributed to the ab initio potential energy surfaces. Hence, providing a comparison of a hypersonic ground test and numerical data anchored to quantum mechanics. The fundamental nature of DMS is leveraged to investigate molecular level mechanisms prevalent in the flow, and comparisons with lower fidelity simulations are presented to highlight the role of these first principles calculations as benchmark solutions.

Science & Technology - Other Topics↗

Star formation in molecular clouds

A general overview is presented of the findings of molecular line and infrared astronomy regarding star formation. Beginning with the statement that stars, and therefore, planets, form from molecular clouds, it is argued that the growing body of data on conditions in molecular clouds should serve as input parameters for theoretical calculations of star and planet formation. The fact that the gravitational potential energy of molecular clouds far exceeds their thermal energy implies that the required rate of star formation is easily explained. Indeed the most important problem in this field is that the rate of star formation predicted from molecular cloud studies is too high. Until this problem is understood, additional modes of inducing star formation, while they may be occurring, only add to the existing difficulties.

Evans, N. J., II↗

Mechanical Properties of Nanostructured Materials Determined Through Molecular Modeling Techniques

The potential for gains in material properties over conventional materials has motivated an effort to develop novel nanostructured materials for aerospace applications. These novel materials typically consist of a polymer matrix reinforced with particles on the nanometer length scale. In this study, molecular modeling is used to construct fully atomistic models of a carbon nanotube embedded in an epoxy polymer matrix. Functionalization of the nanotube which consists of the introduction of direct chemical bonding between the polymer matrix and the nanotube, hence providing a load transfer mechanism, is systematically varied. The relative effectiveness of functionalization in a nanostructured material may depend on a variety of factors related to the details of the chemical bonding and the polymer structure at the nanotube-polymer interface. The objective of this modeling is to determine what influence the details of functionalization of the carbon nanotube with the polymer matrix has on the resulting mechanical properties. By considering a range of degree of functionalization, the structure-property relationships of these materials is examined and mechanical properties of these models are calculated using standard techniques.

Clancy, Thomas C.↗

Sabatier Catalyst Poisoning Investigation

The Carbon Dioxide Reduction Assembly (CRA) on the International Space Station (ISS) has been operational since 2010. The CRA uses a Sabatier reactor to produce water and methane by reaction of the metabolic CO2 scrubbed from the cabin air and the hydrogen byproduct from the water electrolysis system used for metabolic oxygen generation. Incorporating the CRA into the overall air revitalization system has facilitated life support system loop closure on the ISS reducing resupply logistics and thereby enhancing longer term missions. The CRA utilizes CO2 which has been adsorbed in a 5A molecular sieve within the Carbon Dioxide Removal Assembly, CDRA. There is a potential of compounds with molecular dimensions similar to, or less than CO2 to also be adsorbed. In this fashion trace contaminants may be concentrated within the CDRA and subsequently desorbed with the CO2 to the CRA. Currently, there is no provision to remove contaminants prior to entering the Sabatier catalyst bed. The risk associated with this is potential catalyst degradation due to trace organic contaminants in the CRA carbon dioxide feed acting as catalyst poisons. To better understand this risk, United Technologies Aerospace System (UTAS) has teamed with MSFC to investigate the impact of various trace contaminants on the CRA catalyst performance at relative ISS cabin air concentrations and at about 200/400 times of ISS concentrations, representative of the potential concentrating effect of the CDRA molecular sieve. This paper summarizes our initial assessment results.

Nallette, Tim↗

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Draco cloud as a matrix of potential protostars

The Draco molecular cloud is inbound at a z distance of a few hundred parsecs, according to previous observations. IRAS Point Source Catalog data for a 36 sq deg area that includes most of the cloud and a comparison area around it are tabulated, illustrated, and discussed. The sources are classified as 56 unidentified, 38 identified extragalactic objects, and 44 identified stars. The unidentified and unresolved 100 micron sources appear to be globule-like clumps that are candidates for producing protostars, although the data cannot demonstrate currently accreting protostars in the globules. An undetermined percent of 100 micron sources may be artifactual, but only a small percent of them may be galaxies because of the adherence to the Draco cloud. Sources detected also at 60 micron or only at 60 micron do not adhere to the cloud and are probably galaxies. Unless artifactual the unidentified and unresolved 100 micron sources may provide considerable new information about the cloud components and therefore about potential star formation in a new class of sites and origins. Six unidentified sources have 12-25 micron stellar spectral distributions, but neither they nor the identified stars have 12-25 micron spectra characteristic of T Tauri stars.

Johnson, Hugh M.↗

Alterations of p53 in tumorigenic human bronchial epithelial cells correlate with metastatic potential

The cellular and molecular mechanisms of radiation-induced lung cancer are not known. In the present study, alterations of p53 in tumorigenic human papillomavirus-immortalized human bronchial epithelial (BEP2D) cells induced by a single low dose of either alpha-particles or 1 GeV/nucleon (56)Fe were analyzed by PCR-single-stranded conformation polymorphism (SSCP) coupled with sequencing analysis and immunoprecipitation assay. A total of nine primary and four secondary tumor cell lines, three of which were metastatic, together with the parental BEP2D and primary human bronchial epithelial (NHBE) cells were studied. The immunoprecipitation assay showed overexpression of mutant p53 proteins in all the tumor lines but not in NHBE and BEP2D cells. PCR-SSCP and sequencing analysis found band shifts and gene mutations in all four of the secondary tumors. A G-->T transversion in codon 139 in exon 5 that replaced Lys with Asn was detected in two tumor lines. One mutation each, involving a G-->T transversion in codon 215 in exon 6 (Ser-->lle) and a G-->A transition in codon 373 in exon 8 (Arg-->His), was identified in the remaining two secondary tumors. These results suggest that p53 alterations correlate with tumorigenesis in the BEP2D cell model and that mutations in the p53 gene may be indicative of metastatic potential.

Non-NASA Center↗

Development of a deep potential model for F and CF 2 etching of Si and SiO 2

An understanding of plasma-surface interactions at increasingly smaller scales is invaluable for the development of novel technologies and processing techniques. Molecular dynamics (MD) simulations can provide insights into atomic-scale interactions, though they are restricted by the availability of interatomic potentials. Machine learning methods, such as Deep Potential Molecular Dynamics (DeepMD), provide a systematic framework for the development of accurate and flexible ab initio-based models. In this work, we develop DeepMD models for the ion-enhanced etching of Si and SiO 2 by F and CF 2 radicals. We employ an active learning process to expand the data set on which the model is trained and demonstrate its effect on the model accuracy. The DeepMD results are compared to data from classical MD simulations and experiments. Physical sputtering yields of SiO 2 by Ar + ions show good agreement with previous simulation results using conventional interatomic potentials, though the predicted depth profiles are different. Etching yields are calculated as a function of ion energy and neutral to ion flux ratio for the Ar + ion-enhanced etching of SiO 2 and Si by F atoms, as well as for etching of SiO 2 by CF 2 radicals, showing reasonable agreement with experimental data. Finally, an ion-enhanced surface kinetic model is fitted to the DeepMD etch yields, and the fitted parameters are compared to quantities computed directly from DeepMD simulations. This study illustrates how molecular dynamics simulations using machine learning potentials can provide an accurate model of etching processes relevant to device manufacturing.

Kounis-Melas, Andreas [Princeton Univ., NJ (United↗

A cavity-enhanced spectroscopist's lens on molecular polaritons

Polariton chemistry has been hailed as a potential new route to direct molecular processes with electromagnetic fields. To make further strides, the community needs to clarify which unusual polaritonic phenomena are true hallmarks of cavity quantum electrodynamics and which can be rationalized with classical optical physics. Here, in this work, we provide a tutorial perspective on the formation, spectroscopy, and behavior of molecular polaritons using classical optics. Where possible, we draw connections to cavity-enhanced spectroscopy and recast open questions in terms that may be more familiar to the broader community of physical chemists.

Cavity quantum electrodynamics↗

Molecular processes in comets

Potential energy curves for the two lowest 2 sigma- states of OH are computed at the configuration-interaction level using four different basis sets. Electronic transition dipole moments connecting the excited 1 2 sigma- and 2(D) 2 sigma- states with each other and with he ground X2pi state are presented as functions of internuclear distance. The theoretical absorption oscillator strengths for the D 2 sigma-(v prime=0) reverser to X 2pi (v prime prime=0) transition are in good agreement with the empirical value derived from astronomical measurements. The photodissociation cross sections for absorption rom the v prime prime=0,1, and levels of the ground state into the continuum of the 1 2 sigma- state are calculated, and the interstellar and cometary photodissociation rates are derived.

Dalgarno, A.↗

Transition probabilities of the B-prime 1Sigma(u)(+) to X 1Sigma(g)(+) system of molecular hydrogen

From published potential energy curves and transition dipole moments, there are obtained by numerical integration of the equations of nuclear motion the vibrational eigenfunctions of the X 1Sigma(g)(+) and B-prime 1Sigma(u)(+) states of H2. The probabilities of radiative transitions from the discrete vibrational levels of the excited B-prime 1Sigma(u)(+) electronic state of H2 to the discrete and continuum vibrational levels of the ground X 1Sigma(g)(+) electronic state are calculated. The Franck-Condon factors are also presented.

Kwok, T. L.↗

Highly Conjugated Porphyrin Arrays Enable Optical Resolution of Ferromagnetic and Antiferromagnetic Aligned States of the Triplet Exciton and an Incorporated Stable Radical

Well-defined photogenerated molecular spin systems have potential utility in spintronics and quantum information science (QIS). Because molecular magnetic, optical, and electronic properties can be controlled by design, diverse spin systems can be prepared at modest temperatures. Photogenerated molecular spin systems often involve states prepared from the interaction of excitons and charges. Resolving the nature of electron spin alignment in photogenerated spin states described by the coupling of a triplet exciton and a stable radical commonly relies on EPR spectroscopy. Here, we describe ethyne-bridged (porphinato)metal ( PM n ) oligomers that incorporate a macrocycle-bound Cu(II) radical center. Upon photoexcitation of such PM n arrays, a singdoublet ( 2 S 1 ) state is formed; ultrafast internal conversion (IC) then produces a tripdoublet ( 2 T 1 ) state, which undergoes intersystem crossing (ISC) to produce a tripquartet ( 4 T 1 ) state, before relaxation to the ground state ( 2 S 0 ). These highly conjugated Cu(II) radical-containing PMn arrays enable direct observation of copper porphyrin 2 T 1 → 4 T 1 ISC dynamics from the biexponential decay of the near-infrared (NIR) 2,4 T 1 → 2,4 T n transient absorption manifold. Multireference n-electron valence perturbation theory (NEVPT2) computations illuminate how PM n electronic structure controls the relaxation dynamics of these long-lived (>10 ns) electronically excited multiplet states. These studies show that highly conjugated and polarizable porphyrin arrays incorporating stable spin centers provide rare π-delocalized systems where the ferromagnetic and antiferromagnetic alignment between a triplet exciton and a stable radical are both spectrally resolved and addressable using transient optical spectroscopy at wavelengths exceeding 1 μm, providing new opportunities to QIS.

14 SOLAR ENERGY↗

Temperature-Dependent Molecular Diffusional Properties in Deep Eutectic Solvents and Eutectogels

Eutectogels (ETGs) prepared from deep eutectic solvents (DESs), a gelator, and water have many uses in separations, catalysis, and energy storage systems. In these applications, temperature-dependent molecular diffusional properties and intermolecular interactions play a critical role in their function. Diffusional properties of Alexa Fluor 633 and ATTO 647N were measured across a range of temperatures in choline chloride:2glycerol DESs comprised of one molar equivalent of choline chloride and two molar equivalents of glycerol (also known as glyceline) as well as ETGs made from this DES, a xanthan gum gelator, and 10% w/w water or 20% w/w water. Fluorescence recovery after photobleaching (FRAP) was employed to evaluate potential changes in molecular diffusion within DESs and ETGs from 20 to 100 °C. Surprisingly, the ETGs have a larger sample viscosity but also exhibit faster molecular diffusion when compared to the dry DES. This is attributed to macroscopic properties of the ETGs (e.g., pores, three-dimensional structure). The FRAP data also show an irreversible temperature-dependent decrease in the diffusion coefficient of both fluorophores in the ETGs. This is consistent with differential scanning calorimetry data of the ETG, which shows a shift to a more negative glass transition temperature (−36 °C to −47 °C) after the first cooling/heating cycle. Raman data reveal no detectable changes in the intermolecular interactions in either the DESs or ETGs as a function of temperature or water content. Finally, the findings on temperature- and water-dependent diffusional properties of DESs and ETGs provide a foundation for optimizing the use of these materials across various applications, particularly where repeated heating/cooling cycles may be used.

alcohols↗

Hydrogen and related materials at high density: Physics, chemistry and planetary implications

Recent studies of low-Z molecular materials including hydrogen to multimegabar pressures (less than 300 GPa) have uncovered a range of phenomena relevant to understanding the nature of the interiors of the outer planets and their satellites. Synchrotron x ray diffraction measurements (to 42 GPa) have been used to determine the crystal structure of the solid (hexagonal-close packed) and equation of state. Sound velocities in fluid and solid hydrogen (to 24 GPa) have been inverted to obtain elastic constants and aggregate bulk and shear moduli. In addition, an improved intermolecular potential has been determined which fits both static and shock-wave data. Use of the new potential for the molecular envelope of Jupiter suggests the need for major revisions of existing Jovian models or a reanalysis of reported free oscillations for the planet. Studies at higher pressures (greater than 100 GPa) reveal a sequence of pressure-induced symmetry-breaking transitions in molecular hydrogen, giving rise to three high-pressure phases (1, 2, and 3). Phase 1 is the rotationally disordered hcp phase which persists from low pressure to well above 100 GPa at high temperature (e.g., 300 K). Phase 2 is a low-temperature, high-pressure phase (transition at 100 GPa and 77 K in H2) with spectral features indicative of partial rotational ordering and crystallographic distortion. The transition to Phase 3 at 150 GPa is accompanied by a weakening of the molecular bond, gradual changes in orientational ordering, strong enhancement of the infrared intramolecular vibrational absorption, and strong intermolecular interactions similar to those of ambient-pressure network solids. Studies of the phase diagram reveal a triple point near 130 K and 160 GPa. Higher pressure measurements of vibrational spectra place a lower bound of approximately 250 GPa on the predicted transition pressure for dissociation of molecular hydrogen to form a monatomic metal.

Hemley, R. J.↗