Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conjunctions”

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 901 records · Page 50

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗

Effects of Supercritical CO 2 Injection on the Shale Pore Structures and Mass Transport Rates

Characterizing the pore structures and transport properties of low-permeability shales is critical for evaluating these formations as potential seals or storage sites for geological CO 2 sequestration. Here, we use low-pressure gas adsorption in conjunction with nuclear magnetic resonance (NMR) to characterize the pore-size distribution of shales before and after injection of supercritical CO 2 . Nitrogen gas was used as the detecting phase for the adsorption experiments and pentane liquid was used for the NMR experiments. We also performed time-resolved NMR and gravimetric microbalance measurements to observe mass transport during desaturation. We use these data to estimate the self-diffusion coefficient of pentane and changes in the saturation state of the pore network. We analyzed samples with a range of compositions from the Wolfcamp shale before and after exposure to supercritical CO 2 for 3 days. Integrating the gas adsorption and NMR data shows how supercritical CO 2 injection alters the pore-size distribution for pore sizes <1 nm to 1 mm. Finally, our results provide insights on how the pore structure and mass transport properties of different shale lithologies may evolve during storage of supercritical CO 2 .

58 GEOSCIENCES↗

Insights into Constraining Rate Coefficients in Fuel Oxidation Mechanisms Using Genetic Algorithm Optimization

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O 2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. In conclusion, comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

33 ADVANCED PROPULSION SYSTEMS↗

Constraining the Low-Temperature Oxidation Mechanism of n -Hexanol through the Detection and Identification of C 6 Elusive Intermediates

Alcohol-based fuels are currently considered to be viable energy carriers for the transportation sector. Consequently, a comprehensive mechanistic understanding of the low-temperature oxidation of alcohols is essential for application in advanced low-temperature compression engines. Here, in this work, a multidimensional approach involving experimental investigations, kinetic modeling, and theoretical calculations was used to provide new insights into the low-temperature oxidation mechanism of a C 6 alcohol, n -hexanol (CH 3 (CH 2 ) 5 OH), through the detection and identification of elusive C 6 intermediates. The oxidation of n-hexanol was investigated in a jet-stirred reactor under stoichiometric conditions (ϕ = 1.0), an initial fuel concentration of 2%, a residence time of 2 s, a temperature range between 500 and 660 K, and a pressure of 700 Torr. The reactants, intermediates, and final products were detected and identified by means of molecular-beam mass spectrometry coupled with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. Chemical kinetic simulations were performed using a previously published kinetic model (Togbé et al., Energy Fuels 2010, 11, 5859−5875) to predict the reactivity of n-hexanol and elucidate the predominant formation pathways of the observed low-temperature species. Experimental photoionization efficiency curves in conjunction with ab initio calculations, enabled the identification of important low-temperature species, such as C 6 unsaturated alcohols, C 6 olefinic hydroperoxides, C 6 cyclic ethers, C 6 diones, and C 6 ketohydroperoxides. The results of this study provide valuable insight into the mechanism of the low-temperature oxidation chemistry of n -hexanol, contributing to the development of kinetic models for the low-temperature oxidation of n-hexanol and other long-chain linear alcohols.

alcohols↗

Burn Intensity Drives the Alteration of Phenolic Lignin to (Poly) Aromatic Hydrocarbons as Revealed by Pyrolysis Gas Chromatography–Mass Spectrometry (Py-GC/MS)

High-intensity wildfires alter the chemical composition of organic matter, which is expected to be distinctly different from low-intensity prescribed fires. In this study, we used pyrolysis gas chromatography/mass spectrometry (Py-GC/MS), in conjunction with solid-state 13 C nuclear magnetic resonance (NMR) and Fourier transform infrared (FT-IR) spectroscopy, to assess chemical alterations from three wildfires and a long-term frequent prescribed fire site. Our results showed that black ash formed under moderate intensity burns contained less aromatic (ArH), polyaromatic hydrocarbon (PAH), and nitrogen-containing compounds (Ntg) but more lignin (LgC) and phenol compounds (PhC), compared to white ash formed under high intensity burns. Both 13 C NMR and FT-IR confirmed a higher relative percentage of carboxyl carbon in white ash, indicating the potential for higher water solubility and more mobile carbon, relative to black ash. Compared to wildfires, ash from low-intensity prescribed fire contained less ArH, PAH, and Ntg and more LgC and PhC. Controlled laboratory burning trials indicated that organic matter alteration was sensitive to the burn temperature, but not related to the fuel type (pine vs fir) nor oxygen absence/presence at high burn temperatures. This study concludes that higher burn temperatures resulted in higher (poly)aromatic carbon/nitrogen and lower lignin/phenol compounds.

54 ENVIRONMENTAL SCIENCES↗

Supply and Demand Drivers of Global Hydrogen Deployment in the Transition toward a Decarbonized Energy System

The role of hydrogen in energy system decarbonization is being actively examined by the research and policy communities. We evaluate the potential ‘hydrogen economy’ in global climate change mitigation scenarios using the Global Change Analysis Model (GCAM). We consider major hydrogen production methods in conjunction with delivery options to understand how hydrogen infrastructure affects its deployment. We also consider a rich set of hydrogen end-use technologies and vary their costs to understand how demand technologies affect deployment. We find that the availability of hydrogen transmission and distribution infrastructure primarily affects the hydrogen production mix, particularly the share produced centrally versus onsite, whereas assumptions about end-use technology primarily affect the scale of hydrogen deployment. In effect, hydrogen can be a distributed energy vector, enabled by onsite renewable electrolysis and, to a lesser extent, by onsite industrial production using natural gas with carbon capture and storage (CCS). While the share of hydrogen in final energy is small relative to the share of other major energy carriers in our scenarios, hydrogen enables decarbonization in difficult-to-electrify end uses such as industrial high-temperature heat. Its contribution to greenhouse gas mitigation increases as the climate objective is tightened.

54 ENVIRONMENTAL SCIENCES↗

Site Suitability and Air Pollution Impacts of Composting Infrastructure for California’s Organic Waste Diversion Law

The emerging field of polaritonic chemistry explores the behavior of molecules under strong coupling with cavity modes. Despite recent developments in ab initio polaritonic methods for simulating polaritonic chemistry under electronic strong coupling, their capabilities are limited, especially in cases where the molecule also features strong electronic correlation. To bridge this gap, we have developed a novel method for cavity QED calculations utilizing the Density Matrix Renormalization Group (DMRG) algorithm in conjunction with the Pauli–Fierz Hamiltonian. Our approach is applied to investigate the effect of the cavity on the S0–S1 transition of n-oligoacenes, with n ranging from 2 to 5, encompassing 22 fully correlated π orbitals in the largest pentacene molecule. Our findings indicate that the influence of the cavity intensifies with larger acenes. Additionally, we demonstrate that, unlike the full determinantal representation, DMRG efficiently optimizes and eliminates excess photonic degrees of freedom, resulting in an asymptotically constant computational cost as the photonic basis increases.

54 ENVIRONMENTAL SCIENCES↗

Oxidative Uranium Release from Anoxic Sediments under Diffusion-Limited Conditions

Here, uranium (U) contamination occurs as a result of mining and ore processing; often in alluvial aquifers that contain organic-rich, reduced sediments that accumulate tetravalent U, U(IV). Uranium(IV) is sparingly soluble, but may be mobilized upon exposure to nitrate (NO 3 – ) and oxygen (O 2 ), which become elevated in groundwater due to seasonal fluctuations in the water table. The extent to which oxidative U mobilization can occur depends upon the transport properties of the sediments, the rate of U(IV) oxidation, and the availability of inorganic reductants and organic electron donors that consume oxidants. We investigated the processes governing U release upon exposure of reduced sediments to artificial groundwater containing O 2 or NO 3 – under diffusion-limited conditions. Little U was mobilized during the 85-day reaction, despite rapid diffusion of groundwater within the sediments and the presence of nonuraninite U(IV) species. The production of ferrous iron and sulfide in conjunction with rapid oxidant consumption suggested that the sediments harbored large concentrations of bioavailable organic carbon that fueled anaerobic microbial respiration and stabilized U(IV). Our results suggest that seasonal influxes of O 2 and NO 3 – may cause only localized mobilization of U without leading to export of U from the reducing sediments when ample organic carbon is present.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Role of the Relative Humidity-Dependent Physical State of Organic Particulate Matter in the Uptake of Semivolatile Organic Molecules

The uptake of gas-phase dicarboxylic acids to organic particulate matter (PM) was investigated to probe the role of the PM physical state in exchange processes between gas-phase semivolatile organic molecules and organic PM. A homologous series of probe molecules, specifically isotopically labeled 13C-dicarboxylic acids, was used in conjunction with aerosol mass spectrometry to obtain a quantitative characterization of the uptake to organic PM for different relative humidities (RHs). The PM was produced by the dark ozonolysis of unlabeled α-pinene. The uptake of 13C-labeled oxalic, malonic, and α-ketoglutaric acids increased stepwise by 5 to 15 times with increases in RH from 15 to 80%. The enhanced uptake with increasing RH was explained primarily by the higher molecular diffusivity in the particle phase, as associated with changes in the physical state of the organic PM from a nonliquid state to a progressively less-viscous liquid state. At high RH, the partitioning of the probe molecules to the particle phase was more associated with physicochemical interactions with the organic PM than that with the co-absorbed liquid water. Additionally, uptake of the probe molecules also increased with a decrease in volatility along the homologous series. This study quantitatively shows the key roles of the particle physical state in governing the interactions of organic PM with semivolatile organic molecules.

42 ENGINEERING↗

General H2O Outgassing Model Applicable to Silica-Filled Silicones

Moisture-induced agining and degradation are continuing concerns in long-term storage and use of silicone materials. The silica fillers frequently used to provide mechanical reinforcement in silicone formulations have a high affinity for H2O and act as a source of moisture outgassing later when placed in dry or vacuum applications.

36 MATERIALS SCIENCE↗

Prediction of Solution Behavior via Calorimetric Measurements Allows for Detailed Elucidation of Polyoxometalate Transformation

The solution behavior of a polyoxometalate cluster, LiNa-U 24 Pp 12 (Li 24 Na 24 [(UO 2 O 2 ) 24 (P 2 O 7 ) 12 ]) that consists of 24 uranyl ions, peroxide groups, and 12 pyrophosphate linkers, was successfully predicted based on new thermodynamic results using a calorimetric method recently described for uranyl peroxide nanoclusters (UPCs), molybdenum blues, and molybdenum browns. The breakdown of LiNa-U 24 Pp 12 and formation of U 24 (Li 24 [UO 2 O 2 OH] 24 ) was monitored in situ via Raman spectroscopy using a custom heating apparatus. A combination of analytical techniques confirmed the simultaneous existence of U 24 Pp 12 and U 24 midway through the conversion process and U 24 as the single end product. The application of a molecular weight filter resulted in a complete and successful separation of UPCs from solution and, in conjunction with DOSY results, confirmed the presence of large intermediate cluster building blocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin and Orbital Effects on Asymmetric Exchange Interaction in Polar Magnets: M(IO 3 ) 2 (M = Cu and Mn)

We report magnetic polar materials feature an astonishing range of physical properties, such as magnetoelectric coupling, chiral spin textures, and related new spin topology physics. This is primarily attributable to their lack of space inversion symmetry in conjunction with unpaired electrons, potentially facilitating an asymmetric Dzyaloshinskii–Moriya (DM) exchange interaction supported by spin–orbital and electron–lattice coupling. However, engineering the appropriate ensemble of coupled degrees of freedom necessary for enhanced DM exchange has remained elusive for polar magnets. Here, we study how spin and orbital components influence the capability of promoting the magnetic interaction by studying two magnetic polar materials, α-Cu(IO 3 ) 2 ( 2 D) and Mn(IO 3 ) 2 ( 6 S), and connecting their electronic and magnetic properties with their structures. The chemically controlled low-temperature synthesis of these complexes resulted in pure polycrystalline samples, providing a viable pathway to prepare bulk forms of transition-metal iodates. Rietveld refinements of the powder synchrotron X-ray diffraction data reveal that these materials exhibit different crystal structures but crystallize in the same polar and chiral P2 1 space group, giving rise to an electric polarization along the b-axis direction. The presence and absence of an evident phase transition to a possible topologically distinct state observed in α-Cu(IO 3 ) 2 and Mn(IO 3 ) 2 , respectively, imply the important role of spin–orbit coupling. Neutron diffraction experiments reveal helpful insights into the magnetic ground state of these materials. While the long-wavelength incommensurability of α-Cu(IO 3 ) 2 is in harmony with sizable asymmetric DM interaction and low dimensionality of the electronic structure, the commensurate stripe AFM ground state of Mn(IO 3 ) 2 is attributed to negligible DM exchange and isotropic orbital overlapping. The work demonstrates connections between combined spin and orbital effects, magnetic coupling dimensionality, and DM exchange, providing a worthwhile approach for tuning asymmetric interaction, which promotes evolution of topologically distinct spin phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Complexation of NpO 2 + with Amine-Functionalized Diacetamide Ligands in Aqueous Solution: Thermodynamic, Structural, and Computational Studies

Complexation of Np(V) with three structurally related amine-functionalized diacetamide ligands, including 2,2'-azanediylbis( N, N'-dimethylacetamide) (ABDMA), 2,2'-(methylazanediyl)bis( N, N'-dimethylacetamide) (MABDMA), and 2,2'-(benzylazanediyl)bis( N, N'-dimethylacetamide) (BnABDMA), in aqueous solutions was investigated. The stability constants of two successive complexes, namely, NpO 2 L + and NpO 2 L 2 + , where L stands for the ligands, were determined by absorption spectrophotometry. The results suggest that the stability constants of corresponding Np(V) complexes follow the trend: MABDMA > ABDMA ≈ BnABDMA. The data are discussed in terms of the basicity of the ligands and compared with those for the complexation of Np(V) with an ether oxygen-linked diacetamide ligand. Extended X-ray absorption fine structure data indicate that, similar to the complexation with Nd 3+ and UO 2 2+ , the ligands coordinate to NpO 2 + in a tridentate mode through the amine nitrogen and two oxygen atoms of the amide groups. Furthermore, computational results, in conjunction with spectrophotometric data, verify that the 1:2 complexes (NpO 2 (L) 2 + ) in aqueous solutions are highly symmetric with Np at the inversion center, so that the f-f transition of Np(V) is forbidden and NpO 2 (L) 2 + does not display significant absorption in the near-IR region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Homoleptic Imidophosphorane Stabilization of Tetravalent Cerium

We present that the homoleptic complexes of cerium with the tris(piperidinyl)imidophosphorane ligand, [NP(pip) 3 ] - , present the most negative Ce 3+ /4+ redox couple known (<-2.64 V vs Fc/Fc + ). This dramatic stabilization of the cerium tetravalent oxidation state [>4.0 V shift from the Ce 3+/4+ couple in 1 M HClO 4 (aq)] is established through reactivity studies. Spectroscopic studies (UV–vis, electron paramagnetic resonance, and Ce L 3 -edge X-ray absorption spectroscopy), in conjunction with density functional theory studies, reveal the dominant covalent metal–ligand interactions underlying the observed redox chemistry and the dependence of the redox potential on the binding of potassium in the inner coordination sphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of In Silico Multifeature Libraries for Providing Evidence for the Presence of Small Molecules in Synthetic Blinded Samples

The current gold standard for unambiguous molecular identification in metabolomics analysis is comparing two or more orthogonal properties from the analysis of authentic reference materials (standards) to experimental data acquired in the same laboratory with the same analytical methods. This represents a significant limitation for comprehensive chemical identification of small molecules in complex samples. The process is time-consuming and costly, and the majority of molecules are not yet represented by standards. Thus, there is a need to assemble evidence for the presence of small molecules in complex samples through the use of libraries containing calculated chemical properties. To address this need, we developed a Multi-Attribute Matching Engine (MAME) and a library derived in-part from our in silico Chemical Library Engine (ISiCLE). Here, we describe an initial evaluation of these methods in a blinded analysis of synthetic chemical mixtures as part of the U.S. Environmental Protection Agency’s (EPA) Non-Targeted Analysis Collaborative Trial (ENTACT, Phase 1). For molecules in all mixtures, the initial blinded false negative rate (FNR), false discovery rate (FDR), and accuracy were 57%, 77%, and 91%, respectively. For high evidence scores, the FDR was 35%. After unblinding of the sample compositions, we optimized the scoring parameters to better exploit the available evidence, and increased the accuracy for molecules suspected as present. The final FNR, FDR, and accuracy were 67%, 53%, and 96%, respectively. For high evidence scores, the FDR was 10%. This study demonstrates that multi-attribute matching methods in conjunction with in silico libraries may one day enable reduced reliance on experimentally-derived libraries for building evidence for the presence of molecules in complex samples.

Nunez, Jamie↗

Ab Initio Study of Vibronic and Magnetic 5f-to-5f and Dipole-Allowed 5f-to-6d and Charge-Transfer Transitions in [UX 6 ] n – (X = Cl, Br; n = 1, 2)

The absorption spectra of the octahedral [UX 6 ] n– (X = Cl, Br; n = 1, 2) complexes in the near-infrared (NIR) and UV–vis spectral regions were studied theoretically, using a relativistic restricted active space second-order perturbation theory (RASPT2) wavefunction framework, with the spin–orbit (SO) coupling treated by state interaction, in conjunction with Kohn–Sham density functional theory calculations for determining the vibrational normal modes. The electric-dipole-allowed ligand-to-metal charge-transfer (LMCT) and 5f-to-6d transitions, and the electric-dipole-forbidden 5f-to-5f ligand field (LF) transitions, are thereby obtained within the same theoretical framework. For the 5f-to-5f LF transitions, the observed absorption intensity is mostly due to vibronic coupling with low-energy electric-dipole-allowed transitions, but in some cases, the magnetic dipole intensity of the purely electronic transition has comparable intensity to the vibronic transitions. Here, experimental LF spectra of 5f 2 open-shell systems have been reported decades back, but ab initio calculations of their vibronic intensity have not yet been reported in the literature. Although the LF spectra for the 5f 2 systems can be assigned in detail, based on the calculations, the spectra are very complex and the underlying electronic states are strongly multiconfigurational. Therefore, the usefulness of the LF spectra beyond serving as a “fingerprint” of the LF and the metal oxidation state appears to be limited.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron–electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. In conclusion, the transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Importance of Nuclear Quantum Effects on Aqueous Electrolyte Transport under Confinement in Ti 3 C 2 MXenes

Protons display a high chemical activity and strongly affect the charge storage capability in confined interlayer spaces of two-dimensional (2D) materials. As such, an accurate representation of proton dynamics under confinement is important for understanding and predicting charge storage dynamics in these materials. While often ignored in atomistic-scale simulations, nuclear quantum effects (NQEs), e.g., tunneling, can be significant under confinement even at room temperature. Here, using the thermostatted ring polymer molecular dynamics implementation of path integral molecular dynamics (PIMD) in conjunction with the ReaxFF force field, density functional tight binding (DFTB), and NequIP neural network potential simulations, we investigate the role of NQEs on proton and water transport in bulk water and aqueous electrolytes under confinement in Ti 3 C 2 MXenes. Although overall NQEs are relatively small, especially in bulk, we find that they can alter both quantitative values and qualitative trends on both proton transport and water self-diffusion under confinement relative to classical MD predictions. Therefore, our results suggest the need for NQEs to be considered to simulate aqueous systems under confinement for both qualitative and quantitative accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗