Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase shift”

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 973 records · Page 54

Microscopic model for the stacking-fault potential and the exciton wave function in GaAs

Two-dimensional stacking fault defects embedded in a bulk crystal can provide a homogeneous trapping potential for carriers and excitons. Here we utilize state-of-the-art structural imaging cou- pled with density functional and effective-mass theory to build a microscopic model of the stacking- fault exciton. The diamagnetic shift and exciton dipole moment at different magnetic fields are calculated and compared with the experimental photoluminescence of excitons bound to a single stacking fault in GaAs. The model is used to further provide insight into the properties of excitons bound to the double-well potential formed by stacking fault pairs. This microscopic exciton model can be used as an input into models which include exciton-exciton interactions to determine the excitonic phases accessible in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An active, stable cubic molybdenum carbide catalyst for the high-temperature reverse water-gas shift reaction

Although technologically promising, the reduction of carbon dioxide (CO 2 ) to produce carbon monoxide (CO) remains economically challenging owing to the lack of an inexpensive, active, highly selective, and stable catalyst. We show that nanocrystalline cubic molybdenum carbide (α-Mo 2 C), prepared through a facile and scalable route, offers 100% selectivity for CO 2 reduction to CO while maintaining its initial equilibrium conversion at high space velocity after more than 500 hours of exposure to harsh reaction conditions at 600°C. The combination of operando and postreaction characterization of the catalyst revealed that its high activity, selectivity, and stability are attributable to crystallographic phase purity, weak CO-Mo 2 C interactions, and interstitial oxygen atoms, respectively. Mechanistic studies and density functional theory (DFT) calculations provided evidence that the reaction proceeds through an H 2 -aided redox mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental investigation of beam-ion losses induced by magnetic perturbations using the light ion beam probe technique in the ASDEX Upgrade tokamak

The impact of externally applied magnetic perturbations (MPs) on fast-ion losses has been investigated by means of the light ion beam probe (LIBP) technique in the ASDEX Upgrade tokamak. The LIBP technique allows to experimentally infer the fast-ion orbit displacement induced by MPs via first-orbit losses using scintillator based fast-ion loss detector (FILD) measurements. The fast-ion orbit displacement against different applied MP spectra has been studied. These shots were conducted in ELM mitigated H-mode plasmas. A rigid rotation of the MP coils was applied with a frequency of 1 Hz, with an n = 2 configuration and changing the differential phase between the upper and lower set of coils (ΔΦ ul ) on a shot-to-shot basis. Beam sources Q7 (tangential) and Q8 (radial) were used to probe different fast-ion orbits with FILD1. The measured fast-ion orbit displacement ranges from 3 to 20 mm approximately, and no qualitative difference is observed between ions from beam sources Q7 and Q8. The minimum is found for a ΔΦ ul ~ 50°, which is shifted with respect to the minimum of the plasma boundary displacement, found at ΔΦ ul ~ 0°. A first attempt to validate the orbit following code ASCOT—including the plasma response calculated with the MARS-F code—against these experimental measurements is performed. Finally, while the dependence of the first-orbit fast-ion displacement with ΔΦ ul does not match the experimental measurements, these simulations do capture other features such as the order of magnitude of the orbit displacement and the importance of the toroidal spectrum of the applied perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a Multi-Robot System for Autonomous Inspection of Nuclear Waste Tank Pits

This paper introduces the overall design plan, development timeline, and preliminary progress of the Autonomous Pit Exploration System project. This project aims to develop an advanced multi-robot system for the efficient inspection of nuclear waste-storage tank pits. The project is structured into three phases: Phase 1 involves data collection and interface definition in collaboration with Hanford Site experts and university partners, focusing on tank riser geometry and hardware solutions. Phase 2 includes the selection of sensors and robot components, detailed mechanical design, and prototyping. Phase 3 integrates all components into a cohesive system managed by a master control package which also incorporates digital twin and surrogate models, and culminates in comprehensive testing and validation at a simulated tank pit at the Idaho National Laboratory. Additionally, the system’s communication design ensures coordinated operation through shared data, power, and control signals. For transportation and deployment, an electric vehicle (EV) is chosen to support the system for a full 10 h shift with better regulatory compliance for field deployment. A telescopic arm design is selected for its simple configuration and superior reach capability and controllability. Preliminary testing utilizes an educational robot to demonstrate the feasibility of splitting computational tasks between edge and cloud computers. Successful simultaneous localization and mapping (SLAM) tasks validate our distributed computing approach. More design considerations are also discussed, including radiation hardness assurance, SLAM performance, software transferability, and digital twinning strategies.

Nuclear waste management↗

Tackling Disorder in γ‐Ga 2 O 3

Abstract Ga 2 O 3 and its polymorphs are attracting increasing attention. The rich structural space of polymorphic oxide systems such as Ga 2 O 3 offers potential for electronic structure engineering, which is of particular interest for a range of applications, such as power electronics. γ‐Ga 2 O 3 presents a particular challenge across synthesis, characterization, and theory due to its inherent disorder and resulting complex structure–electronic‐structure relationship. Here, density functional theory is used in combination with a machine‐learning approach to screen nearly one million potential structures, thereby developing a robust atomistic model of the γ‐phase. Theoretical results are compared with surface and bulk sensitive soft and hard X‐ray photoelectron spectroscopy, X‐ray absorption spectroscopy, spectroscopic ellipsometry, and photoluminescence excitation spectroscopy experiments representative of the occupied and unoccupied states of γ‐Ga 2 O 3 . The first onset of strong absorption at room temperature is found at 5.1 eV from spectroscopic ellipsometry, which agrees well with the excitation maximum at 5.17 eV obtained by photoluminescence excitation spectroscopy, where the latter shifts to 5.33 eV at 5 K. This work presents a leap forward in the treatment of complex, disordered oxides and is a crucial step toward exploring how their electronic structure can be understood in terms of local coordination and overall structure.

36 MATERIALS SCIENCE↗

Crystal Prediction and Design of Tunable Light Emission in BTB–Based Metal–Organic Frameworks

Metal-organic frameworks (MOFs) have recently been shown to exhibit unique mechanisms of luminescence based on charge transfer between structural units in the framework. These MOFs have the potential to be structural tuned for targeted emission with little or no metal participation. Here a computationally led, material design and synthesis methodology is presented that elucidates the mechanisms of light emission in interpenetrated structures comprised of metal centers (M = In, Ga, InGa, InEu) and BTB (1,3,5-Tris(4-carboxyphenyl)benzene) linkers, forming unique luminescent M-BTB MOF frameworks. Gas phase and periodic electronic structure calculations indicate that the intensity of the emission and the wavelength are overwhelmingly controlled by a combination of the number of interacting stacked linkers and their interatomic spacings, respectively. In the MOF, the ionic radii of the metal centers primarily control the expansion or shrinkage of the linker stacking distances. Experimentally, multiple M-BTB-based MOFs are synthesized and their photoluminescence was tested. Experiments validated the modeling by confirming that shifts in the crystal structure result in variations in light emission. Through this material design method, the mechanisms of tuning luminescence properties in interpenetrated M-BTB MOFs have been identified and applied to the design of MOFs with specific wavelength emission based on their structure.

36 MATERIALS SCIENCE↗

Additive manufacturing of a metastable high entropy alloy: Metastability engineered microstructural control via process variable driven elemental segregation

Compositional paradigm shift in high entropy alloys (HEAs) provided new opportunities for microstructural engineering, whereas process control during laser powder bed fusion (LPBF) additive manufacturing (AM) enable fine microstructure tailoring. Metastability engineering in transformation induced plasticity (TRIP) HEAs by addition of minor alloying elements is an attractive strategy for fine microstructural tuning. This study explored in detail the microstructural evolution during LPBF AM of a metastable Fe 40 Mn 20 Co 20 Cr 15 Si 5 (at.%) dual phase HEA (CS-HEA). LPBF processing window for CS-HEA was established based on quantitative analysis and experiments. Based on melt pool overlap lack of fusion pores were observed at lower energy densities (J) of $\textit{J}$ ≤ 31.25 J/mm 3 and key-hole formation by melt pool destabilization in case of $\textit{J}$ ≥ 75 J/mm 3 . The microstructure of CS-HEA consists of metastable FCC-γ and HCP-ε phases; LPBF process parameters governed the final phase fraction in the alloy which has been correlated to metastability alteration of the high temperature γ phase. Final microstructural engineering was devised by LPBF process control which enabled cooling rate manipulation to guide Mn and Si segregation at the cell boundaries, thereby controlling the matrix metastability and final phase fraction. Additionally, high resolution transmission electron microscopy (TEM) revealed disparity in stacking fault morphology in CS-HEA with LPBF process variable alterations associated with local variations in chemical composition and stacking fault energy. Further, the phase evolution with process parameters also affected the nanomechanical behavior of the alloy.

36 MATERIALS SCIENCE↗

What controls the interannual variation of Hadley cell extent in the Northern Hemisphere: physical mechanism and empirical model for edge variation

Abstract The Hadley circulation is the most prominent atmospheric meridional circulation, reducing the radiatively driven equator-to-pole temperature gradient. While the Hadley cell extent varies by several degrees from year to year, the detailed dynamical mechanisms behind such variations have not been well elucidated. During the expanded phase of the Hadley cell, many regions on the periphery of the subtropics experience unfavorable climatic conditions. In this study, using ERA5 reanalysis data, we examine the physical chain of events responsible for the interannual variation of the Hadley cell edge (HCE) latitude in the Northern Hemisphere. This variation is mainly caused by changing eddy activity and wave breaking from both stationary and transient waves. In particular, we show that transient waves cause the HCE to shift poleward by increasing the eddy momentum flux divergence (EMFD) and reducing the baroclinicity over 20°–40°N, shifting the region of peak baroclinicity poleward. El Niño/La Niña and the Arctic Oscillation (AO) account for a significant portion (60%) of the interannual fluctuation of the HCE latitude. Through the poleward displacement of eddy activity, La Niña and a positive AO state are associated with the poleward shift of the HCE. The analysis of 28 CMIP5 models reveals statistical relationships between EMFD, vertical shear, and HCE latitude similar to those observed.

54 ENVIRONMENTAL SCIENCES↗

Spectroscopic evidence of intra-unit-cell charge redistribution in a charge-neutral magnetic topological insulator

The magnetic topological insulator MnBi 6 Te 10 has emerged as a promising candidate for realizing the quantum anomalous Hall effect (QAHE), owing to its ability to retain ferromagnetism through precise control of anti-site defects. The next important task for realizing the QAHE is to tune the chemical potential into the energy gap formed by the broken time-reversal symmetry. Here we reveal an intra-unit-cell charge redistribution even when the overall doping suggests a near-charge-neutral condition. By performing time- and angle-resolved photoemission spectroscopy (trARPES) on the optimally 18% Sb-doped MnBi 6 Te 10 , we observe transient surface photovoltage (SPV) effects on both the MnBi 2 Te 4 and single-Bi 2 Te 3 terminations. Furthermore, we observe a time-dependent splitting of the band structure indicating multiple SPV shifts with different magnitudes. This observation suggests that adjacent plateaus with nominally the same terminating layer exhibit a strong intra-unit-cell charge redistribution, resulting in spontaneous electrical polarization. This is consistent with static micro-ARPES measurements revealing significant doping deviations from the charge-neutral configuration. Our findings underscore the challenges of engineering the family of Mn-Bi-Te materials to realize QAHE purely through chemical doping. Achieving the desired topological quantum phase requires both a uniform carrier doping and a ferromagnetic ground state. Furthermore, the light-induced polarization within each unit cell of ferromagnetic Mn(Bi 0.82 Sb 0.18 ) 6 Te 10 may open new possibilities for optoelectronic and spintronics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Twist Angle-Dependent Interlayer Exciton Lifetimes in van der Waals Heterostructures

In van der Waals (vdW) heterostructures formed by stacking two monolayers of transition metal dichalcogenides, multiple exciton resonances with highly tunable properties are formed and subject to both vertical and lateral confinement. We investigate how a unique control knob, the twist angle between the two monolayers, can be used to control the exciton dynamics. We observe that the interlayer exciton lifetimes in MoSe 2 /WSe 2 twisted bilayers (TBLs) change by one order of magnitude when the twist angle is varied from 1° to 3.5°. Additionally, using a low-energy continuum model, we theoretically separate two leading mechanisms that influence interlayer exciton radiative lifetimes. The shift to indirect transitions in the momentum space with an increasing twist angle and the energy modulation from the moiré potential both have a significant impact on interlayer exciton lifetimes. We further predict distinct temperature dependence of interlayer exciton lifetimes in TBLs with different twist angles, which is partially validated by experiments. While many recent studies have highlighted how the twist angle in a vdW TBL can be used to engineer the ground states and quantum phases due to many-body interaction, our studies explore its role in controlling the dynamics of optically excited states, thus, expanding the conceptual applications of “twistronics”.

36 MATERIALS SCIENCE↗

Data Integration and Visualization for Enhanced Resilience and Sustainability in Hydropower (DIVERS-H)

U.S. hydropower plants face potential threats from shrinking water supply, rising demands, and warmer stream temperatures from various causes. Power plant owners, operators, and regulators require new tools to take advantage of and interpret the diverse range of scientific data being produced by both observational methods (for example, satellite, radar, stream gauges) and computer modeling methods that evaluate and predict how earth's dynamic systems (atmosphere, oceans, land surface, and sea ice) are changing and interacting. Combining datasets such as these with AI-based analyses introduces a novel decision support system to help users anticipate and address potential impacts on power generation stations. This new technology has been named DIVERS-H for "Data Integration and Visualization for Enhanced Resilience and Sustainability in Hydropower." In Phase I, technical feasibility was established with the development and demonstration of all the new technologies that are required. Most notably, DIVERS-H will use new artificial intelligence (AI) methods to capture the complex dynamics of water availability, demand, and environmental changes. In addition, new data management software was developed, and a prototype user interface was implemented as the precursor to a full scale decision support system. With technical research complete, the project focus now shifts to development of a commercial software product to provide users with actionable insight into water availability and the risk/resilience of critical systems at their locations of interest. Although DIVER-H was originally conceived as a tool for hydroelectric power applications, the same underlying technology can be readily applied to other water-consuming systems including coal, natural gas, oil, and nuclear power plants.

Chaudhary, Aashish [Kitware, Inc., Clifton Park, N↗

An inconsistency in aviation emissions between CMIP5 and CMIP6 and the implications for short-lived species and their radiative forcing

We report on an inconsistency in the latitudinal distribution of aviation emissions between the data products of phases 5 and 6 of the Coupled Model Intercomparison Project (CMIP). Emissions in the CMIP6 data occur at higher latitudes than in the CMIP5 data for all scenarios, years, and emitted species. A comparative simulation with the chemistry–climate model ECHAM/MESSy Atmospheric Chemistry (EMAC) reveals that the difference in nitrogen oxide emission distribution leads to reduced overall ozone changes due to aviation in the CMIP6 scenarios because in those scenarios the distribution of emissions is partly shifted towards the chemically less active higher latitudes. The radiative forcing associated with aviation ozone is 7.6 % higher, and the decrease in methane lifetime is 5.7 % larger for the year 2015 when using the CMIP5 latitudinal distribution of emissions compared to when using the CMIP6 distribution. We do not find a statistically significant difference in the radiative forcing associated with aviation aerosol emissions. In total, future studies investigating the effects of aviation emissions on ozone and climate should consider the inconsistency reported here.

54 ENVIRONMENTAL SCIENCES↗

Sources of Subseasonal-To-Seasonal Predictability of Atmospheric Rivers and Precipitation in the Western United States

Atmospheric rivers (ARs) account for a large portion of winter precipitation in the western US. To evaluate the sources of AR and precipitation predictability at subseasonal-to-seasonal timescales, we examine the relationships between two climate modes, the El Niño/Southern Oscillation (ENSO) and Madden-Julian oscillation (MJO), and winter hydroclimate in the western US. Our analysis uses a large ensemble of the Weather Research and Forecast (WRF) model simulations from 1981 to 2017, facilitating the assessment of uncertainty in climate mode-AR relationship due to climate variability and short length of observations. Over the North Pacific basin, we find ENSO-related latitudinal shifts of ARs, but ENSO has little effect on basin-wide averaged AR frequency. Over the western US, there is some uncertainty in the ENSO-AR connection due to the impacts of ENSO indices and datasets. However, extreme El Niño events defined by the ENSO longitude index are consistently linked to increased landfalling AR activity. The MJO can lead to significantly enhanced or suppressed landfalling AR activity in the large-ensembles, depending on the phase of MJO and time lag, while observations are too short in length to robustly show this signal. ENSO substantially modulates the MJO-AR relationship, triggering variegated responses of landfalling ARs and AR precipitation in La Niña and El Niño years. Our findings highlight the need to evaluate concurrent effects of different climate modes on ARs and precipitation, and may shed light on a path toward more accurate subseasonal-to-seasonal prediction of ARs and precipitation over the western US.

54 ENVIRONMENTAL SCIENCES↗

An industrial policy framework for transforming energy and emissions intensive industries towards zero emissions

The target of zero emissions sets a new standard for industry and industrial policy. Industrial policy in the twenty-first century must aim to achieve zero emissions in the energy and emissions intensive industries. Sectors such as steel, cement, and chemicals have so far largely been sheltered from the effects of climate policy. A major shift is needed, from contemporary industrial policy that mainly protects industry to policy strategies that transform the industry. For this purpose, we draw on a wide range of literatures including engineering, economics, policy, governance, and innovation studies to propose a comprehensive industrial policy framework. The policy framework relies on six pillars: directionality, knowledge creation and innovation, creating and reshaping markets, building capacity for governance and change, international coherence, and sensitivity to socio-economic implications of phase-outs. Complementary solutions relying on technological, organizational, and behavioural change must be pursued in parallel and throughout whole value chains. Current policy is limited to supporting mainly some options, e.g. energy efficiency and recycling, with some regions also adopting carbon pricing, although most often exempting the energy and emissions intensive industries. An extended range of options, such as demand management, materials efficiency, and electrification, must also be pursued to reach zero emissions. New policy research and evaluation approaches are needed to support and assess progress as these industries have hitherto largely been overlooked in domestic climate policy as well as international negotiations.

54 ENVIRONMENTAL SCIENCES↗

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Variational and phase response analysis for limit cycles with hard boundaries, with applications to neuromechanical control problems

Motor systems show an overall robustness, but because they are highly nonlinear, understanding how they achieve robustness is difficult. In many rhythmic systems, robustness against perturbations involves response of both the shape and the timing of the trajectory. This makes the study of robustness even more challenging. To understand how a motor system produces robust behaviors in a variable environment, we consider a neuromechanical model of motor patterns in the feeding apparatus of the marine mollusk Aplysia californica (Shaw et al. in J Comput Neurosci 38(1):25–51, 2015; Lyttle et al. in Biol Cybern 111(1):25–47, 2017). We established in (Wang et al. in SIAM J Appl Dyn Syst 20(2):701–744, 2021. https://doi.org/10.1137/20M1344974) the tools for studying combined shape and timing responses of limit cycle systems under sustained perturbations and here apply them to study robustness of the neuromechanical model against increased mechanical load during swallowing. Interestingly, we discover that nonlinear biomechanical properties confer resilience by immediately increasing resistance to applied loads. In contrast, the effect of changed sensory feedback signal is significantly delayed by the firing rates’ hard boundary properties. Our analysis suggests that sensory feedback contributes to robustness in swallowing primarily by shifting the timing of neural activation involved in the power stroke of the motor cycle (retraction). This effect enables the system to generate stronger retractor muscle forces to compensate for the increased load, and hence achieve strong robustness. The approaches that we are applying to understanding a neuromechanical model in Aplysia, and the results that we have obtained, are likely to provide insights into the function of other motor systems that encounter changing mechanical loads and hard boundaries, both due to mechanical and neuronal firing properties.

59 BASIC BIOLOGICAL SCIENCES↗

Frustration-driven magnetic correlations in the spin-$\frac{5}{2}$ triangular lattice antiferromagnet RbFe⁢(HPO 3 ) 2

Here, a detailed study of the structural and magnetic properties of a spin-5/2 triangular lattice antiferromagnet RbFe⁢(HPO 3 ) 2 is presented using x-ray diffraction, magnetization, heat capacity, and 31 P nuclear magnetic resonance (NMR) experiments on a polycrystalline sample. The crystal structure features an equilateral triangular lattice of Fe 3+ ions. The thermodynamic measurements reveal the onset of a magnetic long-range order at 𝑇 N⁢1 ≃7.8K in zero field, followed by another low-temperature field-induced ordering at 𝑇 N⁢2 in higher fields. The transition at 𝑇 N⁢1 is further confirmed from the NMR spin lattice relaxation measurements. The value of the frustration ratio (𝑓≃7) implies moderate spin frustration in the compound. The 31 P NMR spectra exhibit two distinct spectral lines corresponding to two inequivalent phosphorus sites (P1 and P2), consistent with the crystal structure. The P1 site is strongly coupled with an isotropic hyperfine coupling of 𝐴$^{iso}_{hf}$ = 0.55⁢(2)⁢ T/𝜇 B while the P2 site is weakly coupled with 𝐴$^{iso}_{hf}$ = 0.25⁢(3)⁢ T/𝜇 B with the Fe 3+ ions. The magnetic susceptibility and NMR shift data are described well assuming a spin-5/2 isotropic triangular lattice antiferromagnetic model with an average exchange coupling of 𝐽/𝑘 B = 2.8⁢(2)⁢K. Below 𝑇 N⁢1 , the spectra evolve into a nearly rectangular powder pattern, indicating a commensurate antiferromagnetic type order. The 31 P spin-lattice relaxation rate well below 𝑇N⁢1 follows a 𝑇 3 temperature dependence, implying a two-magnon Raman scattering mechanism in the ordered state. Three well-defined phase regimes are clearly ascertained in the 𝐻−𝑇 phase diagram, reflecting a weak magnetic anisotropy in the compound.

Nagpal, V. [Indian Institute of Science Education ↗

Attosecond Time Delay Trends across the Isoelectronic Noble Gas Sequence

The analysis and measurement of Wigner time delays can provide detailed information about the electronic environment within and around atomic and molecular systems, with one the key differences being the lack of a long-range potential after a halogen ion undergoes photoionization. In this work, we use relativistic random-phase approximation to calculate the average Wigner delay from the highest occupied subshells of the atomic pairings (2p, 2s in Fluorine, Neon), (3p, 3s in Chlorine, Argon), (4p, 4s, 3d, in Bromine, Krypton), and (5p, 5s, 4d in Iodine, Xenon). The qualitative behaviors of the Wigner delays between the isoelectronic pairings were found to be similar in nature, with the only large differences occurring at photoelectron energies less than 20 eV and around Cooper minima. Interestingly, the relative shift in Wigner time delays between negatively charged halogens and noble gases decreases as atomic mass increases. All atomic pairings show large differences at low energies, with noble gas atoms showing large positive Wigner delays, while negatively charged halogen ions show negative delays. The implications for photoionization studies in halide-containing molecules is also discussed.

Grafstrom, Brock↗