Search NASA⌕ Search

SEARCH · Search NASA

Results for “Crystal Field Theory”

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.

69 records · Page 4

A study of the applicability of nucleation theory to quasi-thermodynamic transformations of second and higher Ehrenfest-order

Transient and steady-state phenomena in temperature, stress, and electric, field intensity in ferroelectric polymers were investigated. The application and extension of the theory in the primary stage to the polarization domain nucleation and growth in ferroelectric polymers were developed. The kinetics of this growth were investigated. Expressions describing nucleation under the influence of an electric field were found through the expansion of the Gibbs' free energy in a Maclaurin series. The series was expanded in the electric field strength rather than the degree of undercooling. The resulting expressions were manipulated and applied to the case of nucleation of polarized domains in ferroelectric polymers. The kinetics of the nucleation and growth of polarized domains are also investigated. This was accomplished through the modification of the Johnson-Mehl-Avrami treatment of crystallization kinetics to be applicable to the growth of polarization domains in ferroelectric materials.

Barker, R. E., Jr.↗

Mass flux measurements and correlations in the back flow region of a nozzle plume.

A series of tests were conducted to measure the mass flux in the far field of a nozzle plume in a high vacuum with emphasis on the back flow region. The measurements presented in this report provide fairly accurate data for off-axis angles as large as 140 deg (i.e., the back flow region). The tests utilized five different nozzles. Parameters such as expansion ratio, throat diameter, nozzle lip shape, and plenum (chamber) pressure were varied. Carbon dioxide and nitrogen gases were passed through, and mass flux measurements were taken using quartz crystal microbalances in as many as nine different locations relative to the test nozzle. The tests resulted in a large matrix of data that were correlated and compared to existing prediction theories. Several conclusions with respect to the effect of nozzle and gas parameters on the amount of back flow mass flux are offered, and it was demonstrated that gaseous mass fluxes, which are not predictable by present theories, are encountered in the region behind the nozzle exit plane.

Chirivella, J. E.↗

Dipole-field sums and Lorentz factors for orthorhombic lattices, and implications for polarizable molecules

A method for computing the Lorentz tensor components in single crystals via rapidly convergent sums of Bessels functions is developed using the relationship between dipole-field sums and the tensor components. The Lorentz factors for simple, body-centered, and base-centered orthorhombic lattices are computed using this method, and the derivative Lorentz factors for simple orthorhombic lattices are also determined. Both the Lorentz factors and their derivatives are shown to be very sensitive to a lattice structure. The equivalent of the Clausius-Mossotti relation for general orthorhombic lattices is derived using the Lorentz-factor formalism, and the permanent molecular dipole moment is related to crystal polarization for the case of a ferroelectric of polarizable point dipoles. It is concluded that the polarization enhancement due to self-polarization familiar from classical theory may actually be a reduction in consequences of negative Lorentz factors in one or two lattice directions for noncubic crystals.

Purvis, C. K.↗

Waves in active and passive periodic structures - A review

The theory and recent applications of waves in periodic structures are reviewed. Both the Floquet and coupled waves approach are analyzed in some detail. The theoretical part of the paper includes wave propagation in unbounded and bounded active or passive periodic media, wave scattering from periodic boundaries, source radiation (dipole, Cerenkov, transition, and Smith-Purcell) in periodic media, and pulse transmission through a periodic slab. The applications part covers the recent development in a variety of fields: distributed feedback oscillators, filters, mode converters, couplers, second-harmonic generators, deflectors, modulators, and transducers in the fields of integrated optics and integrated surface acoustics. Work on insect compound eyes, mechanical structures, ocean waves, pulse compressions, temperature waves, and cholesteric liquid crystals, and particles interaction with crystals is briefly reviewed, especially in the case of zeolite crystals and superlattices. Recent advances in fabrication techniques for very fine gratings are also covered.

Elachi, C.↗

A study of the applicability of nucleation theory to quasi-thermodynamic transitions of second and higher Ehrenfest-order, supplement 3

The work includes an investigation of the applicability of the nucleation theory to second and higher order thermodynamic transitions in the Ehrenfest sense, and a number of significant conclusions relevant to first order transitions, as well. The underlying theoretical method consisted of expanding the Gibbs' free energy in a Maclarin or Taylor series and then using fundamental thermodynamic determinable quantities, and interpreting the results. Work was performed on the existence and interpretation of an interfacial energy between phases in a second order transition in addition to an investigation of the solid-liquid interfacial energy for various polymers. Extensive considerations were devoted to various aspects of a particular polymer, polyvinylidene fluoride (PVDF or PVF2), including an experimetal investigation of the effects of an applied electric field on the morphology of melt crystallization and on the nucleation and growth of polarized domains.

Barker, R. E., Jr.↗

Gravitational effects in dendritic growth

The theories of diffusion-controlled dendritic crystallization will be reviewed briefly, along with recently published critical experiments on the kinetics and morphology of dendritic growth in pure substances. The influence of the gravitational body force on dendrite growth kinetics will be shown to be highly dependent on the growth orientation with respect to the gravity vector and on the level of the thermal supercooling. In fact, an abrupt transition occurs at a critical supercooling, above which diffusional transport dominates the growth process and below which convective transport dominates. Our most recent work on binary mixtures shows that dilute solute additions influence the crystallization process indirectly, by altering the interfacial stability, rather than by directly affecting the transport mode. Directions for future studies in this field will also be discussed.

Glicksman, M. E.↗

Radiative Properties of Cirrus Clouds in the Infrared (8-13 microns) Spectral Region

Atmospheric radiation in the infrared (IR) 8-13 microns spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 microns. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 micron and 10000 microns over wavelengths ranging from 8 microns to 13 microns. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 microns. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8-13 microns spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, we have computed the extinction and absorption efficiency for hexagonal ice crystals with sizes ranging from 1 to 10000 microns at 12 wavelengths between 8 and 13 microns Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Parameterization of these bulk scattering properties is carried out by using second-order polynomial functions for the extinction efficiency and the single-scattering albedo and the power law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 microns whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 microns. For effective sizes larger than 100 microns, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing applications involving IR window bands under cirrus cloud conditions.

Yang, Ping↗

Toward Aerosol/Cloud Condensation Nuclei (CCN) Closure during CRYSTAL-FACE

During July 2002, measurements of cloud condensation nuclei were made in the vicinity of southwest Florida as part of the Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) field campaign. These observations, at supersaturations of 0.2 and 0.85%, are presented here. The performance of each of the two CCN counters was validated through laboratory calibration and an in situ intercomparison. The measurements indicate that the aerosol sampled during the campaign was predominantly marine in character: the median concentrations were 233 cm-3 (at S = 0.2%) and 371 cm(sup -3) (at S = 0.85%). Three flights during the experiment differed from this general trend; the aerosol sampled during the two flights on 18 July was more continental in character, and the observations on 28 July indicate high spatial variability and periods of very high aerosol concentrations. This study also includes a simplified aerosol/CCN closure analysis. Aerosol size distributions were measured simultaneously with the CCN observations, and these data are used to predict a CCN concentration using Kohler theory. For the purpose of this analysis, an idealized composition of pure ammonium sulfate was assumed. The analysis indicates that in this case, there was good general agreement between the predicted and observed CCN concentrations: at S = 0.2%, N(sub predicted)/N(sub observed)= 1.047 (R(sup 2)= 0.911)); at S = 0.85%, N(sub predicted)/N(sub observed)=1.201 (R(sup 2)= 0.835)). The impacts of the compositional assumption and of including in-cloud data in the analysis are addressed. The effect of removing the data from the 28 July flight is also examined; doing so improves the result of the closure analysis at S = 0.85%. When omitting that atypical flight, N(sub predicted)/N(sub observed) = 1.085 (R(sup 2) = 0.770) at S = 0.85%.

VanReken, Timothy M.↗

Molecular flux measurements in the back flow region of a nozzle plume

A series of tests were conducted to measure the mass flux in the far field of a nozzle plume in a high vacuum with emphasis on the back flow region. The measurements presented provided fairly accurate data for off-axis angles as large as 140 deg (i.e., in the back flow region). This region, since it is well behind the exit plane, is of paticular interest to those concerned with instrument contamination. Usually sensitive spacecraft surfaces are located in the region affected by the back flow. Parameters such as expansion ratio, throat diameter, nozzle lip shape, and plenum (chamber) pressure were varied, carbon dioxide and nitrogen gases were flowed and mass flux measurements were taken using quartz crystal microbalances in as many as nine different locations relative to the tests nozzle. Several conclusions with respect to the effect of nozzle and gas parameters on the amount of back flow mass flux are offered, and it was demonstrated that gaseous mass fluxes, which are not predictable by present theories, are encountered in the region behind the nozzle exit plane. This knowledge is significant if materials incompatible with the gaseous exhaust products are used in this region.

Chirivella, J. E.↗

Notes on Interface Growth Kinetics 50 Years after Burton, Cabrera and Frank

This overview is devoted to some unresolved basic problems in crystal growth kinetics. The density wave approach to propagation of a spatially diffuse interface between a growing crystal and its simple (e.g., metallic) melt is discussed is Section 2. This approach allows for the calculation of kinetic coefficients and is an alternative to the localized interface concept in which each atom belongs to either a solid or a liquid. Sections 3 and 4 deals mainly with layer growth from solution. Mutual retardation of the growth steps via their bulk and d a c e diffusion fields is the major subject. The influence of solution flow on step bunching (Section 4) suggests the essential influence of bulk diffusion on the surface morphology. The flow within the solution boundary layer enhances step-step interaction, influences the step bunching process and the resulting step pattern morphology on the growing surface. Recent experiments on the rates at which strongly polygonized steps on protein and small molecule crystals propagate during growth from solution are analyzed in Section 5 . We have shown that the step segments may be "singular" and that "one-dimensional nucleation" may be the rate limiting stage for the segments that are shorter or comparable in length to the thermodynamically equilibrium interlink distance. In this case, the reciprocal dependence of the segment propagation rate on the segment length that follow from the Gibbs-Thompson law, should be replaced by an abrupt switch from zero to a finite constant velocity. Until recently, the Kossel crystal remained the only model used in crystal growth theory. In such Kossel Gibbs-Thomson law, should be replaced by an abrupt switch &om zero to a finite constant velocity. crystals, all kinks at the steps are identical and the kink rate is a linear function of the supersaturation. In the non-Kossel crystals, there may be several kink configurations characterized by different geometries and energies. These configurations should appear in a specific sequence when each new lattice unit cell is filled. As a result of such a cooperative interaction within the unit cell, a non-linear dependence of the kink rate on the vapor pressure or solution concentration in excess over the equilibrium value should be expected.

Chernov, A. A.↗

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

For scattering calculations involving nonspherical particles such as ice crystals, we show that the transverse wave condition is not applicable to the refracted electromagnetic wave in the context of geometric optics when absorption is involved. Either the TM wave condition (i.e., where the magnetic field of the refracted wave is transverse with respect to the wave direction) or the TE wave condition (i.e., where the electric field is transverse with respect to the propagating direction of the wave) may be assumed for the refracted wave in an absorbing medium to locally satisfy the electromagnetic boundary condition in the ray tracing calculation. The wave mode assumed for the refracted wave affects both the reflection and refraction coefficients. As a result, a nonunique solution for these coefficients is derived from the electromagnetic boundary condition. In this study we have identified the appropriate solution for the Fresnel reflection/refraction coefficients in light scattering calculation based on the ray tracing technique. We present the 3 x 2 refraction or transmission matrix that completely accounts for the inhomogeneity of the refracted wave in an absorbing medium. Using the Fresnel coefficients for an absorbing medium, we derive an asymptotic solution in an analytical format for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations. The asymptotic theory can produce reasonable accuracy in the phase function calculations in the infrared window region (wavelengths near 10 micron) if the particle size (in diameter) is on the order of 40 micron or larger. However, since strong absorption is assumed in the computation of the single-scattering albedo in the asymptotic theory, the single scattering albedo does not change with variation of the particle size. As a result, the asymptotic theory can lead to substantial errors in the computation of single-scattering albedo for small and moderate particle sizes. However, from comparison of the asymptotic results with the FDTD solution, it is expected that a convergence between the FDTD results and the asymptotic theory results can be reached when the particle size approaches 200 micron. We show that the phase function at side-scattering and backscattering angles is insensitive to particle shape if the random orientation condition is assumed. However, if preferred orientations are assumed for particles, the phase function has a strong dependence on scattering azimuthal angle. The single-scattering albedo also shows very strong dependence on the inclination angle of incident radiation with respect to the rotating axis for the preferred particle orientations.

Yang, Ping↗

A Cryogenic Magnetostrictive Actuator Using a Persistent High Temperature Superconducting Magnet: Concept and Design - Part 1

Cryogenic magnetostrictive materials, such as rare earth zinc crystals, offer high strains and high forces with minimally applied magnetic fields, making the material ideally suited for deformable optics applications. For cryogenic temperature applications the use of superconducting magnets offer the possibility of a persistent mode of operation, i.e., the magnetostrictive material will maintain a strain field without power. High temperature superconductors (HTS) are attractive options if the temperature of operation is higher than 10 degrees Kelvin (K) and below 77 K. However, HTS wires have constraints that limit the minimum radius of winding, and even if good wires can be produced, the technology for joining superconducting wires does not exist. In this paper, the design and capabilities of a rare earth zinc magnetostrictive actuator using bulk HTS is described. Bulk superconductors can be fabricated in the sizes required with excellent superconducting properties. Equivalent permanent magnets, made with this inexpensive material, are persistent, do not require a persistent switch as in HTS wires, and can be made very small. These devices are charged using a technique which is similar to the one used for charging permanent magnets, e.g., by driving them into saturation. A small normal conducting coil can be used for charging or discharging. Because of the magnetic field capability of the superconductor material, a very small amount of superconducting magnet material is needed to actuate the rare earth zinc. In this paper, several designs of actuators using YBCO and BSCCO 2212 superconducting materials are presented. Designs that include magnetic shielding to prevent interaction between adjacent actuators will also be described. Preliminary experimental results and comparison with theory for BSCCO 2212 with a magnetostrictive element will be discussed.

Horner, Garnett↗

A Cryogenic Magnetostrictive Actuator using a Persistent High Temperature Superconducting Magnet, Part 1: Concept and Design: Concept and Design - Part 1

Cryogenic magnetostrictive materials, such as rare earth zinc crystals, offer high strains and high forces with minimally applied magnetic fields, making the material ideally suited for deformable optics applications. For cryogenic temperature applications, such as Next Generation Space Telescope (NGST), the use of superconducting magnets offer the possibility of a persistent mode of operation, i.e., the magnetostrictive material will maintain a strain field without power. High temperature superconductors (HTS) are attractive options if the temperature of operation is higher than 10 degrees Kelvin (K) and below 77 K. However, HTS wires have constraints that limit the minimum radius of winding, and even if good wires can be produced, the technology for joining superconducting wires does not exist. In this paper, the design and capabilities of a rare earth zinc magnetostrictive actuator using bulk HTS is described. Bulk superconductors can be fabricated in the sizes required with excellent superconducting properties. Equivalent permanent magnets, made with this inexpensive material, are persistent, do not require a persistent switch as in HTS wires, and can be made very small. These devices are charged using a technique which is similar to the one used for charging permanent magnets, e.g., by driving them into saturation. A small normal conducting coil can be used for charging or discharging. Very fast charging and discharging of HTS tubes, as short as 100 microseconds, has been demonstrated. Because of the magnetic field capability of the superconductor material, a very small amount of superconducting magnet material is needed to actuate the rare earth zinc. In this paper, several designs of actuators using YBCO and BSCCO 2212 superconducting materials are presented. Designs that include magnetic shielding to prevent interaction between adjacent actuators will also be described. Preliminary experimental results and comparison with theory for BSSCO 2212 with a magnetostrictive element will be discussed.

Horner, Garnett C.↗

Microscopic Foundation and Simulation of Coupled Carrier-Temperature Diffusions in Semiconductor Lasers

A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.

Li, J.↗

Influence of Nucleation Mechanisms on the Radiative Properties of Deep Convective Clouds and Subvisible Cirrus in CRYSTAL/FACE

During the past few years we have conducted work on several different topics, as reflected by our publications. As one of the Co-Project scientists for The Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL FACE) we worked to help design the mission and then conduct it in the field. Another major activity during the past two years has been to pull together various groups to formulate plans for follow on missions to CRYSTAL FACE. We organized a workshop at the University of Colorado during the summer of 2003 to assess the best locations for future missions. Working with a group of about 10 scientists from around the country we prepared a science-planning document (Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4)) that outlined the rationale, locations, strategy to accomplish the goals, and possible payloads for a set of three tropical missions. We also prepared background materials for various NRAs being prepared at NASA Headquarters for missions in Costa Rica, Darwin and Guam. In conjunction with the group at NASA Ames we have helped build a new numerical model for deep convection and have applied that model to simulate the CRYSTAL data. Our goal in particular has been to better understand how convection distributes water vapor isotopes. CRYSTAL observations of water isotopes are very different from those suggested by previous workers who assumed the isotopes would obey Rayleigh fractionation. The water isotope study has several implications. First it is a check on the realism of the deep convection model. Second, the isotopes are a measure of the precipitation removal in the atmosphere. Hence they provide a constraint on a parameter that is difficult to otherwise measure. Finally it has been suggested that isotopes may be the key to unraveling the water transport into the stratosphere and upper troposphere. Such transport is critical both for the radiation balance and for stratospheric chemistry. Ours is the first model that is able to treat this transport. Our initial results have just been submitted to Geophys. Res. Lett (Smith et al., 2005). Essentially we are able to explain the vertical profiles of isotopes in the tropical tropopause transition layer. We are also able to account for stratospheric humidity and isotope abundances with this model. We have also been heavily involved in trying to improve our understanding of nitric acid condensation on ice. Gao et al (2004) have shown that water supersaturations above ice occur when the atmosphere is supersaturated with respect to nitric acid trihydrate. As one of the co-authors of that work, we suggested the mechanism that may explain why this is occurring. Essentially, ice does not like to grow near unit supersaturation, but does so because the water molecules can find sites on the ice surface to attach themselves to before they fly off the ice surface. This phenomena was well known in the 1960s when it was a source of debate about whether condensation and evaporation coefficients for ice would be the same. Evaporation does not require any molecular orientation, while condensation does, so it was possible that the coefficients would differ. They don't differ because the water molecules rapidly move across the surface and find places to attach. Nitric acid may be occupying these preferred sites and therefore the water molecules can't find a desirable place to attach. We anticipate that this research will be the subject of laboratory work during the coming few years. Another possibility that has been suggested is that cubic ice is forming in clouds. We have measured the vapor pressure of cubic ice, and plan to publish that result in the next few months. We have also been working on additional aspects of the condensation of nitric acid on ice. With Y. Kondo we studied the condensation of NOy on ice using the SOLVE data. Gamblin et al. have continued this work. The CRYSTAL NOy and HNO, groups have shown th their data can be fit using standard Langmuir isotherms as suggested in some, but not all, laboratory studies. We have found in the SOLVE data set that this is not the case. Moreover some laboratory studies show there are important kinetic effects that may be occurring in the atmosphere limiting the transfer of nitric acid to the ice. The SOLVE data seem consistent with these studies. We are currently re-analyzing the CRYSTAL data to look for these kinetic effects. There are a number of implications of these studies. One of the more interesting is that the nitric acid coating on ice can be used as a cloud clock to determine how long the cloud parcel has been in existence. We have also been involved with several laboratory studies. We have worked to improve the database on ice optical constants, which are critical for remote sensing. We have also studied the ways in which ice nucleates on clays. We suspect now that the standard theories used for depositional ice nucleation are completely incorrect. Further work will be needed to develop a new theory.

Toon, Owen B.↗