Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

phosaa14SB and phosaa19SB: Updated Amber Force Field Parameters for Phosphorylated Amino Acids

Phosphorylated amino acids are involved in many cell regulatory networks; proteins containing these post-translational modifications are widely studied both experimentally and computationally. Simulations are used to investigate a wide range of structural and dynamic properties of biomolecules, such as ligand binding, enzyme-reaction mechanisms, and protein folding. However, the development of force field parameters for the simulation of proteins containing phosphorylated amino acids using the Amber program has not kept pace with the development of parameters for standard amino acids, and it is challenging to model these modified amino acids with accuracy comparable to proteins containing only standard amino acids. In particular, the popular ff14SB and ff19SB models do not contain parameters for phosphorylated amino acids. Here, the dihedral parameters for the side chains of the most common phosphorylated amino acids are trained against reference data from QM calculations adopting the ff14SB approach, followed by validation against experimental data. Finally, library files and corresponding parameter files are provided, with versions that are compatible with both ff14SB and ff19SB.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Car–Parrinello Monitor for More Robust Born–Oppenheimer Molecular Dynamics

Born–Oppenheimer molecular dynamics (BOMD) is a promising simulation method for exploring the possible reaction pathways of a chemical system, but one significant challenge is the increased difficulty of converging the self-consistent field (SCF) calculation that often accompanies the breaking and forming of chemical bonds. To address this challenge, we developed an enhancement to the BOMD simulation method called the Car–Parrinello monitor (CPMonitor) that uses Car–Parrinello molecular dynamics (CPMD) to recover from SCF convergence failures. CPMonitor works by detecting SCF convergence failures in BOMD and switching to a CPMD Hamiltonian to propagate through the region of configuration space where the SCF calculation is unable to converge, then switching back to BOMD when good convergence behavior is re-established. We present a series of simulation studies that use CPMonitor, including detailed studies of the thermodynamic and dynamical properties of simple systems, as well as ab initio nanoreactor simulations containing transition metal atoms that were previously not possible to simulate using standard BOMD methods. Our studies show that CPMonitor can make BOMD simulations robust to SCF convergence difficulties and improve simulation performance and stability in reaction discovery applications.

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

Quantifying Species Populations in Multivalent Borohydride Electrolytes

Multivalent batteries represent an important beyond Li-ion energy storage concept. The prospect of calcium batteries, in particular, has emerged recently due to novel electrolyte demonstrations, especially that of a ground-breaking combination of the borohydride salt Ca(BH 4 ) 2 dissolved in tetrahydrofuran. Recent analysis of magnesium and calcium versions of this electrolyte led to the identification of divergent speciation pathways for Mg 2+ and Ca 2+ despite identical anions and solvents, owing to differences in cation size and attendant flexibility of coordination. To test these proposed speciation equilibria and develop a more quantitative understanding thereof, we have applied pulsed-field-gradient nuclear magnetic resonance and dielectric relaxation spectroscopy to study these electrolytes. Additionally, concentration-dependent variation in anion diffusivities and solution dipole relaxations, interpreted with the aid of molecular dynamics simulations, confirms these divergent Mg 2+ and Ca 2+ speciation pathways. These results provide a more quantitative description of the electroactive species populations. We find that these species are present in relatively small quantities, even in the highly active Ca(BH 4 ) 2 /tetrahydrofuran electrolyte. This finding helps interpret previous characterizations of metal deposition efficiency and morphology control and thus provides important fundamental insight into the dynamic properties of multivalent electrolytes for next-generation batteries.

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

Accurate Force Field for Carbon Dioxide–Silica Interactions Based on Density Functional Theory

Fluid–silica interfaces are ubiquitous in chemistry, occurring in both natural geochemical environments and practical applications ranging from separations to catalysis. Simulations of these interfaces have been, and continue to be, a significant avenue for understanding their behavior. A constraining factor, however, is the availability of accurate force fields. Most simulations use traditional “mixing rules” to determine nonbonded dispersion interactions, an approach that has not been critically examined. Here, in this study, we present Lennard-Jones parameters for the interaction of carbon dioxide with silica interfaces that are optimized to reproduce density functional theory (DFT)-based binding energies. The modeling is based on the recently developed silica-DDEC force field, whose atomic charges are consistent with DFT calculations. Standard mixing rules are found to predict weaker CO 2 binding to silica than that obtained from DFT, an effect corrected by the optimized parameters given here. This behavior extends to other silica force fields (Clayff and Gulmen-Thompson), and the present Lennard-Jones parameters improve their performance as well. The effects of improved Lennard-Jones parameters on the structural and dynamical properties of condensed CO 2 in silica slit pores are also examined.

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

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗
DOE OSTI2018

The Ideal Crystal Structure of Cristobalite X-I: A Bridge in SiO 2 Densification

On compression of α-cristobalite SiO 2 to pressures above approximately 12 GPa, a new polymorph known as cristobalite X-I forms. The existence of cristobalite X-I has been known for several decades; however, consensus regarding its exact atomic arrangement has not yet been reached. The X-I phase constitutes an important step in the silica densification process, separating low-density tetrahedral framework phases from high-density octahedral polymorphs. It is the only nonquenchable high-density SiO 2 phase, which reverts to the low-density form on decompression at ambient temperature. Recently, an experimental study proposed an octahedral model of SiO 2 X-I with intrinsic structural defects involving partial Si site occupancies. In contrast, our new single-crystal synchrotron X-ray diffraction experiments have shown that the ideal structure of this phase should instead be described by a defect-free model, which does not require partial occupancies. Here, the structure of cristobalite X-I consists of octahedral chains with a 4-60°-2 zigzag chain geometry. This geometry has not been previously considered but is closely related to post-quartz, stishovite, and seifertite. In addition to the ideal, defect-free crystal structure, we also present a description of the defects that are most likely to form within the X-I phase. Density functional theory calculations support our observations, confirming the dynamic stability of the X-I geometry and reasonably reproducing the pressure of the phase transformation. The enthalpy of cristobalite X-I is higher than stishovite and seifertite, but X-I is favored as a high-pressure successor of cristobalite due to a unique transformation pathway. Elastic and lattice dynamical properties of the X-I phase show intermediate values between stable tetrahedral and octahedral polymorphs, confirming the bridge-role of this phase in SiO 2 densification.

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

Modeling the Influence of Correlated Molecular Disorder on the Dynamics of Excitons in Organic Molecular Semiconductors

In this paper, we investigate the role of correlated molecular disorder on the dynamics of excitons in oligothiophene-based organic semiconductors. We simulate exciton dynamics using the Frenkel exciton model and derive parameters for this model so that they reflect the specific characteristics of all-atom molecular systems. By systematically modifying the parameters of the Frenkel exciton model, we isolate the influence of spatial and temporal molecular correlations on the dynamics of excitons in these systems. For the oligothiophene systems considered here, we find that the molecular fluctuations exhibit long-lived memory effects, but these effects do not significantly influence the dynamic properties of excitons. We also find that excitons can be sensitive to the molecular-scale spatial correlations, and that this sensitivity grows with the amount of energetic disorder within the material. We conclude that control over spatial correlations can mitigate the negative influence of disorder on exciton transport.

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

Inverse Design of Ultralow Lattice Thermal Conductivity Materials via Materials Database Screening of Lone Pair Cation Coordination Environment

The presence of lone pair (LP) electrons is strongly associated with the disruption of lattice heat transport, which is a critical component of strategies to achieve efficient thermoelectric energy conversion. By exploiting an empirical relationship between lattice thermal conductivity, κ L , and the bond angles of pnictogen group LP cation coordination environments, we develop an inverse design strategy based on a materials database screening to identify chalcogenide materials with ultralow κ L for thermoelectrics. Screening the ~635 000 real and hypothetical inorganic crystals of the Open Quantum Materials Database based on the constituent elements, nominal electron counting, LP cation coordination environment, and synthesizability, we identify 189 compounds expected to exhibit ultralow κ L . As a validation, we explicitly compute the lattice dynamical properties of two of the compounds (Cu 2 AgBiPbS 4 and MnTl 2 As 2 S 5 ) using first-principles calculations and successfully find both achieve ultralow κ L values at room temperature of ~0.3-0.4 W/(m·K) corresponding to the amorphous limit. Furthermore, our data-driven approach provides promising candidates for thermoelectric materials and opens new avenues for the design of phononic properties of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Polarity of Co-solvents Regulates the Charge Storage Mechanisms in Supercapacitors with Concentrated Electrolytes

Developing better energy storage devices depends on comprehending the underlying mechanisms involved in charge storage. With the continuous conception of new electrolytes, this task becomes progressively more urgent and complex. An example is the utilization of co-solvated concentrated solutions. While these show promising electrochemical responses, their dynamic properties (especially under confinement) and their relationships with performance are not fully understood. Here, we combined modified step potential electrochemical spectroscopy and quasielastic neutron scattering to investigate systems composed of activated mesoporous carbon (AMC) and concentrated solutions of lithium bis(trifluoromethanesulfonyl)imide in acetonitrile co-solvated with either toluene or acetone. Further, we report that acetone does not impair surface-controlled mechanisms, contrary to the case with toluene, which competes with charged species to populate the AMC’s pores without contributing to charge storage. In turn, toluene promotes a greater overall capacitance owing to Faradaic processes, which may be related to changes in the solvation structures under confinement.

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

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen oxygen partial correlation function at the experimental equilibrium density and 311 +/- 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water.

Gaiduk, Alex P.↗
DOE OSTI2021

Quantification of Mono- and Multivalent Counterion-Mediated Bridging in Polyelectrolyte Brushes

Multivalent counterion-induced bridging interactions have been identified as the key mechanism of drastic collapse of the height of polyelectrolyte (PE) brushes. In this article, we employ all-atom molecular dynamics (MD) simulations to quantify the bridging interactions in PE brushes for counterions of different sizes and valences. We identify that unlike the current notion, bridging interactions are not the sole function of the counterion valence. Rather the bridging interactions depend on the fraction of counterions (of a given type) that get physically condensed on the PE backbone as well as the size of the counterion solvation shell. These mechanisms ensure that certain monovalent counterions demonstrate much stronger bridging interactions than that witnessed for certain divalent and trivalent counterions, while certain counterions of identical valences show drastically different bridging. Here, we argue that these counterion-specific bridging interactions eventually enable not only the significant reduction of the PE brush height in presence of certain multivalent screening counterions, but may also give rise to scenarios where the brush height reduction for certain monovalent counterions is larger than certain divalent and trivalent counterions. This latter observation contradicts the experimental findings where the multivalent counterions invariably led to a larger decrease in the height of the PE brushes: we argue that this discrepancy stems from the fact that in our simulations we only consider densely grafted and short (and hence less flexible) PE brushes that hinder the formation of different laterally inhomogeneous structures (like pinned micelles and cylindrical bundles) that would have led to a larger brush height reduction (in experiments, which invariably consider longer and less densely grafted brushes, the formation of such inhomogeneous structures are primarily responsible for larger brush height reduction in presence of multivalent counterions). Finally, we also probe the dynamic properties of the counterions (i.e., their time-dependent displacements) and their bridging interactions (i.e., lifetime of bridging interactions).

42 ENGINEERING↗
DOE OSTI2021

Yielding Behavior of Bottlebrush and Linear Block Copolymers

Block copolymers can exhibit a pronounced yield stress, but the impact of molecular architecture, chemistry, and self-assembly on macroscopic rheology remains poorly understood. Here, we study the linear-viscoelastic and yield-stress fluid behavior of two architectures—bottlebrush copolymers (with statistical or blocky sequences) and linear diblocks—that self-assemble into body-centered cubic (BCC) spheres and hexagonally close-packed cylinders (HEX). The dynamic properties of these polymers were probed by oscillatory frequency and amplitude sweeps at temperatures well below the order–disorder transition (TODT) to furnish insights into the yielding transition. All BCC-forming polymers have a similar signature of yielding: smaller yield strains (γ yBCC ≈ 0.053 < γ y,HEX ≈ 0.18), sharper solid–liquid transitions, and better reversibility than HEX. Statistical bottlebrushes show the most frequency-independent structural modulus (G 0 ) and no signs of defect relaxation. A simple power-law relationship captures the dependence of the normalized structural modulus (G 0 /RT) on the inter-micelle distance (d) across different architectures and morphologies [G 0 /(RT) = 1.31 × 10 4 (nm 2.6 mol/m 3 ) d –2.6 ]. These studies establish quantitative structure–property relationships that are relevant in contemporary applications, for example, extrusion-based 3D printing.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Signatures of Thermoreversible Associations in X-ray and Neutron Scattering from Dilute Polyzwitterion Solutions

In aqueous solutions of polyzwitterions (PZs), an interplay between dipole–dipole interactions and hydration of zwitterionic groups can lead to thermoreversible associations, which have been difficult to detect in experiments. Here, in this study, we investigated dilute aqueous solutions of poly(1-(3-sulphonatopropyl)-2-vinylpyridinium) (P2VPPS) using small-angle X-ray and neutron scattering (SAXS and SANS) to probe the structure and neutron spin-echo (NSE) spectroscopy to probe dynamics. The SAXS and SANS data show that the correlation length increases with an increase in the concentration of P2VPPS for three different molecular weights. The addition of 0.1 M NaCl to one of the solutions led to almost no dependence of the correlation length on the concentration. Such a concentration dependence of the correlation length suggests the formation of clusters driven by thermoreversible associations in the solutions. The NSE measurements show that the solutions with a larger correlation length display slower relaxation, reflecting the reduced mobility of larger clusters. These results should be considered as signatures of thermoreversible associations in the solutions of P2VPPS. To establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS, we combined a thermoreversible gelation theory for the structure of PZ solutions (Li, S.-F.; Muthukumar, M. Theory of Thermoreversible Gelation and Anomalous Concentration Fluctuations in Polyzwitterion Solutions. J. Chem. Phys. 2024, 161, 024903) with a model for the dynamics of the clusters (generalized Zimm model), developed in this work. Using such a theoretical framework, we have predicted the distribution of clusters in the solutions probed with SANS and SAXS. With the distributions, the generalized Zimm model has been used to extract the diffusion constant of the clusters and their characteristic size from the NSE data, where the latter agrees with the values estimated from the SAXS/SANS data. These findings confirm the presence of thermoreversible associations and establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS. Furthermore, with growing interest in technological applications, this work can provide useful insights into the structural and dynamical properties of other PZ solutions.

field theory↗
DOE OSTI2020

Quantification of Dispersion for Weakly and Strongly Correlated Nanofillers in Polymer Nanocomposites

The dispersion of nanoparticles in viscous polymers is dictated by kinetics, interaction potentials between particles, and interfacial compatibility between the matrix and dispersed phases. It was previously proposed that an analogy can be made between thermally dispersed colloids and kinetically dispersed nanoparticles in viscous media when weak interactions exist between particles allowing for a mean-field description under the Ginzburg criterion such as for carbon black dispersed in polybutadiene elastomer. For these cases, the second virial coefficient can be used to quantify the quality of dispersion; additionally, the nanoscale network mesh size can be calculated, which is related to dynamic properties. However, this approach fails for nanoparticles with surface charges or other specific interactions that lead to correlations. Here, these correlated systems are investigated in the context of the mean-field systems in order to gain a comparative description of dispersion using the network mesh size and a derived virial coefficient. The physical origin of the structural parameters from the proposed model for these correlated systems is investigated.

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

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗
DOE OSTI2026

A Model Intercomparison Study to Investigate Mixing Characteristics in Non‐Precipitating Stratocumulus Clouds

Recent aircraft observations of marine stratocumulus clouds consistently showed that cloud microphysical relationships vary with altitude, indicating inhomogeneous mixing characteristics near cloud top and homogeneous mixing characteristics in mid-levels of clouds. Here, we conduct model intercomparison of an idealized, non-precipitating stratocumulus cloud to evaluate model consistency and examine whether simulations can reproduce the observed mixing characteristics. The results show that eleven large-eddy simulations with various dynamics and microphysics schemes show good agreement on the thermodynamical, microphysical, and dynamical properties of the stratocumulus-topped boundary layer in a steady state. The inter-model spread in steady-state liquid water path is significantly reduced compared to previous model intercomparison studies. This improvement might be due to better models and more consistent initial conditions than those used decades ago. In addition, most simulations, including a low-dimensional simulation, capture inhomogeneous mixing characteristics near the cloud top and homogeneous mixing characteristics inside the cloud. Moreover, simulations using Lagrangian microphysics schemes agree better with the observed mixing characteristics compared with those using the bin microphysics schemes. Since most simulations do not fully resolve the entrainment process, the apparent mixing characteristics arise from the variations in the resolved cloud properties. Our results support the vertical circulation mixing hypothesis, which suggests that homogeneous mixing characteristics in mid-levels of clouds are due to the vertical circulation of entrainment-affected and diluted parcels from the cloud top moved to lower levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamic control and quantification of active sites on ceria for CO activation and hydrogenation

Ceria (CeO 2 ) is a widely used oxide catalyst, yet the nature of its active sites remains elusive. This study combines model and powder catalyst studies to elucidate the structure-activity relationships in ceria-catalyzed CO activation and hydrogenation. Well-defined ceria clusters are synthesized on planar CeO 2 (111) and exhibit dynamic and tunable ranges of Ce coordination numbers, which enhance their interaction with CO. Reduced ceria clusters (e.g., Ce 3 O 3 ) bind CO strongly and facilitate its dissociation, while near-stoichiometric clusters (e.g., Ce 3 O 7 ) adsorb CO weakly and promote oxidation via carbonate formation. Unlike planar ceria surfaces, supported ceria clusters exhibit dynamic properties and enhanced catalytic activity, that mimic those of powder ceria catalysts. Insight from model studies provide a method to quantify active sites on powder ceria and guide further optimization of ceria catalysts for syngas conversion. This work marks a leap toward model-guided catalyst design and highlights the importance of site-specific catalysis.

03 NATURAL GAS↗
DOE OSTI2021

Intersystem crossing and exciton–defect coupling of spin defects in hexagonal boron nitride

Despite the recognition of two-dimensional (2D) systems as emerging and scalable host materials of single-photon emitters or spin qubits, the uncontrolled, and undetermined chemical nature of these quantum defects has been a roadblock to further development. Leveraging the design of extrinsic defects can circumvent these persistent issues and provide an ultimate solution. Here, we established a complete theoretical framework to accurately and systematically design quantum defects in wide-bandgap 2D systems. With this approach, essential static and dynamical properties are equally considered for spin qubit discovery. In particular, many-body interactions such as defect–exciton couplings are vital for describing excited state properties of defects in ultrathin 2D systems. Meanwhile, nonradiative processes such as phonon-assisted decay and intersystem crossing rates require careful evaluation, which competes together with radiative processes. From a thorough screening of defects based on first-principles calculations, we identify promising single-photon emitters such as Si VV and spin qubits such as Ti VV and Mo VV in hexagonal boron nitride. This work provided a complete first-principles theoretical framework for defect design in 2D materials.

36 MATERIALS SCIENCE↗