Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic range”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A Practical Framework for Simulating Time-Resolved Spectroscopy Based on a Real-Time Dyson Expansion

Time-resolved spectroscopy is a powerful tool for probing electron dynamics in molecules and solids, revealing transient phenomena on subfemtosecond time scales. The interpretation of experimental results is often enhanced by parallel numerical studies, which can provide insight and validation for experimental hypotheses. However, developing a theoretical framework for simulating time-resolved spectra remains a significant challenge. The most suitable approach involves the many-body nonequilibrium Green's function formalism, which accounts for crucial dynamical many-body correlations during time evolution. While these dynamical correlations are essential for observing emergent behavior in time-resolved spectra, they also render the formalism prohibitively expensive for large-scale simulations. Substantial effort has been devoted to reducing this computational cost─through approximations and numerical techniques─while preserving the key dynamical correlations. The ultimate goal is to enable first-principles simulations of time-dependent systems ranging from small molecules to large, periodic, multidimensional solids. Here, in this perspective, we outline key challenges in developing practical simulations for time-resolved spectroscopy, with a particular focus on Green's function methodologies. We highlight a recent advancement toward a scalable framework: the real-time Dyson expansion (RT-DE) [Phys. Rev. Lett. 2024, 133, 226902]. We introduce the theoretical foundation of RT-DE and discuss strategies for improving scalability, which have already enabled simulations of system sizes beyond the reach of previous fully dynamical approaches. We conclude with an outlook on future directions for extending RT-DE to first-principles studies of dynamically correlated, nonequilibrium systems.

Reeves, Cian C. [Univ. of California, Santa Barbar↗
DOE OSTI2020

Moving beyond the Solvent-Tip Approximation to Determine Site-Specific Variations of Interfacial Water Structure through 3D Force Microscopy

Although interfacial solution structure impacts environmental, biological, and technological phenomena, including colloidal stability, protein assembly, heterogeneous nucleation, and water desalination, its molecular details remain poorly understood. Here, we visualize the three-dimensional (3D) hydration structure at the boehmite(010)-water interface using fast force mapping (FFM). Using a self-consistent scheme to decouple long-range tip-sample interactions from short-range solvation forces, we obtain the solution structure with lattice resolution. The results are benchmarked against molecular dynamics simulations that explicitly include the effects of the tip with different levels of approximation and systematically account for tip size, chemistry, and confinement effects. We find four laterally structured water layers within 1 nm of the surface, with the highest water densities at sites adjacent to hydroxyl groups. The key features beyond the first two layers can only be predicted using a full-scale simulation of the boehmite-water-silica system. Our findings further reveal a complex relationship between site-specific chemistry, water density, and long-range particle interactions; and present important advances toward quantitative data interpretation in 3D FFM.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Two-temperature warm dense hydrogen as a test of quantum protons driven by orbital-free density functional theory electronic forces

We consider a steady-state (but transient) situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature Te, ions at Ti, and Te ≠ Ti. Such states are achievable by pump–probe experiments. For warm dense hydrogen in such a two-temperature situation, we investigate nuclear quantum effects (NQEs) on structure and thermodynamic properties, thereby delineating the limitations of ordinary ab initio molecular dynamics. We use path integral molecular dynamics (PIMD) simulations driven by orbital-free density functional theory (OFDFT) calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence. We calibrate the OFDFT calculations against conventional (explicit orbitals) Kohn–Sham DFT. We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30, the ionic radial distribution function is meaningfully affected by the inclusion of NQEs. Moreover, NQEs induce a substantial increase in both the ionic and electronic pressures. This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen. For Te > 20 kK, increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered. In addition, we confirm that compared with thermostatted ring-polymer molecular dynamics, the primitive PIMD algorithm overestimates electronic pressures, a consequence of the overly localized ionic description from the primitive scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Impact of battery cell imbalance on electric vehicle range

Due to manufacturing variation, battery cells often possess heterogeneous characteristics, leading to battery state-of-charge variation in real-time. Since the lowest cell state-of-charge determines the useful life of battery pack, such variation can negatively impact the battery performance and electric vehicles range. Existing research has been focused on control design to mitigate cell imbalance. However, it is yet unclear how much impacts the cell imbalance can have on electric vehicle range. This paper closes this knowledge gap by using a simulation environment consisting of real-world driving speed data, vehicle longitudinal control, propulsion and vehicle dynamics, and cell level battery modeling. In particular, each battery cell is modeled as an equivalent circuit model, and variations among cell parameters are introduced to assess their impact on electric vehicles range and to identify the most influential parameter variations. Simulation results and analysis can be used to assist balancing control design and to benchmark control performance.

25 ENERGY STORAGE↗
DOE OSTI2024

Comprehensive Dynamics in a Polyelectrolyte Complex Coacervate

The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine its coil radius of gyration, R g , within PECs using small-angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 °C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η 0 , versus chain length, N, was η 0 ~ N 3.1 , in agreement with “sticky reptation” theory. The lifetime and activation energy, E p , of a pair between polyanion and polycation repeat units, Pol + Pol – , were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2E p , showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN– acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol + Pol – pairing efficiently, which decreases effective cross-linking and decreases viscosity. Finally, an equation was derived that quantitatively predicts this effect.

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

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗
DOE OSTI2020

Advanced grazing-emission X-ray fluorescence spectroscopy with small-angle neutron scattering for in-vivo surface composition and defect/morphology surface evolution in tokamak PMI

The dynamic and extreme conditions of thermonuclear fusion tokamak plasmas render material surfaces almost impossible to examine in real-time, surface-sensitive conditions. However, combining particle-probe techniques that are not charge-state dependent open the possibilities of in-vivo diagnosis of the plasma-material interface. The reconstituted surfaces of plasma-facing components ranging from first-wall to divertor-component regions in the device vary depending on the complex plasma transport characteristics of the device. For example, in ITER it is predicted that fluxes of the order 10 22 to 10 24 and incident charged-particle energies of 100’s to 1000’s of eV will reach the surface while in the divertor regions within the private-flux and near the inner and outer-strike points would vary from a few eV to 100’s of eV governed by the complex Chodura magnetic sheath dynamics. This means that the penetration range of most charged-particle fluxes would range between 1-100’s of nm (e.g. for fuel particles of D and T and He).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Non-Adiabatic Excited State Molecular Dynamics Methodologies: comparison and convergence

Direct atomistic simulation of nonadiabatic molecular dynamics is a challenging goal that allows important insights into fundamental physical phenomena. A variety of frameworks, ranging from fully quantum treatment of nuclei to semiclassical and mixed quantum–classical approaches, were developed. These algorithms are then coupled to specific electronic structure techniques. Such diversity and lack of standardized implementation make it difficult to compare the performance of different methodologies when treating realistic systems. Here, we compare three popular methods for large chromophores: Ehrenfest, surface hopping, and multiconfigurational Ehrenfest with ab initio multiple cloning (MCE-AIMC). These approaches are implemented in the NEXMD software, which features a common computational chemistry model. The resulting comparisons reveal the method performance for population relaxation and coherent vibronic dynamics. Finally, we study the numerical convergence of MCE-AIMC algorithms by considering the number of trajectories, cloning thresholds, and Gaussian wavepacket width. Our results provide helpful reference data for selecting an optimal methodology for simulating excited-state molecular dynamics.

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

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Statics and dynamics of skyrmions interacting with disorder and nanostructures

Magnetic skyrmions are topologically stable nanoscale particlelike objects that were discovered in 2009. Since that time, intense research interest in the field has led to the identification of numerous compounds that support skyrmions over a range of conditions spanning from cryogenic to room temperatures. Skyrmions can be set into motion under various types of driving, and the combination of their size, stability, and dynamics makes them ideal candidates for numerous applications. At the same time, skyrmions represent a new class of system in which the energy scales of the skyrmion-skyrmion interactions, sample disorder, temperature, and drive can compete. A growing body of work indicates that the static and dynamic states of skyrmions can be influenced strongly by pinning or disorder in the sample; thus, an understanding of such effects is essential for the eventual use of skyrmions in applications. The current state of knowledge regarding individual skyrmions and skyrmion assemblies interacting with quenched disorder or pinning is reviewed. Additionally, the microscopic mechanisms for skyrmion pinning, including the repulsive and attractive interactions that can arise from impurities, grain boundaries, or nanostructures, are outlined. This is followed by descriptions of depinning phenomena, sliding states over disorder, the effect of pinning on the skyrmion Hall angle, the competition between thermal and pinning effects, the control of skyrmion motion using ordered potential landscapes such as one- or two-dimensional periodic asymmetric substrates, the creation of skyrmion diodes, and skyrmion ratchet effects. Highlighted are the distinctions arising from internal modes and the strong gyrotropic or Magnus forces that cause the dynamical states of skyrmions to differ from those of other systems with pinning, such as vortices in type-II superconductors, charge density waves, or colloidal particles. Throughout this review future directions and open questions related to the pinning and dynamics in skyrmion systems are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Probing correlated states with plasmons

Understanding the nature of strongly correlated states in flat-band materials (such as moiré heterostructures) is at the forefront of both experimental and theoretical pursuits. While magnetotransport, scanning probe, and optical techniques are often very successful in investigating the properties of the underlying order, the exact nature of the ground state often remains unknown. Here, we propose to leverage strong light-matter coupling present in the flat-band systems to gain insight through dynamical dielectric response into the structure of the many-body ground state. We argue that because of the enlargement of the effective lattice of the system arising from correlations, conventional long-range plasmon becomes “folded” to yield a multiband plasmon spectrum. We detail several mechanisms through which the structure of the plasmon spectrum and that of the dynamical dielectric response is susceptible to the underlying order, revealing valued insights such as the interaction-driven band gaps, spin-structure, and the order periodicity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimization of Archimedes Screw for Use in Hydroelectric Projects

In 2017, the Department of Energy’s (DOE) Water Power Technologies Office (WPTO) made federal funding available to several awardees to investigate and provide innovative technological solutions to reduce capital costs of installed hydroelectric facilities at non-powered dam’s (NPD’s) across the United States, while maintaining a high level of efficiency. Awardees were to develop specific technologies that showed promise of reducing installed capital cost (ICC) by at least 20% over a given baseline and maintain a water-to-wire efficiency of at least 80% at head levels of 50 feet or less and flows of 1000 cubic feet per second (cfs) or less. Canyon Hydro (Canyon) was selected as an awardee for this funding and worked to refine the implementation of the Archimedes screw turbine for use at these site conditions. Canyon and its project partners analyzed several dam inventory datasets to determine the most common head and flow range of current NPD’s in order to best select design conditions for a turbine that would perform successfully at common heads and flows. Computational Fluid Dynamic (CFD) analyses were conducted after site conditions had been selected to optimize the turbine geometry for efficiency. Through this effort we were able to estimate water-to-wire efficiency over a range of flows. Equipment design focused on concepts of modularity and pre-fabricated construction. These methodologies were used where possible to reduce costly on-site construction and shorten deployment timelines. Cost modeling was conducted with input from industry experts in general construction, pre-fabricated concrete construction, and composite fabrication. As a hydroelectric turbine manufacturer of 35 years, Canyon drew upon its own knowledge base in steel fabrication, water-to-wire component specification and ancillary systems to develop the remaining cost matrix for the project. As a result of these efforts, Canyon was able to present findings that indicated a modular Archimedes screw turbine system could be designed that reduced ICC by 23.9% and maintain a maximum water-to-wire efficiency of 82.3%, meeting the goals of the project.

13 HYDRO ENERGY↗
DOE OSTI2020

Simulations and symmetries

ABSTRACT We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be efficiently achieved by modern N-body methods with a symmetries-based bias expansion which describes the clustering of any tracer on large scales. We show that at low redshifts, and for moderately biased tracers, the substitution of N-body-determined dynamics improves over an equivalent model using perturbation theory by more than a factor of two in scale, while at high redshifts and for highly biased tracers the gains are more modest. This hybrid approach lends itself well to emulation. By removing the need to identify haloes and subhaloes, and by not requiring any galaxy-formation-related parameters to be included, the emulation task is significantly simplified at the cost of modelling a more limited range in scale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Phase transition dynamics in one-dimensional halide perovskite crystals

Triiodide perovskites CsPbI 3 , CsSnI 3 , and FAPbI 3 (where FA is formamidinium) are highly promising materials for a range of optoelectronic applications in energy conversion. However, they are thermodynamically unstable at room temperature, preferring to form low-temperature (low-T) non-perovskite phases with one-dimensional anisotropic crystal structures. While such thermodynamic behavior represents a major obstacle toward realizing high-performance devices based on their high-temperature (high-T) perovskite phases, the underlying phase transition dynamics are still not well understood. In this work, we use in situ optical micro-spectroscopy to quantitatively study the transition from the low-T to high-T phases in individual CsSnI 3 and FAPbI 3 nanowires. We reveal a large blueshift in the photoluminescence (PL) peak (~38 meV) at the low-T/high-T two-phase interface of partially transitioned FAPbI 3 wire, which may result from the lattice distortion at the phase boundary. Compared to the experimentally derived activation energy of CsSnI 3 (~1.93 eV), the activation energy of FAPbI 3 is relatively small (~0.84 eV), indicating a lower kinetic energy barrier when transitioning from a face-sharing octahedral configuration to a corner-sharing one. Further, the phase propagation rate in CsSnI 3 is observed to be relatively high, which may be attributed to a high concentration of Sn vacancies. Our results could not only facilitate a deeper understanding of phase transition dynamics in halide perovskites with anisotropic crystal structures, but also enable controllable manipulation of optoelectronic properties via local phase engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase transition dynamics in one-dimensional halide perovskite crystals

Triiodide perovskites CsPbI 3 , CsSnI 3 , and FAPbI 3 (where FA is formamidinium) are highly promising materials for a range of optoelectronic applications in energy conversion. However, they are thermodynamically unstable at room temperature, preferring to form low-temperature (low-T) non-perovskite phases with one-dimensional anisotropic crystal structures. While such thermodynamic behavior represents a major obstacle toward realizing high-performance devices based on their high-temperature (high-T) perovskite phases, the underlying phase transition dynamics are still not well understood. Here we use in situ optical micro-spectroscopy to quantitatively study the transition from the low-T to high-T phases in individual CsSnI 3 and FAPbI 3 nanowires. We reveal a large blueshift in the photoluminescence (PL) peak (~38 meV) at the low-T/high-T two-phase interface of partially transitioned FAPbI 3 wire, which may result from the lattice distortion at the phase boundary. Compared to the experimentally derived activation energy of CsSnI 3 (~1.93 eV), the activation energy of FAPbI 3 is relatively small (~0.84 eV), indicating a lower kinetic energy barrier when transitioning from a face-sharing octahedral configuration to a corner-sharing one. Further, the phase propagation rate in CsSnI 3 is observed to be relatively high, which may be attributed to a high concentration of Sn vacancies. Furthermore, our results could not only facilitate a deeper understanding of phase transition dynamics in halide perovskites with anisotropic crystal structures, but also enable controllable manipulation of optoelectronic properties via local phase engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mapping structural and dynamic divergence across the MBOAT family

Membrane-bound O-acyltransferases (MBOATs) are membrane-embedded enzymes that catalyze acyl chain transfer to a diverse group of substrates, including lipids, small molecules, and proteins. MBOATs share a conserved structural core, despite wide-ranging functional specificity across both prokaryotes and eukaryotes. The structural basis of catalytic specificity, regulation and interactions with the surrounding environment remain uncertain. Here, we combine comparative molecular dynamics (MD) simulations with bioinformatics to assess molecular and interactional divergence across the family. In simulations, MBOATs differentially distort the bilayer depending on their substrate type. Additionally, we identify lipid binding sites surrounding reactant gates in the surrounding membrane. Complementary bioinformatic analyses reveal a conserved role for re-entrant loop-2 in MBOAT fold stabilization and a key hydrogen bond bridging DGAT1 dimerization. Finally, we predict differences in MBOAT solvation and water gating properties. These data are pertinent to the design of MBOAT-specific inhibitors that encompass dynamic information within cellular mimetic environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗