SEARCH · Search NASA

Results for “symmetric cells”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Asymmetric ether solvents for high-rate lithium metal batteries

Recent electrolyte solvent design based on weakening lithium-ion solvation have shown promise in enhancing cycling performance of Li-metal batteries. However, they often face slow redox kinetics and poor cycling reversibility at high rate. Here we report using asymmetric solvent molecules substantially accelerates Li redox kinetics. Asymmetric ethers (1-ethoxy-2-methoxyethane, 1-methoxy-2-propoxyethane) showed higher exchange current densities and enhanced high-rate Li 0 plating/stripping reversibility compared to symmetric ethers. Adjusting fluorination levels further improved oxidative stability and Li 0 reversibility. The asymmetric 1-(2,2,2-trifluoro)-ethoxy-2-methoxyethane, with 2 M lithium bis(fluorosulfonyl)imide, exhibited high exchange current density, oxidative stability, compact solid–electrolyte interphase (~10 nm). This electrolyte exhibited superior performance among state-of-the-art electrolytes, enabling over 220 cycles in high-rate Li (50 μm)||LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811, 4.9 mAh cm −2 ) cells and for the first time over 600 cycles in anode-free Cu | |Ni95 pouch cells (200 mAh) under electric vertical take-off and landing cycling protocols. Our findings on asymmetric molecular design strategy points to a new pathway towards achieving fast redox kinetics for high-power Li-metal batteries.

batteries↗

DOE OSTI2022

Spectroscopy of Twisted Bilayer Graphene Correlated Insulators

Here we analytically compute the scanning tunneling microscopy (STM) signatures of integer-filled correlated ground states of the magic angle twisted bilayer graphene (TBG) narrow bands. After experimentally validating the strong-coupling approach at ± 4 electrons/moiré unit cell, we consider the spatial features of the STM signal for 14 different many-body correlated states and assess the possibility of Kekulé distortion (KD) emerging at the graphene lattice scale. Remarkably, we find that coupling the two opposite graphene valleys in the intervalley-coherent (IVC) TBG insulators does not always result in KD. As an example, we show that the Kramers IVC state and its nonchiral U(4) rotations do not exhibit any KD, while the time-reversal-symmetric IVC state does. Our results, obtained over a large range of energies and model parameters, show that the STM signal and Chern number of a state can be used to uniquely determine the nature of the TBG ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Collisionless zonal-flow dynamics in quasisymmetric stellarators

The linear collisionless plasma response to a zonal-density perturbation in quasisymmetric stellarators is studied, including the geodesic-acoustic-mode oscillations and the Rosenbluth–Hinton residual flow. While the geodesic-acoustic-mode oscillations in quasiaxisymmetric configurations are similar to tokamaks, they become non-existent in quasi-helically symmetric configurations when the effective safety factor in helical-angle coordinates is small. Compared with concentric-circular tokamaks, the Rosenbluth–Hinton residual is also found to be multiplied by a geometric factor $\mathcal {C}$ that arises from the flux-surface-averaged classical polarization. Using the near-axis-expansion framework, we derive an analytic expression for $\mathcal {C}$ , which varies significantly among different configurations. These analytic results are compared with numerical simulation results from the global gyrokinetic particle-in-cell code GTC, and good agreement with the theoretical Rosenbluth–Hinton residual level is achieved when the quasisymmetry error is small enough.

Zhu, Hongxuan (ORCID:0000000198446972)↗

DOE OSTI2020

Beyond integration: modeling every pixel to obtain better structure factors from stills

Most crystallographic data processing methods use pixel integration. In serial femtosecond crystallography (SFX), the intricate interaction between the reciprocal lattice point and the Ewald sphere is integrated out by averaging symmetrically equivalent observations recorded across a large number (10 4 -10 6 ) of exposures. Although sufficient for generating biological insights, this approach converges slowly, and using it to accurately measure anomalous differences has proved difficult. This report presents a novel approach for increasing the accuracy of structure factors obtained from SFX data. A physical model describing all observed pixels is defined to a degree of complexity such that it can decouple the various contributions to the pixel intensities. Model dependencies include lattice orientation, unit-cell dimensions, mosaic structure, incident photon spectra and structure factor amplitudes. Maximum likelihood estimation is used to optimize all model parameters. The application of prior knowledge that structure factor amplitudes are positive quantities is included in the form of a reparameterization. The method is tested using a synthesized SFX dataset of ytterbium(III) lysozyme, where each X-ray laser pulse energy is centered at 9034 eV. This energy is 100 eV above the Yb 3+ L -III absorption edge, so the anomalous difference signal is stable at 10 electrons despite the inherent energy jitter of each femtosecond X-ray laser pulse. This work demonstrates that this approach allows the determination of anomalous structure factors with very high accuracy while requiring an order-of-magnitude fewer shots than conventional integration-based methods would require to achieve similar results.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Regulation of Side‐Chain Symmetry for Delaying Triplet Formation and Suppressing Non‐Radiative Loss in Organic Solar Cells

Abstract The structural revolutions of non‐fullerene acceptors (NFAs) have driven continuous efficiency breakthroughs in organic solar cells (OSCs). Rational regulation of NFA structures toward efficient exciton dissociation and mitigated non‐radiative recombination is pivotal for OSCs. The incorporation of asymmetric side chains on NFAs can often achieve these goals by inducing a desirable aggregate state. However, it lacks the studies to directly correlate the side‐chain symmetry of NFAs with the exciton delocalization and triplet dynamics in OSCs. Herein, The influence of structural symmetry on the aggregate properties is systematically investigated and exciton/charge dynamics based on two developed biaxial‐conjugated NFAs with varied side‐chain symmetry. The symmetric NFA having diverse molecular packing orientations can form multiple charge transfer channels in its blend with polymer donor, which cannot be found in that comprising the asymmetric ones. Moreover, a slower rate and lower ratio of the spin‐triplet state are formed in the blend of symmetric NFA, resulting in a much lower non‐radiative voltage loss in corresponding OSCs. This study reveals the distinct advantages of symmetric NFAs in both aggregate properties and exciton/charge dynamics over those of asymmetric ones, paving the way for developing high‐performance OSCs using easier‐to‐prepare, low‐cost symmetric materials.

Chemistry↗

DOE OSTI2022

Lithium Solvation and Mobility in Ionic Liquid Electrolytes with Asymmetric Sulfonyl-Cyano Anion

The solvation structure and transport properties of Li + in ionic liquid (IL) electrolytes based on n-methyl-n-butylpyrrolidinium cyano(trifluoromethanesulfonyl)imide [PYR 14 ][CTFSI] and [Li][CTFSI] (0 ≤ x Li ≤ 0.7) were studied by Raman and Nuclear Magnetic Resonance (NMR) diffusometry, and molecular dynamics (MD) simulations. At x Li < 0.3, Li + coordination is dominated by the cyano group. As x Li is increased, free cyano-sites become limited, resulting in increased coordination via the sulfonyl group. Here, the 1:1 mixture of the symmetric anions bis(trifluoromethanesulfonyl)imide ([TFSI]) and dicyanamide ([DCA]) results in similar physical properties as the IL with [CTFSI]. However, anion asymmetry is shown to increase Li-salt solubility and promote Li+ transference. The lifetimes of Li + -cyano coordination for [CTFSI] are calculated to be shorter than those for [DCA], indicating that the competition from the sulfonyl group weakens its solvation with Li + . This resulted in higher Li + transference for the electrolyte with [CTFSI]. In relation to the utility of these electrolytes in energy storage, the Li–LiFePO 4 half cells assembled with IL electrolyte (x Li = 0.3, 0.5, and 0.7) demonstrated a nominal capacity of 140 mAh/g at 0.1C rate and 90 °C where the cell with x Li = 0.7 IL electrolyte demonstrated 61% capacity retention after 100 cycles and superior rate capability owing to increased electrochemical stability.

25 ENERGY STORAGE↗

DOE OSTI2022

Restored Passivation after Complete Removal of Front Poly-Si Between the Grid in Poly-Si/SiO2 Front/Back Cells

This work explores the use of wet chemical etching to completely remove the front poly-Si layer between grid lines of a front/back poly-Si/SiO2 passivated contact device. The purpose of the work is to provide excellent passivation under the metal grid, but to also maximize Jsc by eliminating optical absorption in the front poly-Si. We show that after a high temperature anneal to distribute carriers in the poly-Si/SiO2/c-Si interface all of the poly-Si can be removed selectively with TMAH, using the SiO2 as an etch-stop layer. The etched surface can be repassivated back to its original level with the addition of a SiNx/Al2O3 stack grown on the preserved tunnelling oxide layer followed by a short forming gas anneal. Using symmetric 100 - 200 nm intrinsically doped poly-Si/SiO2 on textured n-type wafer templates we show absolute device efficiency enhancement from 15% to 21% by increasing the Jsc by 12 mA/cm2 after removing the front poly-Si and repassivating while effectively preserving the Voc, and FF of the device. IQE values are similar to PERC devices in the short wavelength range.

etching↗

DOE OSTI2024

The derivation of CRSS in pure Ti and Ti-Al alloys

The work focuses on the determination of the critical resolved shear stress (CRSS) in titanium (Ti) and titanium-aluminum (Ti-Al) alloys, influenced by an array of factors such as non-symmetric fault energies and minimum energy paths, dislocation core-widths, short-range order (SRO) effects which alter the local atomic environment, and tension-compression (T-C) asymmetry affected by intermittent slip motion. To address these multifaceted complexities, an advanced theory has been developed, offering an in-depth understanding of the mechanisms underlying slip behavior. The active slip systems in these materials are basal, prismatic, and pyramidal planes, with the latter involving both $\langle$a$\rangle$ and $\langle$c + a$\rangle$ dislocations. Each slip system is characterized by distinct Wigner-Seitz cell configurations for misfit energy calculations, varying partial dislocation separation distances, and unique dislocation trajectories—all critical to precise CRSS calculations. The theoretical CRSS results were validated against a comprehensive range of experimental data, demonstrating a strong agreement and underscoring the model's efficacy.

Critical stress↗

DOE OSTI2023

Phase diagram of twisted bilayer graphene at filling factor $v = ±3$

Here, we study the correlated insulating phases of twisted bilayer graphene (TBG) in the absence of lattice strain at integer filling $v = ±3$. Using the self-consistent Hartree-Fock method on a particle-hole symmetric model and allowing translation symmetry breaking terms, we obtain the phase diagram with respect to the ratio of AA interlayer hopping (w 0 ) and AB interlayer hopping (w 1 ). When the interlayer hopping ratio is close to the chiral limit (w 0 /w 1 ≲ 0.5), a quantum anomalous Hall state with Chern number v c = ±1 can be observed consistent with previous studies. Around the realistic value w 0 /w 1 ≈ 0.8, we find a spin and valley polarized, translation symmetry breaking, state with C 2⁢z T symmetry, a charge gap and a doubling of the moiré unit cell, dubbed theC 2⁢z T stripe phase. The real-space total charge distribution of this C 2⁢z T stripe phase in the flat band limit does not have modulation between different moiré unit cells, although the charge density in each layer is modulated, and the translation symmetry is strongly broken. Other symmetries, including C 2⁢z , C 2⁢x , particle-hole symmetry P, and the topology of the C 2⁢z T stripe phase, are also discussed in detail. We observed braiding and annihilation of the Dirac nodes by continuously turning on the order parameter to its fully self-consistent value, and provide a detailed explanation of the mechanism for the charge gap opening despite preserving C 2⁢z T and valley U⁡(1) symmetries. In the transition region between the quantum anomalous Hall phase and the C 2⁢z T stripe phase, we find an additional competing state with comparable energy corresponding to a phase with a tripling of the moiré unit cell.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

De novo design of modular protein hydrogels with programmable intra- and extracellular viscoelasticity

Relating the macroscopic properties of protein-based materials to their underlying component microstructure is an outstanding challenge. Here, we exploit computational design to specify the size, flexibility, and valency of de novo protein building blocks, as well as the interaction dynamics between them, to investigate how molecular parameters govern the macroscopic viscoelasticity of the resultant protein hydrogels. We construct gel systems from pairs of symmetric protein homo-oligomers, each comprising 2, 5, 24, or 120 individual protein components, that are crosslinked either physically or covalently into idealized step-growth biopolymer networks. Through rheological assessment, we find that the covalent linkage of multifunctional precursors yields hydrogels whose viscoelasticity depends on the crosslink length between the constituent building blocks. In contrast, reversibly crosslinking the homo-oligomeric components with a computationally designed heterodimer results in viscoelastic biomaterials exhibiting fluid-like properties under rest and low shear, but solid-like behavior at higher frequencies. Exploiting the unique genetic encodability of these materials, we demonstrate the assembly of protein networks within living mammalian cells and show via fluorescence recovery after photobleaching (FRAP) that mechanical properties can be tuned intracellularly in a manner similar to formulations formed extracellularly. We anticipate that the ability to modularly construct and systematically program the viscoelastic properties of designer protein-based materials could have broad utility in biomedicine, with applications in tissue engineering, therapeutic delivery, and synthetic biology.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Dimer description of the SU(4) antiferromagnet on the triangular lattice

In systems with many local degrees of freedom, high-symmetry points in the phase diagram can provide an important starting point for the investigation of their properties throughout the phase diagram. In systems with both spin and orbital (or valley) degrees of freedom such a starting point gives rise to SU(4)-symmetric models. Here we consider SU(4)-symmetric "spin" models, corresponding to Mott phases at half-filling, i.e. the six-dimensional representation of SU(4). This may be relevant to twisted multilayer graphene. In particular, we study the SU(4) antiferromagnetic "Heisenberg" model on the triangular lattice, both in the classical limit and in the quantum regime. Carrying out a numerical study using the density matrix renormalization group (DMRG), we argue that the ground state is non-magnetic. We then derive a dimer expansion of the SU(4) spin model. An exact diagonalization (ED) study of the effective dimer model suggests that the ground state breaks translation invariance, forming a valence bond solid (VBS) with a 12-site unit cell. Finally, we consider the effect of SU(4)-symmetry breaking interactions due to Hund's coupling, and argue for a possible phase transition between a VBS and a magnetically ordered state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Effect of the J-Shaped Wind Turbine Airfoil Opening Ratio and Thickness on the Performance of Symmetrical Airfoils

Abstract In this work, the effect of the inner opening ratio on the J-shaped airfoils aerodynamic performance was studied and documented for symmetrical airfoils. Three different airfoil thicknesses were investigated: small (NACA0008), medium (NACA0015), and large (NACA0024). For each airfoil thickness, effects of three inner opening ratios were analyzed: one-third, one-half, and two-thirds. The performance of each opening ratio was compared with the performance of the solid airfoil “zero opening ratio” for different angles of attack between 5 deg and 20 deg. All designs were simulated using the computational fluid dynamics (CFD) technology against experimental results for solid NACA4412 airfoil in the University of Wisconsin-Milwaukee (UWM) wind tunnel facility and other published experimental data. It was found that large eddy simulation yields accurate solutions with a smaller number of mesh cells compared to the k–ω turbulence model but with much longer computational time. The lift-to-drag ratio for all studied airfoils has a maximum value for solid airfoils compared to those equipped with openings. For airfoils equipped with 00.00% opening ratio “solid,” NACA0015 airfoil has the maximum lift-to-drag ratio. Furthermore, it was found that NACA0008 equipped with a 33.33% opening ratio has the best performance of all studied J-shaped airfoils.

Energy & Fuels↗

DOE OSTI2020

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Observation of Skewed Electromagnetic Wakefields in an Asymmetric Structure Driven by Flat Electron Bunches

Relativistic charged -particle beams that generate intense longitudinal fields in accelerating structures also inherently couple to transverse modes. The effects of this coupling may lead to beam breakup instability and thus must be countered to preserve beam quality in applications such as linear colliders. Beams with highly asymmetric transverse sizes (flat beams) have been shown to suppress the initial instability in slab -symmetric structures. However, as the coupling to transverse modes remains, this solution serves only to delay instability. In order to understand the hazards of transverse coupling in such a case, we describe here an experiment characterizing the transverse effects on a flat beam, traversing near a planar dielectric lined structure. Further, the measurements reveal the emergence of a previously unobserved skew-quadrupolelike interaction when the beam is canted transversely, which is not present when the flat beam travels parallel to the dielectric surface. We deploy a multipole field fitting algorithm to reconstruct the projected transverse wakefields from the data. We generate the effective kick vector map using a simple two -particle theoretical model, with particle -in -cell simulations used to provide further insight for realistic particle distributions.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Effects of filling, strain, and electric field on the Néel vector in antiferromagnetic CrSb

CrSb is a layered antiferromagnet (AFM) with perpendicular magnetic anisotropy, a high Néel temperature, and large spin-orbit coupling (SOC), which makes it interesting for AFM spintronic applications. To elucidate the various mechanisms of Néel vector control, the effects of strain, band filling, and electric field on the magnetic anisotropy energy (MAE) of bulk and thin-film CrSb are determined and analyzed using density functional theory. The MAE of the bulk crystal is large (1.2 meV per unit cell). Due to the significant ionic nature of the Cr-Sb bond, finite slabs are strongly affected by end termination. Truncation of the bulk crystal to a thin film with one surface terminated with Cr and the other surface terminated with Sb breaks inversion symmetry, creates a large charge dipole and average electric field across the film, and breaks spin degeneracy, such that the thin film becomes a ferrimagnet. The MAE is reduced such that its sign can be switched with realistic strain, and the large SOC gives rise to an intrinsic voltage controlled magnetic anisotropy. A slab terminated on both faces with Cr remains a compensated AFM, but with the compensation occurring nonlocally between mirror-symmetric Cr pairs. In-plane alignment of the moments is preferred, the magnitude of the MAE remains large, similar to that of the bulk, and it is relatively insensitive to filling.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Synthetic Pathways to gamma-Graphyne and Related Allotropes of Carbon

Graphynes, two-dimensional carbon lattices combining sp 1 and sp 2 hybridized atoms, were predicted theoretically more than three decades ago, but few structures have been realized to date. These carbons are believed to possess remarkable mechanical and electronic properties, including high charge carrier mobilities comparable to those in graphene (10 4 to 10 5 cm 2 V -1 s -1 ). Unlike graphene, certain graphynes are predicted to be intrinsic semiconductors. Among these intriguing structures, γ-graphyne stands out as the structurally simplest and most symmetric sp 1 /sp 2 lattice. γ-Graphyne was first theorized in 1987. In contrast with graphene, γ-graphyne will be a semiconductor with a small band gap suitable for fabrication of electronic devices. This solves one of the fundamental problems of carbon-based electronics, the necessity for inducing a band gap in graphene. γ-Graphyne has the potential to form the basis for the next generation of carbon-based electronics operating at speeds unattainable by traditional silicon chips. Unlike silicon, γ-graphyne is a direct band gap semiconductor, and it will feature exceptional strength comparable to that of other 2D carbon allotropes. Such combination of properties may enable a new generation of highly efficient, ultra-light and flexible solar cells. Despite being a potentially “magical” material, γ-graphyne remained synthetically elusive for over three decades. The primary goals of this project were: (1) Synthesis of bulk γ-graphyne phases through solution-phase 2D polymerizations; (2) Experimental exploration of the physical and chemical properties of γ-graphyne; and (3) Mechanistic and theoretical studies of the novel chemical transformations developed in Goal 1. Common pyrolytic and vapor-deposition methodologies used for the synthesis of graphitic allotropes are unsuitable for graphyne and other sp 1 -contaning structures, as acetylenes readily convert to graphene and amorphous carbon at high temperatures. In contrast, this proposal is based on solution-based 2D polymerization. The major advantages of this approach over the traditional high temperature techniques are the potential to adjust the structure of the material with atomic precision, and the possibility of using structurally complex and relatively fragile repeat units. The outcomes of this research can revolutionize carbon nanotechnology, expanding the field’s structural toolbox beyond primarily graphitic and benzenoid structures. Understanding the chemistry of sp 1 carbon allotropes can lead to entirely new classes of structures with unique properties, including graphyne ribbons, nanotubes, quantum dots, and heterostructures with other 2D materials. Furthermore, the development of reliable and robust synthetic pathways towards periodic covalent molecular sheets with atomically precise structures shall have a profound impact on chemistry and materials science.

2D polymerization↗

DOE OSTI2020

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. Additionally, the theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the “accessibility” of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗