Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex scattering”

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

Properties of Backscattered Ultraviolet Reflectances in Cloudy Atmosphere

We use a Plane-Parallel Cloud (PPC) model to illustrate how Mie scattering from cloud particles interacts with Rayleigh scattering in the atmosphere and produces a complex wavelength dependence in the top-of-the-atmosphere (TOA) reflectances measured by satellite instruments that operate in the ultraviolet (UV) part of the spectrum. Comparisons of the PPC model-derived spectral dependence of reflectances with the Total Ozone Mapping Spectrometer (TOMS) measurements show surprisingly good agreement over a wide range of observational conditions. The PPC model results also are compared with the results of two other cloud models: Lambert Equivalent Reflectivity (LER) and Modified Lambert Equivalent Reflectivity (MLER) that have been used to analyze satellite data in the UV. These models assume that clouds are opaque Lambertian reflectors rather than Mie scattering particles. Although one of these models (MLER) agrees reasonably well with the data, the results from this model appear somewhat unphysical and may not be suitable for interpreting satellite data if one desires high accuracy. We also use the PPC model to illustrate how clouds can perturb tropospheric O3 absorption in complex ways that cannot be explained by models that treat them as reflecting surfaces rather than as volume scatterers.

Ahmad, Z.↗

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗

Non-Invasive Biophysical Techniques to Monitor the Structural Plasticity of the Photosynthetic Machinery of Live Diatom Cells

The photosynthetic performance of diatoms depends largely on the organization and structural flexibility of their thylakoid membranes, the densely packed, highly organized membrane vesicles in which light reactions of photosynthesis occur. Different regulatory mechanisms that fine tune the photosynthetic functions affect the organization of the photosynthetic machinery at different levels of structural complexity, from the level of individual protein complexes to the macroarray of membrane proteins and the remodeling of the entire thylakoid membrane system. To monitor these reorganizations, non-invasive techniques are of special value. In this chapter, we focus our attention on three of these techniques, which have been demonstrated to provide unique and useful information on the structure and structural and functional plasticity of live diatom cells: (i) circular dichroism (CD) spectroscopy, which has provided unique information on the chiral (macro-)organization of protein complexes and on their rapid, reversible reorganizations, fine-tuning the light-harvesting processes, as well as on variations in the short-range excitonic interactions in the antenna complexes; (ii) small-angle neutron scattering (SANS), which has been used to determine the periodic organization of the thylakoid membranes and to monitor reversible ultrastructural changes on the time-scale of minutes, induced by variations in the environmental conditions such as changes in temperature or light intensity; and (iii) electrochromic shift absorbance transients (ΔA ECS ), a spectroscopic tool which has been shown to be capable of identifying distinct functional groups of the light-harvesting carotenoid fucoxanthin in different diatoms and in cells exposed to different light intensities. Future use of these techniques will most certainly contribute to the deeper understanding of key regulatory mechanisms of photosynthesis in diatoms.

Szabo, Milan↗

STILGAR End-of-Project Report

The Subsurface Tunnel Imaging LeveraGed by Analysis of Rayleigh wave ellipticity (STILGAR) project demonstrated an integrated geophysical approach for detecting, locating, and characterizing underground structural changes using dense seismic arrays and advanced inversion techniques. Field campaigns were conducted at two operational mines—the Redmond salt mine (Utah) and Graymont Pleasant Gap limestone mine (Pennsylvania)—providing real-world testbeds for monitoring anthropogenic subsurface activity. At the Redmond salt mine, seismic interferometry combined with back-projection inversion successfully identified continuous, low-amplitude signals from mining operations. The approach differentiated stationary from migrating anthropogenic sources, captured daily operational cycles, and validated the potential of passive seismic monitoring for remote detection of underground activity. At the Graymont Pleasant Gap mine, two dense seismic deployments in the spring and fall of 2023 generated over 4 TB of high-resolution data. Key outcomes included the relocation of 199 underground and 8 surface explosions with accuracies within tens of meters and the development of a 3D P-wave velocity model using the triple-difference tomography algorithm (tomoTD) that resolved major structural features such as the mine entrance, low-velocity tunnels, and roof-collapse areas. Ambient noise cross-correlation and back-projection analyses revealed persistent sources linked to ongoing mining activity, whereas horizontal-to-vertical spectral ratio (HVSR) and ellipticity studies confirmed stable site responses across seasons and identified soil thickness trends consistent with regional erosional and depositional processes. Checkerboard and sensitivity tests further validated the robustness of the tomographic results. Overall, the findings emphasize that although significant progress has been made in subsurface imaging, further work is needed to enhance the detection and localization of underground structures. Accurate imaging requires higher frequencies, yet anthropogenic sources tend to dominate the seismic record at those frequencies, and high-frequency surface waves are affected by higher modes that complicate interpretation. The improved detection and localization of human-induced signals enabled detailed temporal and spatial mapping of daily mine operations, demonstrating the feasibility of continuous anthropogenic source monitoring. Sensitivity to signals from nontraditional sources, such as fan operations, highlights the broader applicability of this approach to other industrial environments in which continuous and impulsive signals are present. The field campaigns produced a substantial volume of high-quality seismic data, supporting the development and testing of new methods for seismic source characterization and subsurface imaging. Future deployments should include sensors capable of recording lower frequencies to probe deeper structures, increase bandwidth to enhance resolution and sensitivity to both shallow and deep targets, and collect additional large-scale datasets to refine imaging and source characterization techniques. Moreover, conducting 3D modeling studies of seismic wavefields at higher frequencies will provide a better understanding of wave scattering and cavity–wavefield interactions in complex underground environments. In conclusion, the STILGAR project demonstrated that integrated seismic monitoring can effectively characterize underground operations, capturing both natural and anthropogenic signals. The approaches developed provide a foundation for improved detection, localization, and imaging of subsurface structures and are directly transferable to broader industrial monitoring applications.

58 GEOSCIENCES↗

A resonant state and the ground state of positronium hydride

The lowest-lying resonance occurring in S-wave positronium-hydrogen scattering is reinvestigated, using the complex-rotation method. By employing a generalized Hylleraas-type wave function that includes all six interparticle coordinates, a very accurate value of the resonance position is obtained, along with a good value of the width. The present result for the resonance position (-1.205 plus or minus 0.001 Ry) is lower than the previous result of Drachman and Houston, who omitted the interelectronic coordinate in their trial function. In addition, the lowest ground-state energy of positronium hydride is obtained by using 210 terms in the trial wave function. The effect of the interelectronic coordinate and others on both the resonant energy and the binding energy of PsH is discussed.

Ho, Y. K.↗

Demonstration of a stabilized alumina/ethanol colloidal dispersion technique for seeding high temperature air flows

Laser anemometry enables the measurement of complex flow fields via the light scattered from small particles entrained in the flow. In the study of turbomachinery, refractory seed materials are required for seeding the flow due to the high temperatures encountered. In this work we present a pH stabilization technique commonly employed in ceramic processing to obtain stable dispersions for generating aerosols of refractory seed material. By adding submicron alumina particles to a preadjusted pH solution of ethanol, a stable dispersion is obtained which when atomized, produces a high quality aerosol. Commercial grade alumina powder is used with a moderate size distribution. Other metal oxide powders in various polar solvents could also be used once the point of zero charge (pH(pzc)) of the powder in the solvent has been determined. Laser anemometry measurements obtained using the new seeding technique are compared to measurements obtained using Polystyrene Latex (PSL) spheres as the seed material.

Wernet, Mark P.↗

Retrieval of Aerosol Properties from Multi-Spectral Extinction Measurements

The direct-beam spectral extinction of solar radiation contains information on atmospheric composition in a form that is essentially free from the data analysis complexities that often arise from multiple scattering. Ground based Multi-Filter Shadowband Radiometer (MFRSR) measurements provide such information for the vertical atmospheric column path, while solar occultation measurements from a satellite platform provide horizontal slices through the atmosphere. We describe application of a Multi-Spectral Atmospheric Column Extinction (MACE) analysis technique used to analyze MFRSR data also to occultation measurements made by SAGE II. For analysis, we select the 1985 Nevado del Ruiz volcanic eruption period to retrieve atmospheric profiles of ozone and NO2, and changes in the stratospheric aerosol size and optical depth. The time evolution of volcanic aerosol serves as a passive tracer to study stratospheric dynamics, and changes in particle size put constraints on the sulfur chemistry modeling of volcanic aerosols. Paper presented at The '99 Kyoto Aerosol-Cloud Workshop, held Dec 1-3, 1999, Kyoto, Japan

Lacis, Andrew A.↗

Retrieval of Aerosol Properties from Multi-Spectral Extinction Measurements

The direct-beam spectral extinction of solar radiation contains information on atmospheric composition in a form that is essentially free from the data analysis complexities that often arise from multiple scattering. Ground based Multi-Filter Shadowband Radiometer (MFRSR) measurements provide such information for the vertical atmospheric column path, while solar occultation measurements from a satellite platform provide horizontal slices through the atmosphere. We describe application of a Multi-Spectral Atmospheric Column Extinction (MACE) analysis technique used to analyze MFRSR data also to occultation measurements made by SAGE 11. For analysis, we select the 1985 Nevado del Ruiz volcanic eruption period to retrieve atmospheric profiles of ozone and NO2, and changes in the stratospheric aerosol size and optical depth. The time evolution of volcanic aerosol serves as a passive tracer to study stratospheric dynamics, and changes in particle size put constraints on the sulfur chemistry modeling of volcanic aerosols.

Lacis, Andrew A.↗

Large Scale Finite Element Modeling Using Scalable Parallel Processing

An iterative solver for use with finite element codes was developed for the Cray T3D massively parallel processor at the Jet Propulsion Laboratory. Finite element modeling is useful for simulating scattered or radiated electromagnetic fields from complex three-dimensional objects with geometry variations smaller than an electrical wavelength.

finite element modeling parallel processing iterat↗

Simultaneous Retrieval of Aerosol and Surface Optical Properties from Combined Airborne- and Ground-Based Direct and Diffuse Radiometric Measurements

This paper presents a new method for simultaneously retrieving aerosol and surface reflectance properties from combined airborne and ground-based direct and diffuse radiometric measurements. The method is based on the standard Aerosol Robotic Network (AERONET) method for retrieving aerosol size distribution, complex index of refraction, and single scattering albedo, but modified to retrieve aerosol properties in two layers, below and above the aircraft, and parameters on surface optical properties from combined datasets (Cloud Absorption Radiometer (CAR) and AERONET data). A key advantage of this method is the inversion of all available spectral and angular data at the same time, while accounting for the influence of noise in the inversion procedure using statistical optimization. The wide spectral (0.34-2.30 m) and angular range (180 ) of the CAR instrument, combined with observations from an AERONET sunphotometer, provide sufficient measurement constraints for characterizing aerosol and surface properties with minimal assumptions. The robustness of the method was tested on observations made during four different field campaigns: (a) the Southern African Regional Science Initiative 2000 over Mongu, Zambia, (b) the Intercontinental Transport Experiment-Phase B over Mexico City, Mexico (c) Cloud and Land Surface Interaction Campaign over the Atmospheric Radiation Measurement (ARM) Central Facility, Oklahoma, USA, and (d) the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) over Elson Lagoon in Barrow, Alaska, USA. The four areas are dominated by different surface characteristics and aerosol types, and therefore provide good test cases for the new inversion method.

Gatebe, C. K.↗

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

Continuum, cyclotron line, and absorption variability in the high-mass X-ray binary Vela X-1

Because of its complex clumpy wind, prominent cyclotron resonant scattering features, intrinsic variability, and convenient physical parameters (close distance, high inclination, and small orbital separation), which facilitate the observation and analysis of the system, Vela X-1 is one of the key systems for understanding accretion processes in high-mass X-ray binaries on all scales. We revisit Vela X-1 with two new observations taken with NuSTAR at orbital phases ∼0.68–0.78 and ∼0.36–0.52, which show a plethora of variability and allow us to study the accretion geometry and stellar wind properties of the system. We follow the evolution of spectral parameters down to the pulse period timescale using a partially covered power law continuum with a Fermi-Dirac cutoff to model the continuum and local absorption. We are able to confirm anti-correlations between the photon index and the luminosity and, for low fluxes, between the folding energy and the luminosity, implying a change of properties in the Comptonising plasma. We were not able to confirm a previously seen correlation between the cyclotron line energy and the luminosity of the source in the overall observation, but we observed a drop in the cyclotron line energy following a strong flare. We see strong variability in absorption between the two observations and within one observation (for the ∼0.36–0.52 orbital phases) that can be explained by the presence of a large-scale structure, such as accretion and photoionisation wakes in the system, and our variable line of sight through this structure.

C. M. Diez↗

Freely Suspended Smectic Filaments and the Structure of the B7 Phase of MHOBOW

Our recent discovery of the spontaneous formation of chiral domains in fluid smectic phases of achiral bow-shaped molecules opens up a wide variety of possibilities for new liquid crystal phases and phenomena. The basic, spontaneously chiral layer structure of the highest temperature fluid smectic phases, the B2 and B7, are shown. One of the most intriguing aspects of this structure is the plethora of possible phases coming from different stacking sequences of the polar ordering and tilt directions. The four possibilities of next-nearest neighbor alternation are shown. In the original material studied, NOBOW, the ground states found are antiferroelectric, either the racemic SmC(sub S)P(sub A) or the chiral SmC(sub A)P(sub A). We are currently studying MHOBOW, synthesized by D. Walba which, by virtue of its methyl hexyloxy tail has a tendency to form anticlinic layer interfaces, in the hope of finding a phase with a ferroelectric ground state, either SmC(sub A)P(sub S) or SmC(sub S)P(sub A), which can be obtained in NOBOW only by applying a field. Preliminary observations of MHO-BOW have made its study, from the point of view of understanding novel LC structures, extremely high priority. The following truly remarkable characteristics have been revealed: (i) The smectic phase grows out of the isotropic in the form of helical ribbons. The resulting planar aligned textures of focal conics with layers normal to glass plates exhibit bizarre modulations, including stripes and checker-boards. These have also been seen in other materials suggesting that this is a new phase (tentatively called B7), which is a fluid smectic with some kind of in-layer structure. (ii) It is virtually impossible to make freely suspended films of MHOBOW. Rather it makes the freely suspended filaments which preliminary x-ray scattering experiments reveal to have the nested cylinder layer structure indicated; (iii) The powder x-ray diffraction exhibits four resolution-limited smectic layering peaks, very close in layer spacing, which vary continuously with T. This is further evidence for a more complex three dimensional structure than NOBOW, which has a typical single layering reflection. (iv) The x-ray structure factor of the layering peak of the filaments is extraordinarily complex and rich. Varying in qL (the scattering vector component along the filament axis) from a double slit-like pattern to modulated layer-like patterns, as qH (the scattering vector component normal to the filament axis) is varied over the range where the four powder peaks are located. These results suggest some kind of mosaic structure, perhaps with different layer spacings corresponding to the different stacking sequences. Recent x-ray diffraction experiments show that the peaks are modulated in intensity upon translation along a filament, in domains of several hundred microns dimension. These preliminary experiments suggest that the B7 is a fluid smectic with extremely unusual and fascinating structures. Of all of the many hundreds of fluid smectic materials we have attempted to study in the freely suspended film geometry over the years, only a few have failed to form films, and none showed any great tendency to form filaments, although this clearly should be a possible freely suspended smectic LC morphology. On several occasions in the past we have intentionally tried to make filaments from a variety of smectics without success. Thus the smectic filament formation property makes the B7 phase unique. It seems quite likely that the stability of filaments is related to the in-plane structure. The filaments exhibit other interesting structural and optical features. They are birefringent with a local optic axis which is oblique and which can vary continuously along filament and which can be manipulated with an electric field applied normal to the fiber, as if the field were causing a rotation of the optic axis about the fiber axis. Rapid displacement of the ends of the fiber toward one another causes a macroscopic helixing at low T and causes thick regions to transiently appear at high T, a 1D analog of island formation on a rapidly compressed film.

Clark, N.↗

Electron detachment from negative ions - The effects of isotopic substitution

Absolute total electron-detachment cross sections and relative elastic-scattering cross sections are presented for collisions of H(-) or D(-) with He, Ne, Ar, and N2 for collision energies ranging from about 2 eV up to 100 eV. Special emphasis has been given to the effect of isotopic substitution. It was found that, with the exception of the helium target, the results cannot be in accord with a model using a localized complex potential. On the other hand, the results for H(-) or D(-) with He are in excellent agreement with predictions based on previous (local complex potential) analysis of differential elastic-scattering measurements.

Champion, R. L.↗

Evaluation of Daytime Measurements of Aerosols and Water Vapor made by an Operational Raman Lidar over the Southern Great Plains

Raman lidar water vapor and aerosol extinction profiles acquired during the daytime over the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in northern Oklahoma (36.606 N, 97.50 W, 315 m) are evaluated using profiles measured by in situ and remote sensing instruments deployed during the May 2003 Aerosol Intensive Operations Period (IOP). The automated algorithms used to derive these profiles from the Raman lidar data were first modified to reduce the adverse effects associated with a general loss of sensitivity of the Raman lidar since early 2002. The Raman lidar water vapor measurements, which are calibrated to match precipitable water vapor (PWV) derived from coincident microwave radiometer (MWR) measurements were, on average, 5-10% (0.3-0.6 g/m(exp 3) higher than the other measurements. Some of this difference is due to out-of-date line parameters that were subsequently updated in the MWR PWV retrievals. The Raman lidar aerosol extinction measurements were, on average, about 0.03 km(exp -1) higher than aerosol measurements derived from airborne Sun photometer measurements of aerosol optical thickness and in situ measurements of aerosol scattering and absorption. This bias, which was about 50% of the mean aerosol extinction measured during this IOP, decreased to about 10% when aerosol extinction comparisons were restricted to aerosol extinction values larger than 0.15 km(exp -1). The lidar measurements of the aerosol extinction/backscatter ratio and airborne Sun photometer measurements of the aerosol optical thickness were used along with in situ measurements of the aerosol size distribution to retrieve estimates of the aerosol single scattering albedo (omega(sub o)) and the effective complex refractive index. Retrieved values of omega(sub o) ranged from (0.91-0.98) and were in generally good agreement with omega(sub o) derived from airborne in situ measurements of scattering and absorption. Elevated aerosol layers located between about 2.6 and 3.6 km were observed by the Raman lidar on May 25 and May 27. The airborne measurements and lidar retrievals indicated that these layers, which were likely smoke produced by Siberian forest fires, were primarily composed of relatively large particles (r(sub eff) approximately 0.23 micrometers), and that the layers were relatively nonabsorbing (omega(sub o) approximately 0.96-0.98). Preliminary results show that major modifications that were made to the Raman lidar system during 2004 have dramatically improved the sensitivity in the aerosol and water vapor channels and reduced random errors in the aerosol scattering ratio and water vapor retrievals by an order of magnitude.

Ferrare, Richard↗

Dependent scattering and fractal microstructure determine the transparency of aerogel monoliths

This study reveals how dependent scattering and microstructure significantly affect electromagnetic wave propagation through aerogel monoliths, contributing to their transparency. Light scattering by particle ensembles is considered “dependent” when the scattering properties rely not only on particle size and optical constants but also on their spatial distribution, typically occurring when the average interparticle distance is small in comparison with the wavelength of incident radiation. Addressing dependent scattering requires solving Maxwell’s equations for complex heterogeneous structures, which is computationally demanding and usually limited to sample thicknesses on the same scale as the wavelength. This study combines computer-generated ambigel microstructures of fractal aggregates of polydisperse nanoparticles and the radiative transfer with reciprocal transaction method to predict the transmittance of thick ambigel slabs. Transmittance measurements of ambiently dried aerogel monoliths (ambigels) with porosities from about 50% to 90% closely matched the predicted values for their digital twins. However, ignoring dependent scattering or particle aggregation led to inaccurate predictions. This study validated the computational framework, and its findings offer insights for designing photonic metamaterials and analyzing their interactions with electromagnetic waves.

Yalcin, Refet A. (ORCID:0000000339973494)↗

Stability and bifurcation of quasiparallel Alfven solitons

The inverse scattering transformation (IST) is used to study the one-parameter and two-parameter soliton families of the derivative nonlinear Schroedinger (DNLS) equation. The two-parameter soliton family is determined by the discrete complex eigenvalue spectrum of the Kaup-Newell scattering problem and the one-parameter soliton family corresponds to the discrete real eigenvalue spectrum. The structure of the IST is exploited to discuss the existence of discrete real eigenvalues and to prove their structural stability to perturbations of the initial conditions. Also, though the two-parameter soliton is structurally stable in general, it is shown that a perturbation of the initial conditions may change the two-parameter soliton into a degenerate soliton which, in turn, is structurally unstable. This degenerate, or double pole, soliton may bifurcate due to a perturbation of the initial conditions into a pair of one-parameter solitons. If the initial profile is on compact support, then this pair of one-parameter solitons must be compressive and rarefactive respectively. Finally, the Gelfand-Levitan equations appropriate for the double pole soliton are solved.

Hamilton, R. L.↗

Insights into distorted lamellar phases with small-angle scattering and machine learning

Lamellar phases are essential in various soft matter systems, with topological defects significantly influencing their mechanical properties. In this report, we present a machine-learning approach for quantitatively analyzing the structure and dynamics of distorted lamellar phases using scattering techniques. By leveraging the mathematical framework of Kolmogorov–Arnold networks, we demonstrate that the conformations of these distorted phases – expressed as superpositions of complex waves – can be reconstructed from small-angle scattering intensities. Through the contour analysis of wave field phase singularities, we obtain the statistics of the spatial distribution of topological defects. Furthermore, we establish that the temporal evolution of these defects can be derived from the time-dependent traveling wave field, informed by the dispersion relation of spectral components. This method opens new avenues for investigating the dynamics of distorted lamellar phases using various dynamic scattering techniques such as neutron spin echo and X-ray photon correlation spectroscopy. These findings enhance our microscopic understanding of how defects influence the physical properties of lamellar materials, with implications for both equilibrium and non-equilibrium states in general lamellar systems.

36 MATERIALS SCIENCE↗