Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear structural analysis”

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 163 records · Page 9

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantitative Analysis and Prediction of Thermal Runaway Metrics of High-Nickel Oxide Cathodes by Machine Learning Models

The pursuit of higher energy density in lithium-ion batteries has made high-nickel (Ni) layered oxides leading cathode candidates for next-generation electric vehicles. However, their poor thermal stability, particularly at Ni contents ≥ 90%, increases the risk of cathode-initiated thermal runaway. Furthermore, we present a data-driven framework combining linear and nonlinear machine learning models to predict key thermal runaway descriptors from a high-throughput differential scanning calorimetry database. With cathode composition and state of charge (SOC) as input features, the ensemble model accurately predicts peak temperature, heat release, and peak heat flow. SHAP analysis identifies Ni content and SOC as the dominant factors controlling thermal runaway temperature, while SOC primarily governs heat release and peak heat flow. Al, Mg, and Mn improve thermal stability by strengthening metal–oxygen bonding and delaying structural transformation, whereas B mainly reduces heat release through surface passivation. Validation with a new cathode composition confirms accurate prediction of SOC-dependent thermal runaway behavior and critical SOC.

25 ENERGY STORAGE↗

Probing the Effects of the First Atomic Layer on the Dynamic Behavior of Sub-2 nm MgO/Al 2 O 3 Memristors

As electronic devices continue to scale down from the current sub-5 nm range, atomic-scale control of defects becomes increasingly crucial to suppressing their impact on the physical properties of the devices. Memristors present an excellent example as a nonlinear and dynamic device with high speed and endurance required for electronic applications ranging from neuromorphic computing to nonvolatile memories. Herein we investigate the impact of atomic defects in sub-2 nm thick MgO/Al 2 O 3 atomic layer stack (ALS) memristors that use an M1 (switching layer)/M2 (oxygen vacancy reservoir layer) bilayer structure grown using in vacuo atomic layer deposition (iALD). Intriguingly, we revealed a direct correlation of the atomic defects in the M2 layer with the memristor dynamic behavior using a combined analysis of in situ scanning tunneling spectroscopy (iSTS) on the M2 layer and ex situ characterization on the memristors. Specifically, incomplete coverage of the 1st ALD atomic layer of M2 on the electrode yields defects at the M2/electrode interface. Despite the monotonic increase of ALD coverage, by almost three-fold from ~30% to >90%, at completion of the M2 layer of ~ 0.7 nm in thickness, the impact of the defects on the M2/electrode interface has been found detrimental to both memristor switching speed and endurance. Guided by atomistic simulation, we addressed the issue of interface defects via tuning of the Al surface hydroxylation to increase the first atomic layer ALD coverage to ~75%, leading to improved memristor switching speed and endurance by several orders of magnitude. In conclusion, these findings shed light on the correlation between the atomic defects and the dynamic behavior of sub-2 nm memristors and the importance of minimizing the atomic defects in memristors for future electronic applications.

Atomic Layer Deposition↗

Causal Directions Matter: How Environmental Factors Drive Convective Cloud Detrainment Heights

This study investigates how environmental factors influence the level of maximum detrainment (LMD) in deep convective clouds. Through a novel application of the Linear Non‐Gaussian Acyclic Model (LiNGAM), we discover causal structures between environmental variables and LMD, observed at six tropical sites operated by the Atmospheric Radiation Measurement (ARM) user facility. LiNGAM effectively identifies causal directions among variables of interest, revealing robust relationships such as those among the lifting condensation level (LCL), level of free convection (LFC), and convective inhibition (CIN), aligning with prior knowledge. Relative humidity is shown to directly influence LMD; however, this relationship exhibits strong nonlinearity and becomes difficult to detect when the contrast between oceanic and continental environments is excluded from the analysis. This study highlights the importance of establishing causal relationships before performing statistical inference.

54 ENVIRONMENTAL SCIENCES↗

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗

Tracking the topology of neural manifolds across populations

Neural manifolds summarize the intrinsic structure of the information encoded by a population of neurons. Advances in experimental techniques have made simultaneous recordings from multiple brain regions increasingly commonplace, raising the possibility of studying how these manifolds relate across populations. However, when the manifolds are nonlinear and possibly code for multiple unknown variables, it is challenging to extract robust and falsifiable information about their relationships. We introduce a framework, called the method of analogous cycles, for matching topological features of neural manifolds using only observed dissimilarity matrices within and between neural populations. We demonstrate via analysis of simulations and in vivo experimental data that this method can be used to correctly identify multiple shared circular coordinate systems across both stimuli and inferred neural manifolds. Conversely, the method rejects matching features that are not intrinsic to one of the systems. Further, as this method is deterministic and does not rely on dimensionality reduction or optimization methods, it is amenable to direct mathematical investigation and interpretation in terms of the underlying neural activity. We thus propose the method of analogous cycles as a suitable foundation for a theory of cross-population analysis via neural manifolds.

97 MATHEMATICS AND COMPUTING↗

Stabilization of the Polar Structure and Giant Second-Order Nonlinear Response of Single Crystal γ-NaAs 0.95 Sb 0.05 Se 2

The dearth of suitable materials significantly restricts the practical development of infrared (IR) laser systems with highly efficient and broadband tuning. Recently, γ-NaAsSe 2 is reported, and it exhibits a large nonlinear second-harmonic generation (SHG) coefficient of 590 pm V —1 at 2 µm. However, the crystal growth of γ-NaAsSe 2 is challenging because it undergoes a phase transition to centrosymmetric δ-NaAsSe 2 . Herein, the stabilization of non-centrosymmetric γ-NaAsSe 2 by doping the As site with Sb, which results in γ-NaAs 0.95 Sb 0.05 Se 2 is reported. The congruent melting behavior is confirmed by differential thermal analysis with a melting temperature of 450 °C and crystallization temperature of 415 °C. Single crystals with dimensions of 3 mm × 2 mm are successfully obtained via zone refining and the Bridgman method. The purification of the material plays a significant role in crystal growth and results in a bandgap of 1.78 eV and thermal conductivity of 0.79 Wm —1 K —1 . The single-crystal SHG coefficient of γ-NaAs 0.95 Sb 0.05 Se 2 exhibits an enormous value of |d 11 | = 648 ± 74 pm V —1 , which is comparable to that of γ-NaAsSe 2 and ≈20× larger than that of AgGaSe 2 . Importantly, the bandgap of γ-NaAs 0.95 Sb 0.05 Se 2 (1.78 eV) is similar to that of AgGaSe 2 , thus rendering it highly attractive as a high-performing nonlinear optical material.

36 MATERIALS SCIENCE↗

Resistivity phase diagram of cuprates revisited

The phase diagram of the cuprate superconductors has posed a formidable scientific challenge for more than three decades. This challenge is perhaps best exemplified by the need to understand the normal-state charge transport as the system evolves from Mott insulator to Fermi- liquid metal with doping. Here in this paper, we report a detailed analysis of the temperature (T) and doping (p) dependence of the planar resistivity of simple-tetragonal HgBa 2 CuO 4+δ (Hg1201), the single- CuO 2 -layer cuprate with the highest optimal superconducting transition temperature, T c . The data allow us to test a recently proposed phenomenological model for the cuprate phase diagram that combines a universal transport scattering rate with spatially inhomogeneous (de)localization of the Mott-localized hole. We find that the model provides a good description of the data. We then extend this analysis to prior transport results for several other cuprates, including the Hall number in the overdoped part of the phase diagram, and find little compound-to-compound variation in (de)localization gap scale. The results point to a robust, universal structural origin of the inherent gap inhomogeneity that is unrelated to doping-related disorder. They are inconsistent with the notion that much of the phase diagram is controlled by a quantum critical point, and instead indicate that the unusual electronic properties exhibited by the cuprates are fundamentally related to strong nonlinearities associated with subtle nanoscale inhomogeneity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Develop and verify soil/structure interaction for pile/foundation interaction

Phase II of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was used to verify the implementation of a new soil-structure interaction (SSI) model for use within offshore wind turbine modeling software. The REDWIN Macro-element model implemented and verified in this study enables a computationally efficient way to model the linear and nonlinear SSI problem, including hysteretic damping, of a monopile structure. The modeling approach was integrated into several modeling tools and a series of increasingly complex simulations was conducted using the IEA 10MW reference turbine mounted on a monopile support structure to verify the coupling between the tools and the REDWIN Macro-element SSI model. This campaign includes only numerical verification between various software and modeling approaches so no experimental measurements are available. The load cases (LC) considered include: LC1 – static response of the tower and substructure LC2 – frequency and mode-shape analysis of the tower and substructure LC3 – response of the tower and substructure due to wind-only loading LC4 – response of the tower and substructure due to wave-only loading LC5 – response of the tower and substructure due to wind and wave loading. Detailed properties of the modeled system are found in the following reference, “Bergua, Roger, Amy Robertson, Jason Jonkman, and Andy Platt. 2021. "Specification Document for OC6 Phase II: Verification of an Advanced Soil-Structure Interaction Model for Offshore Wind Turbines.” Golden, CO: National Renewable Energy Laboratory. NREL/TP-5000-79938. https://www.nlr.gov/docs/fy21osti/79938.pdf. Details on the results from the OC6 Phase II project can be found in the following reference, “Bergua R, Robertson A, Jonkman J, et al. OC6 Phase II: Integration and verification of a new soil–structure interaction model for offshore wind design.” Wind Energy. 2022;25(5):793-810. doi:10.1002/we.2698

17 WIND ENERGY↗

Revisiting two thiophosphate compounds constituting d 0 transition metal HfP 2 S 6 and d 10 transition metal α-Ag 4 P 2 S 6 as multifunctional materials for combining second harmonic generation response and photocurrent response

Two acentric thiophosphate compounds, HfP 2 S 6 and α-Ag 4 P 2 S 6 , are revisited and studied as infrared nonlinear optical materials. HfP 2 S 6 and α-Ag 4 P 2 S 6 were structurally characterized without any property measurements. Here, in this work, HfP 2 S 6 and α-Ag 4 P 2 S 6 were synthesized via high temperature salt flux reactions. Low-temperature polymorph acentric α-Ag 4 P 2 S 6 was purified and grown as mm-sized crystals with the aid of AgBr flux. The AgBr flux was revealed to play an important role in stabilizing the acentric α-Ag 4 P 2 S 6 . The acentric α-Ag 4 P 2 S 6 transforms to centrosymmetric β-Ag 4 P 2 S 6 at 850(5) K, which is revealed by differential scanning calorimetry (DSC) analysis and powder X-ray diffraction experiments. HfP 2 S 6 and α-Ag 4 P 2 S 6 are discovered by UV-vis spectrum measurements as indirect bandgap semiconductors with bandgaps of 2.2(1) eV and 2.5(1) eV, respectively, which is supported by DFT calculations and TB-LMTO-ASA calculations. The bonding pictures of α-Ag 4 P 2 S 6 were studied by crystal orbital Hamilton population calculations (COHP) coupled with electron localization function (ELF) analysis. DFT calculations predict that HfP 2 S 6 and α-Ag 4 P 2 S 6 would exhibit different optical performances regardless of being constructed from identical [P 2 S 6 ] motifs. HfP 2 S 6 exhibits a low second harmonic generation (SHG) response, ~0.21 × AGS (for the sample of particle size of 25 μm). α-Ag 4 P 2 S 6 possesses moderate SHG response, ~0.61 × AGS (for the sample of particle size of 225 μm) coupled with a high laser damage threshold (LDT) of ~3.2 × AGS. Characteristics of high ambient stability, moderate bandgap and SHG response, type-I phase-matching capability, and high LDT together with the easy growth of large crystals make α-Ag 4 P 2 S 6 attractive for future infrared nonlinear optical applications. Photocurrent measurements found that α-Ag 4 P 2 S 6 and β-Ag 4 P 2 S 6 have high photocurrent response, 165 nA cm -2 and 135 nA cm -2 , respectively. α-Ag 4 P 2 S 6 is a new multifunctional material of the ternary Ag–P–S system, which combines nonlinear optical (NLO) properties and photocurrent response.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis, Crystal and Electronic Structures, Nonlinear Optical Properties, and Magnetic Properties of Two Thiophosphates: KInP 2 S 7 and KCrP 2 S 7

Two thiophosphates, KInP 2 S 7 and KCrP 2 S 7 , were structurally characterized without investigating any optical properties. Herein in this work, KInP 2 S 7 and KCrP 2 S 7 were revisited to investigate their optical and magnetic properties, respectively. Pure polycrystalline samples and crystals of KInP 2 S 7 and KCrP 2 S 7 were grown by high temperature solid state reactions, where mm-sized crystals of KCrP 2 S 7 were collected. KCrP 2 S 7 is isostructural to KInP 2 S 7 , which features a layered structure. KInP 2 S 7 and KCrP 2 S 7 possess close relationship to the layered thiophosphate M 2 P 2 S 6 (M = Fe, Co, Zn, etc.). The bonding pictures of KInP 2 S 7 were studied using the electron localization function (ELF) coupled with crystal orbital Hamilton population (COHP) calculations. The intrinsically distorted [PS 4 ] tetrahedra and [InS 6 ] octahedra are made by strong covalent P-S interactions and ionic In-S interactions, respectively. Electronic structure analysis confirmed that the optical properties of KInP 2 S 7 are mainly contributed to by [PS 4 ] tetrahedra together with small amounts of the contributions coming from [InS 6 ] octahedra. Magnetic measurement on mm-sized crystals of KCrP 2 S 7 verified that there is an antiferromagnetic transition around 21 K, and the Cr atoms are trivalent. KInP 2 S 7 is predicated to be an indirect bandgap semiconductor of 2.38 eV, which is confirmed by the UV-Vis measurement of 2.4(1) eV. KInP 2 S 7 is not a type-I phase-matching material and exhibits moderate second harmonic generation (SHG) response (0.51 × AgGaS 2 , sample of particle size of 100 µm). The laser damage threshold (LDT) of KInP 2 S 7 is very high of 5.2 × AgGaS 2 . Bandgap engineering were undergone to enhance the SHG response of KInP 2 S 7 .

36 MATERIALS SCIENCE↗

Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation

Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Design and Analysis of High Heat Flux Plasma-Facing Components for NSTX Upgrade

The National Spherical Tokamak Experiment (NSTX) upgrade (NSTX-U) requirements lead to enhanced heat loads on plasma-facing components (PFCs) especially in the divertor regions, where normal heat flux density can reach up to 8 MW/m 2 . For these high heat flux (HHF) regions, the PFCs were redesigned, to use castellations which relieved the thermal stresses associated with high incident heat flux. Another design feature of HHF PFCs is the absence of front face mounting holes which create localized areas of high temperature and thermal stress concentrations. Optimized shaping of the front face of the HHF PFCs eliminates regions where the front face is perpendicular to the heat flux direction during normal helicity operation and spreads the heat load. A unique mechanism was designed to connect PFC tiles to the NSTX-U center stack casing using locking features accessible from the low heat flux regions. Isotropic graphite was selected as the HHF PFC tile material. Initial tile design parameters were accessed using analytical expressions for pulsed heat flux loading. A working prototype of the locking mechanism was created during the initial stages of the design to prove the concept performance. This article presents an overview of the divertor HHF PFC design and the results of the thermal and structural analyses performed using ANSYS software. The results of the analyses cover normal operating conditions and disruptions which impose electromagnetic (EM) loads from eddy and halo currents. The 3-D, transient, nonlinear analyses took into account the temperature-dependent properties of the materials, friction interfaces between the parts, and variable electric properties of the parts and interfaces. The results confirm that the HHF PFC tiles remain within the allowable limits for the loads defined by the NSTX-U Recovery Project. In addition, tolerance stack up analyses were performed to ensure tile performance in the worst possible assembly configuration. Finally, the design process was completed successfully, and the NSTX HHF PFC tiles are currently in the production phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Variationally consistent Maxwell stress in flexoelectric structures under finite deformation and immersed in free space

Maxwell stress refers to the mechanical stress exerted on a dielectric material due to the presence of electric fields. It plays a significant role in the interaction between a dielectric material and the surrounding free space under finite deformation. Previous research on finite deformation of flexoelectricity mainly adopted a modified form of Maxwell stress, potentially not able to correctly capture some physical phenomena, such as the compression of a dielectric droplet in an electric field. In this work, we propose a consistent and complete variational principle for flexoelectricity, in which the Maxwell stress emerges naturally from the derivation, without introducing additional assumptions. An Isogeometric analysis-based numerical framework is developed accordingly and verified by both linear and nonlinear benchmark cases compared with experimental results. The present framework successfully captures and quantifies the behaviors of conductive liquids and soft dielectric solids subjected to an external electric field. Finally, a novel scenario is investigated in which a flexoelectric beam immersed in free space is analyzed, showing the interesting distribution of Maxwell stress-induced tractions at opposing boundaries. The test demonstrates that a higher dielectric constant can effectively enhance the material's stiffness in response to the external electric loading.

36 MATERIALS SCIENCE↗

Implicit Multirate GARK Methods

This article considers multirate generalized-structure additively partitioned Runge–Kutta methods for solving stiff systems of ordinary differential equations with multiple time scales. These methods treat different partitions of the system with different timesteps for a more targeted and efficient solution compared to monolithic single rate approaches. With implicit methods used across all partitions, methods must find a balance between stability and the cost of solving nonlinear equations for the stages. In order to characterize this important trade-off, we explore multirate coupling strategies, problems for assessing linear stability, and techniques to efficiently implement Newton iterations for stage equations. Unlike much of the existing multirate stability analysis which is limited in scope to particular methods, we present general statements on stability and describe fundamental limitations for certain types of multirate schemes. New implicit multirate methods up to fourth order are derived, and their accuracy and efficiency properties are verified with numerical tests.

97 MATHEMATICS AND COMPUTING↗

Insights on the Effects of Magnetic Forces on the Efficiency of Vibration Energy Harvesting Absorbers in Controlling Dynamical Systems

This study investigates the effects of magnetic constraints on a piezoelectric energy harvesting absorber while simultaneously controlling a primary structure and harnessing energy. An accurate forcing representation of the magnetic force is investigated and developed. A reduced-order model is derived using the Euler–Lagrange principle, and the impact of the magnetic force is evaluated on the absorber’s static position and coupled natural frequency of the energy harvesting absorber and the coupled primary absorber system. The results show that attractive magnet configurations cannot improve the system substantially before pull-in occurs. A rigorous eigenvalue problem analysis is performed on the absorber’s substrate thickness and tip mass to effectively design an energy harvesting absorber for multiple initial gap sizes for the repulsive configurations. Then, the effects of the forcing amplitude on the primary structure absorber are studied and characterized by determining an effective design of the system for a simultaneous reduction in the primary structure’s motion and improvement in the harvester’s efficiency.

42 ENGINEERING↗

Structural color from a coupled nanowire pair beyond the bonding and antibonding model

Optical resonances in nanostructures can be harnessed to produce a wide range of structural colors. Conversely, the analysis of structural colors has been used to clarify the nature of optical resonances. Here, we show that silicon nanowire (NW) pairs can display a wide range of structural colors by controlling their radiative coupling. This is accomplished by exciting a series of Fabry–Pérot-like modes where light is repeatedly scattered between two NWs. These modes are beyond the expectation from the conventional chemical bonding model under a quasi-electrostatic approximation, in which only bonding and antibonding modes can be formed in a pair system through modal hybridization. The additional eigenmodes found in a two-resonator system originate from the nonlinear, frequency-dependent coupling strength derived from the radiative nature of low- Q resonators. The Fabry–Pérot modes can be tuned across the entire visible frequency range by varying the distance between two NWs, leading to what we believe is a new type of universal building blocks that can provide structural color within a subwavelength footprint. The presented results pave the way toward the design and usage of highly tunable resonances that exploit the radiative coupling of high-index nanostructures.

Li, Qitong↗