Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular properties”

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 325 records · Page 18

Optical and transport properties of plasma mixtures from ab initio molecular dynamics

Predicting the charged particle transport properties of warm dense matter/hot dense plasma mixtures is a challenge for analytical models. High accuracy ab initio methods are more computationally expensive, but can provide critical insight by explicitly simulating mixtures. In this work, we investigate the transport properties and optical response of warm dense carbon–hydrogen mixtures at varying concentrations under either conserved electronic pressure or mass density at a constant temperature. We compare options for mixing the calculated pure species properties to estimate the results of the mixtures. We find that a combination of the Drude model with the Matthiessen's rule works well for DC electron transport and low-frequency optical response. This breaks down at higher frequencies, where a volumetric mix of pure-species AC conductivities works better.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atomistic simulation of brittle-to-ductile transition in silicon carbide embedded with nano-sized helium bubbles

In this study, the tensile response of cubic silicon carbide (SiC) bulk containing cavities (voids and He bubbles) has been investigated using molecular dynamic simulations. The formation of cavities in SiC leads to a significant degradation in the mechanical properties of SiC with more influence on material fracture than initial elastic deformation. The brittle-to-ductile transition occurs in cavity-embedded SiC as the pressure in He bubbles increases. This is associated with the deformation mechanism that bond breaking at a low He bubble pressure transfers to extensive dislocation activities at a higher He bubble pressure. The cavities can effectively concentrate stress around them in the direction perpendicular to the tension, which leads to preferred cracking in the region with a higher tensile stress. The failure mechanism as revealed by this study improves understanding of property degradation in SiC that may be useful for applications of SiC in advanced nuclear energy systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling the Effect of Hydrogen on Tungsten Grain Boundary Properties

The objective of this study is to assess the effect of hydrogen isotopes on mechanical properties of tungsten that are relevant to first-wall applications in a tokamak. We focus on the effects of hydrogen segregation to grain boundaries and its impact on the grain boundary structure and properties. Through molecular dynamics and other atomistic techniques, we are able to make predictions of grain boundary structure and properties related to it, including grain boundary mobility and impurity diffusivity along the grain boundary. The grain boundary structure can change significantly and abruptly as conditions change, including due to heating and radiation damage. These grain boundary phase transformations can have a pronounced effect on mobility and impurity diffusivity, even when no hydrogen is present. How this behavior changes in the presence of hydrogen is not known. The phase may affect the level of hydrogen segregation and hydrogen storage, and the hydrogen may induce significant grain boundary structure changes. Our objective is to predict whether they do. The results inform thermomechanical modeling of tungsten at tokamak conditions.

36 MATERIALS SCIENCE↗

Associative polymers with controlled sticker placement: How reversible bond distribution and density govern polymer dynamics

Associative polymers with precisely arranged stickers offer opportunities to program material properties with molecular precision. Yet, it remains unclear how the placement and fraction of stickers dictate structure, dynamics, and macroscopic properties. By developing a model unentangled polymer system with hydrogen-bonding stickers, we show that randomly distributed stickers neither form clusters nor change flow properties, whereas stickers placed at chain ends drive nanocluster formation even at low concentrations. Adding more end stickers produces a rubbery plateau spanning eight decades in frequency with two distinct relaxation timescales, in contrast to the single plateau predicted by the classic sticky Rouse model. These results demonstrate that sticker distribution dictates whether associative polymers undergo nanocluster formation or microphase separation, while substantial alterations in dynamics and viscoelasticity require both sticker aggregation and thermomechanical stability of associated domains. Our findings resolve a longstanding debate on associative polymer dynamics and provide molecular design rules for programmable soft materials.

36 MATERIALS SCIENCE↗

Molecular Dynamics Modeling of Carbon Nanotube Composite Fracture Using ReaxFF

Carbon nanotube (CNT) fiber reinforced composites with specific tensile strengths and moduli approaching those of aerospace grade carbon fiber composites have recently been reported. This achievement was enabled by the emerging availability of high N/tex yarns in kilometer-scale quantities. While the production of this yarn is an impressive advance, its strength is still much lower than that of the individual CNTs comprising the yarn. Closing this gap requires understanding load transfer between CNTs at the nanometer dimensional scale. This work uses reactive molecular dynamics simulations to gain an understanding at the nanometer scale of the key factors that determine CNT nanocomposite mechanical performance, and to place more realistic upper bounds on the target properties. While molecular dynamics simulations using conventional force fields can predict elastic properties, the ReaxFF reactive forcefield can also model fracture behavior because of its ability to accurately describe bond breaking and formation during a simulation. The upper and lower bounds of CNT composite properties are investigated by comparing systems composed of CNTs continuously connected across the periodic boundary with systems composed of finite length CNTs. These lengths, effectively infinite for the continuous tubes and an aspect ratio of 13 for the finite length case, result from practical limitations on the number of atoms that can be included in a simulation. Experimentally measured aspect ratios are typically on the order of 100,000, so the calculated results should represent upper and lower limits on experimental mechanical properties. Finally, the effect of various degrees of covalent crosslinking between the CNTs and amorphous carbon matrix is considered to identify the amount of CNT-matrix covalent bonding that maximizes overall composite properties.

Jensen, Benjamin D.↗

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular-scale origins of solution nanostructure and excess thermodynamic properties in a water/amphiphile mixture

The molecular and nanoscale origins of nonideality in excess thermodynamic properties are essential to understanding cosolvent mixtures, yet they remain challenging to determine. Here, we consider a binary mixture of water and an amphiphile, N,N,N',N'-tetramethylmalonamide (TMMA), which is characterized by strong hydrogen bonding between the two components and no hydrogen bonding between amphiphiles. Using molecular dynamics simulation, validated with excess volume measurements and X-ray scattering, we identify three distinct solution regimes across the composition range of the binary mixture and find that the transition between two of these regimes, marked by the water percolation threshold, is closely correlated with minima in the excess volume and excess enthalpy. Structural analysis of the simulations reveals an interplay between local interactions and solution nanostructure, determined by the relative strength of the water–water and water–amphiphile hydrogen bonding interactions. By comparison with other amphiphiles, such as linear alcohols, the relative strength of like and unlike interactions between water and amphiphile affects the relationship between thermodynamics and structural regimes. This provides insight into how molecular forces of mutual solvation interact across length scales and how they manifest in excess thermodynamic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Applying Mechanistic Co-Design Towards Room-Temperature Organic, Optical, and Multiple-Center Molecular Qubits

The next generation of qubits, the building blocks of quantum technologies, are designed using molecules. Molecular qubits have natural advantages over solid-state equivalents in biological and chemical sensing. Chemistry allows the properties of molecular qubits to be fine-tuned, such as the ability of multiple qubits to work together using qubit-qubit coupling. Molecules and structures containing multiple qubits, made from both inorganic and organic elements, are developed in this work to bring quantum technologies closer to room-temperature operation. The goal is achieved by understanding how the core of the qubit, known as spin, interacts with the vibrating parts of the molecule, other spins, and the surrounding environment. A combined approach involving theory, spectroscopy, and synthesis is used to achieve the design goal.

36 MATERIALS SCIENCE↗

Charge Transport Characteristics of Molecular Electronic Junctions Studied by Transition Voltage Spectroscopy

The field of molecular electronics is prompted by tremendous opportunities for using a single-molecule and molecular monolayers as active components in integrated circuits. Until now, a wide range of molecular devices exhibiting characteristic functions, such as diodes, transistors, switches, and memory, have been demonstrated. However, a full understanding of the crucial factors that affect charge transport through molecular electronic junctions should yet be accomplished. Remarkably, recent advances in transition voltage spectroscopy (TVS) elucidate that it can provide key quantities for probing the transport characteristics of the junctions, including, for example, the position of the frontier molecular orbital energy relative to the electrode Fermi level and the strength of the molecule–electrode interactions. These parameters are known to be highly associated with charge transport behaviors in molecular systems and can then be used in the design of molecule-based devices with rationally tuned electronic properties. This article highlights the fundamental principle of TVS and then demonstrates its major applications to study the charge transport properties of molecular electronic junctions.

36 MATERIALS SCIENCE↗

Gas permeation properties of amorphous zeolitic imidazolate framework membranes made by atomic/molecular layer deposition

Amorphous metal-organic framework (MOF) membranes are desirable because they may retain some of the molecular sieving properties of their crystalline counterparts while being free of grain boundary defects, which often hinder the consistent achievement of high membrane performance. However, current methods, like melting and compression, for fabricating amorphous MOF membranes involve multi-step processes that require the formation of a crystalline membrane first, that is then amorphized, and therefore, could be challenging to scale. Here, we utilize atomic/molecular layer deposition (ALD/MLD) of diethylzinc (DEZ) and 2-methylimidazole (2mIm) to directly synthesize ultrathin amorphous zeolitic imidazolate framework (aZIF) deposits on γ-alumina-coated α-alumina supports. As the number of ALD/MLD cycles increased from 10 to 300, gas permeances decreased while ideal selectivities increased. Mixture separation factors for C3H6/C3H8, CO2/N2, and H2/C3H8 as high as 4, 37, and 194, respectively, were obtained. At 200 °C and 2.5 bar equimolar feed of H2 and C3H8, a H2/C3H8 mixture separation factor of 185 is obtained with an H2 permeance of ca. 4.15x10-8 mol/m2-s-Pa (124 GPU). Additionally, the membrane achieves a CO2/N2 mixture separation factor of 37 at 25 °C and 1 bar with a CO2 permeance of ca. 3.96x10-8 mol/m2-s-Pa (118 GPU). Considering the vast array of compositionally distinct aZIFs that can be potentially deposited by this approach, an enormously large parameter space for membrane design is emerging to be explored.

36 MATERIALS SCIENCE↗

Organic Multiferroics (Final Technical Report)

The research objective of this project is to explore the unusual physical properties of molecular electronic crystal materials, comprised of two chemically and electronically distinct donor and acceptor molecules, which combine two or more properties in the same crystal lattice in a manner that is difficult or impossible to achieve in continuous inorganic crystalline solids.

36 MATERIALS SCIENCE↗

Critical role of structural order in bipolar redox-active molecules for organic redox flow batteries

Bipolar redox-active molecules (BRMs) have been suggested as a means to address crossover related issues in all-organic redox flow batteries (RFBs). In such species, electron donors (anolytes) and electron acceptors (catholytes) are linked by a chain, and the same chemical composition is used in both anode and cathode compartments to reduce chemical gradients and unbalanced ion transport. The resultant RFBs resemble the aqueous vanadium RFBs, but offer greater design flexibility and potentially more favorable electrochemical and physicochemical properties. Yet the complex trade-offs in these properties have been a complication. Here we developed a fundamental rubric to uncover the likely origins in the performance metrics of these BRMs that are salient for their use in RFBs. Methylene linked phenothiazine and phthalimide moieties were employed as model BRMs and variations in their properties and molecular conformations were evaluated by systematically varying the methylene linker. Our results revealed that even minimalistic changes in the linker length resulted in dramatic oscillations in the solubility and stability. Using crystallography, quantum chemistry, and molecular dynamics, we clarify that the seemingly fluctuating behavior is due to (1) the inter- and intra-molecular charge transfer between the donor and acceptor through bond or/and through space; (2) the formation of distinctive packing/clustering motifs. Both interactions are strongly dependent on the molecular conformation. As these structural factors modify the electronic structure of both electroactive functional groups and do not change monotonically, the structure–property response in BRMs is often complex. Finally, accordingly, these behaviors need to be taken into account in developing BRMs and knowledge of these structural factors will allow the more rational design of BRMs for grid-scale energy storage.

25 ENERGY STORAGE↗

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leveraging Sequential Doping of Semiconducting Polymers to Enable Functionally Graded Materials for Organic Thermoelectrics

With the ability to modulate electronic properties through molecular doping coupled with ease in processability, semiconducting polymers are at the forefront in enabling organic thermoelectric devices for thermal energy management. In contrast to uniform thermoelectric material properties, an alternative route focuses on functionally graded materials (FGMs) where one spatially controls and optimizes transport properties across the length of a thermoelectric material. While primarily studied in the context of inorganic materials, the concept of FGMs for organic thermoelectrics has not been explored. Herein, we introduce how molecular doping of semiconducting polymers enables spatial compositional control of thin-film FGMs. Specifically, we utilize sequential vapor doping of poly[2,5-bis(3-tetradecylthiophen-2-yl) thieno [3,2-b]thiophene] (PBTTT) with the small molecule acceptor 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ) to fabricate the simplest form of FGMs – double- segmented thin films. The two thin-film segments are of equal length (7.5 mm), but each set to different doping levels. Our study focuses on understanding the thermoelectric properties (Seebeck coefficient, α, and electronic conductivity, σ) and structural properties (through X-ray scattering, UV-vis-NIR spectroscopy, and Raman spectroscopy) within and across the two segments. We observe the presence of a small diffuse interfacial region of 0.5 – 1 mm between the two segments where the doping level and transport properties vary continuously. Despite the diffuse interface, the measured effective Seebeck coefficient (αeff) across the two segments is simply the average of α within each segment. Importantly, this experimental result is consistent with reported mathematical models describing αeff in graded thermoelectric materials. Finally, our results demonstrate the facile fabrication and characterization of functionally graded organic thermoelectric materials, providing guidelines for further development on more complex FGMs.

36 MATERIALS SCIENCE↗

High-throughput screening of tribological properties of monolayer films using molecular dynamics and machine learning

Monolayer films have shown promise as a lubricating layer to reduce friction and wear of mechanical devices with separations on the nanoscale. These films have a vast design space with many tunable properties that can affect their tribological effectiveness. For example, terminal group chemistry, film composition, and backbone chemistry can all lead to films with significantly different tribological properties. This design space, however, is very difficult to explore without a combinatorial approach and an automatable, reproducible, and extensible workflow to screen for promising candidate films. Here, using the Molecular Simulation Design Framework (MoSDeF), a combinatorial screening study was performed to explore 9747 unique monolayer films (116 964 total simulations) and a machine learning (ML) model using a random forest regressor, an ensemble learning technique, to explore the role of terminal group chemistry and its effect on tribological effectiveness. The most promising films were found to contain small terminal groups such as cyano and ethylene. The ML model was subsequently applied to screen terminal group candidates identified from the ChEMBL small molecule library. Approximately 193 131 unique film candidates were screened with approximately a five order of magnitude speed-up in analysis compared to simulation alone. The ML model was thus able to be used as a predictive tool to greatly speed up the initial screening of promising candidate films for future simulation studies, suggesting that computational screening in combination with ML can greatly increase the throughput in combinatorial approaches to generate in silico data and then train ML models in a controlled, self-consistent fashion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗

Understanding the electron-water interaction at the molecular level: Integrating theory and experiment in the cluster regime

In this final technical report, we summarize our major accomplishments resulting from a combined experimental-theoretical approach that exploits the unique properties of molecular clusters to elucidate how local interactions at the molecular level drive the macroscopic behavior of materials and chemical transformations in liquids. The specific targets for our efforts have need to elucidate how charged entities are accommodated by hydrogen-bonding networks and to establish how the behavior of room temperature ionic liquids can be tailored for specific chemical applications. Because of their central importance in aqueous electrolyte chemistry, we have placed particular emphasis understanding the accommodation of excess electrons and protons by water. These surprisingly complex fundamental electrical charges play a fundamental role in a wide range of chemical and biological processes, and yet even a qualitative picture of their speciation in water and at the air-water interface has proven elusive. The main experimental diagnostic in our approach is to measure and analyze vibrational spectra of temperature and composition controlled molecular aggregates that isolate and amplify the spectral markers for local interactions that are often obscured in measurement on bulk or interfacial water. Our combined experimental-theoretical program has led to a greatly enhanced understanding of the nature of excess electrons and protons in water and in the process developed new and powerful methods of chemical analysis that are currently in use across the chemical sciences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A 3D visualization system for molecular structures

The properties of molecules derive in part from their structures. Because of the importance of understanding molecular structures various methodologies, ranging from first principles to empirical technique, were developed for computing the structure of molecules. For large molecules such as polymer model compounds, the structural information is difficult to comprehend by examining tabulated data. Therefore, a molecular graphics display system, called MOLDS, was developed to help interpret the data. MOLDS is a menu-driven program developed to run on the LADC SNS computer systems. This program can read a data file generated by the modeling programs or data can be entered using the keyboard. MOLDS has the following capabilities: draws the 3-D representation of a molecule using stick, ball and ball, or space filled model from Cartesian coordinates, draws different perspective views of the molecule; rotates the molecule on the X, Y, Z axis or about some arbitrary line in space, zooms in on a small area of the molecule in order to obtain a better view of a specific region; and makes hard copy representation of molecules on a graphic printer. In addition, MOLDS can be easily updated and readily adapted to run on most computer systems.

Green, Terry J.↗