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 847 records · Page 47

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dimension-free path-integral molecular dynamics without preconditioning

Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics—including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods—is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann–Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of “dimension-free” numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi–Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

Korol, Roman (ORCID:0000000193076351)↗

Experimental study of mechanistic factors influencing solvent-driven fractional crystallization of calcium sulfate

To advance dimethyl ether-driven fractional crystallization (DME-FC), a more sustainable method of water treatment and mineral recovery, a range of chemical equilibria were measured. These include varying concentrations of miscible organic solvents (MOS) used to experimentally measure the solvent-induced solid-liquid equilibrium (SLE) of calcium sulfate (CaSO 4 ) in water. Seven MOS, including dimethyl ether (DME), acetonitrile (MeCN), 1,4-dioxane, tetrahydrofuran (THF), acetone, ethanol, and diethylamine, were screened to establish trends associated with molecular volume, functional groups, and physical properties. The effect of MOS on CaSO 4 removal differed at concentrations <0.15 mol fraction MOS; MOS with greater molecular volume (THF, 1,4-dioxane, and diethylamine) induced greater CaSO 4 precipitation on a per mole basis. The solvent-induced SLE for all MOS converged between 0.15 and 0.2 mol fraction MOS, reaching a CaSO 4 concentration consistent with a water to MOS hydration ratio of 5:1 to 6:1, which may correspond to the solvent generating a solution-based pseudo-clathrate structure with continuity within the solution. Finally, solution pseudo-clathrate structures provide a mechanistic basis for DME-FC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Vibrational Polaritons Towards Quantum Technologies

In this paper molecular vibrational polaritons (MVPs)—a hybrid molecular-photon quasiparticle—and the development of a proof-of-principle quantum technology platform are discussed to simulate coherence transfer, for use at room temperature. It is shown that MVPs can form qubits, coherence, and have nonlinear interactions, all at room temperature. Some new insights, such as polaritonic nonlinearity dependence on macroscopic properties including cavity thickness and molecular concentrations are also uncovered. It is hoped that these advances can stimulate more research in developing this system into a quantum technology platform free from the constraints imposed by the requirement of cryogenic conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics in Star-forming Cores (DiSCo): project overview and the first look towards the B1 and NGC 1333 regions in Perseus

The internal velocity structure within dense gaseous cores plays a crucial role in providing the initial conditions for star formation in molecular clouds. However, the kinematic properties of dense gas at core scales (~0.01–0.1 pc) has not been extensively characterized because of instrument limitations until the unique capabilities of GBT-Argus became available. The ongoing GBT-Argus Large Program, Dynamics in Star-forming Cores (DiSCo) thus aims to investigate the origin and distribution of angular momentum of star-forming cores. DiSCo will survey all starless cores and Class 0 protostellar cores in the Perseus molecular complex down to ~0.01 pc scales with <0.05 km s –1 velocity resolution using the dense gas tracer N 2 H + . Here, we present the first data sets from DiSCo towards the B1 and NGC 1333 regions in Perseus. Our results suggest that a dense core’s internal velocity structure has little correlation with other core-scale properties, indicating these gas motions may be originated externally from cloud-scale turbulence. These first data sets also reaffirm the ability of GBT-Argus for studying dense core velocity structure and provided an empirical basis for future studies that address the angular momentum problem with a statistically broad sample.

79 ASTRONOMY AND ASTROPHYSICS↗

Combustibility Tests of 1,1,1,2-tetrafluoroethane in a Simulated Compressor Cylinder

The advantages of high-molecular-weight gas (heavy gas) as a wind-tunnel medium have been recognized for some time. The current heavy gas of choice chlorofluorocarbon-12(CFC-12) (refrigerant R12) for the Transonic Dynamics Tunnel(TDT) must be replaced because manufacture of this gas ceased in 1995. An attractive replacement is 1,1,1,2-tetrafluoroethane (refrigerant R134a). Acceptable properties of this gas include molecular weight and speed of sound. Its vapor pressure allows simplified reclamation from mixtures with air. However, it is recognized that R134a is combustible under certain conditions of temperature, pressure, and concentration. A comprehensive study was conducted to identify those conditions and the influence of various parameters on the combustibility of the gas-air mixture.

Babcock, Dale A.↗

Deciphering between enhanced light emission and absorption in multi-mode porphyrin cavity polariton samples

Abstract It remains unclear how the collective strong coupling of cavity-confined photons to the electronic transitions of molecular chromophore leverages the distinct properties of the polaritonic constituents for future technologies. In this study, we design, fabricate, and characterize multiple types of Fabry-Pérot (FP) mirco-resonators containing copper(II) tetraphenyl porphyrin (CuTPP) to show how cavity polariton formation affects radiative relaxation processes in the presence of substantial non-Condon vibronic coupling between two of this molecule’s excited electronic states. Unlike the prototypical enhancement of Q state radiative relaxation of CuTPP in a FP resonator incapable of forming polaritons, we find the light emission processes in multimode cavity polariton samples become enhanced for cavity-exciton energy differences near those of vibrations known to mediate non-Condon vibronic coupling. We propose the value of this detuning is consistent with radiative relaxation of Herzberg-Teller polaritons into collective molecular states coupled to the cavity photon coherently. We contrast the feature stemming from light emission from the HT polariton state with those that occur due to polariton-enhanced light absorption. Our results demonstrate the landscape of molecular and photonic interactions enabled by cavity polariton formation using complex chromophores and how researchers can design resonators to leverage these interactions to characterize and control polaritonic properties.

Odewale, Elizabeth O.↗

Characterizing Transiting Planets with JWST Spectra: Simulations and Retrievals

There are now well over a thousand confirmed exoplanets, ranging from hot to cold and large to small worlds. JWST spectra will provide much more detailed information on the molecular constituents, chemical compositions, and thermal properties of the atmospheres of transiting planets than is now known. We explore this by modeling clear, cloudy,and high mean molecular weight atmospheres of typical hot Jupiter, warm Neptune, warm sub-Neptune, and cool super-Earth planets and then simulating their JWST transmission and emission spectra. These simulations were performed for several JWST instrument modes over 1 - 11 microns and incorporate realistic signal and noise components. We then performed state-of the art retrievals to determine how well temperatures and abundances (CO, CO2, H2O, NH3) will be constrained and over what pressures for these different planet types. Using these results, we appraise what instrument modes will be most useful for determining what properties of the different planets, and we assess how well we can constrain their compositions, CO ratios, and temperature profiles.

extrasolar planets↗

Characterizing Transiting Planets with JWST Spectra: Simulations and Retrievals

There are now well over a thousand confirmed exoplanets, ranging from hot to cold and large to small worlds. JWST spectra will provide much more detailed information on the molecular constituents, chemical compositions, and thermal properties of the atmospheres of transiting planets than is now known. We explore this by modeling clear, cloudy,and high mean molecular weight atmospheres of typical hot Jupiter, warm Neptune, warm sub-Neptune, and cool super-Earth planets and then simulating their JWST transmission and emission spectra. These simulations were performed for several JWST instrument modes over 1 - 11 microns and incorporate realistic signal and noise components. We then performed state-of the art retrievals to determine how well temperatures and abundances (CO, CO2, H2O, NH3) will be constrained and over what pressures for these different planet types. Using these results, we appraise what instrument modes will be most useful for determining what properties of the different planets, and we assess how well we can constrain their compositions, CO ratios, and temperature profiles.

James Webb Space Telescope↗

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗

Structure–property relations of sodium iron phosphate nuclear waste glasses: Effects of iron redox ratio and glass composition

Iron phosphate glasses, known for their exceptional chemical durability and potential applicability in nuclear waste management, have gained significant attention over the years. The structures of these glasses are complicated by the coexistence of Fe 3+ and Fe 2+ , which plays a crucial role in determining their structures and properties. Here, this work uses molecular dynamics simulations to study the structural changes in Na 2 O–Fe 2 O 3 –P 2 O 5 glasses with varying glass composition and Fe 2+ /Fe 3+ redox ratio. It was found that the redox ratio and modifier contents significantly affected the short-range and medium-range orders in the glasses. Significant changes in the local environments around P 5+ and Fe 3+ were observed, as reflected by the bond distances and coordination numbers. Na + cations are found to preferentially associate with Fe 3+ (rather than Fe 2+ ), whereas Fe 2+ has stronger association with P 5+ than Na + , confirming the structural role of Fe 2+ as a glass modifier. The disruptions in P–O–P linkages upon increasing FeO suggest that FeO causes glass depolymerization. These glasses achieved higher connectivity with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, conerting phosphorous Q 2 to Q 3 units and iron Q 5 units to Q 4 units. The decrease of nonbridging oxygen fractions with increasing Fe 3+ / (Fe 3+ + Fe 2+ ) ratios, through creating P–O–Fe linkages, is the main reason of enhanced network connectivity. Quantitative structure–property relationship analyses with different structural descriptors were used to correlate with measured properties. The analyses provided valuable insights into structure–property relationships, emphasizing the importance of choosing relevant energy parameters and defining glass network connectivity, particularly in F net descriptors. It was found the Fe–O–P linkage density exhibits strong correlations to measured dissolution rates, supporting the importance of these linkages in improving the chemical durability in iron phosphate glasses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The FCRAO extragalactic CO survey: Global properties of galaxies

Since stars form in molecular clouds, a critical element in studies of galaxy evolution is knowledge of the molecular content of a large sample of galaxies. To this end, researchers undertook a survey of CO emission from galaxies using the Fire College Radio Astronomy Observatory (FCRAO) 14-m millimeter telescope at 115 GHz. The aim was to better understand the differences found among and within galaxies with regard to the efficiency of star and cloud formation. The galaxies observed as part of the FCRAO Extragalactic CO Survey were selected on the basis of their optical or infrared properties. The galaxies observed thus far are (1) brighter than B sub T (sub o)=12.5 in the blue, or (2) brighter than 20 Jy at 100 microns. From major axis CO observations at 45 seconds resolution and spacing in over 200 galaxies, researchers determined the CO radial distributions, and derived global CO fluxes (cf. Kenney and Young 1988); H2 masses were derived using the conversion factor N(H2)/I sub CO=2.8 times 10 to the 20th power cm(-2)/K km s(-1)) (Bloemen et al. 1986). Here, researchers concentrate on the global galaxy properties within the sample. Neutral hydrogen (HI) masses for the sample galaxies were taken from Huchtmeier et al. (1983), blue luminosities and morphological types were taken from RC2. IR luminosities, colors, dust temperatures and dust masses were determined from coadded Infrared Astronomy Satellite (IRAS) data (Young et al. 1989). They have chosen to first compare absolute luminosities and masses in order to determine the slope and scatter in each correlation; next they investigate luminosity independent ratios in order to intercompare large and small galaxies.

Young, J. S.↗

Computational Chemistry of Adhesive Bonds

This investigation is intended to determine the electrical mechanical, and chemical properties of adhesive bonds at the molecular level. The initial determinations will be followed by investigations of the effects of environmental effects on the chemistry and properties of the bond layer.

Phillips, Donald H.↗

Formulation/cure technology for ultrahigh molecular weight silphenylene-siloxane polymers

Molecular weights above one million were achieved for methylvinylsilphenylene-siloxane terpolymers using a two-stage polymerization technique which was successfully scaled up to 200 grams. The resulting polymer was vulcanized by two different formulations and compared to an identically formulated commercial methylvinyl silicone on the basis of ultimate strength, Young's modulus, percent elongation at failure, and tear strength. Relative thermal/oxidative stabilities of the elastomers were assessed by gradient and isothermal thermogravimetric analyses performed in both air and nitrogen. The experimental elastomer exhibited enhanced thermal/oxidative stability and possed equivalent or superior mechanical properties. The effect of variations in prepolymer molecular weight on mechanical properties was also investigated.

Hundley, N. H.↗

Molecular tail chemistry controls thermal transport in fullerene films

Herein we report on the thermal conductivities of alkyl- and indene-group functionalized fullerene derivative thin films as measured via time domain thermoreflectance and steady state thermoreflectance. The thermal conductivities vary from 0.064 ± 0.007 W m – 1 K – 1 for [6,6]-phenyl C 61 -butyric acid methyl ester (PCBM) to 0.15 ± 0.017 W m – 1 K – 1 for indene- C 60 bisadduct at room temperature and do not exhibit significant temperature dependence from 300 to 375 K. In comparison to the thermal conductivity of PCBM, increasing the length of the alkyl chain, as in the case of [6,6]-phenyl C 61 -butyric acid butyl ester, and [6,6]-phenyl C 61 -butyric acid octyl ester leads to higher thermal conductivities. Likewise, increasing the number of alkyl chains attached to the fullerenes as in the case of bisadduct PCBM leads to a higher thermal conductivity compared to that of PCBM. We present atomistic insights into the role of chemical functionalization on the overall heat transfer in these fullerene derivatives by conducting molecular dynamics simulations and lattice dynamics calculations. The thermal conductivities predicted via our atomistic simulations qualitatively agree with the experimental trends for our fullerene derivatives. Lattice dynamics calculations reveal that one of the main factors dictating the ultralow thermal conductivities in fullerene derivatives is the large reduction in modal diffusivities in the molecular crystals as calculated from the Allen-Feldman model, thus providing an explanation for their largely reduced thermal conductivities as compared to that of C 60 crystals. The low diffusivities result from high degrees of localization of Einstein-like vibrations in fullerene derivatives due to the molecular side chains, providing the ability to dial-in the properties of these low thermal conductivity solids via molecular engineering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Atomistic simulations of nanoindentation on nanoglasses: Effects of grain size and gradient microstructure on the mechanical properties

The use of nanoglass (NG) microstructures has been shown as an effective strategy to improve the ductility of rather brittle metallic glass alloys. To evaluate the effects of grain size and gradient design on the mechanical properties of NGs, we perform molecular dynamics simulations of nanoindentation on Cu 64 Zr 36 NGs. We consider samples with uniform 3 and 7 nm grain sizes as well as gradient microstructures with gradient grain sizes varying from 3 to 7 nm. Here, the results show that the deformation mechanism in NGs with small grain sizes is dominated by the activation and evolution of multiple shear transformation zones. Increasing grain sizes enhances the elastic modulus and hardness at the cost of reduced abrasion resistance. Though the average grain size at the indentation surface plays a crucial role in the deformation behavior of gradient NG models, their strain localization and plastic deformation states are affected by the grain sizes far from the indenter. Desired combinations of mechanical properties can be realized by different NG microstructure designs.

36 MATERIALS SCIENCE↗

Amphiphilic Heterografted Molecular Bottlebrushes with Tertiary Amine-Containing Side Chains as Efficient and Robust pH-Responsive Emulsifiers

Here, by combining the unique characteristics of molecular bottlebrushes (MBBs) and the properties of stimuli-responsive polymers, we show that MBBs with randomly grafted poly(n-butyl acrylate) and pH-responsive poly(2-(N,N-diethylamino)ethyl methacrylate) (PDEAEMA) side chains are efficient and robust pH-responsive emulsifiers. Water-in-toluene emulsions were formed at pH 4.0 and disrupted by increasing the pH to 10.0. The emulsion generation and disruption was reversible over the ten cycles investigated, and the bottlebrushes remained intact. The exceptional emulsion stability stemmed from the high interfacial binding energy of MBBs, imparted by their large molecular size and Janus architecture at the interface, as evidenced by the interfacial jamming and wrinkling of the assemblies upon reducing the interfacial area. At pH 10.0, PDEAEMA became water-insoluble, and the MBBs desorbed from the interface, causing de-emulsification. Consequently, we have shown that the judicious design of MBBs can generate properties of particle emulsifiers from their large size, while the responsiveness of the MBBs enables more potential applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of layer bending on montmorillonite hydration and structure from molecular simulation

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm 2 system that included dehydrated and hydrated regions in the same interlayer.

58 GEOSCIENCES↗