Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electrostatic Potential”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Electron Density and Two-Channel Neutron Emission Measurements in Steady-State Spherical Inertial-Electrostatically Confined Plasmas, with Review of the 1-D Kinetic Model

Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the 1-D Poisson equation for EC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with non-fusing hydrogen-1. The deuterium rates are consistent with predictions from the model.

Dobson, Chris C.↗
DOE OSTI2025

Gating Single Molecules with Counterions

We report atomic-scale gating and visualization of local charge distribution within individual rare-earth-based molecular complexes on a metallic surface. The complexes are formed by a positively charged lanthanum ion coordinated to a (pcam)3 molecule and a negatively charged counterion trapped underneath via electrostatic interactions on a Au(111) surface. Local gating is performed by adding an additional negatively charged counterion to one side of the complex, which results in the redistribution of charges within the complex and a positive shift of the frontier orbitals. This is caused by the internal Stark effect induced by the added counterion. This effect is directly captured using tunneling spectroscopy and spectroscopic mapping at 5 K substrate temperature. The polarizability of the complex is corroborated by density functional theory and analytical calculations based on experimental findings. Furthermore, the influence of charge polarization on nearby complexes is investigated in a cluster purposely assembled using three complexes, which reveals maintaining the charge states as in single complexes. These findings will enable the design of robust charged rare-earth complexes to be tailored for potential solid-state applications.

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

Evidence of Au(111) topological states in a kagome analogue lattice and their robustness beyond ultra-low temperatures and defect-free conditions

The experimental realization of kagome lattices that exhibit the predicted coexistence of topological states with high electron kinetics and non-dispersing quantum states remains challenging. Additionally, the robustness of these states against structural perturbations has rarely been explored. Here, we report on the formation of an analogue kagome structure via the electrostatic self-assembly of 4,7-dibromobenzo[c]-1,2,5-thiadiazole (2Br-BTD) molecules on Au(111). Local spectroscopic measurements, supported by theoretical calculations, reveal that the weak molecular coupling reshapes the topological-induced Shockley surface state of Au(111) by imposing a (7 × 7) periodicity, resulting in new band crossings. The molecular overlayer favours the opening of electron gaps at these positions manifested as sharp peaks in dI/dV spectra and electron localization in either the hexagonal or triangular sublattices of the kagome structure. To explore the robustness of these topological states, we monitored their stability under varied conditions, including different temperatures, the unaltered herringbone reconstruction of the Au(111) surface and local structural relaxation of the molecular assembly. These results demonstrate the degree of topological protection of these states, which holds potential for fundamental and applied research.

Austin, Dave [University of Central Florida, Orlan↗
NASA NTRS2004

Electron Density and Two-Channel Neutron Emission Measurements in Steady-State Spherical Inertial-Electrostatically Confined Plasmas, with Review of the One-Dimensional Kinetic Model

Electron density measurements have been made in steady-state plasmas in a spherical inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with eccentric symmetric minima, and give mean densities of from approx. = 0.4 to 7 x 10(exp 10)/cu cm, the density generally increasing with the neutral gas pressure. Numerical solutions of the one-dimensional Poisson equation for IEC plasmas are reviewed and energy distribution functions are identified which give flat transverse profiles. These functions are used with the plasma approximation to obtain solutions which also give densities consistent with the measurements, and a double potential well solution is obtained which has minima qualitatively similar to those observed. Explicit consideration is given to the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions during operation with nonfusing hydrogen- 1. The deuterium rates are consistent with predictions from the model.

Dobson, Chris C.↗
NASA NTRS2012

Solar Wind Electron Interaction with the Dayside Lunar Surface and Crustal Magnetic Fields: Evidence for Precursor Effects

Electron distributions measured by Lunar Prospector above the dayside lunar surface in the solar wind often have an energy dependent loss cone, inconsistent with adiabatic magnetic reflection. Energy dependent reflection suggests the presence of downward parallel electric fields below the spacecraft, possibly indicating the presence of a standing electrostatic structure. Many electron distributions contain apparent low energy (<100 eV) upwardgoing conics (58% of the time) and beams (12% of the time), primarily in regions with non-zero crustal magnetic fields, implying the presence of parallel electric fields and/or wave-particle interactions below the spacecraft. Some, but not all, of the observed energy dependence comes from the energy gained during reflection from a moving obstacle; correctly characterizing electron reflection requires the use of the proper reference frame. Nonadiabatic reflection may also play a role, but cannot fully explain observations. In cases with upward-going beams, we observe partial isotropization of incoming solar wind electrons, possibly indicating streaming and/or whistler instabilities. The Moon may therefore influence solar wind plasma well upstream from its surface. Magnetic anomaly interactions and/or non-monotonic near surface potentials provide the most likely candidates to produce the observed precursor effects, which may help ensure quasi-neutrality upstream from the Moon.

Solar↗
NASA NTRS1999

Molecular Modeling and Experimental Investigations of Nonlinear Optical Compounds Monosubstituted Derivatives of Dicyanovinylbenzene

A search for potential nonlinear optical compounds was performed using the Cambridge Structure Database and molecular modeling. We investigated a series of monosubstituted derivatives of dicyanovinylbenzene, since the nonlinear optical (NLO) properties of such derivatives (o-methoxy-dicyanovinylbenzene, DIVA) were studied earlier. The molecular geometry of these compounds was investigated with x-ray analysis and discussed along with the results of molecular mechanics and ab initio quantum chemical calculations. The influence of crystal packing on the planarity of the molecules of this series has been revealed. Two new compounds from the series studied, ortho-F and para-Cl-dicyanovinylbenzene, according to powder measurements, were found to be NLO compounds in the crystal state about 10 times more active than urea. The peculiarities of crystal structure formation in the framework of balance between van der Waals and electrostatic interactions have been discussed. The crystal shape of DIVA and two new NLO compounds have been calculated on the basis of the known crystal structure.

Timofeeva, Tatiana V.↗
DOE OSTI2022

Interplay between topological valley and quantum Hall edge transport

An established way of realising topologically protected states in a two-dimensional electron gas is by applying a perpendicular magnetic field thus creating quantum Hall edge channels. In electrostatically gapped bilayer graphene intriguingly, even in the absence of a magnetic field, topologically protected electronic states can emerge at naturally occurring stacking domain walls. While individually both types of topologically protected states have been investigated, their intriguing interplay remains poorly understood. Here, we focus on the interplay between topological domain wall states and quantum Hall edge transport within the eight-fold degenerate zeroth Landau level of high-quality suspended bilayer graphene. We find that the two-terminal conductance remains approximately constant for low magnetic fields throughout the distinct quantum Hall states since the conduction channels are traded between domain wall and device edges. For high magnetic fields, however, we observe evidence of transport suppression at the domain wall, which can be attributed to the emergence of spectral minigaps. This indicates that stacking domain walls potentially do not correspond to a topological domain wall in the order parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Field-aligned auroral electron fluxes

Out of a series of five auroral sounding rocket flights, electron fluxes in three of the flights were observed to be sharply field-aligned. The electrons were detected with an array of electrostatic analyzers making 50-point spectral samples from 0 to 15 keV. The field-aligned fluxes resulted in a peak in the spectrum at a lower energy than that of another monoenergetic peak in the spectrum seen consistently when it was conjugate to bright auroral forms giving multiply peaked spectra. The field-aligned peak in the spectrum is apparently a phenomenon distinct from the isotropic peak and was observed to occur at various times in the energy range from a few hundred electron volts to 5 keV. The field-aligned electrons were highly collimated, having pitch angles less than 20 deg and at times less than 10 deg. The high degree of field alignment and the monoenergetic character support the model that these events can be attributed to relatively cool plasma falling through an electric potential at altitudes not too far removed from the rocket.

Arnoldy, R. L.↗
NASA NTRSDate not supplied

Materials for Power Transmission in Extreme Environments

NASA’s goal for a sustained human presence on the Moon and a future Lunar economy requires the development of novel power transmission materials for extreme environments. As an increasing number of public and private entities set their sights towards a future on the Moon, the demand for surface power will dramatically increase – potentially up to megawatt levels. A Lunar economy is envisioned to include a human habitat, manned and unmanned science exploration assets, excavation rovers and systems, and large-scale infrastructure construction such as solar arrays and power grids. A NASA trade study estimated that 1 to 3 km of power transmission lines (1 to 3 KV ac and dc) will be needed over time. Current high power and voltage SOA cables/wire can easily exceed 1 kg/m where mass is not as critical. Connectors, deployment spools and robots further increase the subsystem mass significantly increasing costs. On the other hand, advances in aeronautics provides path towards space missions. The operating environment on the Moon presents substantial risk to safe and reliable power transmission due to space weathering and degradation. These challenges include cosmic and UV radiation, extreme temperatures swings, the electrostatic Lunar exosphere, and the electrostatic and abrasive lunar regolith. The critical need for lightweight, safe and reliable power cables/wires materials is evident. In this presentation, the extreme high altitude and lunar environments and recent materials development from NASA Glenn Research Center’s High Voltage Materials team towards meeting NASA’s objectives and goals for aeronautics and space missions.

Cu/CNT conductors↗
DOE OSTI2022

2D materials-assisted heterogeneous integration of semiconductor membranes toward functional devices

Heterogeneous integration techniques allow the coupling of highly lattice-mismatched solid-state membranes, including semiconductors, oxides, and two-dimensional materials, to synergistically fuse the functionalities. The formation of heterostructures generally requires two processes: the combination of crystalline growth and a non-destructive lift-off/transfer process enables the formation of high-quality heterostructures. Although direct atomic interaction between the substrate and the target membrane ensures high-quality growth, the strong atomic bonds at the substrate/epitaxial film interface hinder the non-destructive separation of the target membrane from the substrate. Alternatively, a 2D material-coated compound semiconductor substrate can transfer the weakened (but still effective) surface potential field of the surface through the 2D material, allowing both high-quality epitaxial growth and non-destructive lift-off of the grown film. This Perspective reviews 2D/3D heterogeneous integration techniques, along with applications of III–V compound semiconductors and oxides. Here, the advanced heterogeneous integration methods offer an effective method to produce various freestanding membranes for stackable heterostructures with unique functionalities that can be applied to novel electrical, optoelectronic, neuromorphic, and bioelectronic systems.

2D materials↗
NASA NTRS1985

The hidden ion population - Revisited

In an investigation conducted by Olsen (1982) on the basis of particle data taken with an electrostatic analyzer, it was found that a cold plasma population with a density between 10 and 100 per cu cm appeared suddenly when the satellite was eclipsed, but was hidden in sunlight. The present paper has the objective to show further measurements of ordinarily 'hidden' ion populations, in order to resolve some of the questions raised in connection with the Scatha satellite data reported by Olsen. It is found that the retarding ion mass spectrometer (RIMS) detector is capable of measuring the core of the plasma distribution in sunlight and eclipse, though the task is more easily done in eclipse. There are, however, limitations concerning the ability of the detector to measure all the plasma, all the time. It is, therefore, pointed out that continuous effective measurements of the 'hidden' ion population of the magnetosphere still awaits satellites with effective means of potential control.

Olsen, R. C.↗
DOE OSTI2020

Hybrid-Supported Bilayers Formed with Mixed-Charge Surfactants on C 18 -Functionalized Silica Surfaces

Mixtures of cationic-anionic surfactants have been demonstrated to spontaneously form ordered monolayers at hydrophobic-hydrophilic boundaries including air-water and oil-water interfaces. Here, confocal Raman microscopy is used to investigate the structure of hybrid-supported surfactant bilayers (HSSBs) formed by deposition of a distal leaflet of mixed cationic-anionic surfactants onto a proximal leaflet of n-alkane (C 18 ) chains on the interior surfaces of chromatographic silica particles. The surface coverage of the two surfactants in a hybrid bilayer was determined from carbon analysis and the relative Raman scattering of their respective head-groups. Within the measurement uncertainty, the stoichiometric ratio of the two surfactants is one-to-one, equivalent to mixed-charge-surfactant monolayers at air-water and oil-water interfaces and consistent with the role of the head-group electrostatic interactions in their formation. When self-assembled on the hydrophobic surface, pairs of oppositely-charged n-alkyl-chain surfactants resemble a phospholipid (phosphati¬dylcholine) molecule, with its zwitterionic head group and two hydrophobic acyl-chain tails. Indeed, the structure of these hybrid-supported surfactant bilayers on C 18 -modified silica surfaces is similar to that of hybrid-supported lipid bilayers (HSLBs) on the same supports, but with denser and more-ordered n-alkyl chains. Hybrid-supported surfactant bilayers exhibit a melting phase transition (gel to liquid-crystalline phase) with structural and energetic characteristics similar to hybrid-supported bilayers prepared from a zwitterionic phospholipid of the same alkyl chain length. These mixed-charge surfactants on n-alkane modified silica are stable in water over time (months), results that suggest the potential use of these hybrid bilayers for generating supported lipid bilayer-like surfaces or for separation applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2016

Unique Results and Lessons Learned from the TSS Missions

In 1924, Irvin Langmuir and H. M. Mott-Smith published a theoretical model for the complex plasma sheath phenomenon in which they identified some very special cases which greatly simplified the sheath and allowed a closed solution to the problem. The most widely used application is for an electrostatic, or "Langmuir," probe in laboratory plasma. Although the Langmuir probe is physically simple (a biased wire) the theory describing its functional behavior and its current-voltage characteristic is extremely complex and, accordingly, a number of assumptions and approximations are used in the LMS model. These simplifications, correspondingly, place limits on the model's range of application. Adapting the LMS model to real-life conditions is the subject of numerous papers and dissertations. The Orbit-Motion Limited (OML) model that is widely used today is one of these adaptions that is a convenient means of calculating sheath effects. The OML equation for electron current collection by a positively biased body is simply: I is approximately equal to A x j(sub eo) x 2/√𝜋(phi)(exp ½) where A is the area of the body and phi is the electric potential on the body with respect to the plasma. Since the Langmuir probe is a simple biased wire immersed in plasma, it is particularly tempting to use the OML equation in calculating the characteristics of the long, highly biased wires of an Electric Sail in the solar wind plasma. However, in order to arrive at the OML equation, a number of additional simplifying assumptions and approximations (beyond those made by Langmuir-Mott-Smith) are necessary. The OML equation is a good approximation when all conditions are met, but it would appear that the Electric Sail problem lies outside of the limits of applicability.

Stone, Nobie H.↗
DOE OSTI2023

Electronic Structure and Mechanical Properties of Solvated Montmorillonite Clay Using Large-Scale DFT Method

Montmorillonite clay (MMT) has been widely used in engineering and environmental applications as a landfill barrier and toxic waste repository due to its unique property as an expandable clay mineral that can absorb water easily. This absorption process rendered MMT to be highly exothermic due to electrostatic interactions among molecules and hydrogen bonds between surface atoms. A detailed study of a large supercell model of structural clay enables us to predict long-term nuclear waste storage. Herein, a large solvent MMT model with 4071 atoms is studied using ab initio density functional theory. The DFT calculation and analysis clarify the important issues, such as bond strength, solvation effect, elasticity, and seismic wave velocities. These results are compared to our previous study on crystalline MMT (dry). The solvated MMT has reduced shear modulus (G), bulk modulus (K), and Young’s modulus (E). We observe that the conduction band (CB) in the density of states (DOS) of solvated MMT model has a single, conspicuous peak at –8.5 eV. Moreover, the atom-resolved partial density of states (PDOS) summarizes the roles played by each atom in the DOS. These findings illuminate numerous potential sophisticated applications of MMT clay.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Plasma simulation using the massively parallel processor

Two dimensional electrostatic simulation codes using the particle-in-cell model are developed on the Massively Parallel Processor (MPP). The conventional plasma simulation procedure that computes electric fields at particle positions by means of a gridded system is found inefficient on the MPP. The MPP simulation code is thus based on the gridless system in which particles are assigned to processing elements and electric fields are computed directly via Discrete Fourier Transform. Currently, the gridless model on the MPP in two dimensions is about nine times slower that the gridded system on the CRAY X-MP without considering I/O time. However, the gridless system on the MPP can be improved by incorporating a faster I/O between the staging memory and Array Unit and a more efficient procedure for taking floating point sums over processing elements. The initial results suggest that the parallel processors have the potential for performing large scale plasma simulations.

Lin, C. S.↗
DOE OSTI2022

Ultra-selective membrane composed of charge-stabilized fixed carrier and amino acid-based ionic liquid mobile carrier for highly efficient carbon capture

Membrane technology has been extensively studied for CO 2 capture applications, especially for flue gas sources. In the past decades, facilitated transport membranes (FTMs) have made breakthroughs overcoming the permeance-selectivity trade-off upper bound restricting traditional CO 2 separation membranes, but are still facing challenges towards practical applications, including limited performance, such as insufficient CO 2 /N 2 selectivity to achieve 95 % CO 2 dry-base purity by one step separation, and long-term stability. Herein, we designed and fabricated a novel FTM structure containing an ionic liquid (1-ethyl-3-methylimidazolium aminoacetate, [Emim][Gly]) as mobile CO 2 -carrier and a polymeric amine (polyethyleneimine, PEI) as fixed CO 2 -carrier. In this study, the fixed carrier is confined within a carbon nanotube (CNT) framework of 230 nm thickness via electrostatic forces adjusted by a polyelectrolyte (polystyrene sulfonate, PSS), while the mobile carrier diffuses freely within the CNT framework. After optimization of the membrane recipe and spray-coating fabrication procedure, following our previous work, the resulting CNT-PSS-PEI ~ IL membranes demonstrated an ultra-high CO 2 /N 2 selectivity up to 1,000 with CO 2 permeance up to 2,400 GPU (Gas Permeation Unit, 1 GPU = 3.348 × 10 -10 mol·s -1 ·m -2 ·Pa -1 ) under vacuum operation condition. Furthermore, one 100-cm 2 flat sheet membrane sample was prepared and exhibited one-stage CO 2 enrichment from 15 % to 95 % purity (dry-base) for the first time amongst all reported CO 2 separation membranes. The membrane retained a stable performance over 50-h operation period under vacuum condition. The extraordinary CO 2 separation performance illustrates the great potential of the CNT-PSS-PEI ~ IL membranes for flue gas carbon capture application.

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

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗