Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Small systems”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Crystallization kinetics and nanoparticle ordering in semicrystalline polymer nanocomposites

There has been considerable interest in the nucleation and crystallization of polymers in the presence of nanoparticles (NPs, or nanofillers in general, NFs). Most of the extensive work in this area has focused on anisotropic, non-Brownian NFs (e.g., clay sheets, carbon nanotubes) whose spatial dispersion state in these nanocomposites is controlled by the process by which they are formed. Hence, NF spatial dispersion is generally limited and often remains poorly characterized. Thermodynamic handles that can be used to control NF dispersion state in the polymer melt include (a) favorable interactions between the polymer chains and the bare NP surfaces, or (b) the density and length of the chains, with the same chemistry as the matrix, grafted to the NP surface. These relatively large NFs merely act as stationary objects that affect the kinetics of nucleation by providing heterogeneous sites, and the crystallization rate by confining the polymer in the melt state. The dispersion state of the NFs can dramatically affect the nucleation and crystallization of the matrix, but in most cases reported, the NFs increase nucleation efficiency relative to the neat polymer. At higher NF loadings, the effect of polymer confinement by the NFs dominates, leading to a decrease in crystal growth rates. This review describes the most important lessons learned from these commonly studied systems and then extends to polymer composite systems containing small, mobile spherical NPs (typically smaller than 100 nm in size). The role of NP mobility, which provides for dynamic confinement of the polymer melt, on the kinetics of polymer crystallization (nucleation, growth, and overall crystallization) and how this behavior is mostly consistent with the case of immobile NF is a second important focus of this review. In addition to the role of NFs on crystallization kinetics, recently reported nanoparticle ordering phenomena such as the effect of matrix crystallization on the organization of small spherical NPs within the amorphous regions of the semicrystalline morphology are discussed. In conclusion, such phenomena are clearly not observed for large NFs and hence provide a point of departure from past works in this area.

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

InteraChem: Virtual Reality Visualizer for Reactive Interactive Molecular Dynamics

Interactive molecular dynamics in virtual reality (IMD-VR) simulations provide a digital molecular playground for students as an alternative or complement to traditional molecular modeling kits or 2D illustrations. Previous IMD-VR studies have used molecular mechanics to enable simulations of macromolecules such as proteins and nanostructures for the classroom setting with considerable success. Here, we present the InteraChem molecular visualizer, intended for reactive IMD-VR simulation using semiempirical and ab initio methods. InteraChem visualizes not only the molecular geometry but also (1) isosurfaces such as molecular orbitals and electrostatic potentials and (2) two-dimensional graphs of time-varying simulation quantities such as kinetic/potential energy, internal coordinates, and user-applied force. Additionally, InteraChem employs speech recognition to facilitate user interaction and introduces a novel “atom happiness” visualization using emojis to indicate the energetic feasibility of a particular bonding arrangement. Furthermore, we include a set of accompanying exercises that we have used to teach chemical reactivity in small molecular systems.

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

Real-time coupled-cluster approach for the cumulant Green's function

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.

Vila, Fernando D.↗
DOE OSTI2022

GronOR: Scalable and Accelerated Nonorthogonal Configuration Interaction for Molecular Fragment Wave Functions

GronOR is a program package for nonorthogonal configuration interaction calculations. Electronic wave functions are constructed in terms of antisymmetrized products of multiconfiguration molecular fragment wave functions. The computational complexity of the nonorthogonal methodologies implemented in GronOR applied to large molecular assemblies requires a design that takes full advantage of massively parallel supercomputer architectures and accelerator technologies. This work describes the implementation strategy and resulting performance characteristics. In addition to parallelization and acceleration, the software development strategy includes aspects of fault resiliency and heterogeneous computing. The program was designed for large-scale supercomputers but also runs effectively on small clusters and workstations for small molecular systems. GronOR is available as open source to the scientific community.

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

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Scalable Implementation of Mean-Field and Correlation Methods Based on Lie-Algebraic Similarity Transformation of Spin Hamiltonians in the Jordan–Wigner Representation

Recent work has highlighted that the strong correlation inherent in spin Hamiltonians can be effectively reduced by mapping spins to Fermions via the Jordan−Wigner transformation (JW). The Hartree−Fock method is straightforward in the Fermionic domain and may provide a reasonable approximation to the ground state. Correlation with respect to the Fermionic mean field can be recovered based on Lie-algebraic similarity transformation (LAST) with two-body correlators. Specifically, a unitary LAST variant eliminates the dependence on site ordering, while a nonunitary LAST yields size-extensive correlation energies. Whereas the first recent demonstration of such methods was restricted to small spin systems, we present efficient implementations using analytical gradients for the optimization with respect to the mean-field reference and the LAST parameters, thereby enabling the treatment of larger clusters, including systems with local spins s > $\frac{1}{2}$.

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

Improved Treatment of 1–4 Interactions in Force Fields for Molecular Dynamics Simulations

Traditional force fields commonly use a combination of bonded torsional terms and empirically scaled nonbonded interactions to capture 1-4 energies and forces of atoms separated by three bonds in a molecule. While this approach can yield accurate torsional energy barriers, it often leads to inaccurate forces and erroneous geometries and creates an interdependence between dihedral terms and nonbonded interactions, complicating parametrization and reducing transferability. Here, in this paper, we demonstrate that 1-4 interactions can be accurately modeled using only bonded coupling terms, eliminating the need for arbitrarily scaled nonbonded interactions altogether. Furthermore, by leveraging the automated parametrization capabilities of the Q-Force toolkit, we efficiently determine the necessary coupling terms without the need for manual adjustment. Our approach is first validated on a range of small molecule systems, encompassing both flexible and rigid structures, and shows a significant improvement in force field accuracy, obtaining subkcal/mol mean absolute error for every molecule tested. We further extend the bonded-only model for 1-4 interactions to Amber ff14sb, CHARMM36, and OPLS-AA force fields to reproduce ab initio gas and implicit solvent ϕ,ψ surfaces of alanine dipeptide.

Abdullah, Aalim S. [University of California, Berk↗
DOE OSTI2025

Visualizing the Three-Dimensional Arrangement of Hydrogen Atoms in Organic Molecules by Coulomb Explosion Imaging

Structure-sensitive methods based on femtosecond light or electron pulses are now making it possible to measure how molecular structures change during light-induced processes. Despite significant progress, high-fidelity imaging of nuclear positions remains a challenge even for relatively small molecular systems and, notably, regarding the positions of hydrogen atoms. As demonstrated in recent work, X-ray-induced Coulomb explosion imaging (CEI) may overcome this obstacle, as its sensitivity does not depend on the mass of the imaged atoms. The photoinduced ring opening of the heterocyclic molecule 2(5 H )-thiophenone has attracted recent interest. Here, in this work, we show that CEI offers a powerful route to imaging the peripheral H atoms in this molecule and thus, more generally, to tracking detailed nuclear motions (e.g., isomerizations) in organic molecules on ultrafast time scales. Specifically, we record momentum-space Coulomb explosion images that report on the three-dimensional positioning of all nuclei within the molecule, for instance, distinguishing H atoms in C–H bonds that lie within or are directed out of the plane defined by the heavy atoms. The prospect of imaging peripheral H atoms to probe photochemical dynamics is explored by coupling ab initio molecular dynamics with classical Coulomb explosion simulations, thereby differentiating potential photoproduct isomers, including those whose structures primarily differ in the position of the hydrogens.

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

A High-Fidelity Molecular Model of the Cu(111) Repeating Unit

Dynamic processes at surfaces are central to heterogeneous catalysis, but their atomistic mechanism(s) can prove difficult to elucidate due to variations in material structure and the corresponding impact on reactivity. Moreover, disparities between reaction conditions and those employed for spectroscopic characterization at surfaces can inhibit detailed understanding of catalysis-relevant chemistries. Herein, we substantiate the so-called “cluster-surface” analogy by leveraging a low-valent tricopper architecture ( 1 ) as a model system for small molecule activation at Cu(111). Two reaction classes are explored: the adsorption of carbon monoxide (CO) and the dissociative adsorption of dihydrogen (H 2 ). These processes serve as an ideal testbed to compare the reactivity of a molecular cluster ( 1 ) to that of a heterogeneous surface, as both reactions have empirical data from measurements performed on crystalline Cu(111). Cluster 1 reversibly binds CO. Variable temperature NMR analysis with 13 CO reveals a favorable enthalpy but large negative entropy (−5.1 kcal × mol –1 and −22.9 cal × mol –1 × K –1 , respectively) for CO binding, affording a process that is marginally endergonic at room temperature (ΔG ads (298.15 K) = 1.7 ± 0.5 kcal × mol –1 ). Similarly, analogous to a Cu(111) surface, 1 is shown to oxidatively add (chemisorb) H 2 . Kinetic parameters were determined for this process and the activation enthalpy (8.4 ± 0.5 kcal × mol –1 ) closely mirrors that established for H 2 binding at the Cu(111) facet (6.0 to 12.4 kcal × mol –1 ). Together, these results showcase that a trinuclear cluster can reproduce the small molecule binding and activation energetics of a bulk crystalline surface, setting the stage for studying less-defined surface processes in an atomically precise molecular setting.

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

Orientation and Stability of Asymmetric Magnetic Reconnection X Line

Here, the orientation and stability of the reconnection x line in asymmetric geometry is studied using three-dimensional (3-D) particle-in-cell simulations. We initiate reconnection at the center of a large simulation domain to minimize the boundary effect. The resulting x line has sufficient freedom to develop along an optimal orientation, and it remains laminar. Companion 2-D simulations indicate that this x line orientation maximizes the reconnection rate. The divergence of the nongyrotropic pressure tensor breaks the frozen-in condition, consistent with its 2-D counterpart. We then design 3-D simulations with one dimension being short to fix the x line orientation but long enough to allow the growth of the fastest growing oblique tearing modes. This numerical experiment suggests that reconnection tends to radiate secondary oblique tearing modes if it is externally (globally) forced to proceed along an orientation not favored by the local physics. Finally, the development of oblique structure easily leads to turbulence inside small periodic systems.

3-D magnetic reconnection↗
DOE OSTI2022

Transient Storage Model Parameter Optimization Using the Simulated Annealing Method

Hyporheic exchange in streams is critical to ecosystem functions such as nutrient cycling along river corridors, especially for slowly moving or small stream systems. The transient storage model (TSM) has been widely used for modeling of hyporheic exchange. TSM calibration, for hyporheic exchange, is typically used to estimate four parameters, including the mass exchange rate coefficient, the dispersion coefficient, stream cross-sectional area, and hyporheic zone cross-sectional area. Prior studies have raised concerns regarding the non-uniqueness of the inverse problem for the TSM, that is, the occurrence of different parameter vectors resulting in TSM solution that reproduces the observed in-stream tracer break through curve (BTC) with the same error. This leads to practical non-identifiability in determining the unknown parameter vector values even when global-optimal values exist, and the parameter optimization becomes practically non-unique. To address this problem, we applied the simulated annealing method to calibrate the TSM to BTCs, because it is less susceptible to local minima-induced non-identifiability. A hypothetical (or synthetic) tracer test data set with known parameters was developed to demonstrate the capability of the simulated annealing method to find the global minimum parameter vector, and it identified the “hypothetically-true” global minimum parameter vector even with input data that were modified with up to 10% noise without increasing the number of iterations required for convergence. The simulated annealing TSM was then calibrated using two in-stream tracer tests conducted in East Fork Poplar Creek, Tennessee. In conclusion, simulated annealing was determined to be appropriate for quantifying the TSM parameter vector because of its search capability for the global minimum parameter vector.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Monolayer semiconductor nanocavity lasers with ultralow thresholds

Engineering the electromagnetic environment of a nanoscale light emitter by a photonic cavity can significantly enhance its spontaneous emission rate through cavity quantum electrodynamics in the Purcell regime. This effect can greatly reduce the lasing threshold of the emitter 1–5 , providing the ultimate low-threshold laser system with small footprint, low power consumption and ultrafast modulation. A state-of-the-art ultra-low threshold nanolaser has been successfully developed though embedding quantum dots into photonic crystal cavity (PhCC) 6–8 . However, several core challenges impede the practical applications of this architecture, including the random positions and compositional fluctuations of the dots 7 , extreme difficulty in current injection8, and lack of compatibility with electronic circuits 7,8 . Here, we report a new strategy to lase, where atomically thin crystalline semiconductor, i.e., a tungsten-diselenide (WSe 2 ) monolayer, is nondestructively and deterministically introduced as a gain medium at the surface of a pre-fabricated PhCC. A new type of continuous-wave nanolaser operating in the visible regime is achieved with an optical pumping threshold as low as 27 nW at 130 K, similar to the value achieved in quantum dot PhCC lasers 7 . The key to the lasing action lies in the monolayer nature of the gain medium, which confines direct-gap excitons to within 1 nm of the PhCC surface. The surface-gain geometry allows unprecedented accessibilities to multi-functionalize the gain, enabling electrically pumped operation. Our scheme is scalable and compatible with integrated photonics for on-chip optical communication technologies.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Tuning the driving force for exciton dissociation in single-walled carbon nanotube heterojunctions

Understanding the kinetics and energetics of interfacial electron transfer in molecular systems is crucial for the development of a broad array of technologies, including photovoltaics, solar fuel systems and energy storage. The Marcus formulation for electron transfer relates the thermodynamic driving force and reorganization energy for charge transfer between a given donor/acceptor pair to the kinetics and yield of electron transfer. Here we investigated the influence of the thermodynamic driving force for photoinduced electron transfer (PET) between single-walled carbon nanotubes (SWCNTs) and fullerene derivatives by employing time-resolved microwave conductivity as a sensitive probe of interfacial exciton dissociation. For the first time, we observed the Marcus inverted region (in which driving force exceeds reorganization energy) and quantified the reorganization energy for PET for a model SWCNT/acceptor system. The small reorganization energies (about 130 meV, most of which probably arises from the fullerene acceptors) are beneficial in minimizing energy loss in photoconversion schemes.

14 SOLAR ENERGY↗
DOE OSTI2019

Topology, landscapes, and biomolecular energy transport

While ubiquitous, energy redistribution remains a poorly understood facet of the nonequilibrium thermodynamics of biomolecules. At the molecular level, finite-size effects, pronounced nonlinearities, and ballistic processes produce behavior that diverges from the macroscale. Here, we show that transient thermal transport reflects macromolecular energy landscape architecture through the topological characteristics of molecular contacts and the nonlinear processes that mediate dynamics. While the former determines transport pathways via pairwise interactions, the latter reflects frustration within the landscape for local conformational rearrangements. Unlike transport through small-molecule systems, such as alkanes, nonlinearity dominates over coherent processes at even quite short time- and length-scales. Our exhaustive all-atom simulations and novel local-in-time and space analysis, applicable to both theory and experiment, permit dissection of energy migration in biomolecules. The approach demonstrates that vibrational energy transport can probe otherwise inaccessible aspects of macromolecular dynamics and interactions that underly biological function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Universal and versatile morphology engineering via hot fluorous solvent soaking for organic bulk heterojunction

After explosive growth of efficiency in organic solar cells (OSCs), achieving ideal morphology of bulk heterojunction remains crucial and challenging for advancing OSCs into consumer market. Herein, by utilizing the amphiphobic nature and temperature-dependent miscibility of fluorous solvent, hot fluorous solvent soaking method is developed to optimize the morphology with various donor/acceptor combinations including polymer/small-molecule, all-polymer and all-small-molecule systems. By immersing blend film into hot fluorous solvent which is utilized as liquid medium with better thermal conductivity, the molecular reorganization is accelerated. Furthermore, fluorous solvent can be miscible with the residue of chloroform and chloronaphthalene above upper critical solution temperature. This mixed solvent diffuses around inside the active layer and selectively promotes molecular reorganization, leading to optimized morphology. Compared to widely-used thermal annealing, this approach processed under mild conditions achieves superior photovoltaic performance, indicating the practicality and universality for morphological optimization in OSCs as well as other optoelectronic devices.

14 SOLAR ENERGY↗
DOE OSTI2024

Development of a micro-combined heat and power powered by an opposed-piston engine in building applications

Residential homes and light commercial buildings usually require substantial heat and electricity simultaneously. A combined heat and power system enables more efficient and environmentally friendly energy usage than that achieved when heat and electricity are produced in separate processes. However, due to financial and space constraints, residential and light commercial buildings often limit the use of traditional large-scale industrial equipment. Here we develop a micro–combined heat and power system powered by an opposed-piston engine to simultaneously generate electricity and provide heat to residential homes or light commercial buildings. The developed prototype attains the maximum AC electrical efficiency of 35.2%. The electrical efficiency breaks the typical upper boundary of 30% for micro–combined heat and power systems using small internal combustion engines (i.e., <10 kW). Moreover, the developed prototype enables maximum combined electrical and thermal efficiencies greater than 93%. The prototype is optimally designed for natural gas but can also run renewable biogas and hydrogen, supporting the transition from current conventional fossil fuels to zero carbon emissions in the future. The analysis of the unit’s decarbonization and cost-saving potential indicate that, except for specific locations, the developed prototype might excel in achieving decarbonization and cost savings primarily in US northern and middle climate zones.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Electronic structure with direct diagonalization on a D-wave quantum annealer

Quantum chemistry is regarded to be one of the first disciplines that will be revolutionized by quantum computing. Although universal quantum computers of practical scale may be years away, various approaches are currently being pursued to solve quantum chemistry problems on near-term gate-based quantum computers and quantum annealers by developing the appropriate algorithm and software base. This work implements the general Quantum Annealer Eigensolver (QAE) algorithm to solve the molecular electronic Hamiltonian eigenvalue-eigenvector problem on a D-Wave 2000Q quantum annealer. The approach is based on the matrix formulation, efficiently uses qubit resources based on a power-of-two encoding scheme and is hardware-dominant relying on only one classically optimized parameter. We demonstrate the use of D-Wave hardware for obtaining ground and excited electronic states across a variety of small molecular systems. The approach can be adapted for use by a vast majority of electronic structure methods currently implemented in conventional quantum-chemical packages. The results of this work will encourage further development of software such as qbsolv which has promising applications in emerging quantum information processing hardware and has expectation to address large and complex optimization problems intractable for classical computers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Tethered tertiary amines as solid-state n-type dopants for solution-processable organic semiconductors

A scarcity of stable n-type doping strategies compatible with facile processing has been a major impediment to the advancement of organic electronic devices. Localizing dopants near the cores of conductive molecules can lead to improved efficacy of doping. We and others recently showed the effectiveness of tethering dopants covalently to an electron-deficient aromatic molecule using trimethylammonium functionalization with hydroxide counterions linked to a perylene diimide core by alkyl spacers. In this work, we demonstrate that, contrary to previous hypotheses, the main driver responsible for the highly effective doping observed in thin films is the formation of tethered tertiary amine moieties during thin film processing. Furthermore, we demonstrate that tethered tertiary amine groups are powerful and general n-doping motifs for the successful generation of free electron carriers in the solid-state, not only when coupled to the perylene diimide molecular core, but also when linked with other small molecule systems including naphthalene diimide, diketopyrrolopyrrole, and fullerene derivatives. Our findings help expand a promising molecular design strategy for future enhancements of n-type organic electronic materials.

36 MATERIALS SCIENCE↗