Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “AUTOCORRELATION”

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 289 records · Page 16

Anisotropy evolution of elastomeric foams during uniaxial compression measured via in-situ X-ray computed tomography

We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m 3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.

36 MATERIALS SCIENCE↗

Point model equations for neutron correlation counting: Extension of Böhnel's equations to any order

Various methods of autocorrelation neutron analysis may be used to extract information about a measurement item containing spontaneously fissioning material. The two predominant approaches being the time correlation analysis (that make use of a coincidence gate) methods of multiplicity shift register logic and Feynman sampling. The common feature is that the correlated nature of the pulse train can be described by a vector of reduced factorial multiplet rates. We call these singlets, doublets, triplets etc. Within the point reactor model the multiplet rates may be related to the properties of the item, the parameters of the detector, and basic nuclear data constants by a series of coupled algebraic equations – the so called point model equations. Solving, or inverting, the point model equations using experimental calibration model parameters is how assays of unknown items is performed. Currently only the first three multiplets are routinely used. In this work we develop the point model equations to higher order multiplets using the probability generating functions approach combined with the general derivative chain rule, the so called Faà di Bruno Formula. Explicit expression up to 5th order are provided, as well the general iterative formula to calculate any order. This study represents the first necessary step towards determining if higher order multiplets can add value to nondestructive measurement practice for nuclear materials control and accountancy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The effect of deadtime and electronic transients on the predelay bias in neutron coincidence counting

In neutron coincidence counting using the shift register autocorrelation technique, a predelay is inserted before the opening of the (R+A)-gate. Operationally the purpose of the predelay is to ensure that the (R+A)- and A-gates have matched effectiveness, otherwise a bias will result when the difference between the gates is used to calculate the accidentals corrected net reals coincidence rate. The necessity for the predelay was established experimentally in the early practical development and deployment of the coincidence counting method. The choice of predelay for a given detection system is usually made experimentally, but even today long standing traditional values (e.g., 4.5 µs) are often used. This, at least in part, reflects the fact that a deep understanding of why a finite predelay setting is needed and how to control the underlying influences has not been fully worked out. We attempt, in this paper, to gain some insight into the problem. One aspect we consider is the slowing down, thermalization, and diffusion of neutrons in the detector moderator. The other is the influence of deadtime and electronic transients. These may be classified as non-ideal detector behaviors because they are not included in the conventional model used to interpret measurement data. From improved understanding of the effect of deadtime and electronic transients on the predelay bias in neutron coincidence counting, the performance of both future and current coincidence counters may be improved.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimating the Effective System Dead Time Parameter for Correlated Neutron Counting

We present that neutron time correlation analysis is one of the main technical nuclear safeguards techniques used to verify declarations of, or to independently assay, special nuclear materials. Quantitative information is generally extracted from the neutron-event pulse train, collected from moderated assemblies of 3 He proportional counters, in the form of correlated count rates that are derived from event-triggered coincidence gates. These count rates, most commonly referred to as singles, doubles and triples rates etc., when extracted using shift-register autocorrelation logic, are related to the reduced factorial moments of the time correlated clusters of neutrons emerging from the measurement items. Correcting these various rates for dead time losses has received considerable attention recently. The dead time losses for the higher moments in particular, and especially for large mass (high rate and highly multiplying) items, can be significant. Consequently, even in thoughtfully designed systems, accurate dead time treatments are needed if biased mass determinations are to be avoided. In support of this effort, in this paper we discuss a new approach to experimentally estimate the effective system dead time of neutron coincidence counting systems. It involves counting a random neutron source (e.g. AmLi is a good approximation to a source without correlated emission) and relating the second and higher moments of the neutron number distribution recorded in random triggered interrogation coincidence gates to the effective value of dead time parameter. We develop the theoretical basis of the method and apply it to the Oak Ridge Large Volume Active Well Coincidence Counter using sealed AmLi radionuclide neutron sources and standard multiplicity shift register electronics. The method is simple to apply compared to the predominant present approach which involves using a set of 252 Cf sources of wide emission rate, it gives excellent precision in a conveniently short time, and it yields consistent results as a function of the order of the moment used to extract the dead time parameter. In addition, this latter observation is reassuring in that it suggests the assumptions underpinning the theoretical analysis are fit for practical application purposes. However, we found that the effective dead time parameter obtained is not constant, as might be expected for a parameter that in the dead time model is characteristic of the detector system, but rather, varies systematically with gate width.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

On detrending stream velocity time series for robust tidal flow turbulence characterization

We investigated the impact of detrending techniques on turbulence quantities from tidal stream flow data, focusing on the autocorrelation function, $ρ_{uu}$, and velocity spectrum, $Φ(f)$. Standard detrending methods, including high-pass frequency-based and polynomial-based techniques, are examined, alongside a proposed alternative method, the empirical mode decomposition (EMD). Our results highlight that intervals of flow acceleration and deceleration, typical in tidal and riverine flows, significantly affect the estimation of turbulence quantities using high-pass frequency filtering and polynomial detrending of varying orders. These methods can strongly influence $ρ_{uu}$ and $Φ(f)$, thereby affecting the accurate estimation of derived quantities. We examine two variations of detrending data using EMD; the first removes only the EMD residue, and the second removes both the residue and the largest scale intrinsic mode function (IMF). By comparing the detrended spectra with the modeled von Kármán spectra, we demonstrate that the second variation (i.e., removing the residue and the largest scale IMF) successfully removed the large-scale trend of the data while retaining the energy of other scales.

Detrend↗

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗

Accounting for erroneous model structures in biokinetic process models

In engineering practice, model-based design requires not only a good process-based model, but also a good description of stochastic disturbances and measurement errors to learn credible parameter values from observations. However, typical methods use Gaussian error models, which often cannot describe the complex temporal patterns of residuals. Consequently, this results in overconfidence in the identified parameters and, in turn, optimistic reactor designs. Here, we assess the strengths and weaknesses of a method to statistically describe these patterns with autocorrelated error models. This method produces increased widths of the credible prediction intervals following the inclusion of the bias term, in turn leading to more conservative design choices. However, we also show that the augmented error model is not a universal tool, as its application cannot guarantee the desired reliability of the resulting wastewater reactor design.

42 ENGINEERING↗

A Bayesian model for multivariate discrete data using spatial and expert information with application to inferring building attributes

When modeling sparsely observed multivariate data, strong prior information elicited from experts can be used to bolster predictive accuracy and counteract sampling bias. Similarly, modeling autocorrelation in space can help make use of co-occurrence patterns present in many types of spatial data. To make use of both expert prior information and spatial structure, we propose a novel graphical model for a spatial Bayesian network developed specifically to address challenges in inferring the attributes of buildings from geographically sparse observational data. This model is implemented as the sum of a spatial multivariate Gaussian random field and a tabular conditional probability function in real-valued space prior to projection onto the probability simplex. This modeling form is especially suitable for the usage of prior information in the form of sets of atomic rules obtained from experts. To perform inference with missing data, we implement a Markov chain Monte Carlo scheme composed of alternating steps of Gibbs sampling of missing entries and Hamiltonian Monte Carlo for model parameters. A case study in building attribution is presented to highlight the advantages and limitations of this approach.

97 MATHEMATICS AND COMPUTING↗

Understanding electric vehicle ownership using data fusion and spatial modeling

The global shift toward electric vehicles (EVs) for climate sustainability lacks comprehensive insights into the impact of the built environment on EV ownership, especially in varying spatial contexts. This study, focusing on New York State, integrates data fusion techniques across diverse datasets to examine the influence of socioeconomic and built environmental factors on EV ownership. The utilization of spatial regression models reveals consistent coefficient values, highlighting the robustness of the results, with the Spatial Lag model better at capturing spatial autocorrelation. Further, results underscore the significance of charging stations within a 10-mile radius, indicative of a preference for convenient charging options influencing EV ownership decisions. Factors like higher education levels, lower rental populations, and concentrations of older population align with increased EV ownership. Utilizing publicly available data offers a more accessible avenue for understanding EV ownership across regions, complementing traditional survey approaches.

33 ADVANCED PROPULSION SYSTEMS↗

Extended live-tracking and quantitative characterization of wound healing and cell migration with SiR-Hoechst

Highlights: • SiR-Hoechst is a nuclear dye that enables cell tracking without affecting migration. • Single cell analysis provides multiple migration metrics during wound healing. • Cell density-based report of wound healing is less susceptible to subjective bias. • Cell directionality is an important contributor of early wound recovery. Cell migration is essential to many life processes, including immune response, tissue repair, and cancer progression. A reliable quantitative characterization of the cell migration can therefore aid in the high throughput screening of drug efficacy in wound healing and cancer treatments. In this work, we report what we believe is the first use of SiR-Hoechst for extended live tracking and automated analysis of cell migration and wound healing. We showed through rigorous statistical comparisons that this far-red label does not affect migratory behavior. We observed excellent automated tracking of random cell migration, in which the motility parameters (speed, displacement, path length, directionality ratio, persistence time, and direction autocorrelation) obtained closely match those obtained from manual tracking. We also present an analysis framework to characterize the healing of a scratch wound from the perspective of single cells. The use of SiR-Hoechst is advantageous for the crowded environments in wound healing assays because as long as cell nuclei do not overlap, continuous tracking can be maintained even if there is cell-cell contact. In this paper, we report wound recovery based on the number of cells migrating into the wound over time, normalized by the initial cell count prior to the infliction of the wound. This normalized cell count approach is impervious to operator bias during the arbitration of wound edges and is also robust against variability that arises due to differences in the cell density of different samples. Additional wound healing characteristics were also defined based on the evolution of cell speed and directionality during healing. Not unexpected, the wound healing cells exhibited much higher tendency to maintain the same migratory direction in comparison to the randomly migrating cells. The use of SiR-Hoechst thus greatly simplified the automation of single cell and whole population analysis with high spatial and temporal resolution over extended periods of time.

60 APPLIED LIFE SCIENCES↗

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING↗

A method for mapping submicron-scale crystallographic order/disorder applied to human tooth enamel

Tooth enamel, the outermost layer of human teeth, is a complex, hierarchically structured biocomposite. The details of this structure are important in multiple human health contexts, from understanding the progression of dental caries (tooth decay) to understanding the process of amelogenesis and related developmental defects. Enamel is composed primarily of long, nanoscale crystallites of hydroxyapatite that are bundled by the thousands to form micron-scale rods. Studies with transmission electron microscopy show the relationships between small groups of crystallites and X-ray diffraction characterize averages over many rods, but the direct measurement of variations in local crystallographic structure across and between enamel rods has been missing. Here, we describe a synchrotron X-ray-based experimental approach and a novel analysis method developed to address this gap in knowledge. A ~500-nm-wide beam of monochromatic X-rays in conjunction with a sample section only 1 μm in thickness enables 2D diffraction patterns to be collected from small well-separated volumes within the enamel microstructure but still probes enough crystallites (~300 per pattern) to extract population-level statistics on crystallographic features like lattice parameter, crystallite size, and orientation distributions. Furthermore, the development of a quantitative metric to characterize relative order and disorder based on the azimuthal autocorrelation of diffracted intensity enables these crystallographic measurements to be correlated with their location within the enamel microstructure (e.g., between rod and interrod regions). These methods represent a step forward in the characterization of human enamel and will elucidate the variation of the crystallographic structure across and between enamel rods for the first time.

crystallographic order/disorder↗

Localized Phonon Densities of States at Grain Boundaries in Silicon

Since it is now possible to record vibrational spectra at nanometer scales in the electron microscope, it is of interest to explore whether extended defects in crystals such as dislocations or grain boundaries will result in measurable changes of the phonon densities of states (dos) that are reflected in the spectra. Phonon densities of states were calculated for a set of high angle grain boundaries in silicon. The boundaries are modeled by supercells with up to 160 atoms, and the vibrational densities of states were calculated by taking the Fourier transform of the velocity–velocity autocorrelation function from molecular dynamics simulations with larger supercells doubled in all three directions. In selected cases, the results were checked on the original supercells by comparison with the densities of states obtained by diagonalizing the dynamical matrix calculated using density functional theory. Near the core of the grain boundary, the height of the optic phonon peak in the dos at 60 meV was suppressed relative to features due to acoustic phonons that are largely unchanged relative to their bulk values. This can be attributed to the variation in the strength of bonds in grain boundary core regions where there is a range of bond lengths.

Materials Science↗

Strategies and Software for Machine Learning Accelerated Discovery in Transition Metal Chemistry

Machine learning the electronic structure of open shell transition metal complexes presents unique challenges, including robust and automated data set generation. In this report we introduce tools that simplify data acquisition from density functional theory (DFT) and validation of trained machine learning models using the molSimplify automatic design (mAD) workflow. We demonstrate this workflow by training and comparing the performance of LASSO, kernel ridge regression (KRR), and artificial neural network (ANN) models using heuristic, topological revised autocorrelation (RAC) descriptors we have recently introduced for machine learning inorganic chemistry. On a series of open shell transition metal complexes, we evaluate set aside test errors of these models for predicting the HOMO level and HOMO–LUMO gap. The best performing models are ANNs, which show 0.15 and 0.25 eV test set mean absolute errors on the HOMO level and HOMO–LUMO gap, respectively. Poor performing KRR models using the full 153-feature RAC set are improved to nearly the same performance as the ANNs when trained on down-selected subsets of 20–30 features. Analysis of the essential descriptors for HOMO level and HOMO–LUMO gap prediction as well as comparison to subsets previously obtained for other properties reveal the paramount importance of nonlocal, steric properties in determining frontier molecular orbital energetics. We demonstrate our model performance on diverse complexes and in the discovery of molecules with target HOMO–LUMO gaps from a large 15,000 molecule design space in minutes rather than days that full DFT evaluation would require.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MD Simulation of Water Using a Rigid Body Description Requires a Small Time Step to Ensure Equipartition

In simulations of aqueous systems, it is common to freeze the bond vibration and angle bending modes in water to allow for a longer time step δt for integrating the equations of motion. Thus, δt = 2 fs is often used in simulating rigid models of water. We simulate the SPC/E model of water using δt from 0.5 to 3.0 fs and up to 4 fs using hydrogen mass repartitioning. In these simulations, we find that for all but δt = 0.5 fs, equipartition is not obtained between translational and rotational modes, with the rotational modes exhibiting a lower temperature than the translation modes. To probe the reasons for the lack of equipartition, we study the autocorrelation of the translational velocity of the center of mass and the angular velocity of the rigid water molecule, respectively. Here, we find that the rotational relaxation occurs on a timescale comparable to vibrational periods, calling into question the original motivations for freezing the vibrations. Furthermore, a time step with δt ≥ 1 fs is not able to capture accurately the fast rotational relaxation, which reveals its impact as an effective slowing-down of rotational relaxation. The fluctuation–dissipation relation then leads to the conclusion that the rotational temperature should be cooler for δt greater than the reference value of 0.5 fs. Consideration of fluctuation–dissipation in equilibrium molecular dynamics simulations also emphasizes the need to capture the temporal evolution of fluctuations with fidelity and the role of δt in this regard. The time step also influences the solution thermodynamic properties: both the mean system potential energies and the excess entropy of hydration of a soft repulsive cavity are sensitive to δt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine Learning Prediction of the Experimental Transition Temperature of Fe(II) Spin-Crossover Complexes

Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. Here, we leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys. 159, 024120 (2023)] to train three machine learning (ML) models for transition temperature (T 1/2 ) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T 1/2 values. We then compare ML T 1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T 1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T 1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T 1/2 values with a Pearson’s r of 0.82. In contrast, DFA-predicted T 1/2 values have large errors and show no correlation to estimated T 1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adsorption and Solvation Modulate Rhodamine B Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Confinement of solvents and solutes within nanoporous materials frequently leads to the emergence of unique mass transport behaviors that, once fully understood, may lead to improved chemical separations. Here, the diffusion of Rhodamine B (RhB) dye within 10 and 20 nm diameter anodic aluminum oxide (AAO) nanopores filled with binary ethanol/water mixtures is investigated. Mixture compositions spanning from pure ethanol to pure water are employed. The results of confocal fluorescence correlation spectroscopy studies reveal that RhB diffusion occurs by a two-component mechanism comprising composition-dependent fast and slow motions, characterized by diffusion coefficients D f and D s . The results are consistent with those of previous studies performed under more limited conditions [J. Phys. Chem. C, 2023, 127, 411-420]. The fast component scales with mixture viscosity and is assigned to hindered bulk-like diffusion in central pore regions. Slow diffusion likely involves adsorption of RhB to the pore surface and may be described by a desorption mediated mechanism. The occurrence of RhB adsorption to the AAO surface is verified at the single-molecule level by wide-field fluorescence imaging of membrane cross-sectional surfaces. Unique composition dependent trends in the autocorrelation amplitude and in D s that mimic bulk RhB solubility are revealed. D s is found to be smallest in pure ethanol and pure water and largest in intermediate mixtures. These results suggest that RhB surface adsorption is strongest in the pure liquids and weakest in mixtures of intermediate composition, where the dye is least soluble, and most soluble, respectively. As a result, molecular dynamics simulations reveal that a water layer appears on the pore surface under most conditions, while RhB is solvated primarily by ethanol. The composition dependence of RhB diffusion is concluded to reflect its solvation dependent interactions with the pore walls.

36 MATERIALS SCIENCE↗