Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large deviations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Calibration of a calorimeter for measuring the performance of thermal control surfaces

The calibration to evaluate the heat leak in terms of the deviation from a perfect calorimeter is described. A thermal vacuum test was carried out to characterize the performance of the calorimeter. The calorimeter was equipped with a heater to control the sample disc temperature and with specific instrumentation to measure the heat leak. The radiation sink temperature of the black cavity target was varied from -192 to +31C, while the heater power was varied from 0 to 311 milliwatts. A steady state thermal model was developed to correlate the results. The calorimeter performance was characterized in terms of the heat leak as a percentage of the ideal heat flow of the calorimeter disc. Large deviations from ideal performance occur at low sink temperatures. The effect of the use of the heater is discussed. The effects of transient conditions during low Earth orbit are discussed. It is concluded that heat leak calibrations are necessary for a wide range of conditions.

Houseman, John↗
NASA NTRS2001

Coriolis-force-induced trajectory and endpoint deviations in the reaching movements of labyrinthine-defective subjects

When reaching movements are made during passive constant velocity body rotation, inertial Coriolis accelerations are generated that displace both movement paths and endpoints in their direction. These findings directly contradict equilibrium point theories of movement control. However, it has been argued that these movement errors relate to subjects sensing their body rotation through continuing vestibular activity and making corrective movements. In the present study, we evaluated the reaching movements of five labyrinthine-defective subjects (lacking both semicircular canal and otolith function) who cannot sense passive body rotation in the dark and five age-matched, normal control subjects. Each pointed 40 times in complete darkness to the location of a just extinguished visual target before, during, and after constant velocity rotation at 10 rpm in the center of a fully enclosed slow rotation room. All subjects, including the normal controls, always felt completely stationary when making their movements. During rotation, both groups initially showed large deviations of their movement paths and endpoints in the direction of the transient Coriolis forces generated by their movements. With additional per-rotation movements, both groups showed complete adaptation of movement curvature (restoration of straight-line reaches) during rotation. The labyrinthine-defective subjects, however, failed to regain fully accurate movement endpoints after 40 reaches, unlike the control subjects who did so within 11 reaches. Postrotation, both groups' movements initially had mirror image curvatures to their initial per-rotation reaches; the endpoint aftereffects were significantly different from prerotation baseline for the control subjects but not for the labyrinthine-defective subjects reflecting the smaller amount of endpoint adaptation they achieved during rotation. The labyrinthine-defective subjects' movements had significantly lower peak velocity, higher peak elevation, lower terminal velocity, and a more vertical touchdown than those of the control subjects. Thus the way their reaches terminated denied them the somatosensory contact cues necessary for full endpoint adaptation. These findings fully contradict equilibrium point theories of movement control. They emphasize the importance of contact cues in adaptive movement control and indicate that movement errors generated by Coriolis perturbations of limb movements reveal characteristics of motor planning and adaptation in both healthy and clinical populations.

NASA Discipline Neuroscience↗
DOE OSTI2021

Electric field and turbulence in global Braginskii simulations across the ASDEX Upgrade edge and scrape-off layer

Turbulence simulations in diverted geometry across the edge and scrape-off layer (SOL) of ASDEX Upgrade are performed with the GRILLIX code (Stegmeir et al 2019 Phys. Plasmas 26 052517). The underlying global (full-f) drift-reduced Braginskii model allows to concurrently study the self-consistent dynamics of the turbulence and the background as well as the evolution of toroidal and zonal flows. Different contributions to the radial electric field are identified. The dominant contribution on closed flux surfaces comes from the ion pressure gradient, due to the diamagnetic drift in the curved magnetic field. Large deviations can be induced, in particular, by the polarization particle flux, leading to zonal flows. The latter are driven by small-scale eddies, but do not exhibit much impact on the overall transport which is driven by ballooning modes at larger scales. Ion viscosity is found to be important in damping poloidal rotation through adjusting of the parallel velocity profile, but not via direct vorticity damping. The zonal flow drive peaks at the separatrix, where a strong shear layer forms due to the sheath-induced counter-propagating SOL flow, allowing for the formation of a transport barrier. The temperature profile across the separatrix is determined by the competition between cross-field transport and outflow in the SOL, the latter being largely controlled by the parallel heat conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Assessment of Model Hub Height Wind Speed Performance Using DOE Lidar Buoy Data

As the offshore wind industry continues to gain momentum both in the U.S. and globally, reanalysis products remain essential tools for the lifecycle of a wind farm. From project siting and production estimates to determining construction and maintenance windows, reanalysis products play an important role. These regional and global long-term models also serve as inputs to higher resolution models. With such wide-ranging impacts on the offshore wind industry, the need for validation of reanalysis products is strong, particularly at typical turbine hub heights. Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinelTM buoys for the U.S. Department of Energy (DOE) in order to collect meteorological and oceanographic data, including hub height wind speed, in areas of interest for offshore wind development. This work uses DOE buoy observations off the coasts of New Jersey and Virginia to evaluate the performance of Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), the European Center for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5th generation (ERA5), and the National Oceanic and Atmospheric Administration (NOAA) assimilation system Rapid Refresh (RAP). Biases and degrees of correspondence are determined for each reanalysis product in order to provide insights on the performance and uncertainty for long-term wind resource characterization. With the baseline metrics determined, this work proceeds with an investigation into the sources of large deviations between the modeled and observed hub height wind speeds during the East Coast lidar buoy deployments. Consistent sources of reanalysis model error include stably stratified flow conditions, high wind shear, flow parallel to the coastline, and tropical and winter storms. For the near shore New Jersey location, offshore winds (i.e. blowing from land to water) tend to be associated with model overestimation of observed hub height wind speed, while winds parallel to the coastline are correlated with model underestimation of wind speed. Similarly, near the Virginia buoy, during southerly winds (parallel to the coastline) models underestimate the observed wind speed due to coastal upwelling. Reanalysis model wind speeds underestimate the hub height winds with increasing stable atmospheric stratification. Reanalysis model biases are seen when the observed wind speeds approach or exceed typical turbine cut-out speeds (around 20 m s-1). Large model errors during winter and tropical storm events, high wind speed events, and ramp events are often the result of wind speed timing and magnitude offsets. Overall, ERA5 provides the most successful representation of observed offshore hub height wind speeds at the U.S. East Coast buoy locations and is therefore best suited as an input boundary condition for model case studies. ERA5 performs well for a variety of atmospheric phenomena, including storms and sea breezes, however RAP, the next most successful reanalysis model, performs best in capturing the frequency of ramp events.

17 WIND ENERGY↗
DOE OSTI2023

Evidence of Mixed Scaling for Mean Profile Similarity in the Stable Atmospheric Surface Layer

A new mixed scaling parameter Z = z/(Lh)1/2 is proposed for similarity in the stable atmospheric surface layer, where z is the height, L is the Obukhov length, and h is the boundary layer depth. In comparison with the parameter ζ = z/L from Monin–Obukhov similarity theory (MOST), the new parameter Z leads to improved mean profile similarity for wind speed and air temperature in large-eddy simulations. It also yields the same linear similarity relation for CASES-99 field measurements, including in the strongly stable (but still turbulent) regime where large deviations from MOST are observed. Results further suggest that similarity for turbulent energy dissipation rate depends on both Z and ζ. The proposed mixed scaling of Z and relevance of h can be explained by physical arguments related to the limit of z-less stratification that is reached asymptotically above the surface layer. The presented evidence and fitted similarity relations are promising, but the results and arguments are limited to a small sample of idealized stationary stable boundary layers. Here, corroboration is needed from independent datasets and analyses, including for complex and transient conditions not tested here.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

Trajectory Prediction Accuracy and Error Sources for Regional Jet Descents: Results of a 2010 Flight Trial at Denver International Airport using SkyWest Revenue Flights - Part II

The Efficient Descent Advisor (EDA) controller automation tool generates trajectory-based speed, path, and altitude-profile advisories to facilitate efficient, continuous descents into congested terminal airspace. While prior field trials have assessed the trajectory prediction accuracy for large jet (i.e., Boeing and Airbus) types, smaller (i.e., regional and business) jet types present unique challenges involving different descent procedures and Flight Management System (FMS) capabilities. This paper quantifies the trajectory prediction accuracy for small jet revenue flight descents based on SkyWest Canadair Regional Jet 200, 700, and 900 aircraft arrivals to Denver in the fall of 2010. Post flight test data analysis and SkyWest pilot interviews uncovered unexpected variation between flight crews due to different interpretations of (1) which fixed flight path angle (FPA) to fly based on the flight trial procedure and (2) how to fly the descent to achieve the target FPA. Pilot reports were used to select a subset of flights where pilots indicated an FPA according to the flight trial procedure to remove the unexpected variation due to (1) to focus on (2). Results for the subset for en-route descents, from prior to top of descent to the meter fix 30 to 130 nmi downstream, indicate that aircraft arrived to the meter fix six seconds early with about a twelve second standard deviation. Large FPA errors up to one degree relative to the EDA flight trial procedure were detected after the flight trial as a characteristic of the unexpected variation. It is recommended that quantitative validation be performed during future flight trials so that experimental procedures can be adjusted if unexpected results are detected.

Field Test Results↗
DOE OSTI2018

Crystal structure refinements of tetragonal (OH,F)-rich spessartine and henritermierite garnets

Cubic garnet (space group Ia\overline 3 d) has the general formula X 3 Y 2 Z 3 O 12 , where X , Y and Z are cation sites. In the tetragonal garnet (space group I 4 1 / acd ), the corresponding cation sites are X 1 and X 2, Y , and Z 1 and Z 2. In both space groups only the Y site is the same. The crystal chemistry of a tetragonal (OH,F)-rich spessartine sample from Tongbei, near Yunxiao, Fujian Province, China, with composition X (Mn 2.82 Fe^{2+}_{0.14}Ca 0.04 ) Σ3 Y {Al 1.95 Fe^{3+}_{0.05}} Σ2 Z [(SiO 4 ) 2.61 (O 4 H 4 ) 0.28 (F 4 ) 0.11 ] Σ3 (Sps 94 Alm 5 Grs 1 ) was studied with single-crystal X-ray diffraction and space group I 4 1 / acd . The deviation of the unit-cell parameters from cubic symmetry is small [ a = 11.64463(1), c = 11.65481 (2) Å, c / a = 1.0009]. Point analyses and back-scattered electron images, obtained by electron-probe microanalysis, indicate a homogeneous composition. The Z 2 site is fully occupied, but the Z 1 site contains vacancies. The occupied Z 1 and Z 2 sites with Si atoms are surrounded by four O atoms, as in anhydrous cubic garnets. Pairs of split sites are O1 with F11 and O2 with O22. When the Z 1 site is vacant, a larger [(O 2 H 2 )F 2 ] tetrahedron is formed by two OH and two F anions in the O22 and F11 sites, respectively. This [(O 2 H 2 )F 2 ] tetrahedron is similar to the O 4 H 4 tetrahedron in hydrogarnets. These results indicate ^{X}{{\rm Mn}^ {2+}_{3}}\,^{Y}{\rm Al}_{2}^{Z}[({\rm SiO}_{4})_{2}({\rm O}_{2}{\rm H}_{2})_{0.5}({\rm F}_{2})_{0.5}]_{\Sigma3} as a possible end member, which is yet unknown. The H atom that is bonded to the O22 site is not located because of the small number of OH groups. In contrast, tetragonal henritermierite, ideally ^{X}{\rm Ca}_{3}\,^{Y}{\rm Mn}^{3+}_{2}\,^{Z}[({\rm SiO}_{4})_{2}({\rm O}_{4}{\rm H}_{4})_1]_{\Sigma3}, has a vacant Z 2 site that contains the O 4 H 4 tetrahedron. The H atom is bonded to an O3 atom [O3—H3 = 0.73 (2) Å]. Because of O2—Mn 3+ —O2 Jahn–Teller elongation of the Mn 3+ O 6 octahedron, a weak hydrogen bond is formed to the under-bonded O2 atom. This causes a large deviation from cubic symmetry ( c / a = 0.9534).

Antao, Sytle M.↗
DOE OSTI2020

Fisher for complements: extracting cosmology and neutrino mass from the counts-in-cells PDF

In this work, we comprehensively analyse the cosmology dependence of counts-in-cells statistics. Furthermore, we focus on the shape of the one-point probability distribution function (PDF) of the matter density field at mildly non-linear scales. Based on large-deviation statistics, we parametrize the cosmology dependence of the matter PDF in terms of the linear power spectrum, the growth factor, the spherical collapse dynamics, and the non-linear variance. We extend our formalism to include massive neutrinos, finding that the total matter PDF is highly sensitive to the total neutrino mass Mν and can disentangle it from the clustering amplitude σ 8 . Additionally, using more than a million PDFs extracted from the Quijote simulations, we determine the response of the matter PDF to changing parameters in the νΛCDM model and successfully cross-validate the theoretical model and the simulation measurements. We present the first νΛCDM Fisher forecast for the matter PDF at multiple scales and redshifts, and its combination with the matter power spectrum. We establish that the matter PDF and the matter power spectrum are highly complementary at mildly non-linear scales. The matter PDF is particularly powerful for constraining the matter density Ωm, clustering amplitude σ 8 and the total neutrino mass Mν. Adding the mildly non-linear matter PDF to the mildly non-linear matter power spectrum improves constraints on Ω m by a factor of 5 and σ 8 by a factor of 2 when considering the three lowest redshifts. In our joint analysis of the matter PDF and matter power spectrum at three redshifts, the total neutrino mass is constrained to better than 0.01 eV with a total volume of 6 (Gpc h -1 ) 3 . We discuss how density-split statistics can be used to translate those encouraging results for the matter PDF into realistic observables in galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Complete resolution across the neodymium/samarium isotopic envelope with a liquid sampling‐atmospheric pressure glow discharge — Orbitrap mass spectrometer

Rationale Nd and Sm isotope ratios play an important role in geological dating and as nuclear forensic signatures; however, the overlap of the respective 144, 148, 150 Nd/Sm isobars requires prior separations to be performed before analysis on typical MS platforms. The work presented here overcomes these isobaric interferences using ultrahigh‐mass resolution to alleviate interference without prior chemical separations. Methods A liquid sampling‐atmospheric pressure glow discharge ion source was coupled to a standard, QExactive Focus Orbitrap mass spectrometer, providing a mass resolution of ~80 k. A Spectroswiss FTMS booster X2 data acquisition package was used to collect extended transients, providing much higher mass resolution; ~230 k and ~600 k are employed here for Nd and Sm isotopes. Results While the standard Orbitrap resolution is far greater than typical “atomic” MS platforms, it was insufficient to alleviate all isobars. The use of a resolution of ~230 k resulted in baseline separation across the entire isotopic envelope for both Nd and Sm. Isotope ratios obtained from Nd:Sm mixtures using high‐resolution were equivalent to those found for individual‐element solutions, while isotope ratios obtained at a resolution of ~80 k (standard for the OEM data system) showed large deviations. Conclusions Use of ultrahigh‐resolution is an attractive alternative to extensive chemical separations to alleviate severe isobaric interferences. Sufficient mass resolution greatly reduces/eliminates the need for sample manipulations (separations) before analysis while reducing costs and total analysis times.

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

A Characterization of the Reflected Quasipotential

Recent interest in the reflected quasipotential comes from the queueing theory literature, specifically the analysis of so-called (b,A,D) reflected Brownian motion where it is the large deviation rate function for the stationary distribution. Our purpose here is to characterize the reflected quasipotential in terms of a first-order Hamilton–Jacobi equation. Using conventional dynamic programming ideas, along with a complementarity problem formulation of the effect of the Skorokhod map on absolutely continuous paths, we will derive necessary conditions in the form of viscosity-sense boundary conditions. It turns out that even with these boundary conditions solutions are not unique. Thus a unique characterization needs to refer to some additional property of V(⋅). We establish such a characterization in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Asymptotic Properties of the Density of Particles in β-Ensembles

We extend recent results on the asymptotic equipartition property for the density of n particles in β-ensembles, as n tends to infinity. We prove the large deviation principle of the log-density for a general potential and the mod-Gaussian convergence in the classical examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Cluster Tails for Critical Power-Law Inhomogeneous Random Graphs

Recently, the scaling limit of cluster sizes for critical inhomogeneous random graphs of rank-1 type having finite variance but infinite third moment degrees was obtained in Bhamidi et al. (Ann Probab 40:2299–2361, 2012). It was proved that when the degrees obey a power law with exponent τ∈(3,4), the sequence of clusters ordered in decreasing size and multiplied through by n{sup −(τ−2)/(τ−1)} converges as n→∞ to a sequence of decreasing non-degenerate random variables. Here, we study the tails of the limit of the rescaled largest cluster, i.e., the probability that the scaling limit of the largest cluster takes a large value u, as a function of u. This extends a related result of Pittel (J Combin Theory Ser B 82(2):237–269, 2001) for the Erdős–Rényi random graph to the setting of rank-1 inhomogeneous random graphs with infinite third moment degrees. We make use of delicate large deviations and weak convergence arguments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Canonical Structure and Orthogonality of Forces and Currents in Irreversible Markov Chains

We discuss a canonical structure that provides a unifying description of dynamical large deviations for irreversible finite state Markov chains (continuous time), Onsager theory, and Macroscopic Fluctuation Theory (MFT). For Markov chains, this theory involves a non-linear relation between probability currents and their conjugate forces. Within this framework, we show how the forces can be split into two components, which are orthogonal to each other, in a generalised sense. This splitting allows a decomposition of the pathwise rate function into three terms, which have physical interpretations in terms of dissipation and convergence to equilibrium. Similar decompositions hold for rate functions at level 2 and level 2.5. These results clarify how bounds on entropy production and fluctuation theorems emerge from the underlying dynamical rules. We discuss how these results for Markov chains are related to similar structures within MFT, which describes hydrodynamic limits of such microscopic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Front Propagation for Reaction–Diffusion Equations in Composite Structures

We consider asymptotic problems concerning the motion of interface separating the regions of large and small values of the solution of a reaction–diffusion equation in the media consisting of domains with different characteristics (composites). Under certain conditions, the motion can be described by the Huygens principle in the appropriate Finsler (e.g., Riemannian) metric. In general, the motion of the interface has, in a sense, non-local nature. In particular, the interface may move by jumps. We are mostly concerned with the nonlinear term that is of KPP type. The results are based on limit theorems for large deviations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Constitutive Modeling of the High-Temperature Flow Behavior of α-Ti Alloy Tube

In the hot metal gas forming process, the deformation conditions, such as temperature, strain rate and deformation degree, are often prominently changed. The understanding of the flow behavior of α-Ti seamless tubes over a relatively wide range of temperatures and strain rates is important. In this study, the stress–strain curves in the temperature range of 973-1123 K and the initial strain rate range of 0.0004-0.4 s{sup −1} were measured by isothermal tensile tests to conduct a constitutive analysis and a deformation behavior analysis. The results show that the flow stress decreases with the decrease in the strain rate and the increase of the deformation temperature. The Fields–Backofen model and Fields–Backofen–Zhang model were used to describe the stress–strain curves. The Fields–Backofen–Zhang model shows better predictability on the flow stress than the Fields–Backofen model, but there exists a large deviation in the deformation condition of 0.4 s{sup −1}. A modified Fields–Backofen–Zhang model is proposed, in which a strain rate term is introduced. This modified Fields–Backofen–Zhang model gives a more accurate description of the flow stress variation under hot forming conditions with a higher strain rate up to 0.4 s{sup −1}. Accordingly, it is reasonable to adopt the modified Fields–Backofen–Zhang model for the hot forming process which is likely to reach a higher strain rate, such as 0.4 s{sup −1}.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Triple junction drag effects during topological changes in the evolution of polycrystalline microstructures

Experiments, theory and atomistic simulations show that finite triple junction mobility results in non-equilibrium triple junction angles in evolving polycrystalline systems. These angles have been predicted and verified for cases where grain boundary migration is steady-state. Yet, steady-state never occurs during the evolution of polycrystalline microstructures as a result of changing grain size and topological events (e.g., grain face/edge switching - “T 1 ” process, or grain disappearance “T 2 ” or “T 3 ” processes). We examine the non-steady evolution of the triple junction angle in the vicinity of topological events and show that large deviations from equilibrium and/or steady-state angles occur. We analyze $\tau$ the characteristic relaxation time of triple junction angles τ by consideration of a pair of topological events, beginning from steady-state migration. Using numerical results and theoretical analysis we predict how the triple junction angle varies with time and how τ varies with triple junction mobility. We argue that it is precisely those cases where grain boundaries are moving quickly (e.g., topological process in nanocrystalline materials), that the classical steady-state prediction of the triple junction angle about finite triple junction mobility is inapplicable and may only be applied qualitatively.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Exponential number of equilibria and depinning threshold for a directed polymer in a random potential

By extending the Kac–Rice approach to manifolds of finite internal dimension, we show that the mean number ⟨N{sub tot}⟩ of all possible equilibria (i.e. force-free configurations, a.k.a. equilibrium points) of an elastic line (directed polymer), confined in a harmonic well and submitted to a quenched random Gaussian potential in dimension d=1+1, grows exponentially ⟨N{sub tot}⟩ ~exp(r L) with its length L. The growth rate r is found to be directly related to the generalized Lyapunov exponent (GLE) which is a moment-generating function characterizing the large-deviation type fluctuations of the solution to the initial value problem associated with the random Schrödinger operator of the 1D Anderson localization problem. For strong confinement, the rate r is small and given by a non-perturbative (instanton, Lifshitz tail-like) contribution to GLE. For weak confinement, the rate r is found to be proportional to the inverse Larkin length of the pinning theory. As an application, identifying the depinning with a landscape “topology trivialization” phenomenon, we obtain an upper bound for the depinning threshold f{sub c}, in the presence of an applied force, for elastic lines and d-dimensional manifolds, expressed through the mean modulus of the spectral determinant of the Laplace operators with a random potential. We also discuss the question of counting of stable equilibria. Finally, we extend the method to calculate the asymptotic number of equilibria at fixed energy (elastic, potential and total), and obtain the (annealed) distribution of the energy density over these equilibria (i.e. force-free configurations). Some connections with the Larkin model are also established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

On the Mechanism of Bilayer Separation by Extrusion, or Why Your LUVs Are Not Really Unilamellar

Extrusion through porous filters is a widely used method for preparing biomimetic model membranes. Of primary importance in this approach is the efficient production of single bilayer (unilamellar) vesicles that eliminate the influence of interlamellar interactions and strictly define the bilayer surface area available to external reagents such as proteins. Sub-microscopic vesicles produced using extrusion are widely assumed to be unilamellar, and large deviations from this assumption would dramatically impact interpretations from many model membrane experiments. In this work, using three probe-free methods—small-angle X-ray and neutron scattering (SAXS and SANS) and cryogenic electron microscopy (cryoEM)—we report unambiguous evidence of extensive multilamellarity in extruded vesicles composed of neutral phosphatidylcholine lipids, including for the common case of neutral lipids dispersed in physiological buffer and extruded through 100 nm diameter pores. In such preparations, only ~35% of lipids are externally accessible, and this fraction is highly dependent on preparation conditions. Charged lipids promote unilamellarity, as does decreasing solvent ionic strength, indicating the importance of electrostatic interactions in determining the lamellarity of extruded vesicles. Smaller extrusion pore sizes also robustly increase the fraction of unilamellar vesicles, suggesting a role for membrane bending. Taken together, these observations suggest a mechanistic model for extrusion, wherein formation of unilamellar vesicles involves competition between bilayer bending and adhesion energies. The findings presented here have wide-ranging implications for the design and interpretation of model membrane studies, especially ensemble-averaged observations relying on the assumption of unilamellarity.

59 BASIC BIOLOGICAL SCIENCES↗