SEARCH · Engineering Papers

Results for “Boundary condition”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2022

The estimation of Poisson’s ratio by time-averaging and Cornu’s method for isotropic beams

Prior research suggests that direct static (e.g., uniaxial testing) and dynamic (e.g., ultrasonic wave propagation analysis) measurements differ in the estimation of Poisson’s ratio because anisotropies and heterogeneities in the sample material affect the two types of tests differently. Even assuming isotropic and homogeneous material properties, prior research further suggests that discrepancies between static/dynamic test results will exist because the error of the diagnostic techniques for the measurand are inherently different. Finally, thermodynamic effects are not present in static tests but can significantly affect dynamic test results. Given the potential for all these variables to produce discrepancies, it would be helpful to have the measurement of Poisson’s ratio obtainable from the same theory and experimental measurements by either static or dynamic testing methods. Our finite element calculations show that by combining time-averaged scanning digital holography with Cornu’s method, it is theoretically possible to estimate the effective Poisson’s ratio from the anticlastic contours at the antinode of the first out-of-plane bending mode shape. This is true regardless of frequency and therefore applicable for both static and dynamic measurements. Our results show that the estimate of Poisson’s ratio by Cornu’s method using data from simulations of mode shapes approaches the true value of Poisson’s ratio. Additionally, our research suggests that beam geometry and boundary conditions are fundamental factors limiting the convergence of the estimate of Poisson’s ratio to the true value of Poisson’s ratio regardless of performing a static or dynamic test.

42 ENGINEERING↗

DOE OSTI2026

Prediction of dynamic behavior of a system undergoing an exchange of components using a transmission simulator method

In systems with modular design, it is of interest to perform component exchanges to upgrade components and explore new system configurations. If the dynamics of the components are known, then dynamic substructuring can be used to obtain a prediction of system performance for the new configuration. The transmission simulator method is a substructuring technique that utilizes an additional fixture to apply an elastic boundary condition to a test article such that an appropriate dynamic response can be obtained for a given substructuring application, and the technique has been used widely for performing coupling and decoupling operations. Here, in this work, a modified transmission simulator method is proposed that allows a direct exchange of components to be performed. In a numerical demonstration, the proposed method is shown to accurately predict the dynamics of an assembly that has undergone a component exchange in which the exchanged components have different material properties and different geometric properties. This approach can be applied to dynamic qualification frameworks in which structures undergo modifications or upgrades rather than complete redesigns.

Component exchange↗

DOE OSTI2022

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗

DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Interaction of radio frequency waves with cylindrical density filaments: scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas , vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci. , vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

Physics↗

DOE OSTI2026

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Quantifying Impacts of Biomass Pelletization on Fast Pyrolysis Using a Single-Particle Reactor, X-ray Computed Tomography, and Computational Modeling

The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.

09 BIOMASS FUELS↗

DOE OSTI2024

A Reactor Scale-Up Methodology from Lab-Scale to Pilot-Scale Operations: Numerical Modeling of THFA Dehydration to DHP in Packed-Bed Reactors

This manuscript discusses developing a model-based scale-up methodology for a successful technology transfer of gas-phase catalytic reactors from a lab-scale to a pilot-scale operation. The manuscript demonstrates the methodology for gas-phase dehydration of tetrahydrofurfuryl alcohol (THFA) to dihydropyran (DHP) process over commercial Al 2 O 3 catalysts. A two-dimensional reactor model was developed using COMSOL Multiphysics 6.1 software. The model solves heat and mass transport equations in bed-scale and particle scales simultaneously. This powerful feature enables accurate prediction of the heat and mass transfer limitations in pilot-scale reactors, if any exists. Further, the model uses isothermal lab-scale experimental data to derive and validate the reaction chemistry, flow fields and boundary conditions. The model was then scaled-up to project conversion, selectivity, yield and formation rate of DHP in a pilot-scale reactor. The results highlight the complex nature of chemistry, heat, and mass transfer effects in lab-scale and pilot-scale reactors. The model results inform the possible operational limitations of the pilot-scale reactor and design strategies to improve process efficiency. Although the scale-up approach is explained through the THFA dehydration process, the methodology is applicable to any catalytic packed-bed reactor models for a successful process scale-up.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Confining Liquids inside Carbon Nanotubes: Accelerated Molecular Dynamics with Spliced, Soft-Core Potentials and Simulated Annealing

Understanding emergent phenomena of fluids under physical confinement requires the development of advanced tools for rapid and accurate simulation of their physiochemical properties. Simulating liquid molecules commensurate in size with the nanoscale enclosures that confine them is a key challenge. In this work, we demonstrate an accelerated molecular dynamics simulation technique that combines soft-core potentials (SCP) and simulated annealing (SA) to analyze confined liquids. This integrated SCP/SA method relies on a new spliced soft-core potential (SSCP), which enables tunable accuracy with respect to the target hard-core potential (HCP). SCP/SA enables the packing of enclosures with bulk material in a controlled, thermodynamically consistent manner. The enhanced SSCP accuracy is a critical feature of SCP/SA, enabling a smooth transition between the SCP and the HCP at a desired SCP hardness. We applied SCP/SA to the problem of filling a carbon nanotube (CNT) in periodic boundary conditions with a popular ionic liquid (IL), 1- butyl-3-methylimidazolium hexafluorophosphate $[\text{BMIM}^+][\text{PF}^–_6]$. We performed a series of triplicate simulations on systems with varying CNT diameter and charge to demonstrate SCP/SA’s versatility. Beyond this IL/CNT system, the SCP/SA simulation framework has a broad range of potential applications, not limited to nanoscale enclosures and interfaces, including both solid-state and biological systems.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Stochastic Real-Time Second-Order Green’s Function Theory for Neutral Excitations in Molecules and Nanostructures

We present a real-time second-order Green’s function (GF) method for computing excited states in molecules and nanostructures, with a computational scaling of O(N e 3 ), where N e is the number of electrons. The cubic scaling is achieved by adopting the stochastic resolution of the identity to decouple the 4-index electron repulsion integrals. To improve the time propagation and the spectral resolution, we adopt the dynamic mode decomposition technique and assess the accuracy and efficiency of the combined approach for a chain of hydrogen dimer molecules of different lengths. We find that the stochastic implementation accurately reproduces the deterministic results for the electronic dynamics and excitation energies. Furthermore, we provide a detailed analysis of the statistical errors, bias, and long-time extrapolation. Overall, the approach offers an efficient route to investigate excited states in extended systems with open or closed boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

NWPEsSe: an Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au_55, ligated Au_8^(2+), Au_8 supported on graphene oxide and defected rutile, and a large cluster assembly ?[Co?_6 Te_8 (PEt_3 )_6][C_60 ]_n, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems. The work described in this publication was performed at Pacific Northwest National Laboratory (PNNL), which is operated by Battelle for the United States Department of Energy (DOE) under Contract DE-AC05-76RL0180. J. Z. and V.-A. G. acknowledge support from DOE, Office of Science, Office of Basic Energy Sci-ences, Chemical, Geological and Biological Sciences Division and computing resources from PNNL’s Research Computing Facility and the National Energy Research Scientific Computing Center.

Zhang, Jun↗

DOE OSTI2021

Higher Ion Selectivity with Lower Energy Usage Promoted by Electro-osmotic Flow in the Transport through Conical Nanopores

The tradeoff between selectivity and throughput presents fundamental challenges to improve desalination and charge storage, salinity gradient-based energy harvesting, memory device/circuit development, and so forth. The well-known ion-current rectification and the recently resolved time-dependent transport hysteresis in conical nanopores or asymmetric nanointerfaces offer new opportunities for the selective transport of matter. This report shows that electro-osmotic flow (EOF) is an overlooked factor that increases ion selectivity while maintaining enhanced transport throughput in rectified nanoscale electrokinetic transport. The increased selectivity originates primarily from the suppression of anion flux by the fluid flow in the opposite direction under the applied electrical field. By solving the Poisson and Nernst–Planck (PNP) equations without and with coupled Navier–Stokes (PNP–NS) equations, the EOF effects on cation and anion transport are unequivocally revealed in asymmetric nanopipettes. The flux of cations and anions as well as the transference number and flow velocity are elucidated using the models and boundary conditions validated by previous experiments. A dimensionless parameter, radius over the Debye length, reveals optimal ion selectivity and energy cost at intermediate ion concentrations and nanopore sizes, up to hundreds of millimolars and tens of nanometers under time-dependent potential stimulus. Furthermore, the fundamental insights into EOF at nanointerfaces suggest new routes/strategies for better separation, analysis, and energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory

Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe–Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Fractal Model for Drag Reduction on Multiscale Nonwetting Rough Surfaces

Rough surfaces in contact with a flow of fluid exhibit alternating no-slip and free shear boundary conditions at the solid–liquid and air–liquid interfaces, respectively, thereby potentially offering drag reduction benefits. The balance between the dynamic pressure in the flow and the restoring capillary pressure in the interasperity spaces determines the stability of the Cassie state of wettability and is a function of the relative extent of no-slip and free shear regions per unit surface area. In the present study, using a fractal representation of rough surface topography, an analytical model is developed to quantify the stability of the Cassie state of wettability as well as drag reduction and the friction factor for laminar flow in a rectangular channel between nonwetting multiscale rough surfaces. A systematic study is conducted to quantify the effects of fractal parameters of the surfaces and the flow Reynolds number on drag reduction and the friction factor. The studies are used to develop friction factor curves extending the classical Moody diagram to hydrophobic and superhydrophobic surfaces. On the basis of the studies, regime maps are derived for estimating the extent of drag reduction offered by hydrophobic and superhydrophobic surfaces, revealing that superhydrophobic surfaces do not always offer the best drag reduction performance. The application of the fractal model to practical topographies of nonwetting surfaces of copper, aluminum, and zinc oxide fabricated via electrodeposition and etching is also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Design and Preparation of Nematic Colloidal Particles

Colloidal systems are abundant in technology, in biomedical settings, and in our daily life. The so-called “colloidal atoms” paradigm exploits interparticle interactions to self-assemble colloidal analogs of atomic and molecular crystals, liquid crystal glasses, and other types of condensed matter from nanometer- or micrometer-sized colloidal building blocks. Nematic colloids, which comprise colloidal particles dispersed within an anisotropic nematic fluid host medium, provide a particularly rich variety of physical behaviors at the mesoscale, not only matching but even exceeding the diversity of structural and phase behavior in conventional atomic and molecular systems. This feature article, using primarily examples of works from our own group, highlights recent developments in the design, fabrication, and self-assembly of nematic colloidal particles, including the capabilities of preprogramming their behavior by controlling the particle’s surface boundary conditions for liquid crystal molecules at the colloidal surfaces as well as by defining the shape and topology of the colloidal particles. Here recent progress in defining particle-induced defects, elastic multipoles, self-assembly, and dynamics is discussed along with open issues and challenges within this research field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Thermodynamics of Light-Induced Nanoscale Polar Structures in Ferroelectric Superlattices

Here we study the thermodynamics of nanoscale polar structures in PbTiO 3 /SrTiO 3 ferroelectric superlattices induced by above-bandgap optical excitation using a phase-field model explicitly considering both structural and electronic processes. We demonstrate that the light-excited carriers provide the charge compensation of polarization bound charges and the lattice thermal energy, both of which are key to the thermodynamic stabilization of a previously observed supercrystal, a three dimensionally periodic nanostructure, within a window of substrate strains, while different mechanical and electrical boundary conditions can stabilize a number of other nanoscale polar structures by balancing the competing short-range exchange interactions responsible for the domain wall energy and longrange electrostatic and elastic interactions. The insights into the light-induced formation and richness of nanoscale structures from this work offer theoretical guidance for exploring and manipulating the thermodynamic stability of nanoscale polar structures employing a combination of thermal, mechanical, and electrical stimuli as well as light.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Engineering Small HOMO–LUMO Gaps in Polycyclic Aromatic Hydrocarbons with Topologically Protected States

Topological phases in laterally confined low-dimensional nanographenes have emerged as versatile design tools that can imbue otherwise unremarkable materials with exotic band structures ranging from topological semiconductors and quantum dots to intrinsically metallic bands. The periodic boundary conditions that define the topology of a given lattice have thus far prevented the translation of this technology to the quasi-zero-dimensional (0D) domain of small molecular structures. Here, we describe the synthesis of a polycyclic aromatic hydrocarbon (PAH) featuring two localized zero modes (ZMs) formed by the topological junction interface between a trivial and nontrivial phase within a single molecule. First-principles density functional theory calculations predict a strong hybridization between adjacent ZMs that gives rise to an exceptionally small HOMO–LUMO gap. Scanning tunneling microscopy and spectroscopy corroborate the molecular structure of 9/7/9-double quantum dots and reveal an experimental quasiparticle gap of 0.16 eV, corresponding to a carbon-based small molecule long-wavelength infrared (LWIR) absorber.

77 NANOSCIENCE AND NANOTECHNOLOGY↗